BLOWING UP during nesting

Discussion of how to use ROMS on different regional and basin scale applications.

Moderators: arango, robertson

Post Reply
Message
Author
haili
Posts: 6
Joined: Wed Dec 26, 2018 4:37 pm
Location: Brown University

BLOWING UP during nesting

#1 Unread post by haili »

Hi all,

I run a coarse resolution case successfully, I would like to choose a small area of the case area to study the phenomenon of multi-layer nesting in the area.
But now only 2 layers nesting failed.
I do not know the reason.
I created the refine grid file with the Matlab code coarse2fine.
I extracted the refine initial file with the Matlab code d_roms2roms.
My refine tidal forcing files was interpolated from the coarse tidal forcing file.
The ngbay_ngc_2g_ab.nc was generated using the script contact.m.
Other forcing files are the same with the coarse case.
However at the first step the nesting case blow_up :cry: .

What I have already done:
1. check nan of the initial file and the tidal forcing file.
2. decrease the time steps 'dt' from 15 to 1
3. turn off the tidal forcing

Thank you for your help!

Best
Haili

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

7200 ntimes Number of timesteps for 3-D equations.
15.000 dt Timestep size (s) for 3-D equations.
15 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.
480 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
360 nSTA Number of timesteps between the writing of data
the stations file.
T ldefout Switch to create a new output NetCDF file(s).
0 nHIS Number of timesteps between the writing fields
into history file.
0 nQCK Number of timesteps between the writing fields
into quicksave file.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
480 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
5760 ndefAVG Number of timesteps between creation of new
time-averaged 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
2.0000E+00 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
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-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.
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.0000E+01 blk_ZQ Height (m) of surface air humidity measurement.
1.0000E+01 blk_ZT Height (m) of surface air temperature measurement.
1.0000E+01 blk_ZW Height (m) of surface winds measurement.
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
1.0000E+00 theta_s S-coordinate surface control parameter.
1.0000E+00 theta_b S-coordinate bottom control parameter.
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.
0.000 dstart Time-stamp assigned to model initialization (days).
0.000 tide_start Reference time origin for tidal forcing (days).
20060101.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.
1.5000E+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.
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 Aout(idFsur) Write out averaged free-surface.
T Aout(idUbar) Write out averaged 2D U-momentum component.
T Aout(idVbar) Write out averaged 2D V-momentum component.
T Aout(idUvel) Write out averaged 3D U-momentum component.
T Aout(idVvel) Write out averaged 3D V-momentum component.
T Aout(idWvel) Write out averaged W-momentum component.
T Aout(idOvel) Write out averaged omega vertical velocity.
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.
Output/Input Files:

Output Restart File: /users/hwang121/scratch/ngbayhlnestingfile/output/ocean_rst.nc
Prefix for Averages Files: /users/hwang121/scratch/ngbayhlnestingfile/output/ocean_avg
Output Stations File: /users/hwang121/scratch/ngbayhlnestingfile/output/ocean_sta.nc
Input Grid File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc
Nesting grid connectivity File: /users/hwang121/scratch/smallinput/ngbay_ngc_2g_ab.nc
Input Nonlinear Initial File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/osom_grid4_mindep_smlp_mod7_ini_hydro.nc
Tidal Forcing File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/osom_grid4_mindep_smlp_mod6_tidalforcing_2006.nc
Input Forcing File 01: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/frc_wind_field_OSOM_2006_new.nc
Input Forcing File 02: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/frc_bulk8_point_osom_2006.nc
Station positions File: /users/hwang121/data/hwang121/stations_osom_nesting.in

Resolution, Grid 02: 0900x1212x015, Parallel Nodes: 64, Tiling: 008x008


Physical Parameters, Grid: 02
=============================

108000 ntimes Number of timesteps for 3-D equations.
1.000 dt Timestep size (s) for 3-D equations.
15 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.
480 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
360 nSTA Number of timesteps between the writing of data
the stations file.
T ldefout Switch to create a new output NetCDF file(s).
0 nHIS Number of timesteps between the writing fields
into history file.
0 nQCK Number of timesteps between the writing fields
into quicksave file.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
480 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
5760 ndefAVG Number of timesteps between creation of new
time-averaged 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
2.0000E+00 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
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-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.
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.0000E+01 blk_ZQ Height (m) of surface air humidity measurement.
1.0000E+01 blk_ZT Height (m) of surface air temperature measurement.
1.0000E+01 blk_ZW Height (m) of surface winds measurement.
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
1.0000E+00 theta_s S-coordinate surface control parameter.
1.0000E+00 theta_b S-coordinate bottom control parameter.
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.
0.000 dstart Time-stamp assigned to model initialization (days).
0.000 tide_start Reference time origin for tidal forcing (days).
20060101.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.
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 Aout(idFsur) Write out averaged free-surface.
T Aout(idUbar) Write out averaged 2D U-momentum component.
T Aout(idVbar) Write out averaged 2D V-momentum component.
T Aout(idUvel) Write out averaged 3D U-momentum component.
T Aout(idVvel) Write out averaged 3D V-momentum component.
T Aout(idWvel) Write out averaged W-momentum component.
T Aout(idOvel) Write out averaged omega vertical velocity.
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.

Output/Input Files:

Output Restart File: /users/hwang121/scratch/ngbayhlnestingfile/output/oceanhl_rst.nc
Prefix for Averages Files: /users/hwang121/scratch/ngbayhlnestingfile/output/oceanhl_avg
Output Stations File: /users/hwang121/scratch/ngbayhlnestingfile/output/ocean_sta.nc
Input Grid File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc
Nesting grid connectivity File: /users/hwang121/scratch/smallinput/ngbay_ngc_2g_ab.nc
Input Nonlinear Initial File: /users/hwang121/scratch/smallinput/ngbay_grd_ini.nc
Input Forcing File 01: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/frc_wind_field_OSOM_2006_new.nc
Input Forcing File 02: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/frc_bulk8_point_osom_2006.nc
Station positions File: /users/hwang121/data/hwang121/stations_osom_nesting.in

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

GET_GRID - angle between XI-axis and EAST: angler
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -5.73923922E-01 Max = 9.08638297E-01)
GET_GRID - reciprocal XI-grid spacing: pm
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 2.29508154E-03 Max = 2.46946401E-02)
GET_GRID - reciprocal ETA-grid spacing: pn
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 2.57345733E-03 Max = 3.66368034E-02)
GET_GRID - XI-derivative of inverse metric factor pn: dndx
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -6.22353037E+00 Max = 1.29465657E+01)
GET_GRID - ETA-derivative of inverse metric factor pm: dmde
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -1.37933839E+00 Max = 2.00305560E+01)
GET_GRID - Coriolis parameter at RHO-points: f
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 9.47308397E-05 Max = 9.79166886E-05)
GET_GRID - bathymetry at RHO-points: h
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 2.00000000E+00 Max = 1.12145151E+02)
GET_GRID - x-location of RHO-points: x-rho
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -8.08048205E+06 Max = -7.78230926E+06)
GET_GRID - y-location of RHO-points: y_rho
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.93464274E+06 Max = 5.18168211E+06)
GET_GRID - x-location of PSI-points: x_psi
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -8.08034114E+06 Max = -7.78241485E+06)
GET_GRID - y-location of PSI-points: y-psi
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.93490009E+06 Max = 5.18150523E+06)
GET_GRID - x-location of U-points: x_u
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -8.08038301E+06 Max = -7.78236657E+06)
GET_GRID - y-location of U-points: y_u
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.93464257E+06 Max = 5.18157942E+06)
GET_GRID - x-location of V-points: x_v
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -8.08044156E+06 Max = -7.78235475E+06)
GET_GRID - y-location of V-points: y_v
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.93490036E+06 Max = 5.18160159E+06)
GET_GRID - longitude of RHO-points: lon_rho
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.93490036E+06 Max = 5.18160159E+06)
GET_GRID - latitude of RHO-points lat_rho
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.93490036E+06 Max = 5.18160159E+06)
GET_GRID - longitude of PSI-points: lon_psi
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -7.26635975E+01 Max = -6.99866847E+01)
GET_GRID - latitude of PSI-points lat_psi
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.05088007E+01 Max = 4.21724345E+01)
GET_GRID - longitude of U-points: lon_u
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -7.26639628E+01 Max = -6.99862328E+01)
GET_GRID - latitude of U-points: lat_u
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.05071487E+01 Max = 4.21728490E+01)
GET_GRID - longitude of V-points: lon_v
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -7.26650402E+01 Max = -6.99855448E+01)
GET_GRID - latitude of V-points: lat_v
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.05088013E+01 Max = 4.21737681E+01)
GET_GRID - mask on RHO-points: mask_rho
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on PSI-points: mask_psi
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on U-points: mask_u
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on V-points: mask_v
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)

GET_GRID - angle between XI-axis and EAST: angler
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 9.39339461E-02 Max = 4.88527516E-01)
GET_GRID - reciprocal XI-grid spacing: pm
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 2.04210158E-02 Max = 5.53111842E-02)
GET_GRID - reciprocal ETA-grid spacing: pn
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 2.34795268E-02 Max = 4.06173068E-02)
GET_GRID - XI-derivative of inverse metric factor pn: dndx
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -1.49524133E-03 Max = 1.36560555E-02)
GET_GRID - ETA-derivative of inverse metric factor pm: dmde
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -1.72030760E-02 Max = -1.14129901E-02)
GET_GRID - Coriolis parameter at RHO-points: f
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 9.65207017E-05 Max = 9.73545733E-05)
GET_GRID - bathymetry at RHO-points: h
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -1.65801391E+00 Max = 5.26751124E+01)
GET_GRID - x-location of RHO-points: x-rho
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -7.95159584E+06 Max = -7.89327807E+06)
GET_GRID - y-location of RHO-points: y_rho
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 5.08543749E+06 Max = 5.14402907E+06)
GET_GRID - x-location of PSI-points: x_psi
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -7.95157715E+06 Max = -7.89331918E+06)
GET_GRID - y-location of PSI-points: y-psi
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 5.08546640E+06 Max = 5.14400415E+06)
GET_GRID - x-location of U-points: x_u
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -7.95157754E+06 Max = -7.89330716E+06)
GET_GRID - y-location of U-points: y_u
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 5.08544056E+06 Max = 5.14401638E+06)
GET_GRID - x-location of V-points: x_v
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -7.95159545E+06 Max = -7.89329009E+06)
GET_GRID - y-location of V-points: y_v
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 5.08546333E+06 Max = 5.14401684E+06)
GET_GRID - longitude of RHO-points: lon_rho
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 5.08546333E+06 Max = 5.14401684E+06)
GET_GRID - latitude of RHO-points lat_rho
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 5.08546333E+06 Max = 5.14401684E+06)
GET_GRID - longitude of PSI-points: lon_psi
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -7.14921368E+01 Max = -7.10528888E+01)
GET_GRID - latitude of PSI-points lat_psi
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 4.14386890E+01 Max = 4.18767854E+01)
GET_GRID - longitude of U-points: lon_u
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -7.14921482E+01 Max = -7.10527933E+01)
GET_GRID - latitude of U-points: lat_u
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 4.14384986E+01 Max = 4.18768933E+01)
GET_GRID - longitude of V-points: lon_v
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -7.14922510E+01 Max = -7.10526126E+01)
GET_GRID - latitude of V-points: lat_v
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 4.14386717E+01 Max = 4.18768473E+01)
GET_GRID - mask on RHO-points: mask_rho
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on PSI-points: mask_psi
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 0.00000000E+00 Max = 2.00000000E+00)
GET_GRID - mask on U-points: mask_u
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on V-points: mask_v
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)

Vertical S-coordinate System, Grid 01:

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

15 0.0000000 0.0000000 0.000 0.000 0.000 0.000
14 -0.0666667 -0.0064624 -0.085 -0.110 -0.540 -0.901
13 -0.1333333 -0.0257201 -0.181 -0.239 -1.775 -3.199
12 -0.2000000 -0.0573871 -0.286 -0.386 -3.682 -6.852
11 -0.2666667 -0.1008307 -0.401 -0.551 -6.227 -11.792
10 -0.3333333 -0.1551864 -0.524 -0.733 -9.365 -17.924
9 -0.4000000 -0.2193793 -0.656 -0.929 -13.035 -25.130
8 -0.4666667 -0.2921503 -0.794 -1.138 -17.171 -33.273
7 -0.5333333 -0.3720885 -0.938 -1.358 -21.696 -42.199
6 -0.6000000 -0.4576670 -1.086 -1.587 -26.526 -51.741
5 -0.6666667 -0.5472830 -1.238 -1.821 -31.575 -61.724
4 -0.7333333 -0.6393009 -1.391 -2.059 -36.755 -71.969
3 -0.8000000 -0.7320964 -1.546 -2.298 -41.976 -82.299
2 -0.8666667 -0.8241010 -1.699 -2.536 -47.155 -92.543
1 -0.9333333 -0.9138456 -1.851 -2.771 -52.211 -102.540
0 -1.0000000 -1.0000000 -2.000 -3.000 -57.073 -112.145

Vertical S-coordinate System, Grid 02:

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

