too high bottom temperature result in estuary

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
cheng_dada
Posts: 8
Joined: Fri Feb 26, 2021 2:20 pm
Location: Shantou University

too high bottom temperature result in estuary

#1 Unread post by cheng_dada »

Hi.
I encountered high and abnormal bottom temperature results compared to observation data.
The surface temperature of the whole model domain is same like MODIS and HYCOM. But in my study domain, the estuary, its temperature is too high, especially at the bottom, this result is abnormal.
I've tried several advection and mixing (MY-25 and LMD) schemes, but the results didn't get any better.
And both temperature and salinity are given to the river forcing file in my model.

Thanks for your tips.
Fig 1. vertical temperature field in estuary, simulated
Fig 1. vertical temperature field in estuary, simulated
Fig 2. Same to Fig 1. except for observation.
Fig 2. Same to Fig 1. except for observation.
-------------------------------------------------------------------------------------------------------------------
Here is my CPP option
/***************** Model dynamics */
#define SOLVE3D
#define UV_ADV
#define UV_COR
#define UV_QDRAG
#define DJ_GRADPS
#define UV_U3DVECTION
#define TS_C4HADVECTION
#define TS_U3VADVECTION
/************ Grid configuration */
#define MASKING
#define SPHERICAL
#define CURVGRID
/**** Mixing ********************/
#define MIX_GEO_TS
#define MIX_GEO_UV
#define DIFF_GRID
#define VISC_GRID
#define TS_SMAGORINSKY
#define UV_SMAGORINSKY
#define LMD_MIXING

# ifdef MY25_MIXING
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# endif
# ifdef LMD_MIXING
# define LMD_RIMIX
# define LMD_CONVEC
# define LMD_SKPP
# define LMD_BKPP
# define LMD_NONLOCAL
# define LMD_DDMIX
# endif
/**** Equation of State *********/
#define NONLIN_EOS
#define SPLINES
#define SALINITY
/*************** Surface Forcing */
# define SOLAR_SOURCE
# define QCORRECTION
# define DIURNAL_SRFLUX
# define DIURNAL_SRFLUX
/************** Other options ***/
#define RADIATION_2D
#define AVERAGES
#define ANA_BTFLUX
#define ANA_BSFLUX
#define LIMIT_BSTRESS
/******** Tide ************/
# define RAMP_TIDES
# define SSH_TIDES
# define UV_TIDES
# define ADD_M2OBC
# define ADD_FSOBC

#undef OUT_DOUBLE
#define RST_SINGLE
----------------------------------------------

Here is the log

Model Input Parameters: ROMS/TOMS version 3.7
Friday - June 10, 2022 - 10:45:09 AM
-----------------------------------------------------------------------------

Plume

Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /usr/local/bin/mpif90
Compiler flags : -frepack-arrays -O3 -ffast-math -ffree-form -ffree-line-length-none -ffree-form -ffree-line-length-none

Input Script :

SVN Root URL : https://www.myroms.org/svn/src/trunk
Local Root : /home/19xcheng/roms/roms
Header Dir : /home/19xcheng/roms/roms
Header file : stplume.h
Analytical Dir: /home/19xcheng/roms/roms/User/Functionals

Resolution, Grid 01: 0574x0376x030, Parallel Nodes: 72, Tiling: 009x008


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

212544 ntimes Number of timesteps for 3-D equations.
50.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.
1788 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).
72 nHIS Number of timesteps between the writing fields
into history file.
132840 ndefHIS Number of timesteps between creation of new
history files.
0 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
1788 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
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
5.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
5.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)
for tracer 02: salt
5.0000E-06 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
2.5000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
1 lmd_Jwt Jerlov water type.
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.
3.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
44315.000 dstart Time-stamp assigned to model initialization (days).
0.000 tide_start Reference time origin for tidal forcing (days).
19000101.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
3.0000E+01 Znudg Nudging/relaxation time scale (days)
for free-surface.
3.0000E+01 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
3.0000E+01 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
1.0000E+01 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.
35.000 S0 Background salinity (PSU) constant.
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
T LuvSrc Turning ON momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
T LtracerSrc(01) Turning ON point Sources/Sinks on tracer 01: temp
T LtracerSrc(02) Turning ON point Sources/Sinks on tracer 02: salt
F LsshCLM Turning OFF processing of SSH climatology.
F Lm2CLM Turning OFF processing of 2D momentum climatology.
T Lm3CLM Turning ON processing of 3D momentum climatology.
T LtracerCLM(01) Turning ON processing of climatology tracer 01: temp
T LtracerCLM(02) Turning ON processing of climatology tracer 02: salt
F LnudgeM2CLM Turning OFF nudging of 2D momentum climatology.
T LnudgeM3CLM Turning ON nudging of 3D momentum climatology.
T LnudgeTCLM(01) Turning ON nudging of climatology tracer 01: temp
T LnudgeTCLM(02) Turning ON nudging of climatology tracer 02: salt
T Hout(idFsur) Write out free-surface.
T Hout(idUvel) Write out 3D U-momentum component.
T Hout(idVvel) Write out 3D V-momentum component.
T Hout(idTvar) Write out tracer 01: temp
T Hout(idTvar) Write out tracer 02: salt
T Hout(idVvis) Write out vertical viscosity: AKv.
T Hout(idTdif) Write out vertical diffusion: AKt(itemp).
T Hout(idSdif) Write out vertical diffusion: AKt(isalt).

T Aout(idFsur) Write out averaged free-surface.
T Aout(idUvel) Write out averaged 3D U-momentum component.
T Aout(idVvel) Write out averaged 3D V-momentum component.
T Aout(idTvar) Write out averaged tracer 01: temp
T Aout(idTvar) Write out averaged tracer 02: salt
T Aout(idDano) Write out averaged density anomaly.
T Aout(idVvis) Write out averaged vertical viscosity: AKv.
T Aout(idTdif) Write out averaged vertical diffusion: AKt(itemp).
T Aout(idSdif) Write out averaged vertical diffusion: AKt(isalt).
T Aout(idHsbl) Write out averaged depth of surface boundary layer.
T Aout(idHbbl) Write out averaged depth of bottom boundary layer.

