I am simulating the Beaufort Sea with three open boundaries and ICE_MODEL activated, but I found the water velocities terrible.
There is an eddy in the middle and a stronger current along the east boundary, could it be the problem of boundary conditions? I have tried other boundary conditions for 3D velocities like 'RadNud' and 'Gra' but the calculation blows up with 'Diag speed problem'. Do you have any idea about this situation, any help will be appreciated!
The ROMS version is Kate ROMS version.
The calculation grid is curvilinear grid with same spacing in kilometers.
The forcing data are from NCEP(CFsV2), the initial and boundary conditions are from HYcom.
Code: Select all
Physical Parameters, Grid: 01
=============================
434880 ntimes Number of timesteps for 3-D equations.
30.000 dt Timestep size (s) for 3-D equations.
30 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.
2880 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
T ldefout Switch to create a new output NetCDF file(s).
720 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.
2880 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
1 ntsAVG2 Starting timestep for the accumulation of output
secondary time-averaged data.
120 nAVG2 Number of timesteps between the writing of
secondary time-averaged data into averages file.
2880 ndefAVG2 Number of timesteps between creation of new
secondary time-averaged file.
1.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
1.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
5.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
1.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
1.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)
for tracer 02: salt
1.0000E-05 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
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+00 blk_ZQ Height (m) of surface air humidity measurement.
2.0000E+00 blk_ZT Height (m) of surface air temperature measurement.
2.0000E+00 blk_ZW Height (m) of surface winds measurement.
1 lmd_Jwt Jerlov water type.
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
7.0000E+00 theta_s S-coordinate surface control parameter.
2.0000E+00 theta_b S-coordinate bottom control parameter.
250.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).
-1.000 tide_start Reference time origin for tidal forcing (days).
20150802.00 time_ref Reference time for units attribute (yyyymmdd.dd)
3.6000E+02 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
3.6000E+02 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
9.0000E+01 Tnudg_SSS Nudging/relaxation time scale (days)
for sea surface salinity.
3.6000E+02 Znudg Nudging/relaxation time scale (days)
for free-surface.
3.6000E+02 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
3.6000E+02 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
1.2000E+02 obcfac Factor between passive and active
open boundary conditions.
F VolCons(1) NLM western edge boundary volume conservation.
F VolCons(2) NLM southern edge boundary volume conservation.
F VolCons(3) NLM eastern edge boundary volume conservation.
F VolCons(4) NLM northern edge boundary volume conservation.
10.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
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
F LmiCLM Turning OFF processing of ice momentum climatology.
F LaiCLM Turning OFF processing of ice tracer climatology.
F LsiCLM Turning OFF processing of ice stress climatology.
F LnudgeMICLM Turning OFF nudging of ice momentum climatology.
F LnudgeAICLM Turning OFF nudging of ice tracer climatology.
F LnudgeSICLM Turning OFF nudging of ice stress climatology.
T Hout(idFsur) Write out free-surface.
T Hout(idUbar) Write out 2D U-momentum component.
T Hout(idVbar) Write out 2D V-momentum component.
T Hout(idUvel) Write out 3D U-momentum component.
T Hout(idVvel) Write out 3D V-momentum component.
T Hout(idWvel) Write out W-momentum component.
T Hout(idOvel) Write out omega vertical velocity.
T Hout(idTvar) Write out tracer 01: temp
T Hout(idTvar) Write out tracer 02: salt
T Hout(idTsur) Write out surface net heat flux.
T Hout(idTsur) Write out surface net salt flux.
T Hout(idSrad) Write out shortwave radiation flux.
T Hout(idLrad) Write out longwave radiation flux.
T Hout(idLhea) Write out latent heat flux.
T Hout(idShea) Write out sensible heat flux.
T Hout(idEmPf) Write out E-P flux.
T Hout(idevap) Write out evaporation rate.
T Hout(idrain) Write out rain rate.
T Hout(idDano) Write out density anomaly.
T Hout(idVvis) Write out vertical viscosity: AKv.
T Hout(idHsbl) Write out depth of surface boundary layer.
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(idu2dE) Write out averaged 2D U-eastward at RHO-points.
T Aout(idv2dN) Write out averaged 2D V-northward at RHO-points.
T Aout(idUvel) Write out averaged 3D U-momentum component.
T Aout(idVvel) Write out averaged 3D V-momentum component.
T Aout(idu3dE) Write out averaged 3D U-eastward at RHO-points.
T Aout(idv3dN) Write out averaged 3D V-northward at RHO-points.
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(idUsms) Write out averaged surface U-momentum stress.
T Aout(idVsms) Write out averaged surface V-momentum stress.
T Aout(idUbms) Write out averaged bottom U-momentum stress.
T Aout(idVbms) Write out averaged bottom V-momentum stress.
T Aout(idTsur) Write out averaged surface net heat flux.
T Aout(idTsur) Write out averaged surface net salt flux.
T Aout(idSrad) Write out averaged shortwave radiation flux.
T Aout(idLrad) Write out averaged longwave radiation flux.
T Aout(idLhea) Write out averaged latent heat flux.
T Aout(idShea) Write out averaged sensible heat flux.
T Aout(idevap) Write out averaged evaporation rate.
T Aout(idrain) Write out averaged rain rate.
T Aout(idDano) Write out averaged density anomaly.
T Aout(idHsbl) Write out averaged depth of surface boundary layer.
T Aout2(idFsur) Write out averaged free-surface.
T Aout2(idUbar) Write out averaged 2D U-momentum component.
T Aout2(idVbar) Write out averaged 2D V-momentum component.
T Aout2(idu2dE) Write out averaged 2D U-eastward at RHO-points.
T Aout2(idv2dN) Write out averaged 2D V-northward at RHO-points.
T Aout2(idUvel) Write out averaged 3D U-momentum component.
T Aout2(idVvel) Write out averaged 3D V-momentum component.
T Aout2(idu3dE) Write out averaged 3D U-eastward at RHO-points.
T Aout2(idv3dN) Write out averaged 3D V-northward at RHO-points.
T Aout2(idTvar) Write out averaged tracer 01: temp
T Aout2(idTvar) Write out averaged tracer 02: salt
T Aout2(idUsms) Write out averaged surface U-momentum stress.
T Aout2(idVsms) Write out averaged surface V-momentum stress.
T Aout2(idUair) Write out averaged surface U-wind component.
T Aout2(idVair) Write out averaged surface V-wind component.
T Aout2(idUairE) Write out averaged Eastward surface U-wind component.
T Aout2(idVairN) Write out averaged Northward surface V-wind component.
T Aout2(idTsur) Write out averaged surface net heat flux.
T Aout2(idTsur) Write out averaged surface net salt flux.
T Aout2(idSrad) Write out averaged shortwave radiation flux.
T Aout2(idLrad) Write out averaged longwave radiation flux.
T Aout2(idLhea) Write out averaged latent heat flux.
T Aout2(idShea) Write out averaged sensible heat flux.
1 shuffle NetCDF-4/HDF5 file format shuffle filer flag.
1 deflate NetCDF-4/HDF5 file format deflate filer flag.
1 deflate_level NetCDF-4/HDF5 file format deflate level parameter.
Code: Select all
Ice Parameters, Grid: 01
============================
Ice time step = ocean time step
3.0000E+01 DTICE Ice model time step (s).
1 nstrs Number of iterations for nonlinear ice dynamics.
60 nevp Number of elastic steps per plastic step in EVP.
9.1900E+02 rhoice Density of sea ice (kg/m3).
7.5000E-03 cdiw Ice-water drag coefficient (nondimensional).
3.0000E-03 cdai Air-ice drag coefficient (nondimensional).
1.2670E+00 rho_air Air density (kg/m3).
3.3000E+02 rhosnow_dry Dry snow density (kg/m3).
4.5000E+02 rhosnow_wet Wet snow density (kg/m3).
4.5000E+01 alphai Mohr-Coulomb stress angle (degrees).
0.0000E+00 min_h Minimum average ice thickness (m).
1.0000E-30 min_a Minimum ice concentration (nondimensional).
1.0000E+00 max_a Maximum ice concentration (nondimensional).
0.0000E+00 stressang Turning angle for ice-water drag (degrees).
1.0000E-16 tol Numerical tolerance in rheology calculations .
9.7000E-01 ice_emiss Ice emissivity.
1.0040E+03 spec_heat_air Specific heat of air.
1.7500E-03 trans_coeff Transfer coefficient.
2.8340E+06 sublim_latent_Latent_heat of sublimation.
2.7315E+02 t0deg Zero degrees Celsius in degrees Kelvin.
8.0000E+00 lf_k1 Landfast ice tunable k1.
1.5000E+01 lf_k2 Landfast ice tunable k2.
5.0000E-05 lf_u0 Landfast ice minimal velocity.
T Hout(idUice) Write out U-component ice velocity.
T Hout(idVice) Write out V-component ice velocity.
T Hout(idUiceE) Write out East component ice velocity.
T Hout(idViceN) Write out North component ice velocity.
T Hout(idAice) Write out ice concentration.
T Hout(idHice) Write out average ice thickness.
T Hout(idHsno) Write out snow thickness.
T Hout(idTice) Write out ice/snow surface temperature.
T Hout(idTimid) Write out interior ice temperature.
T Hout(idT2ice) Write out top of ice temperature.
T Hout(idAgeice) Write out ice age.
T Hout(idIomflx) Write out ice-ocean mass flux
T Hout(idSig11) Write out internal ice stress component 11.
T Hout(idSig12) Write out internal ice stress component 12.
T Hout(idSig22) Write out internal ice stress component 22.
T Hout(idT0mk) Write out temperature of molecular sublayer under ice.
T Hout(idS0mk) Write out salinity of molecular sublayer under ice.
T Hout(idWfr) Write out frazil ice growth rate.
T Hout(idWai) Write out ice growth/melt rate.
T Hout(idWao) Write out ice growth/melt rate.
T Hout(idWio) Write out ice growth/melt rate.
T Hout(idWro) Write out ice melt runoff rate.
T Hout(idWdiv) Write out ice divergence rate.
T Hout(idAlbe) Write out ocean surface albedo.
T Hout(idAlbI) Write out ice albedo.
T Hout(idLWdn) Write out downward longwave radiation.
T Hout(idSWdn) Write out downward shortwave radiation.
T Hout(idLatI) Write out latent heat flux over ice.
T Hout(idSenI) Write out sensible heat flux over ice.
T Hout(idLWupI) Write out upward longwave over ice.
T Hout(idSWupI) Write out upward shortwave over ice.
T Hout(idSWinI) Write out upward shortwave over ice.
T Hout(idSWthruI) Write out shortwave radiation through ice.
T Hout(idSaltI) Write out salt flux under ice.
T Hout(idSaltL) Write out salt flux in leads.
T Hout(idQao) Write out heat flux in leads.
T Hout(idQai) Write out heat flux on top of ice.
T Hout(idQio) Write out heat flux into ice.
T Hout(idQi2) Write out heat flux in ice.
T Hout(idSnoI) Write out heat flux into ice.
T Aout(idUice) Write out time-averaged U-component ice velocity.
T Aout(idVice) Write out time-averaged V-component ice velocity.
T Aout(idUiceE) Write out time-averaged East component ice velocity.
T Aout(idViceN) Write out time-averaged North component ice velocity.
T Aout(idAice) Write out time-averaged ice concentration.
T Aout(idHice) Write out time-averaged average ice thickness.
T Aout(idHsno) Write out time-averaged snow thickness.
T Aout(idTice) Write out time-averaged ice/snow surface temperature.
T Aout(idTimid) Write out time-averaged interior ice temperature.
T Aout(idT2ice) Write out time-averaged top of ice temperature.
T Aout(idAgeice) Write out time-averaged ice age.
T Aout(idIomflx) Write out time-averaged ice-ocean mass flux
T Aout(idSig11) Write out time-averaged internal ice stress component 11.
T Aout(idSig12) Write out time-averaged internal ice stress component 12.
T Aout(idSig22) Write out time-averaged internal ice stress component 22.
T Aout(idT0mk) Write out time-averaged under ice temperature.
T Aout(idS0mk) Write out time-averaged under ice salinity.
T Aout(idWfr) Write out time-averaged frazil ice growth rate.
T Aout(idWai) Write out time-averaged ice growth/melt rate.
T Aout(idWao) Write out time-averaged ice growth/melt rate.
T Aout(idWio) Write out time-averaged ice growth/melt rate.
T Aout(idWro) Write out time-averaged ice melt runoff rate.
T Aout(idWdiv) Write out time-averaged ice divergence rate.
T Aout(idAlbe) Write out ocean surface albedo.
T Aout(idAlbI) Write out ice albedo.
T Aout(idLWdn) Write out downward longwave radiation.
T Aout(idSWdn) Write out downward shortwave radiation.
T Aout(idLatI) Write out latent heat flux over ice.
T Aout(idSenI) Write out sensible heat flux over ice.
T Aout(idLWupI) Write out upward longwave over ice.
T Aout(idSWupI) Write out upward shortwave over ice.
T Aout(idSWinI) Write out upward shortwave over ice.
T Aout(idSWthruI) Write out shortwave radiation through ice.
T Aout(idSaltI) Write out salt flux under ice.
T Aout(idSaltL) Write out salt flux in leads.
T Aout(idQao) Write out heat flux in leads.
T Aout(idQai) Write out heat flux on top of ice.
T Aout(idQio) Write out heat flux into ice.
T Aout(idQi2) Write out heat flux in ice.
T Aout(idSnoI) Write out heat flux into ice.
T Aout2(idUice) Write out U-component ice velocity.
T Aout2(idVice) Write out V-component ice velocity.
T Aout2(idUiceE) Write out East component ice velocity.
T Aout2(idViceN) Write out North component ice velocity.
T Aout2(idAice) Write out ice concentration.
T Aout2(idHice) Write out average ice thickness.
T Aout2(idHsno) Write out snow thickness.
T Aout2(idTice) Write out ice/snow surface temperature.
T Aout2(idTimid) Write out interior ice temperature.
T Aout2(idAgeice) Write out ice age.
T Aout2(idSig11) Write out internal ice stress component 11.
T Aout2(idSig12) Write out internal ice stress component 12.
T Aout2(idSig22) Write out internal ice stress component 22.
T Aout2(idT0mk) Write out temperature of molecular sublayer under ice.
T Aout2(idS0mk) Write out salinity of molecular sublayer under ice.
T Qout(idUice) Write out U-component ice velocity.
T Qout(idVice) Write out V-component ice velocity.
T Qout(idUiceE) Write out East component ice velocity.
T Qout(idViceN) Write out North component ice velocity.
T Qout(idAice) Write out ice concentration.
T Qout(idHice) Write out average ice thickness.
T Qout(idHsno) Write out snow thickness.
T Qout(idTice) Write out ice/snow surface temperature.
T Qout(idTimid) Write out interior ice temperature.
T Qout(idAgeice) Write out ice age.
Output/Input Files:
Output Restart File: beaufort_rst.nc
Output History File: beaufort_his.nc
Output Averages File: beaufort_avg.nc
Prefix for Averages Files: beaufort_avg2
Input Grid File: /mnt/home/cyue/roms/bohai/dataarc/beaufort2.nc
Input Nonlinear Initial File: /mnt/home/cyue/roms/bohai/dataarc/initial_beaufort3.nc
Tidal Forcing File: /mnt/home/cyue/roms/bohai/dataarc/beauforttide2.nc
Input Forcing File 01: /mnt/home/cyue/roms/bohai/dataarc/gom_Pair_era.nc
Input Forcing File 02: /mnt/home/cyue/roms/bohai/dataarc/gom_Qair_era.nc
Input Forcing File 03: /mnt/home/cyue/roms/bohai/dataarc/gom_Tair_era.nc
Input Forcing File 04: /mnt/home/cyue/roms/bohai/dataarc/gom_rain_era.nc
Input Forcing File 05: /mnt/home/cyue/roms/bohai/dataarc/gom_swrad_era.nc
Input Forcing File 06: /mnt/home/cyue/roms/bohai/dataarc/gom_windU_era01.nc
Input Forcing File 07: /mnt/home/cyue/roms/bohai/dataarc/gom_windV_era01.nc
Input Forcing File 08: /mnt/home/cyue/roms/bohai/dataarc/gom_lwrad_era.nc
Input Forcing File 09: /mnt/home/cyue/roms/bohai/dataarc/albedobea.nc
Input Forcing File 10: /mnt/home/cyue/roms/bohai/dataarc/sss.nc
Input Boundary File 01: /mnt/home/cyue/roms/bohai/dataarc/boundary-beaufort4.nc
Input Boundary File 02: /mnt/home/cyue/roms/bohai/dataarc/boundary-beaufort-ice.nc
Ice Model Parameters File: /mnt/home/cyue/roms/beau/ice.in
Tile partition information for Grid 01: 199x99x30 tiling: 6x6
Code: Select all
Lateral Boundary Conditions: NLM
============================
Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------
zeta 1 Chapman Exp Closed Chapman Exp Chapman Exp
ubar 1 Shchepetkin Closed Shchepetkin Shchepetkin
vbar 1 Shchepetkin Closed Shchepetkin Shchepetkin
u 1 Clamped Closed Clamped Clamped
v 1 Clamped Closed Clamped Clamped
temp 1 Rad + Nud Closed Rad + Nud Rad + Nud
salt 1 Rad + Nud Closed Rad + Nud Rad + Nud
uice 1 Gradient Closed Gradient Gradient
vice 1 Gradient Closed Gradient Gradient
aice 1 Chapman Exp Closed Chapman Exp Chapman Exp
hice 1 Chapman Exp Closed Chapman Exp Chapman Exp
tisrf 1 Gradient Closed Gradient Gradient
snow_thick 1 Gradient Closed Gradient Gradient
apond 1 Closed Closed Closed Closed
hpond 1 Closed Closed Closed Closed
sig11 1 Gradient Closed Gradient Gradient
sig12 1 Gradient Closed Gradient Gradient
sig22 1 Gradient Closed Gradient Gradient
Activated C-preprocessing Options:
BEAUFORT BEAUFORT Seas
ADD_FSOBC Adding tidal elevation to processed OBC data
ADD_M2OBC Adding tidal currents to processed OBC data
ALBEDO_CSIM Ice albedo computed according to CSIM3
ALBEDO_CURVE Ocean shortwave radiation adjusted by albedo curve
ANA_BSFLUX Analytical kinematic bottom salinity flux
ANA_BTFLUX Analytical kinematic bottom temperature flux
ANA_DRAG_GRID Analytical spatially varying quadratic drag coefficient
ASSUMED_SHAPE Using assumed-shape arrays
AVERAGES Writing out time-averaged nonlinear model fields
AVERAGES2 Writing out secondary time-averaged fields
!BOUNDARY_ALLGATHER Using mpi_allreduce in mp_boundary routine
BULK_FLUXES Surface bulk fluxes parameterization
CCSM_FLUXES Surface bulk fluxes parameterization
!COLLECT_ALL... Using mpi_isend/mpi_recv in mp_collect routine
CORE_FORCING Specific humidity from CORE forcing files
CURVGRID Orthogonal curvilinear grid
DEFLATE Setting compression in output NetCDF-4/HDF5 files
DIFF_GRID Horizontal diffusion coefficient scaled by grid size
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002)
DOUBLE_PRECISION Double precision arithmetic numerical kernel
EMINUSP Compute Salt Flux using E-P
FILTERED Writing out time-filtered output
FILTERED_RST Reading and writing time-filtered restart files
HDF5 Creating NetCDF-4/HDF5 format files
ICE_ADVECT Advection of ice tracers
ICE_BULK_FLUXES Ice bulk fluxes from the atmosphere
ICE_CONVSNOW Conversion of flooded snow to ice
ICE_DIAGS Include extra diagnostic output from sea ice
ICE_EVP Elastic-viscous-plastic ice rheology
ICE_I_O Put shortwave heat into ice interior
ICE_LANDFAST Lemieux shorefast ice
ICE_MK Mellor-Kantha ice thermodynamics
ICE_MODEL Include sea ice model
ICE_MOMENTUM Compute ice momentum equations
ICE_SMOLAR Advect ice tracers with MPDATA scheme
ICE_STRENGTH_QUAD Ice strength a quadratic function of thickness
ICE_THERMO Include ice thermodynamics
LIMIT_BSTRESS Limit bottom stress to maintain bottom velocity direction
LMD_BKPP KPP bottom boundary layer mixing
LMD_CONVEC LMD convective mixing due to shear instability
LMD_DDMIX LMD double-diffusive mixing
LMD_MIXING Large/McWilliams/Doney interior mixing
LMD_NONLOCAL LMD convective nonlocal transport
LMD_RIMIX LMD diffusivity due to shear instability
LMD_SHAPIRO Shapiro filtering boundary layer depth
LMD_SKPP KPP surface boundary layer mixing
LONGWAVE_OUT Compute outgoing longwave radiation internally
MASKING Land/Sea masking
MIX_GEO_TS Mixing of tracers along geopotential surfaces
MIX_S_UV Mixing of momentum along constant S-surfaces
MPI MPI distributed-memory configuration
NONLINEAR Nonlinear Model
NONLIN_EOS Nonlinear Equation of State for seawater
NO_WRITE_GRID Not Writing grid arrays into NetCDF ouput files
POWER_LAW Power-law shape time-averaging barotropic filter
PROFILE Time profiling activated
RADIATION_2D Use tangential phase speed in radiation conditions
REDUCE_ALLGATHER Using mpi_allgather in mp_reduce routine
RI_SPLINES Parabolic Spline Reconstruction for Richardson Number
RST_SINGLE Single precision fields in restart NetCDF file
SALINITY Using salinity
SCORRECTION Surface salinity flux correction
SOLAR_SOURCE Solar Radiation Source Term
SOLVE3D Solving 3D Primitive Equations
SPLINES_VDIFF Parabolic Spline Reconstruction for Vertical Diffusion
SPLINES_VVISC Parabolic Spline Reconstruction for Vertical Viscosity
SSH_TIDES Add tidal elevation to SSH climatology
TIDES_ASTRO Compute astronomical (long-period tide) corrections
TS_U3HADVECTION Third-order upstream horizontal advection of tracers
TS_C4VADVECTION Fourth-order centered vertical advection of tracers
TS_DIF2 Harmonic mixing of tracers
TWO_D_TRACER_SOURCE One value per tracer per source
UV_ADV Advection of momentum
UV_COR Coriolis term
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum
UV_C4VADVECTION Fourth-order centered vertical advection of momentum
UV_DRAG_GRID Spatially varying quadratic drag coefficient
UV_QDRAG Quadratic bottom stress
UV_TIDES Add tidal currents to 2D momentum climatologies
UV_VIS2 Harmonic mixing of momentum
VAR_RHO_2D Variable density barotropic mode
VISC_GRID Horizontal viscosity coefficient scaled by grid size
Code: Select all
Vertical S-coordinate System, Grid 01:
level S-coord Cs-curve Z at hmin at hc half way at hmax
30 0.0000000 0.0000000 0.000 0.000 0.000 0.000
29 -0.0333333 -0.0001156 -0.066 -4.181 -7.569 -8.240
28 -0.0666667 -0.0004685 -0.132 -8.392 -15.541 -17.334
27 -0.1000000 -0.0010779 -0.198 -12.635 -23.948 -27.352
26 -0.1333333 -0.0019766 -0.265 -16.914 -32.847 -38.411
25 -0.1666667 -0.0032130 -0.331 -21.235 -42.318 -50.687
24 -0.2000000 -0.0048536 -0.397 -25.607 -52.476 -64.418
23 -0.2333333 -0.0069862 -0.463 -30.040 -63.469 -79.921
22 -0.2666667 -0.0097241 -0.529 -34.549 -75.490 -97.603
21 -0.3000000 -0.0132122 -0.595 -39.152 -88.784 -117.987
20 -0.3333333 -0.0176335 -0.662 -43.871 -103.663 -141.731
19 -0.3666667 -0.0232170 -0.728 -48.735 -120.514 -169.660
18 -0.4000000 -0.0302484 -0.794 -53.781 -139.824 -202.803
17 -0.4333333 -0.0390807 -0.860 -59.052 -162.190 -242.431
16 -0.4666667 -0.0501474 -0.927 -64.602 -188.350 -290.105
15 -0.5000000 -0.0639753 -0.993 -70.497 -219.198 -347.722
14 -0.5333333 -0.0811979 -1.059 -76.816 -255.808 -417.563
13 -0.5666667 -0.1025651 -1.126 -83.654 -299.455 -502.327
12 -0.6000000 -0.1289464 -1.193 -91.118 -351.613 -605.148
11 -0.6333333 -0.1613228 -1.259 -99.332 -413.949 -729.557
10 -0.6666667 -0.2007569 -1.326 -108.428 -488.266 -879.379
9 -0.7000000 -0.2483317 -1.393 -118.541 -576.404 -1058.516
8 -0.7333333 -0.3050421 -1.460 -129.797 -680.050 -1270.550
7 -0.7666667 -0.3716257 -1.527 -142.287 -800.457 -1518.136
6 -0.8000000 -0.4483212 -1.594 -156.040 -938.031 -1802.134
5 -0.8333333 -0.5345589 -1.662 -170.987 -1091.803 -2120.493
4 -0.8666667 -0.6286208 -1.730 -186.911 -1258.859 -2467.027
3 -0.9000000 -0.7273553 -1.797 -203.419 -1433.846 -2830.386
2 -0.9333333 -0.8261019 -1.865 -219.929 -1608.854 -3193.789
1 -0.9666667 -0.9190220 -1.933 -235.711 -1773.971 -3536.211
0 -1.0000000 -1.0000000 -2.000 -250.000 -1918.815 -3835.630
Time Splitting Weights for Grid 01: ndtfast = 30 nfast = 42
==================================
Primary Secondary Accumulated to Current Step
1-0.0008094437383769 0.0333333333333333-0.0008094437383769 0.0333333333333333
2-0.0014053566728197 0.0333603147912792-0.0022148004111966 0.0666936481246126
3-0.0017877524645903 0.0334071600137066-0.0040025528757869 0.1001008081383191
4-0.0019566842408176 0.0334667517625262-0.0059592371166046 0.1335675599008453
5-0.0019122901320372 0.0335319745705535-0.0078715272486418 0.1670995344713988
6-0.0016548570247459 0.0335957175749547-0.0095263842733877 0.2006952520463536
7-0.0011849025289723 0.0336508794757796-0.0107112868023600 0.2343461315221331
8-0.0005032751608631 0.0336903762267453-0.0112145619632232 0.2680365077488784
9 0.0003887272597151 0.0337071520654408-0.0108258347035081 0.3017436598143192
10 0.0014892209965583 0.0336941944901169-0.0093366137069498 0.3354378543044362
11 0.0027955815694920 0.0336445537902317-0.0065410321374578 0.3690824080946679
12 0.0043042707117221 0.0335513677379153-0.0022367614257357 0.4026337758325831
13 0.0060106451121704 0.0334078920475245 0.0037738836864347 0.4360416678801076
14 0.0079087469427945 0.0332075372104522 0.0116826306292293 0.4692492050905598
15 0.0099910761708920 0.0329439123123590 0.0216737068001212 0.5021931174029188
16 0.0122483446563884 0.0326108764399960 0.0339220514565096 0.5348039938429148
17 0.0146692120341107 0.0322025982847830 0.0485912634906203 0.5670065921276978
18 0.0172400033810439 0.0317136245503127 0.0658312668716642 0.5987202166780105
19 0.0199444086685725 0.0311389577709445 0.0857756755402367 0.6298591744489550
20 0.0227631639997064 0.0304741441486588 0.1085388395399431 0.6603333185976138
21 0.0256737146312910 0.0297153720153352 0.1342125541712341 0.6900486906129490
22 0.0286498597812016 0.0288595815276255 0.1628624139524357 0.7189082721405746
23 0.0316613792205220 0.0279045862015855 0.1945237931729577 0.7468128583421600
24 0.0346736416507075 0.0268492068942347 0.2291974348236651 0.7736620652363948
25 0.0376471948657328 0.0256934188392112 0.2668446296893979 0.7993554840756060
26 0.0405373376992232 0.0244385123436867 0.3073819673886211 0.8237939964192927
27 0.0432936737565711 0.0230872677537126 0.3506756411451922 0.8468812641730054
28 0.0458596469320356 0.0216441452951603 0.3965352880772278 0.8685254094681656
29 0.0481720587108284 0.0201154903974257 0.4447073467880562 0.8886408998655914
30 0.0501605672561820 0.0185097551070648 0.4948679140442381 0.9071506549726561
31 0.0517471682814030 0.0168377361985254 0.5466150823256412 0.9239883911711815
32 0.0528456577069106 0.0151128305891453 0.5994607400325517 0.9391012217603267
33 0.0533610761022577 0.0133513086655816 0.6528218161348094 0.9524525304259084
34 0.0531891349131379 0.0115726061288397 0.7060109510479473 0.9640251365547481
35 0.0522156244733761 0.0097996349650684 0.7582265755213234 0.9738247715198165
36 0.0503158038019031 0.0080591141492892 0.8085423793232265 0.9818838856691057
37 0.0473537721847153 0.0063819206892258 0.8558961515079418 0.9882658063583315
38 0.0431818225418188 0.0048034616164019 0.8990779740497606 0.9930692679747335
39 0.0376397765791564 0.0033640675316746 0.9367177506289170 0.9964333355064081
40 0.0305543017255206 0.0021094083123694 0.9672720523544377 0.9985427438187775
41 0.0217382098544505 0.0010909315881854 0.9890102622088881 0.9996336754069629
42 0.0109897377911119 0.0003663245930371 1.0000000000000000 1.0000000000000000
ndtfast, nfast = 30 42 nfast/ndtfast = 1.40000
Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):
1.000000000000 1.047601458608 0.523800729304 1.000000000000 1.000000000000
Power filter parameters, Fgamma, gamma = 0.28400 0.18933
Metrics information for Grid 01:
===============================
Minimum X-grid spacing, DXmin = 3.65231759E+00 km
Maximum X-grid spacing, DXmax = 4.16564098E+00 km
Minimum Y-grid spacing, DYmin = 1.11176924E+01 km
Maximum Y-grid spacing, DYmax = 1.11911327E+01 km
Minimum Z-grid spacing, DZmin = 6.61394005E-02 m
Maximum Z-grid spacing, DZmax = 3.63402678E+02 m
Minimum barotropic Courant Number = 1.77333713E-03
Maximum barotropic Courant Number = 5.51830481E-02
Maximum Coriolis Courant Number = 4.27841684E-03
Horizontal mixing scaled by grid size, GRDMAX = 6.80603442E+00 km
Maximum grid stiffness ratios: rx0 = 7.132708E-01 (Beckmann and Haidvogel)
rx1 = 2.629503E+01 (Haney)
Initial domain volumes: TotVolume = 1.2559256133E+15 m3
MinCellVol = 6.8905412947E+06 m3
MaxCellVol = 1.5705931091E+10 m3
Max/Min = 2.2793464865E+03
NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 000000000001 - 000000434880)
GET_2DFLD - surface u-wind component, 2015-08-01 02:00:00.00
(Grid=01, Rec=0000002, Index=2, File: gom_windU_era01.nc)
(Tmin= -1.0000 Tmax= 152.0000) t = -0.9167
(Min = -4.76425726E-01 Max = 1.14701714E+00) regrid = T
GET_2DFLD - surface v-wind component, 2015-08-01 02:00:00.00
(Grid=01, Rec=0000002, Index=2, File: gom_windV_era01.nc)
(Tmin= -1.0000 Tmax= 152.0000) t = -0.9167
(Min = -9.72492107E-01 Max = 9.75210293E-01) regrid = T
GET_2DFLD - surface air pressure, 2015-08-01 02:00:00.00
(Grid=01, Rec=0000002, Index=2, File: gom_Pair_era.nc)
(Tmin= -1.0000 Tmax= 152.0000) t = -0.9167
(Min = 8.36178548E+02 Max = 1.01244477E+03) regrid = T
GET_2DFLD - surface air temperature, 2015-08-01 02:00:00.00
(Grid=01, Rec=0000002, Index=2, File: gom_Tair_era.nc)
(Tmin= -1.0000 Tmax= 152.0000) t = -0.9167
(Min = -1.28274624E+00 Max = 1.99322102E+01) regrid = T
GET_2DFLD - surface air relative humidity, 2015-08-01 02:00:00.00
(Grid=01, Rec=0000002, Index=2, File: gom_Qair_era.nc)
(Tmin= -1.0000 Tmax= 152.0000) t = -0.9167
(Min = 2.95494827E-03 Max = 7.15484283E-03) regrid = T
GET_2DFLD - solar shortwave radiation flux, 2015-08-01 02:00:00.00
(Grid=01, Rec=0000002, Index=2, File: gom_swrad_era.nc)
(Tmin= -1.0000 Tmax= 152.0000) t = -0.9167
(Min = 5.21664204E-06 Max = 1.33325287E-04) regrid = T
GET_2DFLD - downwelling longwave radiation flux, 2015-08-01 02:00:00.00
(Grid=01, Rec=0000002, Index=2, File: gom_lwrad_era.nc)
(Tmin= -1.0000 Tmax= 152.0000) t = -0.9167
(Min = 5.65535392E-05 Max = 8.76457447E-05) regrid = T
GET_2DFLD - rain fall rate, 2015-08-01 02:00:00.00
(Grid=01, Rec=0000002, Index=2, File: gom_rain_era.nc)
(Tmin= -1.0000 Tmax= 152.0000) t = -0.9167
(Min = 0.00000000E+00 Max = 0.00000000E+00) regrid = T
GET_2DFLD - sea surface salinity climatology, 2015-08-02 00:00:00.00
(Grid=01, Rec=0000002, Index=2, File: sss.nc)
(Tmin= -1.0000 Tmax= 150.3750) t = 0.0000
(Min = 1.79424982E+01 Max = 3.14597013E+01) regrid = T
GET_NGFLD - free-surface western boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -4.68922317E-01 Max = 0.00000000E+00)
GET_NGFLD - free-surface eastern boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -3.05613518E-01 Max = 0.00000000E+00)
GET_NGFLD - free-surface northern boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -3.00482035E-01 Max = -2.36000016E-01)
GET_NGFLD - 2D u-momentum western boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -1.01783112E-01 Max = 7.59038180E-02)
GET_NGFLD - 2D v-momentum western boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -4.07351963E-02 Max = 9.38245356E-02)
GET_NGFLD - 2D u-momentum eastern boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -1.16973929E-01 Max = 1.15304545E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -2.65568998E-02 Max = 2.80445777E-02)
GET_NGFLD - 2D u-momentum northern boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -4.17499617E-02 Max = 4.32500280E-02)
GET_NGFLD - 2D v-momentum northern boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -1.18500076E-01 Max = -8.34892225E-03)
GET_NGFLD - 3D u-momentum western boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -1.73000008E-01 Max = 9.70000029E-02)
GET_NGFLD - 3D v-momentum western boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -6.16486706E-02 Max = 1.14790864E-01)
GET_NGFLD - 3D u-momentum eastern boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -3.04206371E-01 Max = 1.92375675E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -1.07987799E-01 Max = 6.88706189E-02)
GET_NGFLD - 3D u-momentum northern boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -8.26295018E-02 Max = 1.32801592E-01)
GET_NGFLD - 3D v-momentum northern boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -1.41603172E-01 Max = 1.18677234E-02)
GET_NGFLD - temperature western boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -1.54173160E+00 Max = 1.37320738E+01)
GET_NGFLD - salinity western boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = 2.63740005E+01 Max = 3.49451599E+01)
GET_NGFLD - temperature eastern boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -1.58200562E+00 Max = 7.52923489E+00)
GET_NGFLD - salinity eastern boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = 2.59456940E+01 Max = 3.49460564E+01)
GET_NGFLD - temperature northern boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = -1.61400104E+00 Max = -4.66918051E-01)
GET_NGFLD - salinity northern boundary condition, 2015-08-02 00:00:00.00
(Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
(Tmin= -1.0000 Tmax= 151.0000) t = 0.0000
(Min = 2.86680660E+01 Max = 2.99954510E+01)