15 0.0000000 0.0000000 0.000 0.000 0.000 0.000
14 -0.0666667 -0.0064624 0.234 -0.110 -0.326 -0.511
13 -0.1333333 -0.0257201 0.442 -0.239 -0.945 -1.660
12 -0.2000000 -0.0573871 0.624 -0.386 -1.847 -3.428
11 -0.2666667 -0.1008307 0.782 -0.551 -3.017 -5.782
10 -0.3333333 -0.1551864 0.918 -0.733 -4.437 -8.680
9 -0.4000000 -0.2193793 1.033 -0.929 -6.081 -12.068
8 -0.4666667 -0.2921503 1.131 -1.138 -7.921 -15.884
7 -0.5333333 -0.3720885 1.215 -1.358 -9.924 -20.057
6 -0.6000000 -0.4576670 1.286 -1.587 -12.056 -24.512
5 -0.6666667 -0.5472830 1.350 -1.821 -14.281 -29.167
4 -0.7333333 -0.6393009 1.408 -2.059 -16.560 -33.942
3 -0.8000000 -0.7320964 1.466 -2.298 -18.857 -38.756
2 -0.8666667 -0.8241010 1.524 -2.536 -21.136 -43.530
1 -0.9333333 -0.9138456 1.587 -2.771 -23.363 -48.192
0 -1.0000000 -1.0000000 1.658 -3.000 -25.509 -52.675

Time Splitting Weights for Grid 01: ndtfast = 15 nfast = 22
==================================

Primary Secondary Accumulated to Current Step

1-0.0006740978872509 0.0666666666666667-0.0006740978872509 0.0666666666666667
2-0.0001913161153518 0.0667116065258167-0.0008654140026027 0.1333782731924834
3 0.0014472038369421 0.0667243609335069 0.0005817898343395 0.2001026341259903
4 0.0042370228855824 0.0666278806777107 0.0048188127199218 0.2667305148037010
5 0.0081665994120431 0.0663454124853385 0.0129854121319649 0.3330759272890396
6 0.0132119623624627 0.0658009725245357 0.0261973744944276 0.3988768998135752
7 0.0193298623751124 0.0649201750337048 0.0455272368695400 0.4637970748472801
8 0.0264494009361902 0.0636315175420307 0.0719766378057303 0.5274285923893107
9 0.0344621375639418 0.0618682241462847 0.1064387753696720 0.5892968165355954
10 0.0432106750211074 0.0595707483086885 0.1496494503907794 0.6488675648442840
11 0.0524757225556949 0.0566900366406147 0.2021251729464743 0.7055576014848987
12 0.0619616371700794 0.0531916551369017 0.2640868101165537 0.7587492566218004
13 0.0712804429184289 0.0490608793255631 0.3353672530349826 0.8078101359473635
14 0.0799343282324558 0.0443088497976678 0.4153015812674384 0.8521189857450313
15 0.0872966212754959 0.0389798945821708 0.5025982025429343 0.8910988803272021
16 0.0925912433249121 0.0331601198304711 0.5951894458678464 0.9242590001576731
17 0.0948706401828253 0.0269873702754769 0.6900600860506717 0.9512463704331501
18 0.0929921916151718 0.0206626609299552 0.7830522776658435 0.9719090313631054
19 0.0855930988190858 0.0144631814889438 0.8686453764849293 0.9863722128520491
20 0.0710637499186094 0.0087569749010047 0.9397091264035388 0.9951291877530538
21 0.0475195634887281 0.0040193915730974 0.9872286898922669 0.9991485793261512
22 0.0127713101077333 0.0008514206738489 1.0000000000000002 1.0000000000000000

ndtfast, nfast = 15 22 nfast/ndtfast = 1.46667

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

1.000000000000 1.070450639131 0.535225319565 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.09467

Time Splitting Weights for Grid 02: ndtfast = 15 nfast = 22
==================================

Primary Secondary Accumulated to Current Step

1-0.0006740978872509 0.0666666666666667-0.0006740978872509 0.0666666666666667
2-0.0001913161153518 0.0667116065258167-0.0008654140026027 0.1333782731924834
3 0.0014472038369421 0.0667243609335069 0.0005817898343395 0.2001026341259903
4 0.0042370228855824 0.0666278806777107 0.0048188127199218 0.2667305148037010
5 0.0081665994120431 0.0663454124853385 0.0129854121319649 0.3330759272890396
6 0.0132119623624627 0.0658009725245357 0.0261973744944276 0.3988768998135752
7 0.0193298623751124 0.0649201750337048 0.0455272368695400 0.4637970748472801
8 0.0264494009361902 0.0636315175420307 0.0719766378057303 0.5274285923893107
9 0.0344621375639418 0.0618682241462847 0.1064387753696720 0.5892968165355954
10 0.0432106750211074 0.0595707483086885 0.1496494503907794 0.6488675648442840
11 0.0524757225556949 0.0566900366406147 0.2021251729464743 0.7055576014848987
12 0.0619616371700794 0.0531916551369017 0.2640868101165537 0.7587492566218004
13 0.0712804429184289 0.0490608793255631 0.3353672530349826 0.8078101359473635
14 0.0799343282324558 0.0443088497976678 0.4153015812674384 0.8521189857450313
15 0.0872966212754959 0.0389798945821708 0.5025982025429343 0.8910988803272021
16 0.0925912433249121 0.0331601198304711 0.5951894458678464 0.9242590001576731
17 0.0948706401828253 0.0269873702754769 0.6900600860506717 0.9512463704331501
18 0.0929921916151718 0.0206626609299552 0.7830522776658435 0.9719090313631054
19 0.0855930988190858 0.0144631814889438 0.8686453764849293 0.9863722128520491
20 0.0710637499186094 0.0087569749010047 0.9397091264035388 0.9951291877530538
21 0.0475195634887281 0.0040193915730974 0.9872286898922669 0.9991485793261512
22 0.0127713101077333 0.0008514206738489 1.0000000000000002 1.0000000000000000

ndtfast, nfast = 15 22 nfast/ndtfast = 1.46667

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

1.000000000000 1.070450639131 0.535225319565 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.09467

Metrics information for Grid 01:
===============================

Minimum X-grid spacing, DXmin = 4.04946174E-02 km
Maximum X-grid spacing, DXmax = 3.56382138E-01 km
Minimum Y-grid spacing, DYmin = 2.72949578E-02 km
Maximum Y-grid spacing, DYmax = 3.88582312E-01 km
Minimum Z-grid spacing, DZmin = 8.51699402E-02 m
Maximum Z-grid spacing, DZmax = 1.03302185E+01 m

Minimum barotropic Courant Number = 4.20343636E-02
Maximum barotropic Courant Number = 2.97002616E-01
Maximum Coriolis Courant Number = 1.46019263E-03

Horizontal mixing scaled by grid size, GRDMAX = 3.23731774E-01 km

Metrics information for Grid 02:
===============================

Minimum X-grid spacing, DXmin = 1.80523345E-02 km
Maximum X-grid spacing, DXmax = 4.91059154E-02 km
Minimum Y-grid spacing, DYmin = 2.45584637E-02 km
Maximum Y-grid spacing, DYmax = 4.26169539E-02 km
Minimum Z-grid spacing, DZmin = -2.33860479E-01 m
Maximum Z-grid spacing, DZmax = 4.81384947E+00 m

Minimum barotropic Courant Number = 8.72759390E-03
Maximum barotropic Courant Number = 5.95826250E-02
Maximum Coriolis Courant Number = 9.72995226E-05

Horizontal mixing scaled by grid size, GRDMAX = 4.55502426E-02 km

Refined Nested Grid(s) Information:
==================================

Refined Donor Refined Timestep Refined
Grid Grid Scale Ratio Timesteps

02 01 03 15.00000 15

WARNING: Usually the number of Refined Timesteps must be the same
as the Refined Scale for numerical stability.


NLM: GET_STATE - Reading state initial conditions, 2006-01-01 00:00:00.00
(Grid 01, t = 0.0000, File: osom_grid4_mindep_smlp_mod7_ini_hydro.nc, Rec=0001, Index=1)
- free-surface
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- potential temperature
(Min = 4.56679989E+00 Max = 1.23957848E+01)
- salinity
(Min = 3.09659979E-09 Max = 3.36493083E+01)

NLM: GET_STATE - Reading state initial conditions, 2006-01-01 00:00:00.00
(Grid 02, t = 0.0000, File: ngbay_grd_ini.nc, Rec=0001, Index=1)
- free-surface
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- potential temperature
(Min = 5.98502991E+00 Max = 8.41852708E+00)
- salinity
(Min = 0.00000000E+00 Max = 3.17598902E+01)
GET_NGFLD - tidal period
(Grid = 01, Min = 1.49047214E+04 Max = 9.67260838E+04)
GET_2DFLD - tidal elevation amplitude
(Grid = 01, Min = 0.00000000E+00 Max = 7.54251636E-01)
GET_2DFLD - tidal elevation phase angle
(Grid = 01, Min = 6.26773407E-06 Max = 6.28318498E+00)
GET_2DFLD - tidal current inclination angle
(Grid = 01, Min = 0.00000000E+00 Max = 3.14159244E+00)
GET_2DFLD - tidal current phase angle
(Grid = 01, Min = 0.00000000E+00 Max = 6.28318364E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Grid = 01, Min = 0.00000000E+00 Max = 1.49311554E+00)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Grid = 01, Min = -3.37066252E-01 Max = 2.52296461E-01)
GET_2DFLD - surface u-wind component, 2005-12-31 23:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = -0.0417
(Min = -4.81854172E+00 Max = -1.58071958E+00) regrid = T
GET_2DFLD - surface v-wind component, 2005-12-31 23:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = -0.0417
(Min = -2.42937271E+00 Max = 1.13576238E+00) regrid = T
GET_2DFLD - surface air pressure, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 1.00983752E+03 Max = 1.00983752E+03) regrid = F
GET_2DFLD - solar shortwave radiation flux, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 8.57305843E-13 Max = 8.57305843E-13) regrid = F
GET_2DFLD - downwelling longwave radiation flux, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 7.64574591E-05 Max = 7.64574591E-05) regrid = F
GET_2DFLD - surface air temperature, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 3.75000000E-01 Max = 3.75000000E-01) regrid = F
GET_2DFLD - surface air relative humidity, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 9.24499969E-01 Max = 9.24499969E-01) regrid = F
GET_2DFLD - rain fall rate, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 1.17592590E-05 Max = 1.17592590E-05) regrid = F
GET_2DFLD - surface u-wind component, 2005-12-31 23:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = -0.0417
(Min = -3.84400734E+00 Max = -2.83434620E+00) regrid = T
GET_2DFLD - surface v-wind component, 2005-12-31 23:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = -0.0417
(Min = -2.09210197E+00 Max = -4.23612321E-01) regrid = T
GET_2DFLD - surface air pressure, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 1.00983752E+03 Max = 1.00983752E+03) regrid = F
GET_2DFLD - solar shortwave radiation flux, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 8.57305843E-13 Max = 8.57305843E-13) regrid = F
GET_2DFLD - downwelling longwave radiation flux, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 7.64574591E-05 Max = 7.64574591E-05) regrid = F
GET_2DFLD - surface air temperature, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 3.75000000E-01 Max = 3.75000000E-01) regrid = F
GET_2DFLD - surface air relative humidity, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 9.24499969E-01 Max = 9.24499969E-01) regrid = F
GET_2DFLD - rain fall rate, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 1.17592590E-05 Max = 1.17592590E-05) regrid = F

Basin information for Grid 01:

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

Initial domain volumes: TotVolume = 6.2370483442E+11 m3
MinCellVol = 3.4037041132E+02 m3
MaxCellVol = 6.8561716618E+05 m3
Max/Min = 2.0143265789E+03

Basin information for Grid 02:

Maximum grid stiffness ratios: rx0 = 2.592861E-01 (Beckmann and Haidvogel)
rx1 = 6.109278E+00 (Haney)

Initial domain volumes: TotVolume = 3.4468096079E+09 m3
MinCellVol = 4.2938854899E+01 m3
MaxCellVol = 6.6263574929E+03 m3
Max/Min = 1.5432077796E+02

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 000000000001 - 000000007200)
NL ROMS/TOMS: started time-stepping: (Grid: 02 TimeSteps: 000000000001 - 000000108000)

GET_2DFLD - surface u-wind component, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = 0.0000
(Min = -4.81854172E+00 Max = -1.58071958E+00) regrid = T
GET_2DFLD - surface v-wind component, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = 0.0000
(Min = -2.42937271E+00 Max = 1.13576238E+00) regrid = T
GET_2DFLD - surface air pressure, 2006-01-01 03:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 1.00822778E+03 Max = 1.00822778E+03) regrid = F
GET_2DFLD - solar shortwave radiation flux, 2006-01-01 03:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 8.57305843E-13 Max = 8.57305843E-13) regrid = F
GET_2DFLD - downwelling longwave radiation flux, 2006-01-01 03:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 7.66043529E-05 Max = 7.66043529E-05) regrid = F
GET_2DFLD - surface air temperature, 2006-01-01 03:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 8.00000012E-01 Max = 8.00000012E-01) regrid = F
GET_2DFLD - surface air relative humidity, 2006-01-01 03:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 9.20000000E-01 Max = 9.20000000E-01) regrid = F
GET_2DFLD - rain fall rate, 2006-01-01 03:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 1.17592590E-05 Max = 1.17592590E-05) regrid = F

TIME-STEP YYYY-MM-DD hh:mm:ss.ss KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME Grid
C => (i,j,k) Cu Cv Cw Max Speed