Output/Input Files:

Output Restart File: /home/19xcheng/roms/normal/ocean_rst.nc
Prefix for History Files: /home/19xcheng/roms/normal/ocean_his
Output Averages File: /home/19xcheng/roms/normal/ocean_avg.nc
Input Grid File: /home/19xcheng/roms/roms/Data/ROMS/Initial/ST_grid.nc
Input Nonlinear Initial File: /home/19xcheng/roms/roms/Data/ROMS/Initial/Yuedong_ini_ECCO_Y2021M5.nc
Input Sources/Sinks File: /home/19xcheng/roms/roms/Data/ROMS/Initial/Hanjiang_daily.nc
Input Forcing File 01: /home/19xcheng/roms/roms/Data/ROMS/Initial/ST_frc_new.nc
Input Forcing File 02: /home/19xcheng/roms/roms/Data/ROMS/Initial/ST_tide.nc
Input Climatology File: /home/19xcheng/roms/roms/Data/ROMS/Initial/ST_clm.nc
Input Nudge Coefficients File: /home/19xcheng/roms/roms/Data/ROMS/Initial/nudgcoef_b.nc
Input Boundary File: /home/19xcheng/roms/roms/Data/ROMS/Initial/Yuedong_bry_ECCO_Y2021M5.nc
/home/19xcheng/roms/roms/Data/ROMS/Initial/Yuedong_bry_ECCO_Y2021M6.nc
/home/19xcheng/roms/roms/Data/ROMS/Initial/Yuedong_bry_ECCO_Y2021M7.nc
/home/19xcheng/roms/roms/Data/ROMS/Initial/Yuedong_bry_ECCO_Y2021M8.nc

Tile partition information for Grid 01: 0574x0376x0030 tiling: 009x008

tile Istr Iend Jstr Jend Npts

0 1 63 1 47 88830
1 64 127 1 47 90240
2 128 191 1 47 90240
3 192 255 1 47 90240
4 256 319 1 47 90240
5 320 383 1 47 90240
6 384 447 1 47 90240
7 448 511 1 47 90240
8 512 574 1 47 88830
9 1 63 48 94 88830
10 64 127 48 94 90240
11 128 191 48 94 90240
12 192 255 48 94 90240
13 256 319 48 94 90240
14 320 383 48 94 90240
15 384 447 48 94 90240
16 448 511 48 94 90240
17 512 574 48 94 88830
18 1 63 95 141 88830
19 64 127 95 141 90240
20 128 191 95 141 90240
21 192 255 95 141 90240
22 256 319 95 141 90240
23 320 383 95 141 90240
24 384 447 95 141 90240
25 448 511 95 141 90240
26 512 574 95 141 88830
27 1 63 142 188 88830
28 64 127 142 188 90240
29 128 191 142 188 90240
30 192 255 142 188 90240
31 256 319 142 188 90240
32 320 383 142 188 90240
33 384 447 142 188 90240
34 448 511 142 188 90240
35 512 574 142 188 88830
36 1 63 189 235 88830
37 64 127 189 235 90240
38 128 191 189 235 90240
39 192 255 189 235 90240
40 256 319 189 235 90240
41 320 383 189 235 90240
42 384 447 189 235 90240
43 448 511 189 235 90240
44 512 574 189 235 88830
45 1 63 236 282 88830
46 64 127 236 282 90240
47 128 191 236 282 90240
48 192 255 236 282 90240
49 256 319 236 282 90240
50 320 383 236 282 90240
51 384 447 236 282 90240
52 448 511 236 282 90240
53 512 574 236 282 88830
54 1 63 283 329 88830
55 64 127 283 329 90240
56 128 191 283 329 90240
57 192 255 283 329 90240
58 256 319 283 329 90240
59 320 383 283 329 90240
60 384 447 283 329 90240
61 448 511 283 329 90240
62 512 574 283 329 88830
63 1 63 330 376 88830
64 64 127 330 376 90240
65 128 191 330 376 90240
66 192 255 330 376 90240
67 256 319 330 376 90240
68 320 383 330 376 90240
69 384 447 330 376 90240
70 448 511 330 376 90240
71 512 574 330 376 88830

Tile minimum and maximum fractional coordinates for Grid 01:
(interior points only)

tile Xmin Xmax Ymin Ymax grid

