I'm trying to run COAWST (coupling ROMS and SWAN) using free-surface information from HYCOM at the boundary.
In my first test, I've created a file with information just in the East boundary and it worked.
For a second test, I've created a file with information also in the North and South boundary, but in this case the model blows immediately during the first time step. I've set the boundaries to be Chapman and Flather.
I've tried to decrease the time step, but even using 1 second it still blows in the first time step. I've looked at the results for this first time step, trying to figure what is causing the blow, but I can't seem to find anything out of normal.
Do you guys have any idea about what could be causing it or at least how can I try to figure it out?
Below is the log generated
Thanks a million!
Code: Select all
=============================================================
Coupled Input File name = coupling_inlet_test.in
Model Coupling:
Ocean Model MPI nodes: 000 - 063
Waves Model MPI nodes: 064 - 087
Process Information:
Node # 0 (pid= 18176) is active.
Node # 4 (pid= 18180) is active.
Node # 7 (pid= 18183) is active.
Node # 10 (pid= 18186) is active.
Node # 11 (pid= 18187) is active.
Node # 13 (pid= 18189) is active.
Node # 14 (pid= 18190) is active.
Node # 15 (pid= 18191) is active.
Node # 16 (pid= 18192) is active.
Node # 17 (pid= 18193) is active.
Node # 18 (pid= 18194) is active.
Node # 19 (pid= 18195) is active.
Node # 20 (pid= 18196) is active.
Node # 21 (pid= 18197) is active.
Node # 22 (pid= 18198) is active.
Node # 23 (pid= 18199) is active.
Node # 1 (pid= 18177) is active.
Node # 2 (pid= 18178) is active.
Node # 3 (pid= 18179) is active.
Node # 5 (pid= 18181) is active.
Node # 6 (pid= 18182) is active.
Node # 8 (pid= 18184) is active.
Node # 9 (pid= 18185) is active.
Node # 12 (pid= 18188) is active.
Node # 48 (pid= 10489) is active.
Node # 49 (pid= 10490) is active.
Node # 50 (pid= 10491) is active.
Node # 51 (pid= 10492) is active.
Node # 52 (pid= 10493) is active.
Node # 53 (pid= 10494) is active.
Node # 54 (pid= 10495) is active.
Node # 55 (pid= 10496) is active.
Node # 56 (pid= 10497) is active.
Node # 57 (pid= 10498) is active.
Node # 58 (pid= 10499) is active.
Node # 59 (pid= 10500) is active.
Node # 60 (pid= 10501) is active.
Node # 61 (pid= 10502) is active.
Node # 62 (pid= 10503) is active.
Node # 63 (pid= 10504) is active.
Node # 24 (pid= 22731) is active.
Node # 25 (pid= 22732) is active.
Node # 26 (pid= 22733) is active.
Node # 27 (pid= 22734) is active.
Node # 28 (pid= 22735) is active.
Node # 29 (pid= 22736) is active.
Node # 30 (pid= 22737) is active.
Node # 31 (pid= 22738) is active.
Node # 32 (pid= 22739) is active.
Node # 33 (pid= 22740) is active.
Node # 34 (pid= 22741) is active.
Node # 35 (pid= 22742) is active.
Node # 36 (pid= 22743) is active.
Node # 37 (pid= 22744) is active.
Node # 38 (pid= 22745) is active.
Node # 39 (pid= 22746) is active.
Node # 40 (pid= 22747) is active.
Node # 41 (pid= 22748) is active.
Node # 42 (pid= 22749) is active.
Node # 43 (pid= 22750) is active.
Node # 44 (pid= 22751) is active.
Node # 45 (pid= 22752) is active.
Node # 46 (pid= 22753) is active.
Node # 47 (pid= 22754) is active.
Model Input Parameters: ROMS/TOMS version 3.4
Thursday - June 18, 2015 - 6:25:44 PM
-----------------------------------------------------------------------------
SWAN grid 1 is preparing computation
delbayCC2012
Operating system : Linux
CPU/hardware : x86_64
Compiler system : pgi
Compiler command : /opt/shared/openmpi/1.4.4-pgi/bin/mpif90
Compiler flags : -fastsse -Mipa=fast -tp k8-64 -Mfree
Input Script : ocean.in
SVN Root URL : https://www.myroms.org/svn/omlab/branches/jcwarner
SVN Revision : exported
Local Root : /lustre/work/kukulka_lab/androd/COAWST/Rob_Delbay/RobCOAWST
Header Dir : /lustre/work/kukulka_lab/androd/COAWST/Rob_Delbay/RobCOAWST/Projects/delbayCC201
Header file : delbaycc2012.h
Analytical Dir: /lustre/work/kukulka_lab/androd/COAWST/Rob_Delbay/RobCOAWST/Projects/delbayCC201
Resolution, Grid 01: 0148x0298x010, Parallel Nodes: 64, Tiling: 008x008
Physical Parameters, Grid: 01
=============================
192961 ntimes Number of timesteps for 3-D equations.
150.000 dt Timestep size (s) for 3-D equations.
30 ndtfast Number of timesteps for 2-D equations between
each 3D timestep.
1800.000 TI_WAV_OCN Time interval (s) between coupling WAV-OCN models.
12 nOCN_WAV Number of OCN timesteps between coupling to WAV.
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.
576 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).
24 nHIS Number of timesteps between the writing fields
into history file.
17520 ndefHIS Number of timesteps between creation of new
history files.
5.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
5.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
0.0000E+00 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
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.5000E-04 rdrg Linear bottom drag coefficient (m/s).
5.0000E-03 rdrg2 Quadratic bottom drag coefficient.
0.0000E+00 Zob Bottom roughness (m).
1.0000E+00 Dcrit Minimum depth for wetting and drying (m).
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
3.0000E+00 theta_s S-coordinate surface control parameter.
4.0000E-01 theta_b S-coordinate bottom control parameter.
0.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
0.000 dstart Time-stamp assigned to model initialization (days).
-2191.042 tide_start Reference time origin for tidal forcing (days).
20120101.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.
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
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(idTvar) Write out tracer 02: salt
Output/Input Files:
Output Restart File: Output/ocean_rst.nc
Prefix for History Files: Output/ocean_his
Physical parameters File: ocean.in
Input Grid File: delbay_grd.nc
Input Nonlinear Initial File: delbayini2012CCNAM.nc
Input Boundary File: delbay_bry_hycom2012_ENSfsm2.nc
Tile partition information for Grid 01: 0148x0298x0010 tiling: 008x008
tile Istr Iend Jstr Jend Npts
0 1 17 1 35 5950
1 18 36 1 35 6650
2 37 55 1 35 6650
3 56 74 1 35 6650
4 75 93 1 35 6650
5 94 112 1 35 6650
6 113 131 1 35 6650
7 132 148 1 35 5950
8 1 17 36 73 6460
9 18 36 36 73 7220
10 37 55 36 73 7220
11 56 74 36 73 7220
12 75 93 36 73 7220
13 94 112 36 73 7220
14 113 131 36 73 7220
15 132 148 36 73 6460
16 1 17 74 111 6460
17 18 36 74 111 7220
18 37 55 74 111 7220
19 56 74 74 111 7220
20 75 93 74 111 7220
21 94 112 74 111 7220
22 113 131 74 111 7220
23 132 148 74 111 6460
24 1 17 112 149 6460
25 18 36 112 149 7220
26 37 55 112 149 7220
27 56 74 112 149 7220
28 75 93 112 149 7220
29 94 112 112 149 7220
30 113 131 112 149 7220
31 132 148 112 149 6460
32 1 17 150 187 6460
33 18 36 150 187 7220
34 37 55 150 187 7220
35 56 74 150 187 7220
36 75 93 150 187 7220
37 94 112 150 187 7220
38 113 131 150 187 7220
39 132 148 150 187 6460
40 1 17 188 225 6460
41 18 36 188 225 7220
42 37 55 188 225 7220
43 56 74 188 225 7220
44 75 93 188 225 7220
45 94 112 188 225 7220
46 113 131 188 225 7220
47 132 148 188 225 6460
48 1 17 226 263 6460
49 18 36 226 263 7220
50 37 55 226 263 7220
51 56 74 226 263 7220
52 75 93 226 263 7220
53 94 112 226 263 7220
54 113 131 226 263 7220
55 132 148 226 263 6460
56 1 17 264 298 5950
57 18 36 264 298 6650
58 37 55 264 298 6650
59 56 74 264 298 6650
60 75 93 264 298 6650
61 94 112 264 298 6650
62 113 131 264 298 6650
63 132 148 264 298 5950
Tile minimum and maximum fractional grid coordinates:
(interior points only)
tile Xmin Xmax Ymin Ymax grid
0 0.50 17.50 0.50 35.50 RHO-points
1 17.50 36.50 0.50 35.50 RHO-points
2 36.50 55.50 0.50 35.50 RHO-points
3 55.50 74.50 0.50 35.50 RHO-points
4 74.50 93.50 0.50 35.50 RHO-points
5 93.50 112.50 0.50 35.50 RHO-points
6 112.50 131.50 0.50 35.50 RHO-points
7 131.50 148.50 0.50 35.50 RHO-points
8 0.50 17.50 35.50 73.50 RHO-points
9 17.50 36.50 35.50 73.50 RHO-points
10 36.50 55.50 35.50 73.50 RHO-points
11 55.50 74.50 35.50 73.50 RHO-points
12 74.50 93.50 35.50 73.50 RHO-points
13 93.50 112.50 35.50 73.50 RHO-points
14 112.50 131.50 35.50 73.50 RHO-points
15 131.50 148.50 35.50 73.50 RHO-points
16 0.50 17.50 73.50 111.50 RHO-points
17 17.50 36.50 73.50 111.50 RHO-points
18 36.50 55.50 73.50 111.50 RHO-points
19 55.50 74.50 73.50 111.50 RHO-points
20 74.50 93.50 73.50 111.50 RHO-points
21 93.50 112.50 73.50 111.50 RHO-points
22 112.50 131.50 73.50 111.50 RHO-points
23 131.50 148.50 73.50 111.50 RHO-points
24 0.50 17.50 111.50 149.50 RHO-points
25 17.50 36.50 111.50 149.50 RHO-points
26 36.50 55.50 111.50 149.50 RHO-points
27 55.50 74.50 111.50 149.50 RHO-points
28 74.50 93.50 111.50 149.50 RHO-points
29 93.50 112.50 111.50 149.50 RHO-points
30 112.50 131.50 111.50 149.50 RHO-points
31 131.50 148.50 111.50 149.50 RHO-points
32 0.50 17.50 149.50 187.50 RHO-points
33 17.50 36.50 149.50 187.50 RHO-points
34 36.50 55.50 149.50 187.50 RHO-points
35 55.50 74.50 149.50 187.50 RHO-points
36 74.50 93.50 149.50 187.50 RHO-points
37 93.50 112.50 149.50 187.50 RHO-points
38 112.50 131.50 149.50 187.50 RHO-points
39 131.50 148.50 149.50 187.50 RHO-points
40 0.50 17.50 187.50 225.50 RHO-points
41 17.50 36.50 187.50 225.50 RHO-points
42 36.50 55.50 187.50 225.50 RHO-points
43 55.50 74.50 187.50 225.50 RHO-points
44 74.50 93.50 187.50 225.50 RHO-points
45 93.50 112.50 187.50 225.50 RHO-points
46 112.50 131.50 187.50 225.50 RHO-points
47 131.50 148.50 187.50 225.50 RHO-points
48 0.50 17.50 225.50 263.50 RHO-points
49 17.50 36.50 225.50 263.50 RHO-points
50 36.50 55.50 225.50 263.50 RHO-points
51 55.50 74.50 225.50 263.50 RHO-points
52 74.50 93.50 225.50 263.50 RHO-points
53 93.50 112.50 225.50 263.50 RHO-points
54 112.50 131.50 225.50 263.50 RHO-points
55 131.50 148.50 225.50 263.50 RHO-points
56 0.50 17.50 263.50 298.50 RHO-points
57 17.50 36.50 263.50 298.50 RHO-points
58 36.50 55.50 263.50 298.50 RHO-points
59 55.50 74.50 263.50 298.50 RHO-points
60 74.50 93.50 263.50 298.50 RHO-points
61 93.50 112.50 263.50 298.50 RHO-points
62 112.50 131.50 263.50 298.50 RHO-points
63 131.50 148.50 263.50 298.50 RHO-points
0 1.00 17.50 0.50 35.50 U-points
1 17.50 36.50 0.50 35.50 U-points
2 36.50 55.50 0.50 35.50 U-points
3 55.50 74.50 0.50 35.50 U-points
4 74.50 93.50 0.50 35.50 U-points
5 93.50 112.50 0.50 35.50 U-points
6 112.50 131.50 0.50 35.50 U-points
7 131.50 148.00 0.50 35.50 U-points
8 1.00 17.50 35.50 73.50 U-points
9 17.50 36.50 35.50 73.50 U-points
10 36.50 55.50 35.50 73.50 U-points
11 55.50 74.50 35.50 73.50 U-points
12 74.50 93.50 35.50 73.50 U-points
13 93.50 112.50 35.50 73.50 U-points
14 112.50 131.50 35.50 73.50 U-points
15 131.50 148.00 35.50 73.50 U-points
16 1.00 17.50 73.50 111.50 U-points
17 17.50 36.50 73.50 111.50 U-points
18 36.50 55.50 73.50 111.50 U-points
19 55.50 74.50 73.50 111.50 U-points
20 74.50 93.50 73.50 111.50 U-points
21 93.50 112.50 73.50 111.50 U-points
22 112.50 131.50 73.50 111.50 U-points
23 131.50 148.00 73.50 111.50 U-points
24 1.00 17.50 111.50 149.50 U-points
25 17.50 36.50 111.50 149.50 U-points
26 36.50 55.50 111.50 149.50 U-points
27 55.50 74.50 111.50 149.50 U-points
28 74.50 93.50 111.50 149.50 U-points
29 93.50 112.50 111.50 149.50 U-points
30 112.50 131.50 111.50 149.50 U-points
31 131.50 148.00 111.50 149.50 U-points
32 1.00 17.50 149.50 187.50 U-points
33 17.50 36.50 149.50 187.50 U-points
34 36.50 55.50 149.50 187.50 U-points
35 55.50 74.50 149.50 187.50 U-points
36 74.50 93.50 149.50 187.50 U-points
37 93.50 112.50 149.50 187.50 U-points
38 112.50 131.50 149.50 187.50 U-points
39 131.50 148.00 149.50 187.50 U-points
40 1.00 17.50 187.50 225.50 U-points
41 17.50 36.50 187.50 225.50 U-points
42 36.50 55.50 187.50 225.50 U-points
43 55.50 74.50 187.50 225.50 U-points
44 74.50 93.50 187.50 225.50 U-points
45 93.50 112.50 187.50 225.50 U-points
46 112.50 131.50 187.50 225.50 U-points
47 131.50 148.00 187.50 225.50 U-points
48 1.00 17.50 225.50 263.50 U-points
49 17.50 36.50 225.50 263.50 U-points
50 36.50 55.50 225.50 263.50 U-points
51 55.50 74.50 225.50 263.50 U-points
52 74.50 93.50 225.50 263.50 U-points
53 93.50 112.50 225.50 263.50 U-points
54 112.50 131.50 225.50 263.50 U-points
55 131.50 148.00 225.50 263.50 U-points
56 1.00 17.50 263.50 298.50 U-points
57 17.50 36.50 263.50 298.50 U-points
58 36.50 55.50 263.50 298.50 U-points
59 55.50 74.50 263.50 298.50 U-points
60 74.50 93.50 263.50 298.50 U-points
61 93.50 112.50 263.50 298.50 U-points
62 112.50 131.50 263.50 298.50 U-points
63 131.50 148.00 263.50 298.50 U-points
0 0.50 17.50 1.00 35.50 V-points
1 17.50 36.50 1.00 35.50 V-points
2 36.50 55.50 1.00 35.50 V-points
3 55.50 74.50 1.00 35.50 V-points
4 74.50 93.50 1.00 35.50 V-points
5 93.50 112.50 1.00 35.50 V-points
6 112.50 131.50 1.00 35.50 V-points
7 131.50 148.50 1.00 35.50 V-points
8 0.50 17.50 35.50 73.50 V-points
9 17.50 36.50 35.50 73.50 V-points
10 36.50 55.50 35.50 73.50 V-points
11 55.50 74.50 35.50 73.50 V-points
12 74.50 93.50 35.50 73.50 V-points
13 93.50 112.50 35.50 73.50 V-points
14 112.50 131.50 35.50 73.50 V-points
15 131.50 148.50 35.50 73.50 V-points
16 0.50 17.50 73.50 111.50 V-points
17 17.50 36.50 73.50 111.50 V-points
18 36.50 55.50 73.50 111.50 V-points
19 55.50 74.50 73.50 111.50 V-points
20 74.50 93.50 73.50 111.50 V-points
21 93.50 112.50 73.50 111.50 V-points
22 112.50 131.50 73.50 111.50 V-points
23 131.50 148.50 73.50 111.50 V-points
24 0.50 17.50 111.50 149.50 V-points
25 17.50 36.50 111.50 149.50 V-points
26 36.50 55.50 111.50 149.50 V-points
27 55.50 74.50 111.50 149.50 V-points
28 74.50 93.50 111.50 149.50 V-points
29 93.50 112.50 111.50 149.50 V-points
30 112.50 131.50 111.50 149.50 V-points
31 131.50 148.50 111.50 149.50 V-points
32 0.50 17.50 149.50 187.50 V-points
33 17.50 36.50 149.50 187.50 V-points
34 36.50 55.50 149.50 187.50 V-points
35 55.50 74.50 149.50 187.50 V-points
36 74.50 93.50 149.50 187.50 V-points
37 93.50 112.50 149.50 187.50 V-points
38 112.50 131.50 149.50 187.50 V-points
39 131.50 148.50 149.50 187.50 V-points
40 0.50 17.50 187.50 225.50 V-points
41 17.50 36.50 187.50 225.50 V-points
42 36.50 55.50 187.50 225.50 V-points
43 55.50 74.50 187.50 225.50 V-points
44 74.50 93.50 187.50 225.50 V-points
45 93.50 112.50 187.50 225.50 V-points
46 112.50 131.50 187.50 225.50 V-points
47 131.50 148.50 187.50 225.50 V-points
48 0.50 17.50 225.50 263.50 V-points
49 17.50 36.50 225.50 263.50 V-points
50 36.50 55.50 225.50 263.50 V-points
51 55.50 74.50 225.50 263.50 V-points
52 74.50 93.50 225.50 263.50 V-points
53 93.50 112.50 225.50 263.50 V-points
54 112.50 131.50 225.50 263.50 V-points
55 131.50 148.50 225.50 263.50 V-points
56 0.50 17.50 263.50 298.00 V-points
57 17.50 36.50 263.50 298.00 V-points
58 36.50 55.50 263.50 298.00 V-points
59 55.50 74.50 263.50 298.00 V-points
60 74.50 93.50 263.50 298.00 V-points
61 93.50 112.50 263.50 298.00 V-points
62 112.50 131.50 263.50 298.00 V-points
63 131.50 148.50 263.50 298.00 V-points
Maximum halo size in XI and ETA directions:
HaloSizeI(1) = 62
HaloSizeJ(1) = 100
TileSide(1) = 44
TileSize(1) = 1100
Activated C-preprocessing Options:
DELBAYCC2012 delbayCC2012
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_SMFLUX Analytical kinematic surface momentum flux.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ANA_STFLUX Analytical kinematic surface temperature flux.
ASSUMED_SHAPE Using assumed-shape arrays.
CURVGRID Orthogonal curvilinear grid.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
EAST_FSCHAPMAN Eastern edge, free-surface, Chapman condition.
EAST_M2FLATHER Eastern edge, 2D momentum, Flather condition.
EAST_M3RADIATION Eastern edge, 3D momentum, radiation condition.
EAST_TRADIATION Eastern edge, tracers, radiation condition.
FSOBC_REDUCED Using free-surface data in reduced physics condtions
INLINE_2DIO Processing 3D IO level by level to reduce memory needs.
MASKING Land/Sea masking.
MCT_LIB Using Model Coupling Toolkit library.
MIX_GEO_TS Mixing of tracers along geopotential surfaces.
MIX_GEO_UV Mixing of momentum along geopotential surfaces.
MPI MPI distributed-memory configuration.
MY25_MIXING Mellor/Yamada Level-2.5 mixing closure.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
NORTH_FSCHAPMAN Northern edge, free-surface, Chapman condition.
NORTH_M2FLATHER Northern edge, 2D momentum, Flather condition.
NORTH_M3RADIATION Northern edge, 3D momentum, radiation condition.
NORTH_TRADIATION Northern edge, tracers, radiation condition.
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.
SALINITY Using salinity.
SOLVE3D Solving 3D Primitive Equations.
SOUTH_FSCHAPMAN Southern edge, free-surface, Chapman condition.
SOUTH_M2FLATHER Southern edge, 2D momentum, Flather condition.
SOUTH_M3RADIATION Southern edge, 3D momentum, radiation condition.
SOUTH_TRADIATION Southern edge, tracers, radiation condition.
SPLINES Conservative parabolic spline reconstruction.
SWAN_COUPLING SWAN model coupling.
TS_A4HADVECTION Fourth-order Akima horizontal advection of tracers.
TS_A4VADVECTION Fourth-order Akima vertical advection of tracers.
TS_DIF2 Harmonic mixing of tracers.
UV_ADV Advection of momentum.
UV_COR Coriolis term.
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum.
UV_C4VADVECTION Fourth-order centered vertical advection of momentum.
UV_QDRAG Quadratic bottom stress.
UV_TIDES Add tidal currents to 2D momentum climatologies.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
WAVES_OCEAN Two-way wave-ocean models coupling.
WET_DRY Wetting and drying activated.
INITIAL: Configuring and initializing forward nonlinear model ...
Vertical S-coordinate System:
level S-coord Cs-curve Z at hmin at hc half way at hmax
10 0.0000000 0.0000000 0.000 0.000 0.000 0.000
9 -0.1000000 -0.0060605 -0.018 0.000 -0.312 -0.606
8 -0.2000000 -0.0247149 -0.074 0.000 -1.273 -2.471
7 -0.3000000 -0.0573991 -0.172 0.000 -2.956 -5.740
6 -0.4000000 -0.1065532 -0.320 0.000 -5.487 -10.655
5 -0.5000000 -0.1756670 -0.527 0.000 -9.047 -17.567
4 -0.6000000 -0.2692794 -0.808 0.000 -13.868 -26.928
3 -0.7000000 -0.3928500 -1.179 0.000 -20.232 -39.285
2 -0.8000000 -0.5523651 -1.657 0.000 -28.447 -55.237
1 -0.9000000 -0.7534850 -2.260 0.000 -38.804 -75.349
0 -1.0000000 -1.0000000 -3.000 0.000 -51.500 -100.000
Time Splitting Weights: ndtfast = 30 nfast = 42
Primary Secondary Accumulated to Current Step
1-0.0008094437383769 0.0333333333333333-0.0008094437383769 0.0333333333333333
2-0.0014053566728197 0.0333603147912792-0.0022148004111966 0.0666936481246126
3-0.0017877524645903 0.0334071600137066-0.0040025528757869 0.1001008081383191
4-0.0019566842408176 0.0334667517625262-0.0059592371166046 0.1335675599008453
5-0.0019122901320372 0.0335319745705535-0.0078715272486418 0.1670995344713988
6-0.0016548570247459 0.0335957175749547-0.0095263842733877 0.2006952520463536
7-0.0011849025289723 0.0336508794757796-0.0107112868023601 0.2343461315221331
8-0.0005032751608632 0.0336903762267453-0.0112145619632232 0.2680365077488784
9 0.0003887272597151 0.0337071520654408-0.0108258347035082 0.3017436598143192
10 0.0014892209965583 0.0336941944901169-0.0093366137069498 0.3354378543044362
11 0.0027955815694920 0.0336445537902317-0.0065410321374578 0.3690824080946679
12 0.0043042707117221 0.0335513677379153-0.0022367614257357 0.4026337758325831
13 0.0060106451121704 0.0334078920475245 0.0037738836864347 0.4360416678801077
14 0.0079087469427945 0.0332075372104522 0.0116826306292293 0.4692492050905598
15 0.0099910761708919 0.0329439123123590 0.0216737068001212 0.5021931174029188
16 0.0122483446563884 0.0326108764399960 0.0339220514565096 0.5348039938429148
17 0.0146692120341107 0.0322025982847830 0.0485912634906203 0.5670065921276978
18 0.0172400033810439 0.0317136245503127 0.0658312668716642 0.5987202166780105
19 0.0199444086685725 0.0311389577709445 0.0857756755402367 0.6298591744489550
20 0.0227631639997064 0.0304741441486588 0.1085388395399431 0.6603333185976138
21 0.0256737146312911 0.0297153720153352 0.1342125541712341 0.6900486906129490
22 0.0286498597812016 0.0288595815276255 0.1628624139524358 0.7189082721405746
23 0.0316613792205220 0.0279045862015855 0.1945237931729577 0.7468128583421600
24 0.0346736416507075 0.0268492068942347 0.2291974348236652 0.7736620652363948
25 0.0376471948657328 0.0256934188392112 0.2668446296893980 0.7993554840756060
26 0.0405373376992232 0.0244385123436867 0.3073819673886213 0.8237939964192927
27 0.0432936737565710 0.0230872677537126 0.3506756411451923 0.8468812641730054
28 0.0458596469320356 0.0216441452951603 0.3965352880772279 0.8685254094681656
29 0.0481720587108285 0.0201154903974257 0.4447073467880564 0.8886408998655914
30 0.0501605672561820 0.0185097551070648 0.4948679140442384 0.9071506549726561
31 0.0517471682814030 0.0168377361985254 0.5466150823256414 0.9239883911711815
32 0.0528456577069106 0.0151128305891453 0.5994607400325520 0.9391012217603267
33 0.0533610761022577 0.0133513086655816 0.6528218161348097 0.9524525304259084
34 0.0531891349131380 0.0115726061288397 0.7060109510479476 0.9640251365547481
35 0.0522156244733761 0.0097996349650684 0.7582265755213237 0.9738247715198165
36 0.0503158038019030 0.0080591141492892 0.8085423793232267 0.9818838856691057
37 0.0473537721847153 0.0063819206892258 0.8558961515079421 0.9882658063583315
38 0.0431818225418188 0.0048034616164019 0.8990779740497609 0.9930692679747334
39 0.0376397765791564 0.0033640675316746 0.9367177506289173 0.9964333355064080
40 0.0305543017255206 0.0021094083123694 0.9672720523544379 0.9985427438187774
41 0.0217382098544504 0.0010909315881854 0.9890102622088883 0.9996336754069628
42 0.0109897377911118 0.0003663245930371 1.0000000000000002 0.9999999999999999
ndtfast, nfast = 30 42 nfast/ndtfast = 1.40000
Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):
1.000000000000 1.047601458608 0.523800729304 1.000000000000 1.000000000000
Power filter parameters, Fgamma, gamma = 0.28400 0.18933
Minimum X-grid spacing, DXmin = 6.83656115E-01 km
Maximum X-grid spacing, DXmax = 7.98995075E+00 km
Minimum Y-grid spacing, DYmin = 5.84204633E-01 km
Maximum Y-grid spacing, DYmax = 4.61661121E+00 km
Minimum Z-grid spacing, DZmin = 1.81813631E-02 m
Maximum Z-grid spacing, DZmax = 2.46514996E+01 m
Minimum barotropic Courant Number = 3.10445284E-02
Maximum barotropic Courant Number = 2.70482781E-01
Maximum Coriolis Courant Number = 1.39805988E-02
NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(File: delbayini2012CCNAM.nc, Rec=0001, Index=1)
- free-surface
(Min = -2.00000000E+00 Max = 1.17945957E+00)
- vertically integrated u-momentum component
(Min = -6.78342130E-01 Max = 6.58236907E-01)
- vertically integrated v-momentum component
(Min = -5.32789570E-01 Max = 1.49376311E-01)
- u-momentum component
(Min = -8.47913088E-01 Max = 1.19578593E+00)
- v-momentum component
(Min = -9.65962735E-01 Max = 2.50923286E-01)
- potential temperature
(Min = -8.73054692E-02 Max = 9.84094524E+00)
- salinity
(Min = -3.05489670E-01 Max = 3.44343500E+01)
GET_NGFLD - tidal period
(Min = 1.49040000E+04 Max = 9.67320000E+04)
GET_2DFLD - tidal current inclination angle
(Min = 0.00000000E+00 Max = 6.28318531E+00)
GET_2DFLD - tidal current phase angle
(Min = 0.00000000E+00 Max = 6.28317864E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 0.00000000E+00 Max = 1.91539911E+00)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -3.28818765E-01 Max = 2.26390060E-01)
GET_NGFLD - free-surface eastern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -6.24522797E-01 Max = -5.56181716E-01)
GET_NGFLD - free-surface southern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -6.64299099E-01 Max = -6.13464345E-01)
GET_NGFLD - free-surface northern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -5.67338138E-01 Max = -5.10653361E-01)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -1.78650608E-01 Max = 2.38413242E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -4.05335075E-01 Max = 2.60545953E-01)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -1.56434706E-01 Max = 3.81790173E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -2.16383088E-01 Max = 1.65037623E-01)
GET_NGFLD - 2D u-momentum northern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -2.25385881E-02 Max = 1.06998015E-01)
GET_NGFLD - 2D v-momentum northern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -2.31242745E-01 Max = -2.67002905E-02)
Maximum grid stiffness ratios: rx0 = 5.677282E-01 (Beckmann and Haidvogel)
rx1 = 4.038306E+00 (Haney)
Initial basin volumes: TotVolume = 1.5638747369E+12 m3
MinVolume = 9.6644887677E+03 m3
MaxVolume = 8.6750340884E+08 m3
Max/Min = 8.9761955308E+04
== SWAN grid 1 sent wave data to ROMS
** ROMS grid 1 recv data from SWAN
** ROMS grid 1 sent data to SWAN
NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00192961)
GET_NGFLD - free-surface eastern boundary condition, t = 0 03:00:00
(Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -6.27662354E-01 Max = -5.51600357E-01)
GET_NGFLD - free-surface southern boundary condition, t = 0 03:00:00
(Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -6.69472724E-01 Max = -6.10971062E-01)
GET_NGFLD - free-surface northern boundary condition, t = 0 03:00:00
(Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -5.59730947E-01 Max = -5.06058428E-01)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 0 03:00:00
(Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -1.90737257E-01 Max = 1.25957661E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 0 03:00:00
(Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -3.49767328E-01 Max = 1.98288440E-01)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 0 03:00:00
(Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -2.00448120E-01 Max = -1.47098131E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 0 03:00:00
(Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -1.93069435E-01 Max = 2.10398759E-01)
== SWAN grid 1 recvd data from ROMS
GET_NGFLD - 2D u-momentum northern boundary condition, t = 0 03:00:00
(Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -5.45434644E-02 Max = 7.34075623E-02)
GET_NGFLD - 2D v-momentum northern boundary condition, t = 0 03:00:00
(Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
(Tmin= 0.0000 Tmax= 365.8750)
(Min = -1.79090935E-01 Max = -1.54667225E-02)
+time 20120101.001500 , step 1; iteration 1; sweep 1 grid 1
GRID STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME
1 0 0 00:00:00 4.965835E-03 2.557462E+02 2.557512E+02 1.645769E+12
DEF_HIS - creating history file: Output/ocean_his_0001.nc
+time 20120101.001500 , step 1; iteration 1; sweep 2 grid 1
+time 20120101.001500 , step 1; iteration 1; sweep 3 grid 1
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
+time 20120101.001500 , step 1; iteration 1; sweep 4 grid 1
1 1 0 00:02:30 NaN NaN NaN NaN
ROMS/TOMS - Blows up ................ exit_flag: 1
Blowing-up: Saving latest model state into RESTART file
+time 20120101.003000 , step 2; iteration 1; sweep 1 grid 1
+time 20120101.003000 , step 2; iteration 1; sweep 2 grid 1
Node # 23 CPU: 5.568
Node # 1 CPU: 5.568
Node # 2 CPU: 5.568
Node # 12 CPU: 5.568
Node # 3 CPU: 5.568
Node # 5 CPU: 5.568
Node # 4 CPU: 5.568
Node # 6 CPU: 5.568
Node # 8 CPU: 5.568
Node # 7 CPU: 5.568
Node # 10 CPU: 5.568
Node # 9 CPU: 5.568
Node # 11 CPU: 5.568
Node # 13 CPU: 5.568
Node # 14 CPU: 5.568
Node # 15 CPU: 5.568
Node # 16 CPU: 5.569
Node # 17 CPU: 5.568
Node # 18 CPU: 5.568
Node # 19 CPU: 5.568
Node # 20 CPU: 5.568
Node # 21 CPU: 5.568
Node # 22 CPU: 5.568
WRT_RST - wrote re-start fields (Index=1,2) into time record = 0000001
Elapsed CPU time (seconds):
Node # 0 CPU: 5.569
Node # 24 CPU: 5.568
Node # 25 CPU: 5.568
Node # 48 CPU: 5.569
Node # 26 CPU: 5.568
Node # 49 CPU: 5.569
Node # 27 CPU: 5.568
Node # 50 CPU: 5.569
Node # 28 CPU: 5.568
Node # 51 CPU: 5.569
Node # 29 CPU: 5.568
Node # 52 CPU: 5.569
Node # 30 CPU: 5.568
Node # 53 CPU: 5.569
Node # 31 CPU: 5.568
Node # 54 CPU: 5.569
Node # 32 CPU: 5.568
Node # 55 CPU: 5.569
Node # 33 CPU: 5.568
Node # 56 CPU: 5.569
Node # 34 CPU: 5.568
Node # 57 CPU: 5.569
Node # 35 CPU: 5.568
Node # 58 CPU: 5.569
Node # 36 CPU: 5.568
Node # 59 CPU: 5.569
Node # 37 CPU: 5.568
Node # 60 CPU: 5.569
Node # 38 CPU: 5.568
Node # 61 CPU: 5.569
Node # 39 CPU: 5.568
Node # 62 CPU: 5.569
Node # 40 CPU: 5.568
Node # 63 CPU: 5.569
Node # 41 CPU: 5.568
Node # 42 CPU: 5.569
Node # 43 CPU: 5.568
Node # 44 CPU: 5.568
Node # 45 CPU: 5.568
Node # 46 CPU: 5.568
Node # 47 CPU: 5.568
Total: 356.370
Nonlinear model elapsed time profile:
Initialization ................................... 163.414 (45.8552 %)
OI data assimilation ............................. 2.222 ( 0.6236 %)
Reading of input data ............................ 95.292 (26.7396 %)
Processing of input data ......................... 0.983 ( 0.2759 %)
Computation of vertical boundary conditions ...... 0.215 ( 0.0604 %)
Computation of global information integrals ...... 0.760 ( 0.2132 %)
Writing of output data ........................... 53.810 (15.0994 %)
Model 2D kernel .................................. 13.850 ( 3.8863 %)
Tidal forcing .................................... 0.277 ( 0.0776 %)
2D/3D coupling, vertical metrics ................. 0.129 ( 0.0363 %)
Omega vertical velocity .......................... 0.338 ( 0.0949 %)
Equation of state for seawater ................... 0.403 ( 0.1131 %)
My2.5 vertical mixing parameterization ........... 1.133 ( 0.3180 %)
3D equations right-side terms .................... 0.320 ( 0.0897 %)
3D equations predictor step ...................... 2.074 ( 0.5819 %)
Pressure gradient ................................ 1.360 ( 0.3818 %)
Harmonic mixing of tracers, geopotentials ........ 3.138 ( 0.8805 %)
Harmonic stress tensor, geopotentials ............ 3.393 ( 0.9520 %)
Corrector time-step for 3D momentum .............. 0.781 ( 0.2192 %)
Corrector time-step for tracers .................. 0.424 ( 0.1189 %)
Total: 344.316 96.6174
Nonlinear model message Passage profile:
Message Passage: 2D halo exchanges ............... 8.907 ( 2.4995 %)
Message Passage: 3D halo exchanges ............... 1.453 ( 0.4078 %)
Message Passage: 4D halo exchanges ............... 0.513 ( 0.1439 %)
Message Passage: data broadcast .................. 156.046 (43.7875 %)
Message Passage: data reduction .................. 0.148 ( 0.0416 %)
Message Passage: data gathering .................. 18.610 ( 5.2220 %)
Message Passage: data scattering.................. 18.562 ( 5.2085 %)
Message Passage: boundary data gathering ......... 0.239 ( 0.0671 %)
Message Passage: multi-model coupling ............ 22.428 ( 6.2935 %)
Total: 226.906 63.6714
All percentages are with respect to total time = 356.370
ROMS/TOMS - Output NetCDF summary for Grid 01:
number of time records written in HISTORY file = 00000001
number of time records written in RESTART file = 00000001
Analytical header files used:
ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_smflux.h
ROMS/Functionals/ana_stflux.h
ROMS/TOMS - Blows up ................ exit_flag: 1
MAIN: Abnormal termination: BLOWUP.
+time 20120101.003000 , step 2; iteration 1; sweep 3 grid 1
+time 20120101.003000 , step 2; iteration 1; sweep 4 grid 1
== SWAN grid 1 sent wave data to ROMS