0 2006-01-01 00:00:00.00 0.000000E+00 2.254681E+02 2.254681E+02 6.480027E+11 01
(000,0000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_AVG - creating average file, Grid 01: /users/hwang121/scratch/ngbayhlnestingfile/output/ocean_avg_0001.nc
DEF_STATION - creating stations file, Grid 01: /users/hwang121/scratch/ngbayhlnestingfile/output/ocean_sta.nc
GET_2DFLD - surface u-wind component, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = 0.0000
(Min = -3.84400734E+00 Max = -2.83434620E+00) regrid = T
GET_2DFLD - surface v-wind component, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = 0.0000
(Min = -2.09210197E+00 Max = -4.23612321E-01) regrid = T
GET_2DFLD - surface air pressure, 2006-01-01 03:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 1.00822778E+03 Max = 1.00822778E+03) regrid = F
GET_2DFLD - solar shortwave radiation flux, 2006-01-01 03:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 8.57305843E-13 Max = 8.57305843E-13) regrid = F
GET_2DFLD - downwelling longwave radiation flux, 2006-01-01 03:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 7.66043529E-05 Max = 7.66043529E-05) regrid = F
GET_2DFLD - surface air temperature, 2006-01-01 03:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 8.00000012E-01 Max = 8.00000012E-01) regrid = F
GET_2DFLD - surface air relative humidity, 2006-01-01 03:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 9.20000000E-01 Max = 9.20000000E-01) regrid = F
GET_2DFLD - rain fall rate, 2006-01-01 03:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 1.17592590E-05 Max = 1.17592590E-05) regrid = F
0 2006-01-01 00:00:00.00 0.000000E+00 5.118196E+01 5.118196E+01 4.990846E+09 02
(000,0000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_AVG - creating average file, Grid 02: /users/hwang121/scratch/ngbayhlnestingfile/output/oceanhl_avg_0001.nc
DEF_STATION - creating stations file, Grid 02: /users/hwang121/scratch/ngbayhlnestingfile/output/ocean_sta.nc
GET_2DFLD - surface u-wind component, 2006-01-01 01:00:00.00
(Grid=02, Rec=0000003, Index=1, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = 0.0417
(Min = -6.39545228E+00 Max = -4.19266902E+00) regrid = T
GET_2DFLD - surface v-wind component, 2006-01-01 01:00:00.00
(Grid=02, Rec=0000003, Index=1, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = 0.0417
(Min = -1.92643796E+00 Max = -9.67419191E-01) regrid = T
1 2006-01-01 00:00:01.00 NaN NaN NaN NaN 02
(402,1077,01) 4.900151E-05 1.821114E-05 0.000000E+00 1.578565E-03
Found Error: 01 Line: 302 Source: ROMS/Nonlinear/main3d.F
Found Error: 01 Line: 274 Source: ROMS/Drivers/nl_ocean.h

Blowing-up: Saving latest model state into RESTART file

WRT_RST - wrote re-start fields (Index=1,1) in record = 0000001 01

Blowing-up: Saving latest model state into RESTART file
Top
Kosa
Posts: 35
Joined: Mon Jan 12, 2015 4:12 pm
Location: URI GSO

Re: BLOWING UP during nesting

#2 Unread post by Kosa »

1. Is that your full output file? Where are the lateral boundary conditions and C-preprocessing options etc.?
2. What is your grid refinement factor? This should dictate your DT. E.g. if it is 3x refined, it's typical to use a 1/3 smaller DT for the refined grid (so 3 s instead of 15 s). I don't know whether this would help with your NaNs on the first timestep though.
3. I would start with the simplest possible scenario and add complexity one piece at a time. Analytical initial conditions, no atmospheric forcing, no rivers, no tides, etc.
Top
haili
Posts: 6
Joined: Wed Dec 26, 2018 4:37 pm
Location: Brown University

Re: BLOWING UP during nesting

#3 Unread post by haili »

Kosa wrote:1. Is that your full output file? Where are the lateral boundary conditions and C-preprocessing options etc.?
2. What is your grid refinement factor? This should dictate your DT. E.g. if it is 3x refined, it's typical to use a 1/3 smaller DT for the refined grid (so 3 s instead of 15 s). I don't know whether this would help with your NaNs on the first timestep though.
3. I would start with the simplest possible scenario and add complexity one piece at a time. Analytical initial conditions, no atmospheric forcing, no rivers, no tides, etc.
Thank you for your reply.
1. This is not my full output file, because the limit number of the characters.
2. my grid refinement factor is 3, but I have changed the time step to 1/3 smaller than the coarse one. the result did not change.
3. Yeah, we should start from the simplest one , but the boundary and the forcing files have already succeed for the coarse one. I only added the initial condition, the tidal forcing and the contact files for the refined grid.

the following part is the output file:
## SLURM PROLOG ###############################################################
## Job ID : 2931264
## Job Name : laketest
## Nodelist : node[445-447,453]
## CPUs : 64
## Mem/Node : 30720 MB
## Directory : /gpfs/scratch/hwang121/ngbaynestingngc
## Started : Sun Mar 3 22:49:24 EST 2019
###############################################################################
Laketest case started in parallel
--------------------------------------------------------------------------------
Model Input Parameters: ROMS/TOMS version 3.7
Sunday - March 3, 2019 - 10:49:26 PM
--------------------------------------------------------------------------------

NGBAY nesting grids

Operating system : Linux
CPU/hardware : x86_64
Compiler system : ifort
Compiler command : /gpfs/runtime/opt/mpi/mvapich2-2.3a_intel/bin/mpif90
Compiler flags : -heap-arrays -fp-model precise -ip -O3 -free -free

Input Script : ngbay_case56hl.in

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

Local Root : /gpfs/data/bfoxkemp/hwang121/trunk
Header Dir : /users/hwang121/scratch/ngbaynestingngc
Header file : ngbay.h
Analytical Dir: /users/hwang121/scratch/ngbaynestingngc

Resolution, Grid 01: 0998x1098x015, Parallel Nodes: 64, Tiling: 008x008


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

7200 ntimes Number of timesteps for 3-D equations.
15.000 dt Timestep size (s) for 3-D equations.
15 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.
480 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
360 nSTA Number of timesteps between the writing of data
the stations file.
T ldefout Switch to create a new output NetCDF file(s).
0 nHIS Number of timesteps between the writing fields
into history file.
0 nQCK Number of timesteps between the writing fields
into quicksave file.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
480 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
5760 ndefAVG Number of timesteps between creation of new
time-averaged 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
2.0000E+00 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
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-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.
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.0000E+01 blk_ZQ Height (m) of surface air humidity measurement.
1.0000E+01 blk_ZT Height (m) of surface air temperature measurement.
1.0000E+01 blk_ZW Height (m) of surface winds measurement.
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
1.0000E+00 theta_s S-coordinate surface control parameter.
1.0000E+00 theta_b S-coordinate bottom control parameter.
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.
0.000 dstart Time-stamp assigned to model initialization (days).
0.000 tide_start Reference time origin for tidal forcing (days).
20060101.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.
1.5000E+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.
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 Aout(idFsur) Write out averaged free-surface.
T Aout(idUbar) Write out averaged 2D U-momentum component.
T Aout(idVbar) Write out averaged 2D V-momentum component.
T Aout(idUvel) Write out averaged 3D U-momentum component.
T Aout(idVvel) Write out averaged 3D V-momentum component.
T Aout(idWvel) Write out averaged W-momentum component.
T Aout(idOvel) Write out averaged omega vertical velocity.
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.

Output/Input Files:

Output Restart File: /users/hwang121/scratch/ngbayhlnestingfile/output/ocean_rst.nc
Prefix for Averages Files: /users/hwang121/scratch/ngbayhlnestingfile/output/ocean_avg
Output Stations File: /users/hwang121/scratch/ngbayhlnestingfile/output/ocean_sta.nc
Input Grid File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc
Nesting grid connectivity File: /users/hwang121/scratch/smallinput/ngbay_ngc_2g_ab.nc
Input Nonlinear Initial File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/osom_grid4_mindep_smlp_mod7_ini_hydro.nc
Tidal Forcing File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/osom_grid4_mindep_smlp_mod6_tidalforcing_2006.nc
Input Forcing File 01: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/frc_wind_field_OSOM_2006_new.nc
Input Forcing File 02: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/frc_bulk8_point_osom_2006.nc
Station positions File: /users/hwang121/data/hwang121/stations_osom_nesting.in

Resolution, Grid 02: 0900x1212x015, Parallel Nodes: 64, Tiling: 008x008


Physical Parameters, Grid: 02
=============================

108000 ntimes Number of timesteps for 3-D equations.
1.000 dt Timestep size (s) for 3-D equations.
15 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.
480 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
360 nSTA Number of timesteps between the writing of data
the stations file.
T ldefout Switch to create a new output NetCDF file(s).
0 nHIS Number of timesteps between the writing fields
into history file.
0 nQCK Number of timesteps between the writing fields
into quicksave file.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
480 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
5760 ndefAVG Number of timesteps between creation of new
time-averaged 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
2.0000E+00 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
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-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.
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.0000E+01 blk_ZQ Height (m) of surface air humidity measurement.
1.0000E+01 blk_ZT Height (m) of surface air temperature measurement.
1.0000E+01 blk_ZW Height (m) of surface winds measurement.
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
1.0000E+00 theta_s S-coordinate surface control parameter.
1.0000E+00 theta_b S-coordinate bottom control parameter.
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.
0.000 dstart Time-stamp assigned to model initialization (days).
0.000 tide_start Reference time origin for tidal forcing (days).
20060101.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.
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 Aout(idFsur) Write out averaged free-surface.
T Aout(idUbar) Write out averaged 2D U-momentum component.
T Aout(idVbar) Write out averaged 2D V-momentum component.
T Aout(idUvel) Write out averaged 3D U-momentum component.
T Aout(idVvel) Write out averaged 3D V-momentum component.
T Aout(idWvel) Write out averaged W-momentum component.
T Aout(idOvel) Write out averaged omega vertical velocity.
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.

Output/Input Files:

Output Restart File: /users/hwang121/scratch/ngbayhlnestingfile/output/oceanhl_rst.nc
Prefix for Averages Files: /users/hwang121/scratch/ngbayhlnestingfile/output/oceanhl_avg
Output Stations File: /users/hwang121/scratch/ngbayhlnestingfile/output/oceanhl_sta.nc
Input Grid File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc
Nesting grid connectivity File: /users/hwang121/scratch/smallinput/ngbay_ngc_2g_ab.nc
Input Nonlinear Initial File: /users/hwang121/scratch/smallinput/ngbay_ini.nc
Input Forcing File 01: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/frc_wind_field_OSOM_2006_new.nc
Input Forcing File 02: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/frc_bulk8_point_osom_2006.nc
Station positions File: /users/hwang121/data/hwang121/stations_osom_nesting.in

Tile partition information for Grid 01: 0998x1098x0015 tiling: 008x008

tile Istr Iend Jstr Jend Npts

0 1 124 1 135 251100
1 125 249 1 135 253125
2 250 374 1 135 253125
3 375 499 1 135 253125
4 500 624 1 135 253125
5 625 749 1 135 253125
6 750 874 1 135 253125
7 875 998 1 135 251100
8 1 124 136 273 256680
9 125 249 136 273 258750
10 250 374 136 273 258750
11 375 499 136 273 258750
12 500 624 136 273 258750
13 625 749 136 273 258750
14 750 874 136 273 258750
15 875 998 136 273 256680
16 1 124 274 411 256680
17 125 249 274 411 258750
18 250 374 274 411 258750
19 375 499 274 411 258750
20 500 624 274 411 258750
21 625 749 274 411 258750
22 750 874 274 411 258750
23 875 998 274 411 256680
24 1 124 412 549 256680
25 125 249 412 549 258750
26 250 374 412 549 258750
27 375 499 412 549 258750
28 500 624 412 549 258750
29 625 749 412 549 258750
30 750 874 412 549 258750
31 875 998 412 549 256680
32 1 124 550 687 256680
33 125 249 550 687 258750
34 250 374 550 687 258750
35 375 499 550 687 258750
36 500 624 550 687 258750
37 625 749 550 687 258750
38 750 874 550 687 258750
39 875 998 550 687 256680
40 1 124 688 825 256680
41 125 249 688 825 258750
42 250 374 688 825 258750
43 375 499 688 825 258750
44 500 624 688 825 258750
45 625 749 688 825 258750
46 750 874 688 825 258750
47 875 998 688 825 256680
48 1 124 826 963 256680
49 125 249 826 963 258750
50 250 374 826 963 258750
51 375 499 826 963 258750
52 500 624 826 963 258750
53 625 749 826 963 258750
54 750 874 826 963 258750
55 875 998 826 963 256680
56 1 124 964 1098 251100
57 125 249 964 1098 253125
58 250 374 964 1098 253125
59 375 499 964 1098 253125
60 500 624 964 1098 253125
61 625 749 964 1098 253125
62 750 874 964 1098 253125
63 875 998 964 1098 251100

Tile partition information for Grid 02: 0900x1212x0015 tiling: 008x008

tile Istr Iend Jstr Jend Npts

0 1 111 1 150 249750
1 112 224 1 150 254250
2 225 337 1 150 254250
3 338 450 1 150 254250
4 451 563 1 150 254250
5 564 676 1 150 254250
6 677 789 1 150 254250
7 790 900 1 150 249750
8 1 111 151 302 253080
9 112 224 151 302 257640
10 225 337 151 302 257640
11 338 450 151 302 257640
12 451 563 151 302 257640
13 564 676 151 302 257640
14 677 789 151 302 257640
15 790 900 151 302 253080
16 1 111 303 454 253080
17 112 224 303 454 257640
18 225 337 303 454 257640
19 338 450 303 454 257640
20 451 563 303 454 257640
21 564 676 303 454 257640
22 677 789 303 454 257640
23 790 900 303 454 253080
24 1 111 455 606 253080
25 112 224 455 606 257640
26 225 337 455 606 257640
27 338 450 455 606 257640
28 451 563 455 606 257640
29 564 676 455 606 257640
30 677 789 455 606 257640
31 790 900 455 606 253080
32 1 111 607 758 253080
33 112 224 607 758 257640
34 225 337 607 758 257640
35 338 450 607 758 257640
36 451 563 607 758 257640
37 564 676 607 758 257640
38 677 789 607 758 257640
39 790 900 607 758 253080
40 1 111 759 910 253080
41 112 224 759 910 257640
42 225 337 759 910 257640
43 338 450 759 910 257640
44 451 563 759 910 257640
45 564 676 759 910 257640
46 677 789 759 910 257640
47 790 900 759 910 253080
48 1 111 911 1062 253080
49 112 224 911 1062 257640
50 225 337 911 1062 257640
51 338 450 911 1062 257640
52 451 563 911 1062 257640
53 564 676 911 1062 257640
54 677 789 911 1062 257640
55 790 900 911 1062 253080
56 1 111 1063 1212 249750
57 112 224 1063 1212 254250
58 225 337 1063 1212 254250
59 338 450 1063 1212 254250
60 451 563 1063 1212 254250
61 564 676 1063 1212 254250
62 677 789 1063 1212 254250
63 790 900 1063 1212 249750

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 124.50 0.50 135.50 RHO-points
1 124.50 249.50 0.50 135.50 RHO-points
2 249.50 374.50 0.50 135.50 RHO-points
3 374.50 499.50 0.50 135.50 RHO-points
4 499.50 624.50 0.50 135.50 RHO-points
5 624.50 749.50 0.50 135.50 RHO-points
6 749.50 874.50 0.50 135.50 RHO-points
7 874.50 998.50 0.50 135.50 RHO-points
8 0.50 124.50 135.50 273.50 RHO-points
9 124.50 249.50 135.50 273.50 RHO-points
10 249.50 374.50 135.50 273.50 RHO-points
11 374.50 499.50 135.50 273.50 RHO-points
12 499.50 624.50 135.50 273.50 RHO-points
13 624.50 749.50 135.50 273.50 RHO-points
14 749.50 874.50 135.50 273.50 RHO-points
15 874.50 998.50 135.50 273.50 RHO-points
16 0.50 124.50 273.50 411.50 RHO-points
17 124.50 249.50 273.50 411.50 RHO-points
18 249.50 374.50 273.50 411.50 RHO-points
19 374.50 499.50 273.50 411.50 RHO-points
20 499.50 624.50 273.50 411.50 RHO-points
21 624.50 749.50 273.50 411.50 RHO-points
22 749.50 874.50 273.50 411.50 RHO-points
23 874.50 998.50 273.50 411.50 RHO-points
24 0.50 124.50 411.50 549.50 RHO-points
25 124.50 249.50 411.50 549.50 RHO-points
26 249.50 374.50 411.50 549.50 RHO-points
27 374.50 499.50 411.50 549.50 RHO-points
28 499.50 624.50 411.50 549.50 RHO-points
29 624.50 749.50 411.50 549.50 RHO-points
30 749.50 874.50 411.50 549.50 RHO-points
31 874.50 998.50 411.50 549.50 RHO-points
32 0.50 124.50 549.50 687.50 RHO-points
33 124.50 249.50 549.50 687.50 RHO-points
34 249.50 374.50 549.50 687.50 RHO-points
35 374.50 499.50 549.50 687.50 RHO-points
36 499.50 624.50 549.50 687.50 RHO-points
37 624.50 749.50 549.50 687.50 RHO-points
38 749.50 874.50 549.50 687.50 RHO-points
39 874.50 998.50 549.50 687.50 RHO-points
40 0.50 124.50 687.50 825.50 RHO-points
41 124.50 249.50 687.50 825.50 RHO-points
42 249.50 374.50 687.50 825.50 RHO-points
43 374.50 499.50 687.50 825.50 RHO-points
44 499.50 624.50 687.50 825.50 RHO-points
45 624.50 749.50 687.50 825.50 RHO-points
46 749.50 874.50 687.50 825.50 RHO-points
47 874.50 998.50 687.50 825.50 RHO-points
48 0.50 124.50 825.50 963.50 RHO-points
49 124.50 249.50 825.50 963.50 RHO-points
50 249.50 374.50 825.50 963.50 RHO-points
51 374.50 499.50 825.50 963.50 RHO-points
52 499.50 624.50 825.50 963.50 RHO-points
53 624.50 749.50 825.50 963.50 RHO-points
54 749.50 874.50 825.50 963.50 RHO-points
55 874.50 998.50 825.50 963.50 RHO-points
56 0.50 124.50 963.50 1098.50 RHO-points
57 124.50 249.50 963.50 1098.50 RHO-points
58 249.50 374.50 963.50 1098.50 RHO-points
59 374.50 499.50 963.50 1098.50 RHO-points
60 499.50 624.50 963.50 1098.50 RHO-points
61 624.50 749.50 963.50 1098.50 RHO-points
62 749.50 874.50 963.50 1098.50 RHO-points
63 874.50 998.50 963.50 1098.50 RHO-points

0 1.00 124.50 0.50 135.50 U-points
1 124.50 249.50 0.50 135.50 U-points
2 249.50 374.50 0.50 135.50 U-points
3 374.50 499.50 0.50 135.50 U-points
4 499.50 624.50 0.50 135.50 U-points
5 624.50 749.50 0.50 135.50 U-points
6 749.50 874.50 0.50 135.50 U-points
7 874.50 998.00 0.50 135.50 U-points
8 1.00 124.50 135.50 273.50 U-points
9 124.50 249.50 135.50 273.50 U-points
10 249.50 374.50 135.50 273.50 U-points
11 374.50 499.50 135.50 273.50 U-points
12 499.50 624.50 135.50 273.50 U-points
13 624.50 749.50 135.50 273.50 U-points
14 749.50 874.50 135.50 273.50 U-points
15 874.50 998.00 135.50 273.50 U-points
16 1.00 124.50 273.50 411.50 U-points
17 124.50 249.50 273.50 411.50 U-points
18 249.50 374.50 273.50 411.50 U-points
19 374.50 499.50 273.50 411.50 U-points
20 499.50 624.50 273.50 411.50 U-points
21 624.50 749.50 273.50 411.50 U-points
22 749.50 874.50 273.50 411.50 U-points
23 874.50 998.00 273.50 411.50 U-points
24 1.00 124.50 411.50 549.50 U-points
25 124.50 249.50 411.50 549.50 U-points
26 249.50 374.50 411.50 549.50 U-points
27 374.50 499.50 411.50 549.50 U-points
28 499.50 624.50 411.50 549.50 U-points
29 624.50 749.50 411.50 549.50 U-points
30 749.50 874.50 411.50 549.50 U-points
31 874.50 998.00 411.50 549.50 U-points
32 1.00 124.50 549.50 687.50 U-points
33 124.50 249.50 549.50 687.50 U-points
34 249.50 374.50 549.50 687.50 U-points
35 374.50 499.50 549.50 687.50 U-points
36 499.50 624.50 549.50 687.50 U-points
37 624.50 749.50 549.50 687.50 U-points
38 749.50 874.50 549.50 687.50 U-points
39 874.50 998.00 549.50 687.50 U-points
40 1.00 124.50 687.50 825.50 U-points
41 124.50 249.50 687.50 825.50 U-points
42 249.50 374.50 687.50 825.50 U-points
43 374.50 499.50 687.50 825.50 U-points
44 499.50 624.50 687.50 825.50 U-points
45 624.50 749.50 687.50 825.50 U-points
46 749.50 874.50 687.50 825.50 U-points
47 874.50 998.00 687.50 825.50 U-points
48 1.00 124.50 825.50 963.50 U-points
49 124.50 249.50 825.50 963.50 U-points
50 249.50 374.50 825.50 963.50 U-points
51 374.50 499.50 825.50 963.50 U-points
52 499.50 624.50 825.50 963.50 U-points
53 624.50 749.50 825.50 963.50 U-points
54 749.50 874.50 825.50 963.50 U-points
55 874.50 998.00 825.50 963.50 U-points
56 1.00 124.50 963.50 1098.50 U-points
57 124.50 249.50 963.50 1098.50 U-points
58 249.50 374.50 963.50 1098.50 U-points
59 374.50 499.50 963.50 1098.50 U-points
60 499.50 624.50 963.50 1098.50 U-points
61 624.50 749.50 963.50 1098.50 U-points
62 749.50 874.50 963.50 1098.50 U-points
63 874.50 998.00 963.50 1098.50 U-points

0 0.50 124.50 1.00 135.50 V-points
1 124.50 249.50 1.00 135.50 V-points
2 249.50 374.50 1.00 135.50 V-points
3 374.50 499.50 1.00 135.50 V-points
4 499.50 624.50 1.00 135.50 V-points
5 624.50 749.50 1.00 135.50 V-points
6 749.50 874.50 1.00 135.50 V-points
7 874.50 998.50 1.00 135.50 V-points
8 0.50 124.50 135.50 273.50 V-points
9 124.50 249.50 135.50 273.50 V-points
10 249.50 374.50 135.50 273.50 V-points
11 374.50 499.50 135.50 273.50 V-points
12 499.50 624.50 135.50 273.50 V-points
13 624.50 749.50 135.50 273.50 V-points
14 749.50 874.50 135.50 273.50 V-points
15 874.50 998.50 135.50 273.50 V-points
16 0.50 124.50 273.50 411.50 V-points
17 124.50 249.50 273.50 411.50 V-points
18 249.50 374.50 273.50 411.50 V-points
19 374.50 499.50 273.50 411.50 V-points
20 499.50 624.50 273.50 411.50 V-points
21 624.50 749.50 273.50 411.50 V-points
22 749.50 874.50 273.50 411.50 V-points
23 874.50 998.50 273.50 411.50 V-points
24 0.50 124.50 411.50 549.50 V-points
25 124.50 249.50 411.50 549.50 V-points
26 249.50 374.50 411.50 549.50 V-points
27 374.50 499.50 411.50 549.50 V-points
28 499.50 624.50 411.50 549.50 V-points
29 624.50 749.50 411.50 549.50 V-points
30 749.50 874.50 411.50 549.50 V-points
31 874.50 998.50 411.50 549.50 V-points
32 0.50 124.50 549.50 687.50 V-points
33 124.50 249.50 549.50 687.50 V-points
34 249.50 374.50 549.50 687.50 V-points
35 374.50 499.50 549.50 687.50 V-points
36 499.50 624.50 549.50 687.50 V-points
37 624.50 749.50 549.50 687.50 V-points
38 749.50 874.50 549.50 687.50 V-points
39 874.50 998.50 549.50 687.50 V-points
40 0.50 124.50 687.50 825.50 V-points
41 124.50 249.50 687.50 825.50 V-points
42 249.50 374.50 687.50 825.50 V-points
43 374.50 499.50 687.50 825.50 V-points
44 499.50 624.50 687.50 825.50 V-points
45 624.50 749.50 687.50 825.50 V-points
46 749.50 874.50 687.50 825.50 V-points
47 874.50 998.50 687.50 825.50 V-points
48 0.50 124.50 825.50 963.50 V-points
49 124.50 249.50 825.50 963.50 V-points
50 249.50 374.50 825.50 963.50 V-points
51 374.50 499.50 825.50 963.50 V-points
52 499.50 624.50 825.50 963.50 V-points
53 624.50 749.50 825.50 963.50 V-points
54 749.50 874.50 825.50 963.50 V-points
55 874.50 998.50 825.50 963.50 V-points
56 0.50 124.50 963.50 1098.00 V-points
57 124.50 249.50 963.50 1098.00 V-points
58 249.50 374.50 963.50 1098.00 V-points
59 374.50 499.50 963.50 1098.00 V-points
60 499.50 624.50 963.50 1098.00 V-points
61 624.50 749.50 963.50 1098.00 V-points
62 749.50 874.50 963.50 1098.00 V-points
63 874.50 998.50 963.50 1098.00 V-points

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 111.50 0.50 150.50 RHO-points
1 111.50 224.50 0.50 150.50 RHO-points
2 224.50 337.50 0.50 150.50 RHO-points
3 337.50 450.50 0.50 150.50 RHO-points
4 450.50 563.50 0.50 150.50 RHO-points
5 563.50 676.50 0.50 150.50 RHO-points
6 676.50 789.50 0.50 150.50 RHO-points
7 789.50 900.50 0.50 150.50 RHO-points
8 0.50 111.50 150.50 302.50 RHO-points
9 111.50 224.50 150.50 302.50 RHO-points
10 224.50 337.50 150.50 302.50 RHO-points
11 337.50 450.50 150.50 302.50 RHO-points
12 450.50 563.50 150.50 302.50 RHO-points
13 563.50 676.50 150.50 302.50 RHO-points
14 676.50 789.50 150.50 302.50 RHO-points
15 789.50 900.50 150.50 302.50 RHO-points
16 0.50 111.50 302.50 454.50 RHO-points
17 111.50 224.50 302.50 454.50 RHO-points
18 224.50 337.50 302.50 454.50 RHO-points
19 337.50 450.50 302.50 454.50 RHO-points
20 450.50 563.50 302.50 454.50 RHO-points
21 563.50 676.50 302.50 454.50 RHO-points
22 676.50 789.50 302.50 454.50 RHO-points
23 789.50 900.50 302.50 454.50 RHO-points
24 0.50 111.50 454.50 606.50 RHO-points
25 111.50 224.50 454.50 606.50 RHO-points
26 224.50 337.50 454.50 606.50 RHO-points
27 337.50 450.50 454.50 606.50 RHO-points
28 450.50 563.50 454.50 606.50 RHO-points
29 563.50 676.50 454.50 606.50 RHO-points
30 676.50 789.50 454.50 606.50 RHO-points
31 789.50 900.50 454.50 606.50 RHO-points
32 0.50 111.50 606.50 758.50 RHO-points
33 111.50 224.50 606.50 758.50 RHO-points
34 224.50 337.50 606.50 758.50 RHO-points
35 337.50 450.50 606.50 758.50 RHO-points
36 450.50 563.50 606.50 758.50 RHO-points
37 563.50 676.50 606.50 758.50 RHO-points
38 676.50 789.50 606.50 758.50 RHO-points
39 789.50 900.50 606.50 758.50 RHO-points
40 0.50 111.50 758.50 910.50 RHO-points
41 111.50 224.50 758.50 910.50 RHO-points
42 224.50 337.50 758.50 910.50 RHO-points
43 337.50 450.50 758.50 910.50 RHO-points
44 450.50 563.50 758.50 910.50 RHO-points
45 563.50 676.50 758.50 910.50 RHO-points
46 676.50 789.50 758.50 910.50 RHO-points
47 789.50 900.50 758.50 910.50 RHO-points
48 0.50 111.50 910.50 1062.50 RHO-points
49 111.50 224.50 910.50 1062.50 RHO-points
50 224.50 337.50 910.50 1062.50 RHO-points
51 337.50 450.50 910.50 1062.50 RHO-points
52 450.50 563.50 910.50 1062.50 RHO-points
53 563.50 676.50 910.50 1062.50 RHO-points
54 676.50 789.50 910.50 1062.50 RHO-points
55 789.50 900.50 910.50 1062.50 RHO-points
56 0.50 111.50 1062.50 1212.50 RHO-points
57 111.50 224.50 1062.50 1212.50 RHO-points
58 224.50 337.50 1062.50 1212.50 RHO-points
59 337.50 450.50 1062.50 1212.50 RHO-points
60 450.50 563.50 1062.50 1212.50 RHO-points
61 563.50 676.50 1062.50 1212.50 RHO-points
62 676.50 789.50 1062.50 1212.50 RHO-points
63 789.50 900.50 1062.50 1212.50 RHO-points

0 1.00 111.50 0.50 150.50 U-points
1 111.50 224.50 0.50 150.50 U-points
2 224.50 337.50 0.50 150.50 U-points
3 337.50 450.50 0.50 150.50 U-points
4 450.50 563.50 0.50 150.50 U-points
5 563.50 676.50 0.50 150.50 U-points
6 676.50 789.50 0.50 150.50 U-points
7 789.50 900.00 0.50 150.50 U-points
8 1.00 111.50 150.50 302.50 U-points
9 111.50 224.50 150.50 302.50 U-points
10 224.50 337.50 150.50 302.50 U-points
11 337.50 450.50 150.50 302.50 U-points
12 450.50 563.50 150.50 302.50 U-points
13 563.50 676.50 150.50 302.50 U-points
14 676.50 789.50 150.50 302.50 U-points
15 789.50 900.00 150.50 302.50 U-points
16 1.00 111.50 302.50 454.50 U-points
17 111.50 224.50 302.50 454.50 U-points
18 224.50 337.50 302.50 454.50 U-points
19 337.50 450.50 302.50 454.50 U-points
20 450.50 563.50 302.50 454.50 U-points
21 563.50 676.50 302.50 454.50 U-points
22 676.50 789.50 302.50 454.50 U-points
23 789.50 900.00 302.50 454.50 U-points
24 1.00 111.50 454.50 606.50 U-points
25 111.50 224.50 454.50 606.50 U-points
26 224.50 337.50 454.50 606.50 U-points
27 337.50 450.50 454.50 606.50 U-points
28 450.50 563.50 454.50 606.50 U-points
29 563.50 676.50 454.50 606.50 U-points
30 676.50 789.50 454.50 606.50 U-points
31 789.50 900.00 454.50 606.50 U-points
32 1.00 111.50 606.50 758.50 U-points
33 111.50 224.50 606.50 758.50 U-points
34 224.50 337.50 606.50 758.50 U-points
35 337.50 450.50 606.50 758.50 U-points
36 450.50 563.50 606.50 758.50 U-points
37 563.50 676.50 606.50 758.50 U-points
38 676.50 789.50 606.50 758.50 U-points
39 789.50 900.00 606.50 758.50 U-points
40 1.00 111.50 758.50 910.50 U-points
41 111.50 224.50 758.50 910.50 U-points
42 224.50 337.50 758.50 910.50 U-points
43 337.50 450.50 758.50 910.50 U-points
44 450.50 563.50 758.50 910.50 U-points
45 563.50 676.50 758.50 910.50 U-points
46 676.50 789.50 758.50 910.50 U-points
47 789.50 900.00 758.50 910.50 U-points
48 1.00 111.50 910.50 1062.50 U-points
49 111.50 224.50 910.50 1062.50 U-points
50 224.50 337.50 910.50 1062.50 U-points
51 337.50 450.50 910.50 1062.50 U-points
52 450.50 563.50 910.50 1062.50 U-points
53 563.50 676.50 910.50 1062.50 U-points
54 676.50 789.50 910.50 1062.50 U-points
55 789.50 900.00 910.50 1062.50 U-points
56 1.00 111.50 1062.50 1212.50 U-points
57 111.50 224.50 1062.50 1212.50 U-points
58 224.50 337.50 1062.50 1212.50 U-points
59 337.50 450.50 1062.50 1212.50 U-points
60 450.50 563.50 1062.50 1212.50 U-points
61 563.50 676.50 1062.50 1212.50 U-points
62 676.50 789.50 1062.50 1212.50 U-points
63 789.50 900.00 1062.50 1212.50 U-points

0 0.50 111.50 1.00 150.50 V-points
1 111.50 224.50 1.00 150.50 V-points
2 224.50 337.50 1.00 150.50 V-points
3 337.50 450.50 1.00 150.50 V-points
4 450.50 563.50 1.00 150.50 V-points
5 563.50 676.50 1.00 150.50 V-points
6 676.50 789.50 1.00 150.50 V-points
7 789.50 900.50 1.00 150.50 V-points
8 0.50 111.50 150.50 302.50 V-points
9 111.50 224.50 150.50 302.50 V-points
10 224.50 337.50 150.50 302.50 V-points
11 337.50 450.50 150.50 302.50 V-points
12 450.50 563.50 150.50 302.50 V-points
13 563.50 676.50 150.50 302.50 V-points
14 676.50 789.50 150.50 302.50 V-points
15 789.50 900.50 150.50 302.50 V-points
16 0.50 111.50 302.50 454.50 V-points
17 111.50 224.50 302.50 454.50 V-points
18 224.50 337.50 302.50 454.50 V-points
19 337.50 450.50 302.50 454.50 V-points
20 450.50 563.50 302.50 454.50 V-points
21 563.50 676.50 302.50 454.50 V-points
22 676.50 789.50 302.50 454.50 V-points
23 789.50 900.50 302.50 454.50 V-points
24 0.50 111.50 454.50 606.50 V-points
25 111.50 224.50 454.50 606.50 V-points
26 224.50 337.50 454.50 606.50 V-points
27 337.50 450.50 454.50 606.50 V-points
28 450.50 563.50 454.50 606.50 V-points
29 563.50 676.50 454.50 606.50 V-points
30 676.50 789.50 454.50 606.50 V-points
31 789.50 900.50 454.50 606.50 V-points
32 0.50 111.50 606.50 758.50 V-points
33 111.50 224.50 606.50 758.50 V-points
34 224.50 337.50 606.50 758.50 V-points
35 337.50 450.50 606.50 758.50 V-points
36 450.50 563.50 606.50 758.50 V-points
37 563.50 676.50 606.50 758.50 V-points
38 676.50 789.50 606.50 758.50 V-points
39 789.50 900.50 606.50 758.50 V-points
40 0.50 111.50 758.50 910.50 V-points
41 111.50 224.50 758.50 910.50 V-points
42 224.50 337.50 758.50 910.50 V-points
43 337.50 450.50 758.50 910.50 V-points
44 450.50 563.50 758.50 910.50 V-points
45 563.50 676.50 758.50 910.50 V-points
46 676.50 789.50 758.50 910.50 V-points
47 789.50 900.50 758.50 910.50 V-points
48 0.50 111.50 910.50 1062.50 V-points
49 111.50 224.50 910.50 1062.50 V-points
50 224.50 337.50 910.50 1062.50 V-points
51 337.50 450.50 910.50 1062.50 V-points
52 450.50 563.50 910.50 1062.50 V-points
53 563.50 676.50 910.50 1062.50 V-points
54 676.50 789.50 910.50 1062.50 V-points
55 789.50 900.50 910.50 1062.50 V-points
56 0.50 111.50 1062.50 1212.00 V-points
57 111.50 224.50 1062.50 1212.00 V-points
58 224.50 337.50 1062.50 1212.00 V-points
59 337.50 450.50 1062.50 1212.00 V-points
60 450.50 563.50 1062.50 1212.00 V-points
61 563.50 676.50 1062.50 1212.00 V-points
62 676.50 789.50 1062.50 1212.00 V-points
63 789.50 900.50 1062.50 1212.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 556
HaloSizeJ(1) = 608
TileSide(1) = 146
TileSize(1) = 19418


Maximum halo size in XI and ETA directions:

HaloSizeI(2) = 508
HaloSizeJ(2) = 664
TileSide(2) = 160
TileSize(2) = 19360
Top
haili
Posts: 6
Joined: Wed Dec 26, 2018 4:37 pm
Location: Brown University

Re: BLOWING UP during nesting

#4 Unread post by haili »

Kosa wrote:1. Is that your full output file? Where are the lateral boundary conditions and C-preprocessing options etc.?
2. What is your grid refinement factor? This should dictate your DT. E.g. if it is 3x refined, it's typical to use a 1/3 smaller DT for the refined grid (so 3 s instead of 15 s). I don't know whether this would help with your NaNs on the first timestep though.
3. I would start with the simplest possible scenario and add complexity one piece at a time. Analytical initial conditions, no atmospheric forcing, no rivers, no tides, etc.
The rest part：



Lateral Boundary Conditions: NLM
============================

Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------

zeta 1 Periodic Closed Periodic Closed
2 Nested Nested Nested Nested

ubar 1 Periodic Closed Periodic Closed
2 Nested Nested Nested Nested

vbar 1 Periodic Closed Periodic Closed
2 Nested Nested Nested Nested

u 1 Periodic Closed Periodic Closed
2 Nested Nested Nested Nested

v 1 Periodic Closed Periodic Closed
2 Nested Nested Nested Nested

temp 1 Periodic Closed Periodic Closed
2 Nested Nested Nested Nested

salt 1 Periodic Closed Periodic Closed
2 Nested Nested Nested Nested

tke 1 Periodic Closed Periodic Closed
2 Nested Nested Nested Nested

Activated C-preprocessing Options:

NGBAY NGBAY nesting grids
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.
!ASSEMBLE_ALL... Using mpi_isend/mpi_recv in mp_assemble routine.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged nonlinear model fields.
!BOUNDARY_ALLGATHER Using mpi_allreduce in mp_boundary routine.
BULK_FLUXES Surface bulk fluxes parameterization.
!COLLECT_ALL... Using mpi_isend/mpi_recv in mp_collect routine.
CURVGRID Orthogonal curvilinear grid.
DIFF_GRID Horizontal diffusion coefficient scaled by grid size.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
EMINUSP Compute Salt Flux using E-P.
GLS_MIXING Generic Length-Scale turbulence closure.
LONGWAVE_OUT Compute outgoing longwave radiation internally.
MASKING Land/Sea masking.
MIX_S_TS Mixing of tracers along constant S-surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
MPI MPI distributed-memory configuration.
NESTING Nesting grids: Composite and Refinement.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
NO_WRITE_GRID Not Writing grid arrays into NetCDF ouput files.
ONE_WAY One-way nesting in refinement grids.
PERFECT_RESTART Processing perfect restart variables.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
K_GSCHEME Third-order upstream advection of TKE fields.
RADIATION_2D Use tangential phase speed in radiation conditions.
REDUCE_ALLGATHER Using mpi_allgather in mp_reduce routine.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLVE3D Solving 3D Primitive Equations.
SSH_TIDES Add tidal elevation to SSH climatology.
STATIONS Writing out station data.
TS_MPDATA Recursive flux corrected MPDATA 3D 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_SADVECTION Parabolic splines vertical advection of momentum.
UV_LOGDRAG Logarithmic bottom stress.
UV_TIDES Add tidal currents to 2D momentum climatologies.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
VISC_GRID Horizontal viscosity coefficient scaled by grid size.

Process Information:

Node # 0 (pid= 30031) is active.
Node # 3 (pid= 30034) is active.
Node # 7 (pid= 30038) is active.
Node # 9 (pid= 30040) is active.
Node # 14 (pid= 30045) is active.
Node # 1 (pid= 30032) is active.
Node # 2 (pid= 30033) is active.
Node # 5 (pid= 30036) is active.
Node # 6 (pid= 30037) is active.
Node # 8 (pid= 30039) is active.
Node # 10 (pid= 30041) is active.
Node # 11 (pid= 30042) is active.
Node # 12 (pid= 30043) is active.
Node # 13 (pid= 30044) is active.
Node # 15 (pid= 30046) is active.
Node # 4 (pid= 30035) is active.
Node # 20 (pid= 20243) is active.
Node # 33 (pid= 31531) is active.
Node # 57 (pid= 12425) is active.
Node # 47 (pid= 31545) is active.
Node # 22 (pid= 20245) is active.
Node # 46 (pid= 31544) is active.
Node # 29 (pid= 20252) is active.
Node # 34 (pid= 31532) is active.
Node # 26 (pid= 20249) is active.
Node # 35 (pid= 31533) is active.
Node # 30 (pid= 20253) is active.
Node # 37 (pid= 31535) is active.
Node # 21 (pid= 20244) is active.
Node # 40 (pid= 31538) is active.
Node # 23 (pid= 20246) is active.
Node # 42 (pid= 31540) is active.
Node # 18 (pid= 20241) is active.
Node # 38 (pid= 31536) is active.
Node # 31 (pid= 20254) is active.
Node # 43 (pid= 31541) is active.
Node # 28 (pid= 20251) is active.
Node # 45 (pid= 31543) is active.
Node # 24 (pid= 20247) is active.
Node # 44 (pid= 31542) is active.
Node # 17 (pid= 20240) is active.
Node # 41 (pid= 31539) is active.
Node # 25 (pid= 20248) is active.
Node # 36 (pid= 31534) is active.
Node # 16 (pid= 20239) is active.
Node # 32 (pid= 31530) is active.
Node # 27 (pid= 20250) is active.
Node # 39 (pid= 31537) is active.
Node # 19 (pid= 20242) is active.
Node # 62 (pid= 12430) is active.
Node # 51 (pid= 12419) is active.
Node # 52 (pid= 12420) is active.
Node # 56 (pid= 12424) is active.
Node # 49 (pid= 12417) is active.
Node # 58 (pid= 12426) is active.
Node # 54 (pid= 12422) is active.
Node # 55 (pid= 12423) is active.
Node # 60 (pid= 12428) is active.
Node # 61 (pid= 12429) is active.
Node # 59 (pid= 12427) is active.
Node # 53 (pid= 12421) is active.
Node # 50 (pid= 12418) is active.
Node # 48 (pid= 12416) is active.
Node # 63 (pid= 12431) is active.

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

GET_GRID - angle between XI-axis and EAST: angler
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -5.73923922E-01 Max = 9.08638297E-01)
GET_GRID - reciprocal XI-grid spacing: pm
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 2.29508154E-03 Max = 2.46946401E-02)
GET_GRID - reciprocal ETA-grid spacing: pn
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 2.57345733E-03 Max = 3.66368034E-02)
GET_GRID - XI-derivative of inverse metric factor pn: dndx
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -6.22353037E+00 Max = 1.29465657E+01)
GET_GRID - ETA-derivative of inverse metric factor pm: dmde
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -1.37933839E+00 Max = 2.00305560E+01)
GET_GRID - Coriolis parameter at RHO-points: f
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 9.47308397E-05 Max = 9.79166886E-05)
GET_GRID - bathymetry at RHO-points: h
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 2.00000000E+00 Max = 1.12145151E+02)
GET_GRID - x-location of RHO-points: x-rho
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -8.08048205E+06 Max = -7.78230926E+06)
GET_GRID - y-location of RHO-points: y_rho
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.93464274E+06 Max = 5.18168211E+06)
GET_GRID - x-location of PSI-points: x_psi
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -8.08034114E+06 Max = -7.78241485E+06)
GET_GRID - y-location of PSI-points: y-psi
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.93490009E+06 Max = 5.18150523E+06)
GET_GRID - x-location of U-points: x_u
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -8.08038301E+06 Max = -7.78236657E+06)
GET_GRID - y-location of U-points: y_u
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.93464257E+06 Max = 5.18157942E+06)
GET_GRID - x-location of V-points: x_v
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -8.08044156E+06 Max = -7.78235475E+06)
GET_GRID - y-location of V-points: y_v
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.93490036E+06 Max = 5.18160159E+06)
GET_GRID - longitude of RHO-points: lon_rho
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.93490036E+06 Max = 5.18160159E+06)
GET_GRID - latitude of RHO-points lat_rho
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.93490036E+06 Max = 5.18160159E+06)
GET_GRID - longitude of PSI-points: lon_psi
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -7.26635975E+01 Max = -6.99866847E+01)
GET_GRID - latitude of PSI-points lat_psi
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.05088007E+01 Max = 4.21724345E+01)
GET_GRID - longitude of U-points: lon_u
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -7.26639628E+01 Max = -6.99862328E+01)
GET_GRID - latitude of U-points: lat_u
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.05071487E+01 Max = 4.21728490E+01)
GET_GRID - longitude of V-points: lon_v
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = -7.26650402E+01 Max = -6.99855448E+01)
GET_GRID - latitude of V-points: lat_v
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 4.05088013E+01 Max = 4.21737681E+01)
GET_GRID - mask on RHO-points: mask_rho
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on PSI-points: mask_psi
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on U-points: mask_u
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on V-points: mask_v
(Grid = 01, File: /users/hwang121/roms/ngbay/Projects/ngbayfiles/coarseinput/ngbay_grd.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)