0 0.50 63.50 0.50 47.50 RHO-points
1 63.50 127.50 0.50 47.50 RHO-points
2 127.50 191.50 0.50 47.50 RHO-points
3 191.50 255.50 0.50 47.50 RHO-points
4 255.50 319.50 0.50 47.50 RHO-points
5 319.50 383.50 0.50 47.50 RHO-points
6 383.50 447.50 0.50 47.50 RHO-points
7 447.50 511.50 0.50 47.50 RHO-points
8 511.50 574.50 0.50 47.50 RHO-points
9 0.50 63.50 47.50 94.50 RHO-points
10 63.50 127.50 47.50 94.50 RHO-points
11 127.50 191.50 47.50 94.50 RHO-points
12 191.50 255.50 47.50 94.50 RHO-points
13 255.50 319.50 47.50 94.50 RHO-points
14 319.50 383.50 47.50 94.50 RHO-points
15 383.50 447.50 47.50 94.50 RHO-points
16 447.50 511.50 47.50 94.50 RHO-points
17 511.50 574.50 47.50 94.50 RHO-points
18 0.50 63.50 94.50 141.50 RHO-points
19 63.50 127.50 94.50 141.50 RHO-points
20 127.50 191.50 94.50 141.50 RHO-points
21 191.50 255.50 94.50 141.50 RHO-points
22 255.50 319.50 94.50 141.50 RHO-points
23 319.50 383.50 94.50 141.50 RHO-points
24 383.50 447.50 94.50 141.50 RHO-points
25 447.50 511.50 94.50 141.50 RHO-points
26 511.50 574.50 94.50 141.50 RHO-points
27 0.50 63.50 141.50 188.50 RHO-points
28 63.50 127.50 141.50 188.50 RHO-points
29 127.50 191.50 141.50 188.50 RHO-points
30 191.50 255.50 141.50 188.50 RHO-points
31 255.50 319.50 141.50 188.50 RHO-points
32 319.50 383.50 141.50 188.50 RHO-points
33 383.50 447.50 141.50 188.50 RHO-points
34 447.50 511.50 141.50 188.50 RHO-points
35 511.50 574.50 141.50 188.50 RHO-points
36 0.50 63.50 188.50 235.50 RHO-points
37 63.50 127.50 188.50 235.50 RHO-points
38 127.50 191.50 188.50 235.50 RHO-points
39 191.50 255.50 188.50 235.50 RHO-points
40 255.50 319.50 188.50 235.50 RHO-points
41 319.50 383.50 188.50 235.50 RHO-points
42 383.50 447.50 188.50 235.50 RHO-points
43 447.50 511.50 188.50 235.50 RHO-points
44 511.50 574.50 188.50 235.50 RHO-points
45 0.50 63.50 235.50 282.50 RHO-points
46 63.50 127.50 235.50 282.50 RHO-points
47 127.50 191.50 235.50 282.50 RHO-points
48 191.50 255.50 235.50 282.50 RHO-points
49 255.50 319.50 235.50 282.50 RHO-points
50 319.50 383.50 235.50 282.50 RHO-points
51 383.50 447.50 235.50 282.50 RHO-points
52 447.50 511.50 235.50 282.50 RHO-points
53 511.50 574.50 235.50 282.50 RHO-points
54 0.50 63.50 282.50 329.50 RHO-points
55 63.50 127.50 282.50 329.50 RHO-points
56 127.50 191.50 282.50 329.50 RHO-points
57 191.50 255.50 282.50 329.50 RHO-points
58 255.50 319.50 282.50 329.50 RHO-points
59 319.50 383.50 282.50 329.50 RHO-points
60 383.50 447.50 282.50 329.50 RHO-points
61 447.50 511.50 282.50 329.50 RHO-points
62 511.50 574.50 282.50 329.50 RHO-points
63 0.50 63.50 329.50 376.50 RHO-points
64 63.50 127.50 329.50 376.50 RHO-points
65 127.50 191.50 329.50 376.50 RHO-points
66 191.50 255.50 329.50 376.50 RHO-points
67 255.50 319.50 329.50 376.50 RHO-points
68 319.50 383.50 329.50 376.50 RHO-points
69 383.50 447.50 329.50 376.50 RHO-points
70 447.50 511.50 329.50 376.50 RHO-points
71 511.50 574.50 329.50 376.50 RHO-points

0 1.00 63.50 0.50 47.50 U-points
1 63.50 127.50 0.50 47.50 U-points
2 127.50 191.50 0.50 47.50 U-points
3 191.50 255.50 0.50 47.50 U-points
4 255.50 319.50 0.50 47.50 U-points
5 319.50 383.50 0.50 47.50 U-points
6 383.50 447.50 0.50 47.50 U-points
7 447.50 511.50 0.50 47.50 U-points
8 511.50 574.00 0.50 47.50 U-points
9 1.00 63.50 47.50 94.50 U-points
10 63.50 127.50 47.50 94.50 U-points
11 127.50 191.50 47.50 94.50 U-points
12 191.50 255.50 47.50 94.50 U-points
13 255.50 319.50 47.50 94.50 U-points
14 319.50 383.50 47.50 94.50 U-points
15 383.50 447.50 47.50 94.50 U-points
16 447.50 511.50 47.50 94.50 U-points
17 511.50 574.00 47.50 94.50 U-points
18 1.00 63.50 94.50 141.50 U-points
19 63.50 127.50 94.50 141.50 U-points
20 127.50 191.50 94.50 141.50 U-points
21 191.50 255.50 94.50 141.50 U-points
22 255.50 319.50 94.50 141.50 U-points
23 319.50 383.50 94.50 141.50 U-points
24 383.50 447.50 94.50 141.50 U-points
25 447.50 511.50 94.50 141.50 U-points
26 511.50 574.00 94.50 141.50 U-points
27 1.00 63.50 141.50 188.50 U-points
28 63.50 127.50 141.50 188.50 U-points
29 127.50 191.50 141.50 188.50 U-points
30 191.50 255.50 141.50 188.50 U-points
31 255.50 319.50 141.50 188.50 U-points
32 319.50 383.50 141.50 188.50 U-points
33 383.50 447.50 141.50 188.50 U-points
34 447.50 511.50 141.50 188.50 U-points
35 511.50 574.00 141.50 188.50 U-points
36 1.00 63.50 188.50 235.50 U-points
37 63.50 127.50 188.50 235.50 U-points
38 127.50 191.50 188.50 235.50 U-points
39 191.50 255.50 188.50 235.50 U-points
40 255.50 319.50 188.50 235.50 U-points
41 319.50 383.50 188.50 235.50 U-points
42 383.50 447.50 188.50 235.50 U-points
43 447.50 511.50 188.50 235.50 U-points
44 511.50 574.00 188.50 235.50 U-points
45 1.00 63.50 235.50 282.50 U-points
46 63.50 127.50 235.50 282.50 U-points
47 127.50 191.50 235.50 282.50 U-points
48 191.50 255.50 235.50 282.50 U-points
49 255.50 319.50 235.50 282.50 U-points
50 319.50 383.50 235.50 282.50 U-points
51 383.50 447.50 235.50 282.50 U-points
52 447.50 511.50 235.50 282.50 U-points
53 511.50 574.00 235.50 282.50 U-points
54 1.00 63.50 282.50 329.50 U-points
55 63.50 127.50 282.50 329.50 U-points
56 127.50 191.50 282.50 329.50 U-points
57 191.50 255.50 282.50 329.50 U-points
58 255.50 319.50 282.50 329.50 U-points
59 319.50 383.50 282.50 329.50 U-points
60 383.50 447.50 282.50 329.50 U-points
61 447.50 511.50 282.50 329.50 U-points
62 511.50 574.00 282.50 329.50 U-points
63 1.00 63.50 329.50 376.50 U-points
64 63.50 127.50 329.50 376.50 U-points
65 127.50 191.50 329.50 376.50 U-points
66 191.50 255.50 329.50 376.50 U-points
67 255.50 319.50 329.50 376.50 U-points
68 319.50 383.50 329.50 376.50 U-points
69 383.50 447.50 329.50 376.50 U-points
70 447.50 511.50 329.50 376.50 U-points
71 511.50 574.00 329.50 376.50 U-points

