ROMS-SWAN problem

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
mistral
Posts: 23
Joined: Wed Jan 16, 2008 11:17 pm
Location: Privat

ROMS-SWAN problem

#1 Unread post by mistral »

Hello everyone
some days I'm trying to work with 2 models coupled.
test cases (inlet_test , test_head) are executed with a good outcome.
There are problems when using a realistic case.
The following attached file and error.
I hope very much that can help me.



NTIMES == 1440
DT == 30.0d0
NDTFAST == 20
.
.
.
DSTART = 0.0d0 ! days
TIDE_START = 0.0d0 ! days
TIME_REF = 0.0d0 ! yyyymmdd.dd

...
TimeInterval = 60.0d0
...

swan.in


PROJECT 'TEST HEAD' ''
'ocean_test_head test case to drive ROMS'
''
''

MODE NONSTATIONARY TWODIMENSIONAL
SET PWTAIL 5 INRHOG 1 NAUTICAL
COORDINATES SPHERICAL

&& KEYWORDS TO CREATE AND READ COMPUTATIONAL GRID &&
CGRID CURVILINEAR 113 155 EXC 9.999000e+003 9.999000e+003 &
CIRCLE 36 0.04 1.0 30
READGRID COORDINATES 1 'Waves/SWAN/External/maratea_grid_coord.grd' 4 0 0 FREE

&& KEYWORDS TO CREATE AND READ BATHYMETRY GRID &&
INPGRID BOTTOM CURVILINEAR 0 0 113 155 EXC 9.999000e+003
READINP BOTTOM 1 'Waves/SWAN/External/maratea_bathy.bot' 4 0 FREE

&& KEYWORD TO CREATE CURRENT GRID &&
INPGRID CURRENT CURVILINEAR 0 0 113 155 EXC 9.999000e+003 &
NONSTATIONARY 20050401.0000 1 DAY 20050402.0000

&& KEYWORD TO CREATE WATER LEVEL GRID &&
INPGRID WLEV CURVILINEAR 0 0 113 155 EXC 9.999000e+003 &
NONSTATIONARY 20050401.0000 1 DAY 20050402.0000

&& KEYWORD TO CREATE FRICTION LEVEL GRID &&
INPGRID FRIC CURVILINEAR 0 0 113 155 EXC 9.999000e+003 &
NONSTATIONARY 20050401.0000 1 DAY 20050402.0000


BOUNDPAR1 SHAPESPEC JONSWAP 3.3 PEAK DSPR DEGREES
BOUNDPAR2 SEGMENT IJ 0 27 57 27 CONSTANT PAR 2.0 10.0 0. 20.

& Restart name **********************************
&INIT HOTSTART 'swan_restart.dat'

& PHYSICS ****************************************
&BREAKING CONSTANT 1.0 0.73
FRICTION MADSEN 0.05
OFF QUAD
GEN3

PROP BSBT

BLOCK 'COMPGRID' NOHEADER 'depth.mat' LAY 4 DEPTH 1. OUTPUT 20050401.000000 2 HR
BLOCK 'COMPGRID' NOHEADER 'dissip.mat' LAY 4 DISSIP 1. OUTPUT 20050401.000000 2 HR
BLOCK 'COMPGRID' NOHEADER 'force.mat' LAY 4 FORCE 1. OUTPUT 20050401.000000 2 HR
BLOCK 'COMPGRID' NOHEADER 'fric.mat' LAY 4 FRCOEEF 1. OUTPUT 20050401.000000 2 HR
BLOCK 'COMPGRID' NOHEADER 'hsig.mat' LAY 4 HSIGN 1. OUTPUT 20050401.000000 2 HR
BLOCK 'COMPGRID' NOHEADER 'rtp.mat' LAY 4 RTP 1. OUTPUT 20050401.000000 2 HR
BLOCK 'COMPGRID' NOHEADER 'tm01.mat' LAY 4 TM01 1. OUTPUT 20050401.000000 2 HR
BLOCK 'COMPGRID' NOHEADER 'setup.mat' LAY 4 SETUP 1. OUTPUT 20050401.000000 2 HR
BLOCK 'COMPGRID' NOHEADER 'tmbot.mat' LAY 4 TMBOT 1. OUTPUT 20050401.000000 2 HR
BLOCK 'COMPGRID' NOHEADER 'ubot.mat' LAY 4 UBOT 1. OUTPUT 20050401.000000 2 HR
BLOCK 'COMPGRID' NOHEADER 'vel.mat' LAY 4 VEL 1. OUTPUT 20050401.000000 2 HR
BLOCK 'COMPGRID' NOHEADER 'wdir.mat' LAY 4 DIR 1. OUTPUT 20050401.000000 2 HR
BLOCK 'COMPGRID' NOHEADER 'wlen.mat' LAY 4 WLEN 1. OUTPUT 20050401.000000 2 HR
BLOCK 'COMPGRID' NOHEADER 'xp.mat' LAY 4 XP 1. OUTPUT 20050401.000000 2 HR
BLOCK 'COMPGRID' NOHEADER 'yp.mat' LAY 4 YP 1. OUTPUT 20050401.000000 2 HR