GET_GRID - angle between XI-axis and EAST: angler
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 9.39339461E-02 Max = 4.88527516E-01)
GET_GRID - reciprocal XI-grid spacing: pm
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 2.04210158E-02 Max = 5.53111842E-02)
GET_GRID - reciprocal ETA-grid spacing: pn
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 2.34795268E-02 Max = 4.06173068E-02)
GET_GRID - XI-derivative of inverse metric factor pn: dndx
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -1.49524133E-03 Max = 1.36560555E-02)
GET_GRID - ETA-derivative of inverse metric factor pm: dmde
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -1.72030760E-02 Max = -1.14129901E-02)
GET_GRID - Coriolis parameter at RHO-points: f
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 9.65207017E-05 Max = 9.73545733E-05)
GET_GRID - bathymetry at RHO-points: h
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -1.65801391E+00 Max = 5.26751124E+01)
GET_GRID - x-location of RHO-points: x-rho
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -7.95159584E+06 Max = -7.89327807E+06)
GET_GRID - y-location of RHO-points: y_rho
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 5.08543749E+06 Max = 5.14402907E+06)
GET_GRID - x-location of PSI-points: x_psi
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -7.95157715E+06 Max = -7.89331918E+06)
GET_GRID - y-location of PSI-points: y-psi
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 5.08546640E+06 Max = 5.14400415E+06)
GET_GRID - x-location of U-points: x_u
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -7.95157754E+06 Max = -7.89330716E+06)
GET_GRID - y-location of U-points: y_u
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 5.08544056E+06 Max = 5.14401638E+06)
GET_GRID - x-location of V-points: x_v
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -7.95159545E+06 Max = -7.89329009E+06)
GET_GRID - y-location of V-points: y_v
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 5.08546333E+06 Max = 5.14401684E+06)
GET_GRID - longitude of RHO-points: lon_rho
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 5.08546333E+06 Max = 5.14401684E+06)
GET_GRID - latitude of RHO-points lat_rho
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 5.08546333E+06 Max = 5.14401684E+06)
GET_GRID - longitude of PSI-points: lon_psi
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -7.14921368E+01 Max = -7.10528888E+01)
GET_GRID - latitude of PSI-points lat_psi
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 4.14386890E+01 Max = 4.18767854E+01)
GET_GRID - longitude of U-points: lon_u
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -7.14921482E+01 Max = -7.10527933E+01)
GET_GRID - latitude of U-points: lat_u
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 4.14384986E+01 Max = 4.18768933E+01)
GET_GRID - longitude of V-points: lon_v
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -7.14922510E+01 Max = -7.10526126E+01)
GET_GRID - latitude of V-points: lat_v
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 4.14386717E+01 Max = 4.18768473E+01)
GET_GRID - mask on RHO-points: mask_rho
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on PSI-points: mask_psi
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 0.00000000E+00 Max = 2.00000000E+00)
GET_GRID - mask on U-points: mask_u
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID - mask on V-points: mask_v
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)