0 0.50 63.50 1.00 47.50 V-points
1 63.50 127.50 1.00 47.50 V-points
2 127.50 191.50 1.00 47.50 V-points
3 191.50 255.50 1.00 47.50 V-points
4 255.50 319.50 1.00 47.50 V-points
5 319.50 383.50 1.00 47.50 V-points
6 383.50 447.50 1.00 47.50 V-points
7 447.50 511.50 1.00 47.50 V-points
8 511.50 574.50 1.00 47.50 V-points
9 0.50 63.50 47.50 94.50 V-points
10 63.50 127.50 47.50 94.50 V-points
11 127.50 191.50 47.50 94.50 V-points
12 191.50 255.50 47.50 94.50 V-points
13 255.50 319.50 47.50 94.50 V-points
14 319.50 383.50 47.50 94.50 V-points
15 383.50 447.50 47.50 94.50 V-points
16 447.50 511.50 47.50 94.50 V-points
17 511.50 574.50 47.50 94.50 V-points
18 0.50 63.50 94.50 141.50 V-points
19 63.50 127.50 94.50 141.50 V-points
20 127.50 191.50 94.50 141.50 V-points
21 191.50 255.50 94.50 141.50 V-points
22 255.50 319.50 94.50 141.50 V-points
23 319.50 383.50 94.50 141.50 V-points
24 383.50 447.50 94.50 141.50 V-points
25 447.50 511.50 94.50 141.50 V-points
26 511.50 574.50 94.50 141.50 V-points
27 0.50 63.50 141.50 188.50 V-points
28 63.50 127.50 141.50 188.50 V-points
29 127.50 191.50 141.50 188.50 V-points
30 191.50 255.50 141.50 188.50 V-points
31 255.50 319.50 141.50 188.50 V-points
32 319.50 383.50 141.50 188.50 V-points
33 383.50 447.50 141.50 188.50 V-points
34 447.50 511.50 141.50 188.50 V-points
35 511.50 574.50 141.50 188.50 V-points
36 0.50 63.50 188.50 235.50 V-points
37 63.50 127.50 188.50 235.50 V-points
38 127.50 191.50 188.50 235.50 V-points
39 191.50 255.50 188.50 235.50 V-points
40 255.50 319.50 188.50 235.50 V-points
41 319.50 383.50 188.50 235.50 V-points
42 383.50 447.50 188.50 235.50 V-points
43 447.50 511.50 188.50 235.50 V-points
44 511.50 574.50 188.50 235.50 V-points
45 0.50 63.50 235.50 282.50 V-points
46 63.50 127.50 235.50 282.50 V-points
47 127.50 191.50 235.50 282.50 V-points
48 191.50 255.50 235.50 282.50 V-points
49 255.50 319.50 235.50 282.50 V-points
50 319.50 383.50 235.50 282.50 V-points
51 383.50 447.50 235.50 282.50 V-points
52 447.50 511.50 235.50 282.50 V-points
53 511.50 574.50 235.50 282.50 V-points
54 0.50 63.50 282.50 329.50 V-points
55 63.50 127.50 282.50 329.50 V-points
56 127.50 191.50 282.50 329.50 V-points
57 191.50 255.50 282.50 329.50 V-points
58 255.50 319.50 282.50 329.50 V-points
59 319.50 383.50 282.50 329.50 V-points
60 383.50 447.50 282.50 329.50 V-points
61 447.50 511.50 282.50 329.50 V-points
62 511.50 574.50 282.50 329.50 V-points
63 0.50 63.50 329.50 376.00 V-points
64 63.50 127.50 329.50 376.00 V-points
65 127.50 191.50 329.50 376.00 V-points
66 191.50 255.50 329.50 376.00 V-points
67 255.50 319.50 329.50 376.00 V-points
68 319.50 383.50 329.50 376.00 V-points
69 383.50 447.50 329.50 376.00 V-points
70 447.50 511.50 329.50 376.00 V-points
71 511.50 574.50 329.50 376.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 152
HaloSizeJ(1) = 118
TileSide(1) = 70
TileSize(1) = 3710


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 Clamped Rad + Nud Rad + Nud Closed

temp 1 Closed Rad + Nud Rad + Nud Closed

salt 1 Closed Rad + Nud Rad + Nud Closed

Activated C-preprocessing Options:

RIVERPLUME Plume
ADD_FSOBC Adding tidal elevation to processed OBC data.
ADD_M2OBC Adding tidal currents to processed OBC data.
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged nonlinear model fields.
CURVGRID Orthogonal curvilinear grid.
DIFF_3DCOEF Horizontal, time-dependent 3D diffusion coefficient.
DIURNAL_SRFLUX Modulate shortwave radiation by the local diurnal cycle.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
LIMIT_BSTRESS Limit bottom stress to maintain bottom velocity direction.
LMD_BKPP KPP bottom boundary layer mixing.
LMD_CONVEC LMD convective mixing due to shear instability.
LMD_DDMIX LMD double-diffusive mixing.
LMD_MIXING Large/McWilliams/Doney interior mixing.
LMD_NONLOCAL LMD convective nonlocal transport.
LMD_RIMIX LMD diffusivity due to shear instability.
LMD_SKPP KPP surface boundary layer mixing.
MASKING Land/Sea masking.
MPI MPI distributed-memory configuration.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
QCORRECTION Surface net heat flux correction.
RADIATION_2D Use tangential phase speed in radiation conditions.
RAMP_TIDES Ramping tidal forcing for one day.
RST_SINGLE Single precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLAR_SOURCE Solar Radiation Source Term.
SOLVE3D Solving 3D Primitive Equations.
SPLINES Conservative parabolic spline reconstruction.
SPHERICAL Spherical grid configuration.
SSH_TIDES Add tidal elevation to SSH climatology.
TS_C4HADVECTION Fourth-order centered horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection 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.
VAR_RHO_2D Variable density barotropic mode.
VISC_3DCOEF Horizontal, time-dependent 3D viscosity coefficient.