COMPUTE NONSTAT NONSTATIONARY 20050401.0000 60 SEC 20050402.0000

STOP
anto@aaa:~/aaa/ROMSanto> mpirun -np 2 ./oceanM ROMS/External/coupling_test_head.in
WARNING: Can't read mpd.hosts for list of hosts, start only on current
Coupled Input File name = ROMS/External/coupling_test_head.in

Model Coupling Parallel Threads:

Ocean Model MPI nodes: 000 - 000
Waves Model MPI nodes: 001 - 001


Ocean Export: bath:SSH:Ubar:Vbar:ZO
Waves Export: Wdir:Wamp:Wlen:Wptop:Wpbot:Wubot


Process Information:

Node # 0 (pid= 21797) is active.

Model Input Parameters: ROMS/TOMS version 3.4
Thursday - February 25, 2010 - 8:03:53 PM
-----------------------------------------------------------------------------

Sediment Test Headland Case

Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /opt/intel/ict/3.0/mpi/3.0/bin64/mpif90
Compiler flags : -frepack-arrays -O3 -ffast-math -I/home/anto/aaa/ROMSanto/MCT_1/mct -I/home/anto/aaa/ROMSanto/MCT_1/mpeu -ffree-form -ffree-line-length-none

Input Script : ROMS/External/ocean_test_head.in

SVN Root URL : https://www.myroms.org/svn/src/trunk
SVN Revision : exportado

Local Root : /home/anto/aaa/ROMSanto
Header Dir : /home/anto/aaa/ROMSanto/ROMS/Include
Header file : test_head.h
Analytical Dir: /home/anto/aaa/ROMSanto/ROMS/Functionals

Resolution, Grid 01: 0112x0154x030, Parallel Nodes: 1, Tiling: 001x001


Physical Parameters, Grid: 01
=============================