Vertical S-coordinate System, Grid 01:

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

15 0.0000000 0.0000000 0.000 0.000 0.000 0.000
14 -0.0666667 -0.0064624 -0.085 -0.110 -0.540 -0.901
13 -0.1333333 -0.0257201 -0.181 -0.239 -1.775 -3.199
12 -0.2000000 -0.0573871 -0.286 -0.386 -3.682 -6.852
11 -0.2666667 -0.1008307 -0.401 -0.551 -6.227 -11.792
10 -0.3333333 -0.1551864 -0.524 -0.733 -9.365 -17.924
9 -0.4000000 -0.2193793 -0.656 -0.929 -13.035 -25.130
8 -0.4666667 -0.2921503 -0.794 -1.138 -17.171 -33.273
7 -0.5333333 -0.3720885 -0.938 -1.358 -21.696 -42.199
6 -0.6000000 -0.4576670 -1.086 -1.587 -26.526 -51.741
5 -0.6666667 -0.5472830 -1.238 -1.821 -31.575 -61.724
4 -0.7333333 -0.6393009 -1.391 -2.059 -36.755 -71.969
3 -0.8000000 -0.7320964 -1.546 -2.298 -41.976 -82.299
2 -0.8666667 -0.8241010 -1.699 -2.536 -47.155 -92.543
1 -0.9333333 -0.9138456 -1.851 -2.771 -52.211 -102.540
0 -1.0000000 -1.0000000 -2.000 -3.000 -57.073 -112.145