Process Information:

Node # 0 (pid= 188606) is active.
Node # 1 (pid= 188607) is active.
Node # 7 (pid= 188613) is active.
Node # 22 (pid= 188628) is active.
Node # 25 (pid= 188631) is active.
Node # 41 (pid= 188647) is active.
Node # 44 (pid= 188650) is active.
Node # 45 (pid= 188651) is active.
Node # 2 (pid= 188608) is active.
Node # 4 (pid= 188610) is active.
Node # 6 (pid= 188612) is active.
Node # 8 (pid= 188614) is active.
Node # 18 (pid= 188624) is active.
Node # 23 (pid= 188629) is active.
Node # 26 (pid= 188632) is active.
Node # 31 (pid= 188637) is active.
Node # 32 (pid= 188638) is active.
Node # 39 (pid= 188645) is active.
Node # 46 (pid= 188652) is active.
Node # 66 (pid= 188672) is active.
Node # 67 (pid= 188673) is active.
Node # 68 (pid= 188674) is active.
Node # 69 (pid= 188675) is active.
Node # 70 (pid= 188676) is active.
Node # 3 (pid= 188609) is active.
Node # 5 (pid= 188611) is active.
Node # 11 (pid= 188617) is active.
Node # 14 (pid= 188620) is active.
Node # 17 (pid= 188623) is active.
Node # 28 (pid= 188634) is active.
Node # 29 (pid= 188635) is active.
Node # 30 (pid= 188636) is active.
Node # 33 (pid= 188639) is active.
Node # 38 (pid= 188644) is active.
Node # 40 (pid= 188646) is active.
Node # 43 (pid= 188649) is active.
Node # 47 (pid= 188653) is active.
Node # 52 (pid= 188658) is active.
Node # 55 (pid= 188661) is active.
Node # 57 (pid= 188663) is active.
Node # 58 (pid= 188664) is active.
Node # 59 (pid= 188665) is active.
Node # 60 (pid= 188666) is active.
Node # 61 (pid= 188667) is active.
Node # 63 (pid= 188669) is active.
Node # 10 (pid= 188616) is active.
Node # 12 (pid= 188618) is active.
Node # 13 (pid= 188619) is active.
Node # 15 (pid= 188621) is active.
Node # 16 (pid= 188622) is active.
Node # 20 (pid= 188626) is active.
Node # 21 (pid= 188627) is active.
Node # 24 (pid= 188630) is active.
Node # 27 (pid= 188633) is active.
Node # 35 (pid= 188641) is active.
Node # 37 (pid= 188643) is active.
Node # 51 (pid= 188657) is active.
Node # 56 (pid= 188662) is active.
Node # 62 (pid= 188668) is active.
Node # 65 (pid= 188671) is active.
Node # 71 (pid= 188677) is active.
Node # 9 (pid= 188615) is active.
Node # 19 (pid= 188625) is active.
Node # 34 (pid= 188640) is active.
Node # 36 (pid= 188642) is active.
Node # 42 (pid= 188648) is active.
Node # 48 (pid= 188654) is active.
Node # 49 (pid= 188655) is active.
Node # 50 (pid= 188656) is active.
Node # 53 (pid= 188659) is active.
Node # 54 (pid= 188660) is active.
Node # 64 (pid= 188670) is active.

INITIAL: Configuring and initializing forward nonlinear model ...
*******

Vertical S-coordinate System, Grid 01:

level S-coord Cs-curve Z at hmin at hc half way at hmax

