i have encountered a problem when running ROMS version 3.7 that appears to be related to the
use of the WRITE_WATER option in conjunction with the specification of river inflow (using LuvSrc). i am using the
MASKING option and i specify the river inflow at a v-point that is masked as land (which is adjacent to
a water rho-point). when i run the model with WRITE_WATER undefined, all is well. i see the river inflow in
the vbar array at the right location. but when i define WRITE_WATER, there is an error that occurs when
the history file is first written: "error while writing variable: vbar into history NetCDF file...".
at the bottom of the output file is the following:
ERROR: Abnormal termination: NetCDF OUTPUT.
REASON: NetCDF: Start+count exceeds dimension bound
when i move the river point so that it is at a water v-point, this error does not occur (but the model
eventually blows up at that point as might be expected).
so it seems to me that there is something incompatible with the specification of rivers using LuvSrc and the
WRITE_WATER option.
has anyone else encountered this problem?
any suggestions?
will attach the std output below.
thanks,
dave
Code: Select all
Model Input Parameters: ROMS/TOMS version 3.7
Wednesday - September 7, 2016 - 10:33:58 AM
-----------------------------------------------------------------------------
OSOM-grid4
Operating system : Linux
CPU/hardware : x86_64
Compiler system : ifort
Compiler command : /gpfs/runtime/opt/mvapich2/2.2a-intel/bin/mpif90
Compiler flags : -heap-arrays -fp-model precise -ip -O3 -free -free -free
Input Script :
SVN Root URL : https://www.myroms.org/svn/src/trunk
SVN Revision :
Local Root : /gpfs_home/dullman1/ROMS/trunk
Header Dir : /users/dullman1/projects/OSOM
Header file : osom.h
Analytical Dir: /gpfs_home/dullman1/ROMS/trunk/ROMS/Functionals
Resolution, Grid 01: 0998x1098x015, Parallel Nodes: 100, Tiling: 010x010
Physical Parameters, Grid: 01
=============================
129600 ntimes Number of timesteps for 3-D equations.
20.000 dt Timestep size (s) for 3-D equations.
10 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.
4320 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
45 nSTA Number of timesteps between the writing of data
the stations file.
T ldefout Switch to create a new output NetCDF file(s).
180 nHIS Number of timesteps between the writing fields
into history file.
4320 ndefHIS Number of timesteps between creation of new
history files.
4.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
4.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
4.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-06 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.
1.0000E-03 Zob Bottom roughness (m).
1.0000E-03 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)
1.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
1.0000E+00 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
1.0000E+00 Znudg Nudging/relaxation time scale (days)
for free-surface.
1.0000E+00 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
1.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).
T LuvSrc Turning ON momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
T LtracerSrc(01) Turning ON point Sources/Sinks on tracer 01: temp
T LtracerSrc(02) Turning ON point Sources/Sinks on tracer 02: salt
F LsshCLM Turning OFF processing of SSH climatology.
F Lm2CLM Turning OFF processing of 2D momentum climatology.
F Lm3CLM Turning OFF processing of 3D momentum climatology.
F LtracerCLM(01) Turning OFF processing of climatology tracer 01: temp
F LtracerCLM(02) Turning OFF processing of climatology tracer 02: salt
F LnudgeM2CLM Turning OFF nudging of 2D momentum climatology.
F LnudgeM3CLM Turning OFF nudging of 3D momentum climatology.
F LnudgeTCLM(01) Turning OFF nudging of climatology tracer 01: temp
F LnudgeTCLM(02) Turning OFF nudging of climatology tracer 02: salt
T Hout(idFsur) Write out free-surface.
T Hout(idUbar) Write out 2D U-momentum component.
T Hout(idVbar) Write out 2D V-momentum component.
T Hout(idUvel) Write out 3D U-momentum component.
T Hout(idVvel) Write out 3D V-momentum component.
T Hout(idTvar) Write out tracer 01: temp
T Hout(idTvar) Write out tracer 02: salt
T Hout(idMtke) Write out turbulent kinetic energy.
T Hout(idMtls) Write out turbulent generic length-scale.
Output/Input Files:
Output Restart File: /users/dullman1/data/dullman1/OSOM/test/ocean_rst.nc
Prefix for History Files: /users/dullman1/data/dullman1/OSOM/test/ocean_his
Output Stations File: ocean_sta.nc
Input Grid File: osom_grid4_mindep_smlp_mod6.nc
Input Nonlinear Initial File: osom_grid4_mindep_smlp_mod6_ini_hydro.nc
Input Sources/Sinks File: frc_rivers_osom_2006-2010.nc
Input Forcing File 01: osom_grid4_mindep_smlp_mod6_tidalforcing_2006.nc
Input Forcing File 02: frc_wind_field_OSOM_2006_new.nc
Input Forcing File 03: frc_bulk8_point_osom_2006.nc
Input Boundary File: bry_osom_grid4_mindep_smlp_mod6_2006.nc
Station positions File: stations_osom.in
Lateral Boundary Conditions: NLM
============================
Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------
zeta 1 Chapman Imp Chapman Imp Chapman Imp Closed
ubar 1 Flather Flather Flather Closed
vbar 1 Flather Flather Flather Closed
u 1 Rad + Nud Rad + Nud Rad + Nud Closed
v 1 Rad + Nud Rad + Nud Rad + Nud Closed
temp 1 Rad + Nud Rad + Nud Rad + Nud Closed
salt 1 Rad + Nud Rad + Nud Rad + Nud Closed
tke 1 Radiation Radiation Radiation Closed
Activated C-preprocessing Options:
OSOM OSOM-grid4
ADD_FSOBC Adding tidal elevation to processed OBC data.
ADD_M2OBC Adding tidal currents to processed OBC data.
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ASSUMED_SHAPE Using assumed-shape arrays.
BULK_FLUXES Surface bulk fluxes parameterization.
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.
INLINE_2DIO Processing 3D IO level by level to reduce memory needs.
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.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
NO_WRITE_GRID Not Writing grid arrays into NetCDF ouput files.
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.
RAMP_TIDES Ramping tidal forcing for one day.
!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.
WRITE_WATER Writing data at water points only.
INITIAL: Configuring and initializing forward nonlinear model ...
*******
Vertical S-coordinate System, Grid 01:
level S-coord Cs-curve Z at hmin at hc half way at hmax
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
Time Splitting Weights for Grid 01: ndtfast = 10 nfast = 15
==================================
Primary Secondary Accumulated to Current Step
1 0.0014850648765542 0.1000000000000000 0.0014850648765542 0.1000000000000000
2 0.0059386031399402 0.0998514935123446 0.0074236680164944 0.1998514935123445
3 0.0133457074936657 0.0992576331983506 0.0207693755101601 0.2991091267106951
4 0.0236484051180389 0.0979230624489840 0.0444177806281990 0.3970321891596791
5 0.0366959666818608 0.0955582219371801 0.0811137473100598 0.4925904110968592
6 0.0521753389587947 0.0918886252689940 0.1332890862688545 0.5844790363658532
7 0.0695217010484137 0.0866710913731145 0.2028107873172682 0.6711501277389678
8 0.0878091442019233 0.0797189212682732 0.2906199315191915 0.7508690490072409
9 0.1056214752525633 0.0709380068480808 0.3962414067717548 0.8218070558553218
10 0.1209031436506851 0.0603758593228245 0.5171445504224399 0.8821829151781463
11 0.1307902921035067 0.0482855449577560 0.6479348425259466 0.9304684601359023
12 0.1314219308195445 0.0352065157474053 0.7793567733454911 0.9656749758833076
13 0.1177312353577226 0.0220643226654509 0.8970880087032138 0.9877392985487585
14 0.0832169680811588 0.0102911991296786 0.9803049767843727 0.9980304976784371
15 0.0196950232156274 0.0019695023215627 1.0000000000000000 0.9999999999999999
ndtfast, nfast = 10 15 nfast/ndtfast = 1.50000
Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):
1.000000000000 1.083803072626 0.541901536313 1.000000000000 1.000000000000
Power filter parameters, Fgamma, gamma = 0.28400 0.00000
Metrics information for Grid 01:
===============================
Minimum X-grid spacing, DXmin = 4.04946174E-02 km
Maximum X-grid spacing, DXmax = 4.35714366E-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 = 8.40687272E-02
Maximum barotropic Courant Number = 5.94005231E-01
Maximum Coriolis Courant Number = 1.94692351E-03
Horizontal mixing scaled by grid size, GRDMAX = 3.23731774E-01 km
NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(Grid 01, File: osom_grid4_mindep_smlp_mod6_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.09658034E-09 Max = 3.36493083E+01)
GET_NGFLD - tidal period
(Min = 1.49047214E+04 Max = 9.67260838E+04)
GET_2DFLD - tidal elevation amplitude
(Min = 0.00000000E+00 Max = 7.54251636E-01)
GET_2DFLD - tidal elevation phase angle
(Min = 6.26773407E-06 Max = 6.28318498E+00)
GET_2DFLD - tidal current inclination angle
(Min = 0.00000000E+00 Max = 3.14159244E+00)
GET_2DFLD - tidal current phase angle
(Min = 0.00000000E+00 Max = 6.28318364E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 0.00000000E+00 Max = 1.49311554E+00)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -3.37066252E-01 Max = 2.52296461E-01)
GET_NGFLD - river runoff XI-positions at RHO-points
(Min = 6.36000000E+02 Max = 6.37000000E+02)
GET_NGFLD - river runoff ETA-positions at RHO-points
(Min = 1.02400000E+03 Max = 1.02400000E+03)
GET_NGFLD - river runoff direction
(Min = 1.00000000E+00 Max = 1.00000000E+00)
GET_NGFLD - river runoff mass transport vertical profile
(Min = 6.66666667E-02 Max = 6.66666667E-02)
GET_NGFLD - river runoff mass transport, t = 0 **:**:**
(Rec=0000001, Index=2, File: frc_rivers_osom_2006-2010.nc)
(Tmin= -0.2917 Tmax= 1826.7083)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - river runoff potential temperature, t = 0 **:**:**
(Rec=0000001, Index=2, File: frc_rivers_osom_2006-2010.nc)
(Tmin= -0.2917 Tmax= 1826.7083)
(Min = 4.27218930E-01 Max = 4.27218930E-01)
GET_NGFLD - river runoff salinity, t = 0 **:**:**
(Rec=0000001, Index=2, File: frc_rivers_osom_2006-2010.nc)
(Tmin= -0.2917 Tmax= 1826.7083)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface u-wind component, t = 0 00:00:00
(Rec=0000002, Index=1, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000)
(Min = -4.81854172E+00 Max = -1.58071958E+00)
GET_2DFLD - surface v-wind component, t = 0 00:00:00
(Rec=0000002, Index=1, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000)
(Min = -2.42937271E+00 Max = 1.13576238E+00)
GET_2DFLD - surface air pressure, t = 0 00:00:00
(Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 1.00983752E+03 Max = 1.00983752E+03)
GET_2DFLD - solar shortwave radiation flux, t = 0 00:00:00
(Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 8.57305843E-13 Max = 8.57305843E-13)
GET_2DFLD - downwelling longwave radiation flux, t = 0 00:00:00
(Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 7.64574591E-05 Max = 7.64574591E-05)
GET_2DFLD - surface air temperature, t = 0 00:00:00
(Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 3.75000000E-01 Max = 3.75000000E-01)
GET_2DFLD - surface air relative humidity, t = 0 00:00:00
(Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 9.24499969E-01 Max = 9.24499969E-01)
GET_2DFLD - rain fall rate, t = 0 00:00:00
(Rec=0000001, Index=1, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 1.17592590E-05 Max = 1.17592590E-05)
GET_NGFLD - free-surface western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 3D u-momentum western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 3D v-momentum western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - temperature western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 6.42663527E+00)
GET_NGFLD - salinity western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 2.96997375E+01)
GET_NGFLD - temperature eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 5.38652849E+00)
GET_NGFLD - salinity eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 3.13541832E+01)
GET_NGFLD - temperature southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 1.24026794E+01)
GET_NGFLD - salinity southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 3.36533089E+01)
Basin information for Grid 01:
Maximum grid stiffness ratios: rx0 = 2.000008E-01 (Beckmann and Haidvogel)
rx1 = 5.076469E+00 (Haney)
Initial basin volumes: TotVolume = 6.2370337146E+11 m3
MinVolume = 3.0294087096E+02 m3
MaxVolume = 6.8561716618E+05 m3
Max/Min = 2.2632045785E+03
NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00129600)
GET_NGFLD - river runoff mass transport, t = 0 17:00:00
(Rec=0000002, Index=1, File: frc_rivers_osom_2006-2010.nc)
(Tmin= -0.2917 Tmax= 1826.7083)
(Min = -2.80837539E+00 Max = -2.80837539E+00)
GET_NGFLD - river runoff potential temperature, t = 0 17:00:00
(Rec=0000002, Index=1, File: frc_rivers_osom_2006-2010.nc)
(Tmin= -0.2917 Tmax= 1826.7083)
(Min = 3.59214325E-01 Max = 3.59214325E-01)
GET_NGFLD - river runoff salinity, t = 0 17:00:00
(Rec=0000002, Index=1, File: frc_rivers_osom_2006-2010.nc)
(Tmin= -0.2917 Tmax= 1826.7083)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface u-wind component, t = 0 00:59:59
(Rec=0000003, Index=2, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000)
(Min = -7.11108624E+00 Max = -3.16260144E+00)
GET_2DFLD - surface v-wind component, t = 0 00:59:59
(Rec=0000003, Index=2, File: frc_wind_field_OSOM_2006_new.nc)
(Tmin= -0.0417 Tmax= 365.0000)
(Min = -1.93170730E+00 Max = 2.58984371E+00)
GET_2DFLD - surface air pressure, t = 0 03:00:00
(Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 1.00822778E+03 Max = 1.00822778E+03)
GET_2DFLD - solar shortwave radiation flux, t = 0 03:00:00
(Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 8.57305843E-13 Max = 8.57305843E-13)
GET_2DFLD - downwelling longwave radiation flux, t = 0 03:00:00
(Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 7.66043529E-05 Max = 7.66043529E-05)
GET_2DFLD - surface air temperature, t = 0 03:00:00
(Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 8.00000012E-01 Max = 8.00000012E-01)
GET_2DFLD - surface air relative humidity, t = 0 03:00:00
(Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 9.20000000E-01 Max = 9.20000000E-01)
GET_2DFLD - rain fall rate, t = 0 03:00:00
(Rec=0000002, Index=2, File: frc_bulk8_point_osom_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 1.17592590E-05 Max = 1.17592590E-05)
GET_NGFLD - free-surface western boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -4.23068777E-02 Max = 0.00000000E+00)
GET_NGFLD - free-surface eastern boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -2.07332429E-02 Max = 0.00000000E+00)
GET_NGFLD - free-surface southern boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -1.39146959E-02 Max = 3.52518074E-02)
GET_NGFLD - 2D u-momentum western boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -4.16707098E-02 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum western boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -3.41589265E-02 Max = 4.35377558E-04)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -9.58518870E-03 Max = 4.02900055E-02)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -3.39524075E-02 Max = 5.07995533E-03)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -8.91573448E-03 Max = 5.82007766E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -6.17771409E-02 Max = 2.03768909E-02)
GET_NGFLD - 3D u-momentum western boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -6.19846620E-02 Max = 0.00000000E+00)
GET_NGFLD - 3D v-momentum western boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -5.00817895E-02 Max = 1.04351330E-03)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -1.17071811E-02 Max = 5.05722351E-02)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -4.30298895E-02 Max = 5.63606620E-03)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -2.29379460E-02 Max = 6.88057020E-02)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = -7.85559937E-02 Max = 3.22754867E-02)
GET_NGFLD - temperature western boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 6.45185947E+00)
GET_NGFLD - salinity western boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 2.98226089E+01)
GET_NGFLD - temperature eastern boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 4.65854502E+00)
GET_NGFLD - salinity eastern boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 3.14130192E+01)
GET_NGFLD - temperature southern boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 1.23872366E+01)
GET_NGFLD - salinity southern boundary condition, t = 0 06:00:00
(Rec=0000002, Index=1, File: bry_osom_grid4_mindep_smlp_mod6_2006.nc)
(Tmin= 0.0000 Tmax= 365.0000)
(Min = 0.00000000E+00 Max = 3.36308708E+01)
STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME
C => (i,j,k) Cu Cv Cw Max Speed
0 0 00:00:00 0.000000E+00 2.254689E+02 2.254689E+02 6.480012E+11
(000,0000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS - creating history file, Grid 01: /users/dullman1/data/dullman1/OSOM/test/ocean_his_0001.nc
WRT_HIS - error while writing variable: vbar
into history NetCDF file for time record: 1
Elapsed CPU time (seconds):
Node # 0 CPU: 9.521
Node # 1 CPU: 17.072
Node # 8 CPU: 17.200
Node # 7 CPU: 17.182
Node # 2 CPU: 17.159
Node # 10 CPU: 17.060
Node # 68 CPU: 17.537
Node # 66 CPU: 17.551
Node # 88 CPU: 17.492
Node # 73 CPU: 17.544
Node # 99 CPU: 17.481
Node # 18 CPU: 17.143
Node # 76 CPU: 17.535
Node # 21 CPU: 17.088
Node # 85 CPU: 17.473
Node # 72 CPU: 17.546
Node # 90 CPU: 17.477
Node # 61 CPU: 17.529
Node # 36 CPU: 17.084
Node # 86 CPU: 17.482
Node # 69 CPU: 17.535
Node # 98 CPU: 17.480
Node # 40 CPU: 17.543
Node # 78 CPU: 17.544
Node # 81 CPU: 17.481
Node # 47 CPU: 17.563
Node # 67 CPU: 17.546
Node # 24 CPU: 17.107
Node # 71 CPU: 17.524
Node # 87 CPU: 17.483
Node # 57 CPU: 17.522
Node # 14 CPU: 17.253
Node # 93 CPU: 17.484
Node # 42 CPU: 17.545
Node # 27 CPU: 17.082
Node # 19 CPU: 17.214
Node # 49 CPU: 17.556
Node # 91 CPU: 17.494
Node # 29 CPU: 17.093
Node # 52 CPU: 17.535
Node # 80 CPU: 17.465
Node # 20 CPU: 17.084
Node # 79 CPU: 17.545
Node # 4 CPU: 17.177
Node # 12 CPU: 17.161
Node # 26 CPU: 17.092
Node # 65 CPU: 17.540
Node # 38 CPU: 17.109
Node # 44 CPU: 17.548
Node # 17 CPU: 17.198
Node # 31 CPU: 17.092
Node # 16 CPU: 17.168
Node # 83 CPU: 17.478
Node # 45 CPU: 17.547
Node # 35 CPU: 17.100
Node # 3 CPU: 17.087
Node # 39 CPU: 17.097
Node # 92 CPU: 17.462
Node # 11 CPU: 17.090
Node # 96 CPU: 17.482
Node # 63 CPU: 17.546
Node # 51 CPU: 17.547
Node # 5 CPU: 17.166
Node # 13 CPU: 17.169
Node # 9 CPU: 17.181
Node # 15 CPU: 17.194
Node # 75 CPU: 17.545
Node # 6 CPU: 17.122
Node # 54 CPU: 17.546
Node # 95 CPU: 17.480
Node # 37 CPU: 17.093
Node # 64 CPU: 17.545
Node # 50 CPU: 17.543
Node # 41 CPU: 17.554
Node # 55 CPU: 17.545
Node # 25 CPU: 17.110
Node # 77 CPU: 17.545
Node # 53 CPU: 17.547
Node # 62 CPU: 17.541
Node # 82 CPU: 17.474
Node # 59 CPU: 17.545
Node # 89 CPU: 17.477
Node # 28 CPU: 17.115
Node # 74 CPU: 17.540
Node # 70 CPU: 17.532
Node # 48 CPU: 17.547
Node # 94 CPU: 17.480
Node # 56 CPU: 17.535
Node # 60 CPU: 17.535
Node # 97 CPU: 17.489
Node # 46 CPU: 17.544
Node # 58 CPU: 17.546
Node # 43 CPU: 17.546
Node # 23 CPU: 17.092
Node # 84 CPU: 17.490
Node # 22 CPU: 17.104
Node # 32 CPU: 17.098
Node # 30 CPU: 17.088
Node # 33 CPU: 17.086
Node # 34 CPU: 17.087
Total: 1728.766
Nonlinear model elapsed time profile, Grid: 01
Allocation and array initialization .............. 26.320 ( 1.5225 %)
Ocean state initialization ....................... 9.014 ( 0.5214 %)
Reading of input data ............................ 1156.506 (66.8978 %)
Processing of input data ......................... 1.576 ( 0.0911 %)
Computation of vertical boundary conditions ...... 1.050 ( 0.0607 %)
Computation of global information integrals ...... 2.225 ( 0.1287 %)
Tidal forcing .................................... 5.464 ( 0.3161 %)
2D/3D coupling, vertical metrics ................. 5.089 ( 0.2944 %)
Omega vertical velocity .......................... 1.305 ( 0.0755 %)
Equation of state for seawater ................... 7.016 ( 0.4058 %)
Atmosphere-Ocean bulk flux parameterization ...... 6.634 ( 0.3837 %)
Total: 1222.198 70.6977
Nonlinear model message Passage profile, Grid: 01
Message Passage: 2D halo exchanges ............... 13.104 ( 0.7580 %)
Message Passage: 3D halo exchanges ............... 11.133 ( 0.6440 %)
Message Passage: 4D halo exchanges ............... 1.783 ( 0.1031 %)
Message Passage: data broadcast .................. 1339.021 (77.4553 %)
Message Passage: data reduction .................. 5.131 ( 0.2968 %)
Message Passage: data gathering .................. 4.691 ( 0.2714 %)
Message Passage: data scattering.................. 270.364 (15.6391 %)
Message Passage: boundary data gathering ......... 3.732 ( 0.2159 %)
Message Passage: point data gathering ............ 0.681 ( 0.0394 %)
Total: 1649.641 95.4230
All percentages are with respect to total time = 1728.766
ROMS/TOMS - Output NetCDF summary for Grid 01:
number of time records written in HISTORY file = 00000001
Analytical header files used:
ROMS/Functionals/ana_btflux.h
ROMS/TOMS - Output error ............ exit_flag: 3
ERROR: Abnormal termination: NetCDF OUTPUT.
REASON: NetCDF: Start+count exceeds dimension bound