Vertical S-coordinate System, Grid 02:

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

15 0.0000000 0.0000000 0.000 0.000 0.000 0.000
14 -0.0666667 -0.0064624 0.234 -0.110 -0.326 -0.511
13 -0.1333333 -0.0257201 0.442 -0.239 -0.945 -1.660
12 -0.2000000 -0.0573871 0.624 -0.386 -1.847 -3.428
11 -0.2666667 -0.1008307 0.782 -0.551 -3.017 -5.782
10 -0.3333333 -0.1551864 0.918 -0.733 -4.437 -8.680
9 -0.4000000 -0.2193793 1.033 -0.929 -6.081 -12.068
8 -0.4666667 -0.2921503 1.131 -1.138 -7.921 -15.884
7 -0.5333333 -0.3720885 1.215 -1.358 -9.924 -20.057
6 -0.6000000 -0.4576670 1.286 -1.587 -12.056 -24.512
5 -0.6666667 -0.5472830 1.350 -1.821 -14.281 -29.167
4 -0.7333333 -0.6393009 1.408 -2.059 -16.560 -33.942
3 -0.8000000 -0.7320964 1.466 -2.298 -18.857 -38.756
2 -0.8666667 -0.8241010 1.524 -2.536 -21.136 -43.530
1 -0.9333333 -0.9138456 1.587 -2.771 -23.363 -48.192
0 -1.0000000 -1.0000000 1.658 -3.000 -25.509 -52.675

Time Splitting Weights for Grid 01: ndtfast = 15 nfast = 22
==================================

Primary Secondary Accumulated to Current Step

1-0.0006740978872509 0.0666666666666667-0.0006740978872509 0.0666666666666667
2-0.0001913161153518 0.0667116065258167-0.0008654140026027 0.1333782731924834
3 0.0014472038369421 0.0667243609335069 0.0005817898343395 0.2001026341259903
4 0.0042370228855824 0.0666278806777107 0.0048188127199218 0.2667305148037010
5 0.0081665994120431 0.0663454124853385 0.0129854121319649 0.3330759272890396
6 0.0132119623624627 0.0658009725245357 0.0261973744944276 0.3988768998135752
7 0.0193298623751124 0.0649201750337048 0.0455272368695400 0.4637970748472801
8 0.0264494009361902 0.0636315175420307 0.0719766378057303 0.5274285923893107
9 0.0344621375639418 0.0618682241462847 0.1064387753696720 0.5892968165355954
10 0.0432106750211074 0.0595707483086885 0.1496494503907794 0.6488675648442840
11 0.0524757225556949 0.0566900366406147 0.2021251729464743 0.7055576014848987
12 0.0619616371700794 0.0531916551369017 0.2640868101165537 0.7587492566218004
13 0.0712804429184289 0.0490608793255631 0.3353672530349826 0.8078101359473635
14 0.0799343282324558 0.0443088497976678 0.4153015812674384 0.8521189857450313
15 0.0872966212754959 0.0389798945821708 0.5025982025429343 0.8910988803272021
16 0.0925912433249121 0.0331601198304711 0.5951894458678464 0.9242590001576731
17 0.0948706401828253 0.0269873702754769 0.6900600860506717 0.9512463704331501
18 0.0929921916151718 0.0206626609299552 0.7830522776658435 0.9719090313631054
19 0.0855930988190858 0.0144631814889438 0.8686453764849293 0.9863722128520491
20 0.0710637499186094 0.0087569749010047 0.9397091264035388 0.9951291877530538
21 0.0475195634887281 0.0040193915730974 0.9872286898922669 0.9991485793261512
22 0.0127713101077333 0.0008514206738489 1.0000000000000002 1.0000000000000000

ndtfast, nfast = 15 22 nfast/ndtfast = 1.46667

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

1.000000000000 1.070450639131 0.535225319565 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.09467

Time Splitting Weights for Grid 02: ndtfast = 15 nfast = 22
==================================