30 0.0000000 0.0000000 0.000 0.000 0.000 0.000
29 -0.0333333 -0.0006695 -0.051 -0.051 -0.768 -1.437
28 -0.0666667 -0.0026838 -0.104 -0.104 -2.879 -5.559
27 -0.1000000 -0.0060605 -0.159 -0.159 -6.351 -12.402
26 -0.1333333 -0.0108285 -0.216 -0.216 -11.211 -22.024
25 -0.1666667 -0.0170290 -0.276 -0.276 -17.502 -34.506
24 -0.2000000 -0.0247149 -0.337 -0.337 -25.276 -49.955
23 -0.2333333 -0.0339514 -0.401 -0.401 -34.599 -68.500
22 -0.2666667 -0.0448161 -0.467 -0.467 -45.547 -90.297
21 -0.3000000 -0.0573991 -0.536 -0.536 -58.211 -115.525
20 -0.3333333 -0.0718033 -0.608 -0.608 -72.693 -144.390
19 -0.3666667 -0.0881450 -0.682 -0.682 -89.110 -177.124
18 -0.4000000 -0.1065532 -0.760 -0.760 -107.591 -213.985
17 -0.4333333 -0.1271703 -0.841 -0.841 -128.277 -255.258
16 -0.4666667 -0.1501521 -0.925 -0.925 -151.324 -301.252
15 -0.5000000 -0.1756670 -1.014 -1.014 -176.901 -352.305
14 -0.5333333 -0.2038961 -1.106 -1.106 -205.187 -408.779
13 -0.5666667 -0.2350323 -1.203 -1.203 -236.377 -471.058
12 -0.6000000 -0.2692794 -1.304 -1.304 -270.673 -539.549
11 -0.6333333 -0.3068504 -1.410 -1.410 -308.287 -614.679
10 -0.6666667 -0.3479655 -1.522 -1.522 -349.441 -696.886
9 -0.7000000 -0.3928500 -1.639 -1.639 -394.358 -786.620
8 -0.7333333 -0.4417302 -1.763 -1.763 -443.265 -884.334
7 -0.7666667 -0.4948299 -1.892 -1.892 -496.385 -990.474
6 -0.8000000 -0.5523651 -2.029 -2.029 -553.934 -1105.472
5 -0.8333333 -0.6145374 -2.172 -2.172 -616.114 -1229.730
4 -0.8666667 -0.6815273 -2.322 -2.322 -683.103 -1363.609
3 -0.9000000 -0.7534850 -2.480 -2.480 -755.053 -1507.409
2 -0.9333333 -0.8305209 -2.646 -2.646 -832.074 -1661.350
1 -0.9666667 -0.9126940 -2.819 -2.819 -914.224 -1825.550
0 -1.0000000 -1.0000000 -3.000 -3.000 -1001.500 -2000.000

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.0011514592651644 0.0501156998674241-0.0034654566136460 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.0000580677400949 0.3006719858911027
7 0.0054132615310267 0.0499970966129953 0.0054713292711215 0.3506690825040980
8 0.0084687837865133 0.0497264335364439 0.0139401130576348 0.4003955160405419
9 0.0120633394191050 0.0493029943471183 0.0260034524767398 0.4496985103876601
10 0.0161716623600090 0.0486998273761630 0.0421751148367488 0.4983983377638231
11 0.0207585511322367 0.0478912442581626 0.0629336659689855 0.5462895820219857
12 0.0257765478740990 0.0468533167015507 0.0887102138430846 0.5931428987235364
13 0.0311633730493853 0.0455644893078458 0.1198735868924699 0.6387073880313822
14 0.0368391158442262 0.0440063206553765 0.1567127027366961 0.6827137086867586
15 0.0427031802506397 0.0421643648631652 0.1994158829873358 0.7248780735499238
16 0.0486309868367616 0.0400292058506332 0.2480468698240974 0.7649072794005570
17 0.0544704302037592 0.0375976565087951 0.3025173000278565 0.8025049359093521
18 0.0600380921294286 0.0348741349986072 0.3625553921572851 0.8373790709079593
19 0.0651152103984763 0.0318722303921357 0.4276706025557614 0.8692513013000950
20 0.0694434033194839 0.0286164698722119 0.4971140058752453 0.8978677711723069
21 0.0727201499285570 0.0251442997062377 0.5698341558038023 0.9230120708785446
22 0.0745940258796570 0.0215082922098099 0.6444281816834592 0.9445203630883545
23 0.0746596950216179 0.0177785909158270 0.7190878767050771 0.9622989540041815
24 0.0724526566618460 0.0140456061647461 0.7915405333669231 0.9763445601689277
25 0.0674437485167025 0.0104229733316538 0.8589842818836255 0.9867675335005816
26 0.0590334053485720 0.0070507859058187 0.9180176872321975 0.9938183194064003
27 0.0465456732896125 0.0040991156383901 0.9645633605218099 0.9979174350447905
28 0.0292219798521905 0.0017718319739095 0.9937853403740003 0.9996892670187000
29 0.0062146596259993 0.0003107329813000 0.9999999999999997 1.0000000000000000

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.53347320E-01 km
Maximum X-grid spacing, DXmax = 3.57200684E+01 km
Minimum Y-grid spacing, DYmin = 3.13846812E-01 km
Maximum Y-grid spacing, DYmax = 1.83601711E+01 km
Minimum Z-grid spacing, DZmin = 5.10042627E-02 m
Maximum Z-grid spacing, DZmax = 1.74450342E+02 m

Minimum barotropic Courant Number = 1.82372013E-03
Maximum barotropic Courant Number = 9.68670119E-01
Maximum Coriolis Courant Number = 3.28529472E-03

Horizontal mixing scaled by grid size, GRDMAX = 2.52100647E+01 km