1440 ntimes Number of timesteps for 3-D equations.
30.000 dt Timestep size (s) for 3-D equations.
20 ndtfast Number of timesteps for 2-D equations between
each 3D timestep.
1 ERstr Starting ensemble/perturbation run number.
1 ERend Ending ensemble/perturbation run number.
0 nrrec Number of restart records to read from disk.
T LcycleRST Switch to recycle time-records in restart file.
120 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
T ldefout Switch to create a new output NetCDF file(s).
120 nHIS Number of timesteps between the writing fields
into history file.
1.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
1.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)
for tracer 02: salt
1.0000E-05 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
5.0000E-06 Akk_bak Background vertical mixing coefficient (m2/s)
for turbulent energy.
5.0000E-06 Akp_bak Background vertical mixing coefficient (m2/s)
for turbulent generic statistical field.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
3.0000E-03 Zob Bottom roughness (m).
1 Vtransform S-coordinate transformation equation.
1 Vstretching S-coordinate stretching function.
1.0000E+00 theta_s S-coordinate surface control parameter.
1.0000E+00 theta_b S-coordinate bottom control parameter.
20.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1024.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
0.000 dstart Time-stamp assigned to model initialization (days).
20050401.00 time_ref Reference time for units attribute (yyyymmdd.dd)
0.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
0.0000E+00 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
0.0000E+00 Znudg Nudging/relaxation time scale (days)
for free-surface.
0.0000E+00 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
0.0000E+00 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
0.0000E+00 obcfac Factor between passive and active
open boundary conditions.
10.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
1024.000 R0 Background density (kg/m3) used in linear Equation
of State.
1.7000E-04 Tcoef Thermal expansion coefficient (1/Celsius).
7.6000E-04 Scoef Saline contraction coefficient (1/PSU).
-1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
T Hout(idBath) Write out time-dependent bathymetry.
T Hout(idFsur) Write out free-surface.
T Hout(idUbar) Write out 2D U-momentum component.
T Hout(idVbar) Write out 2D V-momentum component.
T Hout(idUvel) Write out 3D U-momentum component.
T Hout(idVvel) Write out 3D V-momentum component.
T Hout(idWvel) Write out W-momentum component.
T Hout(idOvel) Write out omega vertical velocity.
T Hout(idTvar) Write out tracer 01: temp
T Hout(idTvar) Write out tracer 02: salt
T Hout(idUsms) Write out surface U-momentum stress.
T Hout(idVsms) Write out surface V-momentum stress.
T Hout(idUbms) Write out bottom U-momentum stress.
T Hout(idVbms) Write out bottom V-momentum stress.
T Hout(idUbrs) Write out bottom U-current stress.
T Hout(idVbrs) Write out bottom V-current stress.
T Hout(idUbws) Write out wind-induced, bottom U-wave stress.
T Hout(idVbws) Write out wind-induced, bottom V-wave stress.
T Hout(idUbcs) Write out max wind + current, bottom U-wave stress.
T Hout(idVbcs) Write out max wind + current, bottom V-wave stress.
T Hout(idUbot) Write out bed wave orbital U-velocity.
T Hout(idVbot) Write out bed wave orbital V-velocity.
T Hout(idUbur) Write out bottom U-momentum above bed.
T Hout(idVbvr) Write out bottom V-momentum above bed.
T Hout(idW2xx) Write out 2D radiation stress, Sxx.
T Hout(idW2xy) Write out 2D radiation stress, Sxy.
T Hout(idW2yy) Write out 2D radiation stress, Syy.
T Hout(idU2rs) Write out total 2D u-radiation stress.
T Hout(idV2rs) Write out total 2D v-radiation stress.
T Hout(idU2Sd) Write out 2D u-momentum stokes velocity.
T Hout(idV2Sd) Write out 2D v-momentum stokes velocity.
T Hout(idW3xx) Write out 3D horizonrtal radiation stress, Sxx.
T Hout(idW3xy) Write out 3D horizonrtal radiation stress, Sxy.
T Hout(idW3yy) Write out 3D horizonrtal radiation stress, Syy.
T Hout(idW3zx) Write out 3D vertical radiation stress, Szx.
T Hout(idW3zy) Write out 3D vertical radiation stress, Szy.
T Hout(idU3rs) Write out total 3D u-radiation stress.
T Hout(idV3rs) Write out total 3D v-radiation stress.
T Hout(idU3Sd) Write out 3D u-momentum stokes velocity.
T Hout(idV3Sd) Write out 3D v-momentum stokes velocity.
T Hout(idWamp) Write out wave height.
T Hout(idWlen) Write out wavelength.
T Hout(idWdir) Write out wave direction.
T Hout(idWptp) Write out wave surface period.
T Hout(idWpbt) Write out wave bottom period.
T Hout(idWorb) Write out wave bottom orbital velocity.
T Hout(idWdis) Write out wave dissipation.
T Hout(idBott) Write out bottom property 01: grain_diameter
T Hout(idBott) Write out bottom property 02: grain_density
T Hout(idBott) Write out bottom property 03: settling_vel
T Hout(idBott) Write out bottom property 04: erosion_stress
T Hout(idBott) Write out bottom property 05: ripple_length
T Hout(idBott) Write out bottom property 06: ripple_height
T Hout(idBott) Write out bottom property 07: bed_wave_amp
T Hout(idBott) Write out bottom property 09: Zo_app
T Hout(idBott) Write out bottom property 10: Zo_Nik
T Hout(idBott) Write out bottom property 11: Zo_bio
T Hout(idBott) Write out bottom property 15: active_layer_thickness
T Hout(idTsur) Write out surface net heat flux.
T Hout(idDano) Write out density anomaly.
T Hout(idVvis) Write out vertical viscosity coefficient.
T Hout(idTdif) Write out vertical T-diffusion coefficient.