Primary Secondary Accumulated to Current Step

1-0.0006740978872509 0.0666666666666667-0.0006740978872509 0.0666666666666667
2-0.0001913161153518 0.0667116065258167-0.0008654140026027 0.1333782731924834
3 0.0014472038369421 0.0667243609335069 0.0005817898343395 0.2001026341259903
4 0.0042370228855824 0.0666278806777107 0.0048188127199218 0.2667305148037010
5 0.0081665994120431 0.0663454124853385 0.0129854121319649 0.3330759272890396
6 0.0132119623624627 0.0658009725245357 0.0261973744944276 0.3988768998135752
7 0.0193298623751124 0.0649201750337048 0.0455272368695400 0.4637970748472801
8 0.0264494009361902 0.0636315175420307 0.0719766378057303 0.5274285923893107
9 0.0344621375639418 0.0618682241462847 0.1064387753696720 0.5892968165355954
10 0.0432106750211074 0.0595707483086885 0.1496494503907794 0.6488675648442840
11 0.0524757225556949 0.0566900366406147 0.2021251729464743 0.7055576014848987
12 0.0619616371700794 0.0531916551369017 0.2640868101165537 0.7587492566218004
13 0.0712804429184289 0.0490608793255631 0.3353672530349826 0.8078101359473635
14 0.0799343282324558 0.0443088497976678 0.4153015812674384 0.8521189857450313
15 0.0872966212754959 0.0389798945821708 0.5025982025429343 0.8910988803272021
16 0.0925912433249121 0.0331601198304711 0.5951894458678464 0.9242590001576731
17 0.0948706401828253 0.0269873702754769 0.6900600860506717 0.9512463704331501
18 0.0929921916151718 0.0206626609299552 0.7830522776658435 0.9719090313631054
19 0.0855930988190858 0.0144631814889438 0.8686453764849293 0.9863722128520491
20 0.0710637499186094 0.0087569749010047 0.9397091264035388 0.9951291877530538
21 0.0475195634887281 0.0040193915730974 0.9872286898922669 0.9991485793261512
22 0.0127713101077333 0.0008514206738489 1.0000000000000002 1.0000000000000000

ndtfast, nfast = 15 22 nfast/ndtfast = 1.46667

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

1.000000000000 1.070450639131 0.535225319565 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.09467

Metrics information for Grid 01:
===============================

Minimum X-grid spacing, DXmin = 4.04946174E-02 km
Maximum X-grid spacing, DXmax = 3.56382138E-01 km
Minimum Y-grid spacing, DYmin = 2.72949578E-02 km
Maximum Y-grid spacing, DYmax = 3.88582312E-01 km
Minimum Z-grid spacing, DZmin = 8.51699402E-02 m
Maximum Z-grid spacing, DZmax = 1.03302185E+01 m

Minimum barotropic Courant Number = 4.20343636E-02
Maximum barotropic Courant Number = 2.97002616E-01
Maximum Coriolis Courant Number = 1.46019263E-03

Horizontal mixing scaled by grid size, GRDMAX = 3.23731774E-01 km

Metrics information for Grid 02:
===============================

Minimum X-grid spacing, DXmin = 1.00000000E-40 km
Maximum X-grid spacing, DXmax = 4.89691606E-02 km
Minimum Y-grid spacing, DYmin = 1.00000000E-40 km
Maximum Y-grid spacing, DYmax = 4.26169539E-02 km
Minimum Z-grid spacing, DZmin = -2.33860479E-01 m
Maximum Z-grid spacing, DZmax = 9.27954758E+35 m

Minimum barotropic Courant Number = 8.72759390E-03
Maximum barotropic Courant Number = 9.33809402E+54
Maximum Coriolis Courant Number = 1.00000000E+37

Horizontal mixing scaled by grid size, GRDMAX = 4.54778126E-02 km

Refined Nested Grid(s) Information:
==================================

Refined Donor Refined Timestep Refined
Grid Grid Scale Ratio Timesteps

02 01 03 15.00000 15

WARNING: Usually the number of Refined Timesteps must be the same
as the Refined Scale for numerical stability.


NLM: GET_STATE - Reading state initial conditions, 2006-01-01 00:00:00.00
(Grid 01, t = 0.0000, File: osom_grid4_mindep_smlp_mod7_ini_hydro.nc, Rec=0001, Index=1)
- free-surface
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- potential temperature
(Min = 4.56679989E+00 Max = 1.23957848E+01)
- salinity
(Min = 3.09659979E-09 Max = 3.36493083E+01)

NLM: GET_STATE - Reading state initial conditions, 2006-01-01 00:00:00.00
(Grid 02, t = 0.0000, File: ngbay_ini.nc, Rec=0001, Index=1)
- free-surface
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- potential temperature
(Min = 4.55619148E+00 Max = 1.23957848E+01)
- salinity
(Min = 1.04572727E-21 Max = 3.36493083E+01)
GET_NGFLD - tidal period
(Grid = 01, Min = 1.49047214E+04 Max = 9.67260838E+04)
GET_2DFLD - tidal elevation amplitude
(Grid = 01, Min = 0.00000000E+00 Max = 7.54251636E-01)
GET_2DFLD - tidal elevation phase angle
(Grid = 01, Min = 6.26773407E-06 Max = 6.28318498E+00)
GET_2DFLD - tidal current inclination angle
(Grid = 01, Min = 0.00000000E+00 Max = 3.14159244E+00)
GET_2DFLD - tidal current phase angle
(Grid = 01, Min = 0.00000000E+00 Max = 6.28318364E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Grid = 01, Min = 0.00000000E+00 Max = 1.49311554E+00)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Grid = 01, Min = -3.37066252E-01 Max = 2.52296461E-01)
GET_2DFLD - surface u-wind component, 2005-12-31 23:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = -0.0417
(Min = -4.81854172E+00 Max = -1.58071958E+00) regrid = T
GET_2DFLD - surface v-wind component, 2005-12-31 23:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = -0.0417
(Min = -2.42937271E+00 Max = 1.13576238E+00) regrid = T
GET_2DFLD - surface air pressure, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 1.00983752E+03 Max = 1.00983752E+03) regrid = F
GET_2DFLD - solar shortwave radiation flux, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 8.57305843E-13 Max = 8.57305843E-13) regrid = F
GET_2DFLD - downwelling longwave radiation flux, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 7.64574591E-05 Max = 7.64574591E-05) regrid = F
GET_2DFLD - surface air temperature, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 3.75000000E-01 Max = 3.75000000E-01) regrid = F
GET_2DFLD - surface air relative humidity, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 9.24499969E-01 Max = 9.24499969E-01) regrid = F
GET_2DFLD - rain fall rate, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 1.17592590E-05 Max = 1.17592590E-05) regrid = F
GET_2DFLD - surface u-wind component, 2005-12-31 23:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = -0.0417
(Min = -3.84400734E+00 Max = -2.83434620E+00) regrid = T
GET_2DFLD - surface v-wind component, 2005-12-31 23:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = -0.0417
(Min = -2.09210197E+00 Max = -4.23612321E-01) regrid = T
GET_2DFLD - surface air pressure, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 1.00983752E+03 Max = 1.00983752E+03) regrid = F
GET_2DFLD - solar shortwave radiation flux, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 8.57305843E-13 Max = 8.57305843E-13) regrid = F
GET_2DFLD - downwelling longwave radiation flux, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 7.64574591E-05 Max = 7.64574591E-05) regrid = F
GET_2DFLD - surface air temperature, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 3.75000000E-01 Max = 3.75000000E-01) regrid = F
GET_2DFLD - surface air relative humidity, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 9.24499969E-01 Max = 9.24499969E-01) regrid = F
GET_2DFLD - rain fall rate, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.0000
(Min = 1.17592590E-05 Max = 1.17592590E-05) regrid = F

Basin information for Grid 01:

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

Initial domain volumes: TotVolume = 6.2370483442E+11 m3
MinCellVol = 3.4037041132E+02 m3
MaxCellVol = 6.8561716618E+05 m3
Max/Min = 2.0143265789E+03

Basin information for Grid 02:

Maximum grid stiffness ratios: rx0 = 2.592861E-01 (Beckmann and Haidvogel)
rx1 = 6.109278E+00 (Haney)

Initial domain volumes: TotVolume = 3.4468096079E+09 m3
MinCellVol = 4.2938854899E+01 m3
MaxCellVol = 6.6263574929E+03 m3
Max/Min = 1.5432077796E+02

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 000000000001 - 000000007200)
NL ROMS/TOMS: started time-stepping: (Grid: 02 TimeSteps: 000000000001 - 000000108000)

GET_2DFLD - surface u-wind component, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = 0.0000
(Min = -4.81854172E+00 Max = -1.58071958E+00) regrid = T
GET_2DFLD - surface v-wind component, 2006-01-01 00:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = 0.0000
(Min = -2.42937271E+00 Max = 1.13576238E+00) regrid = T
GET_2DFLD - surface air pressure, 2006-01-01 03:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 1.00822778E+03 Max = 1.00822778E+03) regrid = F
GET_2DFLD - solar shortwave radiation flux, 2006-01-01 03:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 8.57305843E-13 Max = 8.57305843E-13) regrid = F
GET_2DFLD - downwelling longwave radiation flux, 2006-01-01 03:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 7.66043529E-05 Max = 7.66043529E-05) regrid = F
GET_2DFLD - surface air temperature, 2006-01-01 03:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 8.00000012E-01 Max = 8.00000012E-01) regrid = F
GET_2DFLD - surface air relative humidity, 2006-01-01 03:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 9.20000000E-01 Max = 9.20000000E-01) regrid = F
GET_2DFLD - rain fall rate, 2006-01-01 03:00:00.00
(Grid=01, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 1.17592590E-05 Max = 1.17592590E-05) regrid = F

TIME-STEP YYYY-MM-DD hh:mm:ss.ss KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME Grid
C => (i,j,k) Cu Cv Cw Max Speed