Minimum horizontal diffusion coefficient = 0.00000000E+00 墉+? Maximum horizontal diffusion coefficient = 0.00000000E+00 墉+?
Minimum horizontal viscosity coefficient = 0.00000000E+00 墉+? Maximum horizontal viscosity coefficient = 0.00000000E+00 墉+?
NLM: GET_STATE - Read state initial conditions, t = 44315 00:00:00
(Grid 01, File: Yuedong_ini_ECCO_Y2021M5.nc, Rec=0001, Index=1)
- free-surface
(Min = 4.81336866E-01 Max = 1.04336820E+00)
- vertically integrated u-momentum component
(Min = -5.85548353E-02 Max = 4.78658900E-02)
- vertically integrated v-momentum component
(Min = -2.66768481E-02 Max = 3.97719391E-02)
- u-momentum component
(Min = -3.63861133E-01 Max = 4.76323310E-01)
- v-momentum component
(Min = -3.95165058E-01 Max = 3.74744901E-01)
- potential temperature
(Min = 2.27813267E+00 Max = 2.88403156E+01)
- salinity
(Min = 3.15444638E+01 Max = 3.46058391E+01)
GET_NGFLD - tidal period
(Min = 4.30819210E+04 Max = 2.38071599E+06)
GET_2DFLD - tidal elevation amplitude
(Min = 2.12440005E-04 Max = 1.84178955E+01)
GET_2DFLD - tidal elevation phase angle
(Min = 1.48785271E-05 Max = 6.28317460E+00)
GET_2DFLD - tidal current inclination angle
(Min = 1.69949498E-06 Max = 6.28318475E+00)
GET_2DFLD - tidal current phase angle
(Min = 5.10562512E-06 Max = 6.28318513E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 1.70600766E-06 Max = 1.51971823E+00)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -1.50056927E-01 Max = 3.51568002E-01)
GET_NGFLD - river runoff XI-positions at RHO-points
(Min = 6.30000000E+01 Max = 3.50000000E+02)
GET_NGFLD - river runoff ETA-positions at RHO-points
(Min = 3.28000000E+02 Max = 3.70000000E+02)
GET_NGFLD - river runoff direction
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - river runoff mass transport vertical profile
(Min = 3.33333333E-02 Max = 3.33333333E-02)
GET_NGFLD - river runoff mass transport, t = 44315 00:00:00
(Rec=0000121, Index=2, File: Hanjiang_daily.nc)
(Tmin= 44195.0000 Tmax= 44467.0000)
(Min = 6.53333333E+01 Max = 2.86500000E+03)
GET_NGFLD - river runoff potential temperature, t = 44315 00:00:00
(Rec=0000121, Index=2, File: Hanjiang_daily.nc)
(Tmin= 44195.0000 Tmax= 44467.0000)
(Min = 2.11298064E+01 Max = 2.11298064E+01)
GET_NGFLD - river runoff salinity, t = 44315 00:00:00
(Rec=0000121, Index=2, File: Hanjiang_daily.nc)
(Tmin= 44195.0000 Tmax= 44467.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface u-momentum stress, t = 44315 00:00:00
(Rec=0000005, Index=1, File: ST_frc_new.nc)
(Tmin= 44314.0000 Tmax= 44438.7500)
(Min = -2.15780504E-05 Max = 8.66460633E-05)
GET_2DFLD - surface v-momentum stress, t = 44315 00:00:00
(Rec=0000005, Index=1, File: ST_frc_new.nc)
(Tmin= 44314.0000 Tmax= 44438.7500)
(Min = -9.40035117E-06 Max = 3.82513170E-05)
GET_2DFLD - solar shortwave radiation flux, t = 106 12:45:04
(Rec=0000004, Index=1, File: ST_frc_new.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = 4.74603928E-05 Max = 6.46473377E-05)
GET_2DFLD - surface net heat flux, t = 106 12:45:04
(Rec=0000004, Index=1, File: ST_frc_new.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = -1.15152480E-05 Max = 3.33501383E-05)
GET_2DFLD - sea surface temperature climatology, t = 106 12:45:04
(Rec=0000004, Index=1, File: ST_frc_new.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = 1.74614755E+01 Max = 2.70000000E+01)
GET_2DFLD - surface net heat flux sensitivity to SST, t = 106 12:45:04
(Rec=0000004, Index=1, File: ST_frc_new.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = -1.05280444E-05 Max = -6.93547198E-06)
GET_2DFLD - surface net freswater flux, (E-P), t = 106 12:45:04
(Rec=0000004, Index=1, File: ST_frc_new.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = -2.62930894E-08 Max = 3.23239649E-08)
GET_NGFLD - free-surface eastern boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = 4.81336866E-01 Max = 9.22173492E-01)
GET_NGFLD - free-surface southern boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = 4.85756595E-01 Max = 1.04336820E+00)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -1.91002008E-03 Max = 1.65378829E-02)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -8.33784421E-03 Max = 3.97719391E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -2.97307725E-02 Max = 4.78658900E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -2.66768481E-02 Max = 3.97719391E-02)
GET_NGFLD - 3D v-momentum western boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -1.52677455E-01 Max = 9.86929809E-02)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -1.61544932E-01 Max = 1.02712584E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -3.04256552E-01 Max = 4.45800260E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -3.42674777E-01 Max = 3.52863544E-01)
GET_NGFLD - temperature eastern boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = 2.16154427E+01 Max = 2.71410352E+01)
GET_NGFLD - salinity eastern boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = 3.28669468E+01 Max = 3.36892556E+01)
GET_NGFLD - temperature southern boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = 2.27873434E+00 Max = 2.88403156E+01)
GET_NGFLD - salinity southern boundary condition, t = 44314 00:00:00
(Rec=0000001, Index=2, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = 3.27062456E+01 Max = 3.46058391E+01)
GET_3DFLD - u-momentum component climatology, t = 15 05:15:04
(Rec=0000001, Index=1, File: ST_clm.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = -4.80002956E-01 Max = 5.63855070E-01)
GET_3DFLD - v-momentum component climatology, t = 15 05:15:04
(Rec=0000001, Index=1, File: ST_clm.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = -3.15815396E-01 Max = 2.73184611E-01)
GET_3DFLD - potential temperature climatology, t = 15 05:15:04
(Rec=0000001, Index=1, File: ST_clm.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = 0.00000000E+00 Max = 2.43247781E+01)
GET_3DFLD - salinity climatology, t = 15 05:15:04
(Rec=0000001, Index=1, File: ST_clm.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = 3.27365877E+01 Max = 3.50000000E+01)

Basin information for Grid 01:

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

Initial basin volumes: TotVolume = 1.0636556266E+14 m3
MinVolume = 5.8101952466E+03 m3
MaxVolume = 1.2398779567E+10 m3
Max/Min = 2.1339695210E+06

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00212544)