Output/Input Files:

Output Restart File: ocean_rst.nc
Output History File: ocean_his.nc
Physical parameters File: ROMS/External/ocean_test_head.in
Input Grid File: Data/ROMS/Grid/maratea-grid.nc
Input Nonlinear Initial File: Data/ROMS/Initial/maratea_init.nc

Tile partition information for Grid 01: 0112x0154x0030 tiling: 001x001

tile Istr Iend Jstr Jend Npts

0 1 112 1 154 517440

Tile minimum and maximum fractional grid coordinates:
(interior points only)

tile Xmin Xmax Ymin Ymax grid

0 0.50 112.50 0.50 154.50 RHO-points

0 1.00 112.00 0.50 154.50 U-points

0 0.50 112.50 1.00 154.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 369
HaloSizeJ(1) = 495
TileSide(1) = 159
TileSize(1) = 18603



Sediment Parameters, Grid: 01
=============================


Size Sd50 Csed Srho Wsed Erate poros
Class (mm) (kg/m3) (kg/m3) (mm/s) (kg/m2/s) (nondim)

1 2.5000E-01 0.0000E+00 2.6500E+03 3.0000E+01 5.0000E-04 9.0000E-01

tau_ce tau_cd nl_tnu2 nl_tnu4 Akt_bak Tnudg
(N/m2) (N/m2) (m2/s) (m4/s) (m2/s) (day)

1 1.9000E-01 6.0000E-02 0.0000E+00 0.0000E+00 5.0000E-06 0.0000E+00

morph_fac
(nondim)

1 1.0000E+01

New bed layer formed when deposition exceeds 0.10000E-01 (m).
Two first layers are combined when 2nd layer smaller than 0.00000E+00 (m).
Rate coefficient for bed load transport = 0.50000E-01

T Hout(idTvar) Write out sediment01: mud_01
T Hout(idfrac) Write out bed fraction, sediment 01: mudfrac_01
T Hout(idfrac) Write out mass, sediment 01: mudmass_01
T Hout(idUbld) Write out bed load at U-points, sediment 01: bedload_Umud_01
T Hout(idVbld) Write out bed load at V-points, sediment 01: bedload_Vmud_01
T Hout(idSbed) Write out BED property 01: bed_thickness
T Hout(idSbed) Write out BED property 02: bed_age
T Hout(idSbed) Write out BED property 03: bed_porosity
T Hout(idSbed) Write out BED property 04: bed_biodiff

Activated C-preprocessing Options:

TEST_HEAD Sediment Test Headland Case
ANA_BPFLUX Analytical bottom passive tracers fluxes.
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_SMFLUX Analytical kinematic surface momentum flux.
ANA_SPFLUX Analytical surface passive tracer fluxes.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ANA_STFLUX Analytical kinematic surface temperature flux.
ASSUMED_SHAPE Using assumed-shape arrays.
BEDLOAD_SOULSBY Activate bed load sediment transport Soulsby formula.
CURVGRID Orthogonal curvilinear grid.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
EASTERN_WALL Wall boundary at Eastern edge.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
MASKING Land/Sea masking.
MCT_LIB Using Model Coupling Toolkit library.
MPI MPI distributed-memory configuration.
MY25_MIXING Mellor/Yamada Level-2.5 mixing closure.
NEARSHORE_MELLOR Nearshore RAdiation Stress Terms.
NONLINEAR Nonlinear Model.
!NONLIN_EOS Linear Equation of State for seawater.
NORTHERN_WALL Wall boundary at Northern edge.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
OUT_DOUBLE Double precision output fields in NetCDF files.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
K_GSCHEME Third-order upstream advection of TKE fields.
!RST_SINGLE Double precision fields in restart NetCDF file.
SEDIMENT Cohesive and noncohesive sediments.
SED_MORPH Allow bottom model elevation to evolve.
SUSPLOAD Activate suspended sediment transport.
SOLVE3D Solving 3D Primitive Equations.
SOUTHERN_WALL Wall boundary at Southern edge.
SPLINES Conservative parabolic spline reconstruction.
SSW_BBL Styles and Glenn Bottom Boundary Layer - modified.
SSW_CALC_ZNOT Internal computation of bottom roughness.
SWAN_COUPLING Two-way SWAN/ROMS coupling.
THREE_GHOST Using three Ghost Points in halo regions.
TS_MPDATA Recursive flux corrected MPDATA 3D advection of tracers.
UV_ADV Advection of momentum.
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum.
UV_C4VADVECTION Fourth-order centered vertical advection of momentum.
VAR_RHO_2D Variable density barotropic mode.
WAVES_OCEAN Two-way wave-ocean models coupling.
WESTERN_WALL Wall boundary at Western edge.

