Initial and boundary file is HYCOM + NCODA Global 1/12° Analysis (GLBa0.08/expt_91.1) data. Initial field of surface temperature seemed to similar with satellite observation field, but simulated surface temperature was increased. As a result, simulated surface temperature was about 5 celsius degree higher compared to observed.
Thanks to your tips.
Fig 1. Initial surface temperature field, simulated.
Fig 2. Same to Fig 1. except for satellite observation.
Fig 3. Surface temperature time plot.
Here is my cpp option.
-------------------------------------------------------
#define ROMS_MODEL
#undef NESTING
#ifdef ROMS_MODEL
/* Physics + numerics */
# define UV_ADV
# define UV_COR
# define UV_VIS2
# define MIX_S_UV
# undef TS_FIXED
# define TS_U3HADVECTION
# define TS_C4VADVECTION
# undef TS_MPDATA
# undef SSW_BBL
# ifdef SSW_BBL
# define SSW_CALC_ZNOT
# /*define ANA_SEDIMENT*/
# else
# define UV_LOGDRAG
# endif
# define DJ_GRADPS
# define TS_DIF2
# define MIX_GEO_TS
# define CURVGRID
# define SALINITY
# define SOLVE3D
# define SPLINES_VDIFF
# define SPLINES_VVISC
# define AVERAGES
# define NONLIN_EOS
/* Grid and Initial */
# define MASKING
/* Forcing */
# ifdef WRF_MODEL
# undef BULK_FLUXES
# define ATM2OCN_FLUXES
# define ANA_SSFLUX
# undef LONGWAVE_OUT
# define ATM_PRESS
# else
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_SSFLUX
# endif
# define ANA_BTFLUX
# define ANA_BSFLUX
# define EMINUSP
# undef SOLAR_SOURCE
/* Turbulence closure */
# define GLS_MIXING
# undef MY25_MIXING
# define AKLIMIT
# if defined GLS_MIXING || defined MY25_MIXING
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# endif
/* Output */
# undef DIAGNOSTICS_UV
# undef DIAGNOSTICS_TS
#endif
-----------------------------------------------------------
Here is log
-----------------------------------------------------------
Model Input Parameters: ROMS/TOMS version 3.7
Thursday - February 9, 2017 - 10:22:28 PM
-----------------------------------------------------------------------------
Yo-Hwan Fog case
Operating system : Linux
CPU/hardware : x86_64
Compiler system : pgi
Compiler command : /usr/local/mpi/pgi/pgi-133/openmpi-1.8.3/bin/mpif90
Compiler flags : -g -C -Mfree
Input Script :
SVN Root URL : https:://myroms.org/svn/src
SVN Revision : 1058M
Local Root : /home/PSuniv/COAWST
Header Dir : /home/PSuniv/COAWST/ROMS/Include
Header file : one.h
Analytical Dir: /home/PSuniv/COAWST/ROMS/Functionals
Resolution, Grid 01: 0286x0304x016, Parallel Nodes: 64, Tiling: 008x008
Physical Parameters, Grid: 01
=============================
20160 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.
120 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
120 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
0.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
0.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
1.0000E-01 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
F LuvSponge Turning OFF sponge on horizontal momentum.
F LtracerSponge(01) Turning OFF sponge on tracer 01: temp
F LtracerSponge(02) Turning OFF sponge on tracer 02: salt
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.000 gls_p GLS stability exponent.
1.500 gls_m GLS turbulent kinetic energy exponent.
-1.000 gls_n GLS turbulent length scale exponent.
7.6000E-06 gls_Kmin GLS minimum value of turbulent kinetic energy.
1.0000E-12 gls_Pmin GLS minimum value of dissipation.
5.4770E-01 gls_cmu0 GLS stability coefficient.
1.4400E+00 gls_c1 GLS shear production coefficient.
1.9200E+00 gls_c2 GLS dissipation coefficient.
-4.0000E-01 gls_c3m GLS stable buoyancy production coefficient.
1.0000E+00 gls_c3p GLS unstable buoyancy production coefficient.
1.0000E+00 gls_sigk GLS constant Schmidt number for TKE.
1.3000E+00 gls_sigp GLS constant Schmidt number for PSI.
1400.000 charnok_alpha Charnok factor for Zos calculation.
0.500 zos_hsig_alpha Factor for Zos calculation using Hsig(Awave).
0.250 sz_alpha Factor for Wave dissipation surface tke flux .
100.000 crgban_cw Factor for Craig/Banner surface tke flux.
0.000 wec_alpha WEC factor for roller/breaking energy distribution.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
2.0000E-02 Zos Surface roughness (m).
1 Vtransform S-coordinate transformation equation.
1 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.
5.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
57061.000 dstart Time-stamp assigned to model initialization (days).
18581117.00 time_ref Reference time for units attribute (yyyymmdd.dd)
1.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
1.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.
1.0000E+00 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
0.0000E+00 obcfac Factor between passive and active
open boundary conditions.
F VolCons(1) NLM western edge boundary volume conservation.
F VolCons(2) NLM southern edge boundary volume conservation.
F VolCons(3) NLM eastern edge boundary volume conservation.
F VolCons(4) NLM northern edge boundary volume conservation.
10.000 T0 Background potential temperature (C) constant.
30.000 S0 Background salinity (PSU) constant.
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
F LuvSrc Turning OFF momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
F LtracerSrc(01) Turning OFF point Sources/Sinks on tracer 01: temp
F LtracerSrc(02) Turning OFF point Sources/Sinks on tracer 02: salt
F LsshCLM Turning OFF processing of SSH climatology.
F Lm2CLM Turning OFF processing of 2D momentum climatology.
F Lm3CLM Turning OFF processing of 3D momentum climatology.
F LtracerCLM(01) Turning OFF processing of climatology tracer 01: temp
F LtracerCLM(02) Turning OFF processing of climatology tracer 02: salt
F LnudgeM2CLM Turning OFF nudging of 2D momentum climatology.
F LnudgeM3CLM Turning OFF nudging of 3D momentum climatology.
F LnudgeTCLM(01) Turning OFF nudging of climatology tracer 01: temp
F LnudgeTCLM(02) Turning OFF nudging of climatology tracer 02: salt
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(idVvis) Write out vertical viscosity: AKv.
T Hout(idMtke) Write out turbulent kinetic energy.
T Hout(idMtls) Write out turbulent generic length-scale.
T Aout(idSdif) Write out averaged vertical diffusion: AKt(isalt).
Output/Input Files:
Output Restart File: YH_roms_rst.nc
Output History File: YH_roms_his.nc
Output Averages File: YH_roms_avg.nc
Input Grid File: Projects/tidal_mixing/YH_test_grid.nc
Input Nonlinear Initial File: Projects/tidal_mixing/20150208-20150214/coawst_ini.nc
Input Boundary File 01: Projects/tidal_mixing/20150208-20150214/merged_coawst_bdy.nc
Tile partition information for Grid 01: 0286x0304x0016 tiling: 008x008
tile Istr Iend Jstr Jend Npts
Number of tracers: 2
0 1 35 1 38 21280
1 36 71 1 38 21888
2 72 107 1 38 21888
3 108 143 1 38 21888
4 144 179 1 38 21888
5 180 215 1 38 21888
6 216 251 1 38 21888
7 252 286 1 38 21280
8 1 35 39 76 21280
9 36 71 39 76 21888
10 72 107 39 76 21888
11 108 143 39 76 21888
12 144 179 39 76 21888
13 180 215 39 76 21888
14 216 251 39 76 21888
15 252 286 39 76 21280
16 1 35 77 114 21280
17 36 71 77 114 21888
18 72 107 77 114 21888
19 108 143 77 114 21888
20 144 179 77 114 21888
21 180 215 77 114 21888
22 216 251 77 114 21888
23 252 286 77 114 21280
24 1 35 115 152 21280
25 36 71 115 152 21888
26 72 107 115 152 21888
27 108 143 115 152 21888
28 144 179 115 152 21888
29 180 215 115 152 21888
30 216 251 115 152 21888
31 252 286 115 152 21280
32 1 35 153 190 21280
33 36 71 153 190 21888
34 72 107 153 190 21888
35 108 143 153 190 21888
36 144 179 153 190 21888
37 180 215 153 190 21888
38 216 251 153 190 21888
39 252 286 153 190 21280
40 1 35 191 228 21280
41 36 71 191 228 21888
42 72 107 191 228 21888
43 108 143 191 228 21888
44 144 179 191 228 21888
45 180 215 191 228 21888
46 216 251 191 228 21888
47 252 286 191 228 21280
48 1 35 229 266 21280
49 36 71 229 266 21888
50 72 107 229 266 21888
51 108 143 229 266 21888
52 144 179 229 266 21888
53 180 215 229 266 21888
54 216 251 229 266 21888
55 252 286 229 266 21280
56 1 35 267 304 21280
57 36 71 267 304 21888
58 72 107 267 304 21888
59 108 143 267 304 21888
60 144 179 267 304 21888
61 180 215 267 304 21888
62 216 251 267 304 21888
63 252 286 267 304 21280
Tile minimum and maximum fractional coordinates for Grid 01:
(interior points only)
tile Xmin Xmax Ymin Ymax grid
0 0.50 35.50 0.50 38.50 RHO-points
1 35.50 71.50 0.50 38.50 RHO-points
2 71.50 107.50 0.50 38.50 RHO-points
3 107.50 143.50 0.50 38.50 RHO-points
4 143.50 179.50 0.50 38.50 RHO-points
5 179.50 215.50 0.50 38.50 RHO-points
6 215.50 251.50 0.50 38.50 RHO-points
7 251.50 286.50 0.50 38.50 RHO-points
8 0.50 35.50 38.50 76.50 RHO-points
9 35.50 71.50 38.50 76.50 RHO-points
10 71.50 107.50 38.50 76.50 RHO-points
11 107.50 143.50 38.50 76.50 RHO-points
12 143.50 179.50 38.50 76.50 RHO-points
13 179.50 215.50 38.50 76.50 RHO-points
14 215.50 251.50 38.50 76.50 RHO-points
15 251.50 286.50 38.50 76.50 RHO-points
16 0.50 35.50 76.50 114.50 RHO-points
17 35.50 71.50 76.50 114.50 RHO-points
18 71.50 107.50 76.50 114.50 RHO-points
19 107.50 143.50 76.50 114.50 RHO-points
20 143.50 179.50 76.50 114.50 RHO-points
21 179.50 215.50 76.50 114.50 RHO-points
22 215.50 251.50 76.50 114.50 RHO-points
23 251.50 286.50 76.50 114.50 RHO-points
24 0.50 35.50 114.50 152.50 RHO-points
25 35.50 71.50 114.50 152.50 RHO-points
26 71.50 107.50 114.50 152.50 RHO-points
27 107.50 143.50 114.50 152.50 RHO-points
28 143.50 179.50 114.50 152.50 RHO-points
29 179.50 215.50 114.50 152.50 RHO-points
30 215.50 251.50 114.50 152.50 RHO-points
31 251.50 286.50 114.50 152.50 RHO-points
32 0.50 35.50 152.50 190.50 RHO-points
33 35.50 71.50 152.50 190.50 RHO-points
34 71.50 107.50 152.50 190.50 RHO-points
35 107.50 143.50 152.50 190.50 RHO-points
36 143.50 179.50 152.50 190.50 RHO-points
37 179.50 215.50 152.50 190.50 RHO-points
38 215.50 251.50 152.50 190.50 RHO-points
39 251.50 286.50 152.50 190.50 RHO-points
40 0.50 35.50 190.50 228.50 RHO-points
41 35.50 71.50 190.50 228.50 RHO-points
42 71.50 107.50 190.50 228.50 RHO-points
43 107.50 143.50 190.50 228.50 RHO-points
44 143.50 179.50 190.50 228.50 RHO-points
45 179.50 215.50 190.50 228.50 RHO-points
46 215.50 251.50 190.50 228.50 RHO-points
47 251.50 286.50 190.50 228.50 RHO-points
48 0.50 35.50 228.50 266.50 RHO-points
49 35.50 71.50 228.50 266.50 RHO-points
50 71.50 107.50 228.50 266.50 RHO-points
51 107.50 143.50 228.50 266.50 RHO-points
52 143.50 179.50 228.50 266.50 RHO-points
53 179.50 215.50 228.50 266.50 RHO-points
54 215.50 251.50 228.50 266.50 RHO-points
55 251.50 286.50 228.50 266.50 RHO-points
56 0.50 35.50 266.50 304.50 RHO-points
57 35.50 71.50 266.50 304.50 RHO-points
58 71.50 107.50 266.50 304.50 RHO-points
59 107.50 143.50 266.50 304.50 RHO-points
60 143.50 179.50 266.50 304.50 RHO-points
61 179.50 215.50 266.50 304.50 RHO-points
62 215.50 251.50 266.50 304.50 RHO-points
63 251.50 286.50 266.50 304.50 RHO-points
0 1.00 35.50 0.50 38.50 U-points
1 35.50 71.50 0.50 38.50 U-points
2 71.50 107.50 0.50 38.50 U-points
3 107.50 143.50 0.50 38.50 U-points
4 143.50 179.50 0.50 38.50 U-points
5 179.50 215.50 0.50 38.50 U-points
6 215.50 251.50 0.50 38.50 U-points
7 251.50 286.00 0.50 38.50 U-points
8 1.00 35.50 38.50 76.50 U-points
9 35.50 71.50 38.50 76.50 U-points
10 71.50 107.50 38.50 76.50 U-points
11 107.50 143.50 38.50 76.50 U-points
12 143.50 179.50 38.50 76.50 U-points
13 179.50 215.50 38.50 76.50 U-points
14 215.50 251.50 38.50 76.50 U-points
15 251.50 286.00 38.50 76.50 U-points
16 1.00 35.50 76.50 114.50 U-points
17 35.50 71.50 76.50 114.50 U-points
18 71.50 107.50 76.50 114.50 U-points
19 107.50 143.50 76.50 114.50 U-points
20 143.50 179.50 76.50 114.50 U-points
21 179.50 215.50 76.50 114.50 U-points
22 215.50 251.50 76.50 114.50 U-points
23 251.50 286.00 76.50 114.50 U-points
24 1.00 35.50 114.50 152.50 U-points
25 35.50 71.50 114.50 152.50 U-points
26 71.50 107.50 114.50 152.50 U-points
27 107.50 143.50 114.50 152.50 U-points
28 143.50 179.50 114.50 152.50 U-points
29 179.50 215.50 114.50 152.50 U-points
30 215.50 251.50 114.50 152.50 U-points
31 251.50 286.00 114.50 152.50 U-points
32 1.00 35.50 152.50 190.50 U-points
33 35.50 71.50 152.50 190.50 U-points
34 71.50 107.50 152.50 190.50 U-points
35 107.50 143.50 152.50 190.50 U-points
36 143.50 179.50 152.50 190.50 U-points
37 179.50 215.50 152.50 190.50 U-points
38 215.50 251.50 152.50 190.50 U-points
39 251.50 286.00 152.50 190.50 U-points
40 1.00 35.50 190.50 228.50 U-points
41 35.50 71.50 190.50 228.50 U-points
42 71.50 107.50 190.50 228.50 U-points
43 107.50 143.50 190.50 228.50 U-points
44 143.50 179.50 190.50 228.50 U-points
45 179.50 215.50 190.50 228.50 U-points
46 215.50 251.50 190.50 228.50 U-points
47 251.50 286.00 190.50 228.50 U-points
48 1.00 35.50 228.50 266.50 U-points
49 35.50 71.50 228.50 266.50 U-points
50 71.50 107.50 228.50 266.50 U-points
51 107.50 143.50 228.50 266.50 U-points
52 143.50 179.50 228.50 266.50 U-points
53 179.50 215.50 228.50 266.50 U-points
54 215.50 251.50 228.50 266.50 U-points
55 251.50 286.00 228.50 266.50 U-points
56 1.00 35.50 266.50 304.50 U-points
57 35.50 71.50 266.50 304.50 U-points
58 71.50 107.50 266.50 304.50 U-points
59 107.50 143.50 266.50 304.50 U-points
60 143.50 179.50 266.50 304.50 U-points
61 179.50 215.50 266.50 304.50 U-points
62 215.50 251.50 266.50 304.50 U-points
63 251.50 286.00 266.50 304.50 U-points
0 0.50 35.50 1.00 38.50 V-points
1 35.50 71.50 1.00 38.50 V-points
2 71.50 107.50 1.00 38.50 V-points
3 107.50 143.50 1.00 38.50 V-points
4 143.50 179.50 1.00 38.50 V-points
5 179.50 215.50 1.00 38.50 V-points
6 215.50 251.50 1.00 38.50 V-points
7 251.50 286.50 1.00 38.50 V-points
8 0.50 35.50 38.50 76.50 V-points
9 35.50 71.50 38.50 76.50 V-points
10 71.50 107.50 38.50 76.50 V-points
11 107.50 143.50 38.50 76.50 V-points
12 143.50 179.50 38.50 76.50 V-points
13 179.50 215.50 38.50 76.50 V-points
14 215.50 251.50 38.50 76.50 V-points
15 251.50 286.50 38.50 76.50 V-points
16 0.50 35.50 76.50 114.50 V-points
17 35.50 71.50 76.50 114.50 V-points
18 71.50 107.50 76.50 114.50 V-points
19 107.50 143.50 76.50 114.50 V-points
20 143.50 179.50 76.50 114.50 V-points
21 179.50 215.50 76.50 114.50 V-points
22 215.50 251.50 76.50 114.50 V-points
23 251.50 286.50 76.50 114.50 V-points
24 0.50 35.50 114.50 152.50 V-points
25 35.50 71.50 114.50 152.50 V-points
26 71.50 107.50 114.50 152.50 V-points
27 107.50 143.50 114.50 152.50 V-points
28 143.50 179.50 114.50 152.50 V-points
29 179.50 215.50 114.50 152.50 V-points
30 215.50 251.50 114.50 152.50 V-points
31 251.50 286.50 114.50 152.50 V-points
32 0.50 35.50 152.50 190.50 V-points
33 35.50 71.50 152.50 190.50 V-points
34 71.50 107.50 152.50 190.50 V-points
35 107.50 143.50 152.50 190.50 V-points
36 143.50 179.50 152.50 190.50 V-points
37 179.50 215.50 152.50 190.50 V-points
38 215.50 251.50 152.50 190.50 V-points
39 251.50 286.50 152.50 190.50 V-points
40 0.50 35.50 190.50 228.50 V-points
41 35.50 71.50 190.50 228.50 V-points
42 71.50 107.50 190.50 228.50 V-points
43 107.50 143.50 190.50 228.50 V-points
44 143.50 179.50 190.50 228.50 V-points
45 179.50 215.50 190.50 228.50 V-points
46 215.50 251.50 190.50 228.50 V-points
47 251.50 286.50 190.50 228.50 V-points
48 0.50 35.50 228.50 266.50 V-points
49 35.50 71.50 228.50 266.50 V-points
50 71.50 107.50 228.50 266.50 V-points
51 107.50 143.50 228.50 266.50 V-points
52 143.50 179.50 228.50 266.50 V-points
53 179.50 215.50 228.50 266.50 V-points
54 215.50 251.50 228.50 266.50 V-points
55 251.50 286.50 228.50 266.50 V-points
56 0.50 35.50 266.50 304.00 V-points
57 35.50 71.50 266.50 304.00 V-points
58 71.50 107.50 266.50 304.00 V-points
59 107.50 143.50 266.50 304.00 V-points
60 143.50 179.50 266.50 304.00 V-points
61 179.50 215.50 266.50 304.00 V-points
62 215.50 251.50 266.50 304.00 V-points
63 251.50 286.50 266.50 304.00 V-points
Maximum halo size in XI and ETA directions:
HaloSizeI(1) = 96
HaloSizeJ(1) = 100
TileSide(1) = 44
TileSize(1) = 1848
Lateral Boundary Conditions: NLM
============================
Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------
zeta 1 Closed Chapman Imp Chapman Imp Closed
ubar 1 Closed Flather Flather Closed
vbar 1 Closed Flather Flather Closed
u 1 Closed Rad + Nud Rad + Nud Closed
v 1 Closed Rad + Nud Rad + Nud Closed
temp 1 Closed Rad + Nud Rad + Nud Closed
salt 1 Closed Rad + Nud Rad + Nud Closed
tke 1 Closed Closed Closed Closed
Activated C-preprocessing Options:
ONE Yo-Hwan Fog case
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.
AVERAGES Writing out time-averaged nonlinear model fields.
CURVGRID Orthogonal curvilinear grid.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
GLS_MIXING Generic Length-Scale turbulence closure.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
MASKING Land/Sea masking.
MIX_GEO_TS Mixing of tracers along geopotential surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
MPI MPI distributed-memory configuration.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
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.
TS_U3HADVECTION Third-order upstream horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered 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_LOGDRAG Logarithmic bottom stress.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
Process Information:
Node # 0 (pid= 13453) is active.
Node # 32 (pid= 25569) is active.
Node # 12 (pid= 13475) is active.
Node # 44 (pid= 25591) is active.
Node # 5 (pid= 13461) is active.
Node # 15 (pid= 13481) is active.
Node # 3 (pid= 13458) is active.
Node # 6 (pid= 13463) is active.
Node # 13 (pid= 13477) is active.
Node # 37 (pid= 25577) is active.
Node # 38 (pid= 25579) is active.
Node # 45 (pid= 25593) is active.
Node # 47 (pid= 25597) is active.
Node # 35 (pid= 25572) is active.
Node # 11 (pid= 13473) is active.
Node # 2 (pid= 13455) is active.
Node # 34 (pid= 25571) is active.
Node # 43 (pid= 25589) is active.
Node # 16 (pid= 31885) is active.
Node # 48 (pid= 24553) is active.
Node # 21 (pid= 31893) is active.
Node # 28 (pid= 31907) is active.
Node # 53 (pid= 24561) is active.
Node # 60 (pid= 24575) is active.
Node # 31 (pid= 31913) is active.
Node # 22 (pid= 31895) is active.
Node # 29 (pid= 31909) is active.
Node # 54 (pid= 24563) is active.
Node # 63 (pid= 24581) is active.
Node # 61 (pid= 24577) is active.
Node # 19 (pid= 31888) is active.
Node # 27 (pid= 31905) is active.
Node # 59 (pid= 24573) is active.
Node # 51 (pid= 24556) is active.
Node # 18 (pid= 31887) is active.
Node # 50 (pid= 24555) is active.
Node # 8 (pid= 13467) is active.
Node # 40 (pid= 25583) is active.
Node # 4 (pid= 13459) is active.
Node # 14 (pid= 13479) is active.
Node # 36 (pid= 25575) is active.
Node # 10 (pid= 13471) is active.
Node # 46 (pid= 25595) is active.
Node # 9 (pid= 13469) is active.
Node # 41 (pid= 25585) is active.
Node # 42 (pid= 25587) is active.
Node # 7 (pid= 13465) is active.
Node # 1 (pid= 13454) is active.
Node # 39 (pid= 25581) is active.
Node # 33 (pid= 25570) is active.
INITIAL: Configuring and initializing forward nonlinear model ...
*******
Node # 24 (pid= 31899) is active.
Node # 20 (pid= 31891) is active.
Node # 56 (pid= 24567) is active.
Node # 52 (pid= 24559) is active.
Node # 30 (pid= 31911) is active.
Node # 62 (pid= 24579) is active.
Node # 25 (pid= 31901) is active.
Node # 26 (pid= 31903) is active.
Node # 57 (pid= 24569) is active.
Node # 58 (pid= 24571) is active.
Node # 23 (pid= 31897) is active.
Node # 55 (pid= 24565) is active.
Node # 17 (pid= 31886) is active.
Node # 49 (pid= 24554) is active.
Vertical S-coordinate System, Grid 01:
level S-coord Cs-curve Z at hmin at hc half way at hmax
16 0.0000000 0.0000000 0.000 0.000 0.000 0.000
15 -0.0625000 -0.0201876 -0.313 -0.313 -1.834 -3.356
14 -0.1250000 -0.0441682 -0.625 -0.625 -3.955 -7.285
13 -0.1875000 -0.0732955 -0.938 -0.938 -6.463 -11.989
12 -0.2500000 -0.1089095 -1.250 -1.250 -9.461 -17.672
11 -0.3125000 -0.1520926 -1.563 -1.563 -13.029 -24.495
10 -0.3750000 -0.2033374 -1.875 -1.875 -17.205 -32.534
9 -0.4375000 -0.2622537 -2.188 -2.188 -21.959 -41.730
8 -0.5000000 -0.3275288 -2.500 -2.500 -27.193 -51.885
7 -0.5625000 -0.3973005 -2.813 -2.813 -32.765 -62.718
6 -0.6250000 -0.4698651 -3.125 -3.125 -38.548 -73.972
5 -0.6875000 -0.5443911 -3.438 -3.438 -44.480 -85.522
4 -0.7500000 -0.6213093 -3.750 -3.750 -50.591 -97.432
3 -0.8125000 -0.7023017 -4.063 -4.063 -57.010 -109.957
2 -0.8750000 -0.7900507 -4.375 -4.375 -63.937 -123.500
1 -0.9375000 -0.8879634 -4.688 -4.688 -71.632 -138.576
0 -1.0000000 -1.0000000 -5.000 -5.000 -80.391 -155.781
Time Splitting Weights for Grid 01: 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.1501639558342179
4-0.0003735756740661 0.0501732728306823-0.0038390322877122 0.2003372286649002
5 0.0009829200513762 0.0501919516143856-0.0028561122363360 0.2505291802792858
6 0.0029141799764308 0.0501428056118168 0.0000580677400948 0.3006719858911026
7 0.0054132615310267 0.0499970966129953 0.0054713292711215 0.3506690825040979
8 0.0084687837865132 0.0497264335364439 0.0139401130576347 0.4003955160405418
9 0.0120633394191050 0.0493029943471183 0.0260034524767397 0.4496985103876601
10 0.0161716623600090 0.0486998273761630 0.0421751148367487 0.4983983377638231
11 0.0207585511322367 0.0478912442581626 0.0629336659689854 0.5462895820219856
12 0.0257765478740990 0.0468533167015507 0.0887102138430845 0.5931428987235363
13 0.0311633730493854 0.0455644893078458 0.1198735868924699 0.6387073880313821
14 0.0368391158442262 0.0440063206553765 0.1567127027366961 0.6827137086867585
15 0.0427031802506397 0.0421643648631652 0.1994158829873358 0.7248780735499237
16 0.0486309868367616 0.0400292058506332 0.2480468698240974 0.7649072794005569
17 0.0544704302037592 0.0375976565087951 0.3025173000278566 0.8025049359093520
18 0.0600380921294286 0.0348741349986072 0.3625553921572852 0.8373790709079592
19 0.0651152103984763 0.0318722303921357 0.4276706025557616 0.8692513013000949
20 0.0694434033194840 0.0286164698722119 0.4971140058752456 0.8978677711723068
21 0.0727201499285570 0.0251442997062377 0.5698341558038026 0.9230120708785445
22 0.0745940258796570 0.0215082922098099 0.6444281816834596 0.9445203630883544
23 0.0746596950216180 0.0177785909158270 0.7190878767050777 0.9622989540041814
24 0.0724526566618460 0.0140456061647461 0.7915405333669236 0.9763445601689276
25 0.0674437485167026 0.0104229733316538 0.8589842818836262 0.9867675335005814
26 0.0590334053485720 0.0070507859058187 0.9180176872321982 0.9938183194064001
27 0.0465456732896125 0.0040991156383901 0.9645633605218107 0.9979174350447901
28 0.0292219798521905 0.0017718319739095 0.9937853403740011 0.9996892670186996
29 0.0062146596259991 0.0003107329813000 1.0000000000000002 0.9999999999999996
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
Metrics information for Grid 01:
===============================
Minimum X-grid spacing, DXmin = 3.35152087E+00 km
Maximum X-grid spacing, DXmax = 3.44191975E+00 km
Minimum Y-grid spacing, DYmin = 3.35242046E+00 km
Maximum Y-grid spacing, DYmax = 3.44114726E+00 km
Minimum Z-grid spacing, DZmin = 3.12500000E-01 m
Maximum Z-grid spacing, DZmax = 1.72055392E+01 m
Minimum barotropic Courant Number = 4.31998051E-03
Maximum barotropic Courant Number = 2.45862776E-02
Maximum Coriolis Courant Number = 2.86662245E-03
NLM: GET_STATE - Read state initial conditions, t = 57061 00:00:00
(Grid 01, File: coawst_ini.nc, Rec=0001, Index=1)
- free-surface
(Min = -8.82477820E-01 Max = 8.41171682E-01)
- vertically integrated u-momentum component
(Min = -3.71012151E-01 Max = 7.83710659E-01)
- vertically integrated v-momentum component
(Min = -6.31825864E-01 Max = 3.07241231E-01)
- u-momentum component
(Min = -5.74320376E-01 Max = 9.57133770E-01)
- v-momentum component
(Min = -1.30848122E+00 Max = 3.31409514E-01)
- potential temperature
(Min = -9.13685441E-01 Max = 1.81215267E+01)
- salinity
(Min = 2.72288475E+01 Max = 3.80091476E+01)
CHECK_MULTIFILE - Error while checking input Lateral Boundary file:
Projects/tidal_mixing/20150208-20150214/merged_coawst_bdy.nc
last data time record available is for day: 57067 00:00:00
but data is needed to finish run until day: 57068 00:00:00
GET_NGFLD - free-surface eastern boundary condition, t = 57061 00:00:00
(Rec=0000001, Index=2, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -6.80020043E-02 Max = 2.25579692E-01)
GET_NGFLD - free-surface southern boundary condition, t = 57061 00:00:00
(Rec=0000001, Index=2, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -6.80019909E-02 Max = 8.32440454E-01)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 57061 00:00:00
(Rec=0000001, Index=2, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -2.70189606E-02 Max = 1.62398726E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 57061 00:00:00
(Rec=0000001, Index=2, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -4.56406921E-02 Max = 6.02269657E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 57061 00:00:00
(Rec=0000001, Index=2, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -6.13243841E-02 Max = 2.54670888E-01)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 57061 00:00:00
(Rec=0000001, Index=2, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -6.09676242E-01 Max = 2.83141583E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 57061 00:00:00
(Rec=0000001, Index=2, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -3.48781720E-02 Max = 3.63639772E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 57061 00:00:00
(Rec=0000001, Index=2, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -4.57415313E-01 Max = 1.42103016E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 57061 00:00:00
(Rec=0000001, Index=2, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -8.29695240E-02 Max = 2.83121318E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 57061 00:00:00
(Rec=0000001, Index=2, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -1.01872611E+00 Max = 8.68435875E-02)
GET_NGFLD - temperature eastern boundary condition, t = 57061 00:00:00
(Rec=0000001, Index=2, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = 6.45003271E+00 Max = 1.81215267E+01)
GET_NGFLD - salinity eastern boundary condition, t = 57061 00:00:00
(Rec=0000001, Index=2, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = 3.21437225E+01 Max = 3.78672295E+01)
GET_NGFLD - temperature southern boundary condition, t = 57061 00:00:00
(Rec=0000001, Index=2, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = 5.78649187E+00 Max = 1.78717747E+01)
GET_NGFLD - salinity southern boundary condition, t = 57061 00:00:00
(Rec=0000001, Index=2, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = 3.19170570E+01 Max = 3.79794960E+01)
Basin information for Grid 01:
Maximum grid stiffness ratios: rx0 = 2.689877E-01 (Beckmann and Haidvogel)
rx1 = 6.321894E+00 (Haney)
Initial basin volumes: TotVolume = 2.1387919345E+13 m3
MinVolume = 3.5224366674E+06 m3
MaxVolume = 1.8307785795E+08 m3
Max/Min = 5.1974776339E+01
NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00020160)
GET_NGFLD - free-surface eastern boundary condition, t = 57062 00:00:00
(Rec=0000002, Index=1, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -3.08899307E-02 Max = 2.93030547E-01)
GET_NGFLD - free-surface southern boundary condition, t = 57062 00:00:00
(Rec=0000002, Index=1, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -3.32379120E-02 Max = 2.93030547E-01)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 57062 00:00:00
(Rec=0000002, Index=1, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -5.01045547E-02 Max = 1.49097964E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 57062 00:00:00
(Rec=0000002, Index=1, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -2.60220729E-02 Max = 1.11977592E-01)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 57062 00:00:00
(Rec=0000002, Index=1, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -7.91303590E-02 Max = 7.60832429E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 57062 00:00:00
(Rec=0000002, Index=1, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -1.25897735E-01 Max = 1.52920187E-01)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 57062 00:00:00
(Rec=0000002, Index=1, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -7.39130154E-02 Max = 2.44209051E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 57062 00:00:00
(Rec=0000002, Index=1, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -3.06345820E-01 Max = 1.33361116E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 57062 00:00:00
(Rec=0000002, Index=1, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -1.22631989E-01 Max = 9.40430388E-02)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 57062 00:00:00
(Rec=0000002, Index=1, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = -2.63097614E-01 Max = 1.91671550E-01)
GET_NGFLD - temperature eastern boundary condition, t = 57062 00:00:00
(Rec=0000002, Index=1, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = 6.51260519E+00 Max = 1.84249763E+01)
GET_NGFLD - salinity eastern boundary condition, t = 57062 00:00:00
(Rec=0000002, Index=1, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = 3.21496544E+01 Max = 3.78586769E+01)
GET_NGFLD - temperature southern boundary condition, t = 57062 00:00:00
(Rec=0000002, Index=1, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = 5.63879538E+00 Max = 1.81941452E+01)
GET_NGFLD - salinity southern boundary condition, t = 57062 00:00:00
(Rec=0000002, Index=1, File: merged_coawst_bdy.nc)
(Tmin= 57061.0000 Tmax= 57067.0000)
(Min = 3.19435997E+01 Max = 3.79206848E+01)