GET_NGFLD - river runoff mass transport, t = 44316 00:00:00
(Rec=0000122, Index=1, File: Hanjiang_daily.nc)
(Tmin= 44195.0000 Tmax= 44467.0000)
(Min = 4.00000000E+01 Max = 2.84500000E+03)
GET_NGFLD - river runoff potential temperature, t = 44316 00:00:00
(Rec=0000122, Index=1, File: Hanjiang_daily.nc)
(Tmin= 44195.0000 Tmax= 44467.0000)
(Min = 2.12231645E+01 Max = 2.12231645E+01)
GET_NGFLD - river runoff salinity, t = 44316 00:00:00
(Rec=0000122, Index=1, File: Hanjiang_daily.nc)
(Tmin= 44195.0000 Tmax= 44467.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface u-momentum stress, t = 44315 06:00:00
(Rec=0000006, Index=2, File: ST_frc_new.nc)
(Tmin= 44314.0000 Tmax= 44438.7500)
(Min = -2.19630449E-05 Max = 8.25646255E-05)
GET_2DFLD - surface v-momentum stress, t = 44315 06:00:00
(Rec=0000006, Index=2, File: ST_frc_new.nc)
(Tmin= 44314.0000 Tmax= 44438.7500)
(Min = -1.72191299E-05 Max = 6.50793943E-05)
GET_2DFLD - solar shortwave radiation flux, t = 136 23:15:04
(Rec=0000005, Index=2, File: ST_frc_new.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = 5.43239961E-05 Max = 6.88190550E-05)
GET_2DFLD - surface net heat flux, t = 136 23:15:04
(Rec=0000005, Index=2, File: ST_frc_new.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = 1.09375809E-05 Max = 3.70641525E-05)
GET_2DFLD - sea surface temperature climatology, t = 136 23:15:04
(Rec=0000005, Index=2, File: ST_frc_new.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = 2.12222387E+01 Max = 2.90939773E+01)
GET_2DFLD - surface net heat flux sensitivity to SST, t = 136 23:15:04
(Rec=0000005, Index=2, File: ST_frc_new.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = -1.02865180E-05 Max = -7.38687187E-06)
GET_2DFLD - surface net freswater flux, (E-P), t = 136 23:15:04
(Rec=0000005, Index=2, File: ST_frc_new.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = -6.72115299E-08 Max = 2.44504539E-08)
GET_NGFLD - free-surface eastern boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = 4.81336866E-01 Max = 9.22173492E-01)
GET_NGFLD - free-surface southern boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = 4.85756595E-01 Max = 1.04336820E+00)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -1.91002008E-03 Max = 1.65378829E-02)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -8.33784421E-03 Max = 3.97719391E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -2.97307725E-02 Max = 4.78658900E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -2.66768481E-02 Max = 3.97719391E-02)
GET_NGFLD - 3D v-momentum western boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -1.52677455E-01 Max = 9.86929809E-02)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -1.61544932E-01 Max = 1.02712584E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -3.04256552E-01 Max = 4.45800260E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = -3.42674777E-01 Max = 3.52863544E-01)
GET_NGFLD - temperature eastern boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = 2.16154427E+01 Max = 2.71410352E+01)
GET_NGFLD - salinity eastern boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = 3.28669468E+01 Max = 3.36892556E+01)
GET_NGFLD - temperature southern boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = 2.27873434E+00 Max = 2.88403156E+01)
GET_NGFLD - salinity southern boundary condition, t = 44317 00:00:00
(Rec=0000002, Index=1, File: Yuedong_bry_ECCO_Y2021M5.nc)
(Tmin= 44314.0000 Tmax= 44347.0000)
(Min = 3.27062456E+01 Max = 3.46058391E+01)
GET_3DFLD - u-momentum component climatology, t = 45 15:45:04
(Rec=0000002, Index=2, File: ST_clm.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = -3.90936258E-01 Max = 6.78131241E-01)
GET_3DFLD - v-momentum component climatology, t = 45 15:45:04
(Rec=0000002, Index=2, File: ST_clm.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = -2.83800600E-01 Max = 3.73359857E-01)
GET_3DFLD - potential temperature climatology, t = 45 15:45:04
(Rec=0000002, Index=2, File: ST_clm.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = 0.00000000E+00 Max = 2.66074829E+01)
GET_3DFLD - salinity climatology, t = 45 15:45:04
(Rec=0000002, Index=2, File: ST_clm.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = 3.27471674E+01 Max = 3.50000000E+01)

jcwarner
Posts: 1172
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: too high bottom temperature result in estuary

#2 Unread post by jcwarner »

you need to look at how the temperature is evolving over time.
does it get warmer from the surface net fluxes and then mix downwards?
does it get warmer from the offshore boundary and then propagate in?

i see:
GET_3DFLD - potential temperature climatology, t = 45 15:45:04
(Rec=0000002, Index=2, File: ST_clm.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = 0.00000000E+00 Max = 2.66074829E+01)

GET_3DFLD - salinity climatology, t = 45 15:45:04
(Rec=0000002, Index=2, File: ST_clm.nc)
(Tmin= 15.2188 Tmax= 350.0313)
(Min = 3.27471674E+01 Max = 3.50000000E+01)

That Temp max look rather high (350)
and the Saltmax is large also.
-j

User avatar
wilkin
Posts: 875
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: too high bottom temperature result in estuary

#3 Unread post by wilkin »

Tmax is time.
The max value 26.6 is OK.

Consider your jerlov water type WTYPE. If it's 1 (clear ocean water) then in a shallow estuary a lot of sunlight reaches the seafloor and is trapped there because the bottom is effectively insulated. Use type 3 or 4 and absorption is more rapid keeping heating near surface.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

User avatar
susonic
Posts: 167
Joined: Tue Aug 21, 2007 5:44 pm
Location: UST21 / Korea
Contact:

Re: too high bottom temperature result in estuary

#4 Unread post by susonic »

#define TS_C4HADVECTION
#define TS_U3VADVECTION
needs to be written as
#define TS_U3HADVECTION
#define TS_C4VADVECTION


Seems like you provide 6 hourly wind stress while the heat flux is coming from mean value.

I would not recommend diurnal_srflux if this is the case. From what I experienced, it increases the SST.

Changing the water type will influence the structure of the temperature profile.
Joonho Lee

User avatar
arango
Site Admin
Posts: 1347
Joined: Wed Feb 26, 2003 4:41 pm
Location: DMCS, Rutgers University
Contact:

Re: too high bottom temperature result in estuary

#5 Unread post by arango »

WOW, all the tracer advection CPP options were eliminated two years ago! You are using a very old code, please update. Check the following :arrow: track ticket for details. ROMS is a living entity, still evolving and improving.

cheng_dada
Posts: 8
Joined: Fri Feb 26, 2021 2:20 pm
Location: Shantou University

Re: too high bottom temperature result in estuary

#6 Unread post by cheng_dada »

wilkin wrote: Mon Jun 13, 2022 3:58 pm Tmax is time.
The max value 26.6 is OK.

Consider your jerlov water type WTYPE. If it's 1 (clear ocean water) then in a shallow estuary a lot of sunlight reaches the seafloor and is trapped there because the bottom is effectively insulated. Use type 3 or 4 and absorption is more rapid keeping heating near surface.
Hi wilkin,
Thanks for your advice! There is a marked decline in temperature when I use type 5(Figure 1). But now there's a new problem, there has been a significant increase in temperatures in non-estuarine areas. Is there any tips to solve this problem?
Figure 1
Figure 1
3786cf2a61be7b3dd930a8f4927f2e8.png (62 KiB) Viewed 10011 times

User avatar
wilkin
Posts: 875
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: too high bottom temperature result in estuary

#7 Unread post by wilkin »

If you now have too-warm conditions at the surface where water is clearer (because heating is too surface intensified with wtype=5), then you can try a regionally varying water type. See https://www.myroms.org/projects/src/ticket/609. This must be static - not time varying.

Other users have had some success using salinity as a proxy for waters with increased absorption due to river sediment loads, or CDOM in terrestrial run-off, but that requires some hacking to the code. See the light absorption methodology in the fennel biology model for guidance if you want to go that route.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

Post Reply