INITIAL: Configuring and initializing forward nonlinear model ...


Vertical S-coordinate System:

level S-coord Cs-curve at_hmin over_slope at_hmax

30 0.0000000 0.0000000 0.000 0.000 0.000
29 -0.0333333 -0.0287969 -0.033 -1.098 -2.162
28 -0.0666667 -0.0584418 -0.067 -2.227 -4.386
27 -0.1000000 -0.0889041 -0.100 -3.386 -6.671
26 -0.1333333 -0.1201476 -0.133 -4.574 -9.014
25 -0.1666667 -0.1521307 -0.167 -5.789 -11.412
24 -0.2000000 -0.1848065 -0.200 -7.030 -13.860
23 -0.2333333 -0.2181229 -0.233 -8.295 -16.356
22 -0.2666667 -0.2520228 -0.267 -9.581 -18.895
21 -0.3000000 -0.2864445 -0.300 -10.887 -21.473
20 -0.3333333 -0.3213219 -0.333 -12.209 -24.084
19 -0.3666667 -0.3565853 -0.367 -13.545 -26.724
18 -0.4000000 -0.3921616 -0.400 -14.894 -29.387
17 -0.4333333 -0.4279749 -0.433 -16.251 -32.068
16 -0.4666667 -0.4639475 -0.467 -17.613 -34.760
15 -0.5000000 -0.5000000 -0.500 -18.979 -37.458
14 -0.5333333 -0.5360525 -0.533 -20.345 -40.157
13 -0.5666667 -0.5720251 -0.567 -21.708 -42.849
12 -0.6000000 -0.6078384 -0.600 -23.065 -45.529
11 -0.6333333 -0.6434147 -0.633 -24.413 -48.192
10 -0.6666667 -0.6786781 -0.667 -25.749 -50.832
9 -0.7000000 -0.7135555 -0.700 -27.072 -53.444
8 -0.7333333 -0.7479772 -0.733 -28.377 -56.021
7 -0.7666667 -0.7818771 -0.767 -29.664 -58.560
6 -0.8000000 -0.8151935 -0.800 -30.928 -61.056
5 -0.8333333 -0.8478693 -0.833 -32.169 -63.505
4 -0.8666667 -0.8798524 -0.867 -33.385 -65.902
3 -0.9000000 -0.9110959 -0.900 -34.573 -68.245
2 -0.9333333 -0.9415582 -0.933 -35.732 -70.530
1 -0.9666667 -0.9712031 -0.967 -36.861 -72.755
0 -1.0000000 -1.0000000 -1.000 -37.958 -74.917

Time Splitting Weights: ndtfast = 20 nfast = 29

Primary Secondary Accumulated to Current Step