0 2006-01-01 00:00:00.00 0.000000E+00 2.254681E+02 2.254681E+02 6.480027E+11 01
(000,0000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_AVG - creating average file, Grid 01: /users/hwang121/scratch/ngbayhlnestingfile/output/ocean_avg_0001.nc
DEF_STATION - creating stations file, Grid 01: /users/hwang121/scratch/ngbayhlnestingfile/output/ocean_sta.nc
GET_2DFLD - surface u-wind component, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = 0.0000
(Min = -3.84400734E+00 Max = -2.83434620E+00) regrid = T
GET_2DFLD - surface v-wind component, 2006-01-01 00:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = 0.0000
(Min = -2.09210197E+00 Max = -4.23612321E-01) regrid = T
GET_2DFLD - surface air pressure, 2006-01-01 03:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 1.00822778E+03 Max = 1.00822778E+03) regrid = F
GET_2DFLD - solar shortwave radiation flux, 2006-01-01 03:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 8.57305843E-13 Max = 8.57305843E-13) regrid = F
GET_2DFLD - downwelling longwave radiation flux, 2006-01-01 03:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 7.66043529E-05 Max = 7.66043529E-05) regrid = F
GET_2DFLD - surface air temperature, 2006-01-01 03:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 8.00000012E-01 Max = 8.00000012E-01) regrid = F
GET_2DFLD - surface air relative humidity, 2006-01-01 03:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 9.20000000E-01 Max = 9.20000000E-01) regrid = F
GET_2DFLD - rain fall rate, 2006-01-01 03:00:00.00
(Grid=02, Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000) t = 0.1250
(Min = 1.17592590E-05 Max = 1.17592590E-05) regrid = F
0 2006-01-01 00:00:00.00 0.000000E+00 5.120693E+01 5.120693E+01 4.990846E+09 02
(000,0000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_AVG - creating average file, Grid 02: /users/hwang121/scratch/ngbayhlnestingfile/output/oceanhl_avg_0001.nc
DEF_STATION - creating stations file, Grid 02: /users/hwang121/scratch/ngbayhlnestingfile/output/oceanhl_sta.nc
GET_2DFLD - surface u-wind component, 2006-01-01 01:00:00.00
(Grid=02, Rec=0000003, Index=1, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = 0.0417
(Min = -6.39545228E+00 Max = -4.19266902E+00) regrid = T
GET_2DFLD - surface v-wind component, 2006-01-01 01:00:00.00
(Grid=02, Rec=0000003, Index=1, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000) t = 0.0417
(Min = -1.92643796E+00 Max = -9.67419191E-01) regrid = T
1 2006-01-01 00:00:01.00 NaN NaN NaN NaN 02
(310,0595,01) 2.617588E-03 1.370244E-03 0.000000E+00 NaN
Found Error: 01 Line: 302 Source: ROMS/Nonlinear/main3d.F
Found Error: 01 Line: 274 Source: ROMS/Drivers/nl_ocean.h

Blowing-up: Saving latest model state into RESTART file

WRT_RST - wrote re-start fields (Index=1,1) in record = 0000001 01

Blowing-up: Saving latest model state into RESTART file

Node # 48 CPU: 54.900
Node # 32 CPU: 55.032
Node # 16 CPU: 55.014
Node # 49 CPU: 54.874
Node # 33 CPU: 55.116
Node # 17 CPU: 55.050
WRT_RST - wrote re-start fields (Index=1,2) in record = 0000001 02

Elapsed CPU time (seconds):

Node # 0 CPU: 43.720
Total: 3509.985
Node # 50 CPU: 55.066
Node # 34 CPU: 55.123
Node # 18 CPU: 55.050
Node # 1 CPU: 55.114
Node # 51 CPU: 55.057
Node # 35 CPU: 55.108
Node # 19 CPU: 54.905
Node # 2 CPU: 54.913
Node # 52 CPU: 55.081
Node # 36 CPU: 55.115
Node # 20 CPU: 55.107
Node # 3 CPU: 55.006
Node # 53 CPU: 55.057
Node # 37 CPU: 55.030
Node # 21 CPU: 54.779
Node # 4 CPU: 55.121
Node # 54 CPU: 54.888
Node # 38 CPU: 55.100
Node # 22 CPU: 54.683
Node # 5 CPU: 55.003
Node # 55 CPU: 54.844
Node # 39 CPU: 55.112
Node # 23 CPU: 55.113
Node # 6 CPU: 55.110
Node # 56 CPU: 55.019
Node # 40 CPU: 55.091
Node # 24 CPU: 54.822
Node # 7 CPU: 54.974
Node # 57 CPU: 55.085
Node # 41 CPU: 55.046
Node # 25 CPU: 54.997
Node # 8 CPU: 54.970
Node # 58 CPU: 55.096
Node # 42 CPU: 55.060
Node # 26 CPU: 55.079
Node # 9 CPU: 55.052
Node # 59 CPU: 55.044
Node # 43 CPU: 55.114
Node # 27 CPU: 54.923
Node # 10 CPU: 55.166
Node # 60 CPU: 54.970
Node # 44 CPU: 55.081
Node # 28 CPU: 54.986
Node # 11 CPU: 55.125
Node # 61 CPU: 55.067
Node # 45 CPU: 55.002
Node # 29 CPU: 55.057
Node # 12 CPU: 55.069
Node # 62 CPU: 54.950
Node # 46 CPU: 55.084
Node # 30 CPU: 54.994
Node # 13 CPU: 54.928
Node # 63 CPU: 54.955
Node # 47 CPU: 55.108
Node # 31 CPU: 55.109
Node # 14 CPU: 54.848
Node # 15 CPU: 54.923

Nonlinear model elapsed CPU time profile, Grid: 01

Allocation and array initialization .............. 32.764 ( 0.9334 %)
Ocean state initialization ....................... 553.223 (15.7614 %)
Reading of input data ............................ 611.396 (17.4188 %)
Processing of input data ......................... 3.559 ( 0.1014 %)
Processing of output time averaged data .......... 0.003 ( 0.0001 %)
Computation of vertical boundary conditions ...... 0.170 ( 0.0048 %)
Computation of global information integrals ...... 0.567 ( 0.0161 %)
Writing of output data ........................... 34.367 ( 0.9791 %)
Model 2D kernel .................................. 13.601 ( 0.3875 %)
Tidal forcing .................................... 0.813 ( 0.0232 %)
2D/3D coupling, vertical metrics ................. 5.117 ( 0.1458 %)
Omega vertical velocity .......................... 1.034 ( 0.0295 %)
Equation of state for seawater ................... 4.133 ( 0.1178 %)
Atmosphere-Ocean bulk flux parameterization ...... 3.019 ( 0.0860 %)
GLS vertical mixing parameterization ............. 12.478 ( 0.3555 %)
3D equations right-side terms .................... 2.062 ( 0.0587 %)
3D equations predictor step ...................... 1.942 ( 0.0553 %)
Pressure gradient ................................ 0.839 ( 0.0239 %)
Harmonic mixing of tracers, S-surfaces ........... 0.299 ( 0.0085 %)
Harmonic stress tensor, S-surfaces ............... 0.623 ( 0.0177 %)
Corrector time-step for 3D momentum .............. 2.789 ( 0.0795 %)
Corrector time-step for tracers .................. 16.648 ( 0.4743 %)
Nesting algorithm ................................ 43.169 ( 1.2299 %)
Reading model state vector ....................... 144.918 ( 4.1287 %)
Total: 1489.532 42.4370

Nonlinear model message Passage profile, Grid: 01

Message Passage: 2D halo exchanges ............... 16.447 ( 0.4686 %)
Message Passage: 3D halo exchanges ............... 7.907 ( 0.2253 %)
Message Passage: 4D halo exchanges ............... 8.309 ( 0.2367 %)
Message Passage: data broadcast .................. 967.364 (27.5604 %)
Message Passage: data reduction .................. 4.707 ( 0.1341 %)
Message Passage: data gathering .................. 203.935 ( 5.8101 %)
Message Passage: data scattering.................. 145.144 ( 4.1352 %)
Message Passage: point data gathering ............ 7.135 ( 0.2033 %)
Message Passage: nesting point data gathering .... 462.264 (13.1700 %)
Message Passage: synchronization barrier ......... 0.795 ( 0.0227 %)
Total: 1824.007 51.9662

Nonlinear model elapsed CPU time profile, Grid: 02

Allocation and array initialization .............. 32.764 ( 0.9335 %)
Ocean state initialization ....................... 551.690 (15.7177 %)
Reading of input data ............................ 10.013 ( 0.2853 %)
Processing of input data ......................... 0.341 ( 0.0097 %)
Processing of output time averaged data .......... 0.001 ( 0.0000 %)
Computation of vertical boundary conditions ...... 0.116 ( 0.0033 %)
Computation of global information integrals ...... 0.871 ( 0.0248 %)
Writing of output data ........................... 14.395 ( 0.4101 %)
Model 2D kernel .................................. 13.747 ( 0.3917 %)
Tidal forcing .................................... 0.000 ( 0.0000 %)
2D/3D coupling, vertical metrics ................. 5.327 ( 0.1518 %)
Omega vertical velocity .......................... 1.228 ( 0.0350 %)
Equation of state for seawater ................... 5.938 ( 0.1692 %)
Atmosphere-Ocean bulk flux parameterization ...... 2.577 ( 0.0734 %)
GLS vertical mixing parameterization ............. 12.301 ( 0.3504 %)
3D equations right-side terms .................... 2.082 ( 0.0593 %)
3D equations predictor step ...................... 1.800 ( 0.0513 %)
Pressure gradient ................................ 0.858 ( 0.0244 %)
Harmonic mixing of tracers, S-surfaces ........... 0.305 ( 0.0087 %)
Harmonic stress tensor, S-surfaces ............... 0.646 ( 0.0184 %)
Corrector time-step for 3D momentum .............. 2.408 ( 0.0686 %)
Corrector time-step for tracers .................. 14.788 ( 0.4213 %)
Nesting algorithm ................................ 16.886 ( 0.4811 %)
Reading model state vector ....................... 110.424 ( 3.1460 %)
Total: 801.506 22.8350

Nonlinear model message Passage profile, Grid: 02

Message Passage: 2D halo exchanges ............... 11.532 ( 0.3285 %)
Message Passage: 3D halo exchanges ............... 7.760 ( 0.2211 %)
Message Passage: 4D halo exchanges ............... 6.648 ( 0.1894 %)
Message Passage: data broadcast .................. 358.380 (10.2103 %)
Message Passage: data reduction .................. 1.822 ( 0.0519 %)
Message Passage: data gathering .................. 209.274 ( 5.9623 %)
Message Passage: data scattering.................. 93.113 ( 2.6528 %)
Message Passage: point data gathering ............ 3.059 ( 0.0872 %)
Message Passage: nesting point data gathering .... 648.264 (18.4691 %)
Total: 1339.852 38.1726

Unique code regions profiled ..................... 2291.037 65.2720 %
Residual, non-profiled code ...................... 1218.947 34.7280 %


All percentages are with respect to total time = 3509.985


ROMS/TOMS - Output NetCDF summary for Grid 01:
number of time records written in RESTART file = 00000001
number of time records written in STATION file = 00000001

ROMS/TOMS - Output NetCDF summary for Grid 02:
number of time records written in RESTART file = 00000001
number of time records written in STATION file = 00000001

Analytical header files used:

ROMS/Functionals/ana_btflux.h

ROMS/TOMS: DONE... Sunday - March 3, 2019 - 10:50:22 PM
Laketest case complete
Top
Kosa
Posts: 35
Joined: Mon Jan 12, 2015 4:12 pm
Location: URI GSO

Re: BLOWING UP during nesting

#5 Unread post by Kosa »

haili wrote:I only added the initial condition, the tidal forcing and the contact files for the refined grid.
Can you clarify this sentence? You don't have a tidal forcing file for the refined grid so are you saying that you haven't run the coarse grid standalone with its initial conditions file or its tidal forcing file yet? If so, I would recommend running the coarse grid on its own with all of your initial/forcing files first.
Top
haili
Posts: 6
Joined: Wed Dec 26, 2018 4:37 pm
Location: Brown University

Re: BLOWING UP during nesting

#6 Unread post by haili »

Kosa wrote:
haili wrote:I only added the initial condition, the tidal forcing and the contact files for the refined grid.
Can you clarify this sentence? You don't have a tidal forcing file for the refined grid so are you saying that you haven't run the coarse grid standalone with its initial conditions file or its tidal forcing file yet? If so, I would recommend running the coarse grid on its own with all of your initial/forcing files first.
I have already ran the coarse grid with its own initial/forcing files successfully, and then I add the initial condition, the tidal forcing and the contact files of the refined grid to do the nesting.
Top
jcwarner
Posts: 1200
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: BLOWING UP during nesting

#7 Unread post by jcwarner »

are you using wet/dry for the nest? how shallow do the water depths get?
-j
Top
User avatar
wilkin
Posts: 922
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:
Contact wilkin

Re: BLOWING UP during nesting

#8 Unread post by wilkin »

Metrics information for Grid 02:
===============================

...
Minimum Z-grid spacing, DZmin = -2.33860479E-01 m
Maximum Z-grid spacing, DZmax = 9.27954758E+35 m
You have negative depths.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu
Top
wavewatch
Posts: 12
Joined: Thu Nov 14, 2013 1:51 am
Location: institute of oceanology

Re: BLOWING UP during nesting

#9 Unread post by wavewatch »

I looked through your output file and found that ROMS did not read the boundary conditions. There was no such problem in the old version of ROMS. I would like to ask if anyone has encountered this problem. :?
Top
User avatar
wilkin
Posts: 922
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:
Contact wilkin

Re: BLOWING UP during nesting

#10 Unread post by wilkin »

When updating your source code remember you may need to update the roms.in (ocean.in). New options may have been introduced to manage file input/output that your old set-up does not have. It is not difficult to diff your input script with the trunk version it was based on to highlight possible missing new options. Also, generally speaking the release notes accompanying a major update will alert you to new options you may need to add. I recommend reading those.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu
Top
wavewatch
Posts: 12
Joined: Thu Nov 14, 2013 1:51 am
Location: institute of oceanology

Re: BLOWING UP during nesting

#11 Unread post by wavewatch »

wilkin wrote: ↑Mon Jun 15, 2020 1:09 pm
Dear Dr. Wilkin:
According to your suggestion, I modified the build script file and *.in file of the new version, and recompiled and ran the Two-way nesting examples.The result is the following error:

Code: Select all

 TIME-STEP YYYY-MM-DD hh:mm:ss.ss  KINETIC_ENRG   POTEN_ENRG    TOTAL_ENRG    NET_VOLUME  Grid
                     C => (i,j,k)       Cu            Cv            Cw         Max Speed

         0 2000-01-02 00:00:00.00  1.314949E-03  2.184196E+04  2.184196E+04  1.005743E+17  01
                     (214,052,04)  3.291530E-02  3.888322E-02  0.000000E+00  9.450913E-01
      DEF_HIS     - creating  history      file, Grid 01: ./output/SP_s1_his_0001.nc
      WRT_HIS     - wrote history     fields (Index=1,1) in record = 0000001               01
      DEF_AVG     - creating  average      file, Grid 01: ./output/SP_s1_avg_0001.nc

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
--------------------------------------------------------------------------
A process has executed an operation involving a call to the
"fork()" system call to create a child process.  Open MPI is currently
operating in a condition that could result in memory corruption or
other system errors; your job may hang, crash, or produce silent
data corruption.  The use of fork() (or system() or other calls that
create child processes) is strongly discouraged.

The process that invoked fork was:

  Local host:          [[64185,1],0] (PID 13572)

If you are *absolutely sure* that your application will successfully
and correctly survive a call to fork(), you may disable this warning
by setting the mpi_warn_on_fork MCA parameter to 0.
--------------------------------------------------------------------------
#0  0x2B687FAD66F7
#1  0x2B687FAD6D3E
#2  0x2B688078426F
#3  0x5255BB in __nesting_mod_MOD_nesting
#4  0x46D4ED in main3d_
#5  0x403C5D in __ocean_control_mod_MOD_roms_run
#6  0x403868 in MAIN__ at master.f90:?
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 0 on node comput43 exited on signal 11 (Segmentation fault).
--------------------------------------------------------------------------
But when I define ONE_WAY in the .h file, the mode worked, but blowing up quickly. In both cases, the boundary conditions can be found in Output/Input Files, however, no boundary values are read.
Sorry, There is no way to upload the attachment. I reposted a post with details on it.Looking forward to reply.Thank you
viewtopic.php?f=14&t=5529
Top
johnluick

Re: BLOWING UP during nesting

#12 Unread post by johnluick »

Just wondered if this was what John Wilkin meant to paste in about negative depths.
GET_GRID - bathymetry at RHO-points: h
(Grid = 02, File: /users/hwang121/scratch/smallinput/ngbay_grd_fine1.nc)
(Min = -1.65801391E+00 Max = 5.26751124E+01)
Top
Post Reply

Return to “ROMS Usage”