1-0.0009651193358779 0.0500000000000000-0.0009651193358779 0.0500000000000000
2-0.0013488780126037 0.0500482559667939-0.0023139973484816 0.1000482559667939
3-0.0011514592651645 0.0501156998674241-0.0034654566136461 0.1501639558342180
4-0.0003735756740661 0.0501732728306823-0.0038390322877122 0.2003372286649003
5 0.0009829200513762 0.0501919516143856-0.0028561122363360 0.2505291802792859
6 0.0029141799764308 0.0501428056118168 0.0000580677400948 0.3006719858911027
7 0.0054132615310267 0.0499970966129953 0.0054713292711214 0.3506690825040980
8 0.0084687837865132 0.0497264335364439 0.0139401130576347 0.4003955160405419
9 0.0120633394191050 0.0493029943471183 0.0260034524767396 0.4496985103876602
10 0.0161716623600090 0.0486998273761630 0.0421751148367486 0.4983983377638232
11 0.0207585511322367 0.0478912442581626 0.0629336659689853 0.5462895820219857
12 0.0257765478740990 0.0468533167015507 0.0887102138430844 0.5931428987235364
13 0.0311633730493853 0.0455644893078458 0.1198735868924697 0.6387073880313822
14 0.0368391158442262 0.0440063206553765 0.1567127027366959 0.6827137086867586
15 0.0427031802506397 0.0421643648631652 0.1994158829873356 0.7248780735499238
16 0.0486309868367616 0.0400292058506332 0.2480468698240972 0.7649072794005570
17 0.0544704302037592 0.0375976565087951 0.3025173000278564 0.8025049359093521
18 0.0600380921294286 0.0348741349986072 0.3625553921572851 0.8373790709079593
19 0.0651152103984763 0.0318722303921357 0.4276706025557614 0.8692513013000950
20 0.0694434033194840 0.0286164698722119 0.4971140058752453 0.8978677711723069
21 0.0727201499285570 0.0251442997062377 0.5698341558038024 0.9230120708785446
22 0.0745940258796570 0.0215082922098099 0.6444281816834594 0.9445203630883545
23 0.0746596950216180 0.0177785909158270 0.7190878767050775 0.9622989540041815
24 0.0724526566618460 0.0140456061647461 0.7915405333669234 0.9763445601689277
25 0.0674437485167026 0.0104229733316538 0.8589842818836260 0.9867675335005816
26 0.0590334053485720 0.0070507859058187 0.9180176872321980 0.9938183194064002
27 0.0465456732896125 0.0040991156383901 0.9645633605218105 0.9979174350447902
28 0.0292219798521904 0.0017718319739095 0.9937853403740009 0.9996892670186998
29 0.0062146596259990 0.0003107329812999 0.9999999999999999 0.9999999999999997

ndtfast, nfast = 20 29 nfast/ndtfast = 1.45000

Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):

1.000000000000 1.060707743385 0.530353871693 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.14200

Minimum X-grid spacing, DXmin = 5.95474669E-03 km
Maximum X-grid spacing, DXmax = 3.27720523E-02 km
Minimum Y-grid spacing, DYmin = 7.34701695E-03 km
Maximum Y-grid spacing, DYmax = 3.80203002E-02 km
Minimum Z-grid spacing, DZmin = 3.33333333E-02 m
Maximum Z-grid spacing, DZmax = 2.69821760E+00 m

Minimum barotropic Courant Number = 4.10503999E-01
Maximum barotropic Courant Number = 2.91159817E+00
Maximum Coriolis Courant Number = 2.81394738E-03


NLM: GET_STATE - Read state initial conditions, t = 13462 00:00:00
(File: maratea_init.nc, Rec=0001, Index=1)
- free-surface
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- potential temperature
(Min = 1.00000000E+01 Max = 1.00000000E+01)
- salinity
(Min = 3.00000000E+01 Max = 3.00000000E+01)
- suspended cohesive sediment, size class 01
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- cohesive sediment fraction, size class 01
(Min = 1.00000000E+00 Max = 1.00000000E+00)
- cohesive sediment mass, size class 01
(Min = 1.02300000E+03 Max = 1.02300000E+03)
- sediment bed layer thickness
(Min = 1.00000000E+00 Max = 1.00000000E+00)
- sediment layer age
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- sediment layer porosity
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- biodiffusivity at bottom of each layer
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- sediment median grain diameter size
(Min = 2.00000000E-04 Max = 2.00000000E-04)
- sediment median grain density
(Min = 1.02300000E+03 Max = 1.02300000E+03)
- sediment median grain settling velocity
(Min = 1.00000000E-02 Max = 1.00000000E-02)
- sediment median critical erosion stress
(Min = 1.95503421E-04 Max = 1.95503421E-04)

Maximum grid stiffness ratios: rx0 = 3.637517E-01 (Beckmann and Haidvogel)
rx1 = 2.422988E+01 (Haney)


Initial basin volumes: TotVolume = 1.6590147525E+08 m3
MinVolume = 1.2431839364E+01 m3
MaxVolume = 1.8561442070E+03 m3
Max/Min = 1.4930567817E+02

OCN2WAV - (06) imported and (05) exported fields, t = 13462 00:00:00
- ROMS coupling exchanges wait clock (s):
(Recv= 4.60428700E+00 Send= 4.00000000E-03)
- ROMS Import: wave direction
(Min= 0.00000000E+00 Max= 1.14873150E-07)
- ROMS Import: significant wave height
(Min= 0.00000000E+00 Max= 1.99365467E-02)
- ROMS Import: average wave length
(Min= 0.00000000E+00 Max= 3.80385232E+00)
- ROMS Import: surface wave relative peak period
(Min= 0.00000000E+00 Max= 1.90365243E+00)
- ROMS Import: bottom wave period
(Min= 0.00000000E+00 Max= 0.00000000E+00)
- ROMS Import: wave bottom orbital velocity
(Min= 0.00000000E+00 Max= 0.00000000E+00)
- ROMS Export: bathymetry
(Min= 1.00000000E+00 Max= 7.49166890E+01)
- ROMS Export: free-surface
(Min= 0.00000000E+00 Max= 0.00000000E+00)
- ROMS Export: vertically integrated u-momentum component
(Min= 0.00000000E+00 Max= 0.00000000E+00)
- ROMS Export: vertically integrated v-momentum component
(Min= 0.00000000E+00 Max= 0.00000000E+00)
- ROMS Export: bottom roughness
(Min= 1.00000000E-04 Max= 1.00000000E-04)


NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 38770561 - 00001440)

Elapsed CPU time (seconds):

Node # 0 CPU: 5.688
Total: 5.688

Nonlinear model elapsed time profile:

Initialization ................................... 4.892 (86.0056 %)
2D/3D coupling, vertical metrics ................. 0.024 ( 0.4219 %)
Omega vertical velocity ..........................


SWAN is preparing computation ...

WAV2OCN - (05) imported and (06) exported fields, t = 20050401.000000
- SWAN coupling exchanges wait clock (s):
(Recv= 8.00200000E-03 Send= 0.00000000E+00)

I have the impression that this is an error on the start time of the two models.

What do you think

Best regards

hyunmin.eom
Posts: 7
Joined: Tue Jan 06, 2009 3:58 pm
Location: Inha University

Re: ROMS-SWAN problem

#2 Unread post by hyunmin.eom »

Dear mistral

Now timestep of ROMS is from 38770561 to 00001440.
...NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 38770561 - 00001440)...
Because initial condition time is 13462 00:00:00.
NLM: GET_STATE - Read state initial conditions, t = 13462 00:00:00
(File: maratea_init.nc, Rec=0001, Index=1)
Maybe this is the cause.

Good luck to you.
Hyunmin.

mistral
Posts: 23
Joined: Wed Jan 16, 2008 11:17 pm
Location: Privat

Re: ROMS-SWAN problem

#3 Unread post by mistral »

ok...thanks for your advice
I changed the start time and now it works,
but I have a persistent error, this error occurs even with the test case ...
depends on what you know??


...
86397 13462 23:59:57 7.221366E-04 2.485874E+02 2.485881E+02 1.711143E+08
86398 13462 23:59:58 7.221378E-04 2.485874E+02 2.485881E+02 1.711143E+08
86399 13462 23:59:59 7.221389E-04 2.485874E+02 2.485881E+02 1.711143E+08
86400 13463 00:00:00 7.221401E-04 2.485874E+02 2.485881E+02 1.711143E+08
WAV2OCN - (05) imported and (06) exported fields, t = 20050402.000000
- SWAN coupling exchanges wait clock (s):
(Recv= 2.04396773E+02 Send= 4.00000000E-03)




At line 3791 of file swanparll.f90 (unit = 454, file = 'vel.mat-001')
Fortran runtime error: Bad value during floating point read
rank 26 in job 1 aaa_33055 caused collective abort of all ranks
exit status of rank 26: return code 2

Post Reply