History file output fails using WRITE_WATER and river inflow

Bug reports, work arounds and fixes

Moderators: arango, robertson

Post Reply
Message
Author
dullman
Posts: 3
Joined: Fri Feb 05, 2016 2:51 pm
Location: Graduate School of Oceanography, URI

History file output fails using WRITE_WATER and river inflow

#1 Post by dullman »

Hi,
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

dullman
Posts: 3
Joined: Fri Feb 05, 2016 2:51 pm
Location: Graduate School of Oceanography, URI

Re: History file output fails using WRITE_WATER and river in

#2 Post by dullman »

I delved a bit further into this problem and I think that I have found the source of the problem. I also have what I think is the solution.

The problem:
The run would fail when writing ubar because the program was attempting to write more ubar values than the history netcdf file was dimensioned for. the dimension of ubar is set in def_his.F and makes use of the variable IOBOUNDS(ng)%xy_u which contains the number of water points (at u-points). this value is determined in routine wpoints.F based on the umask array. but when ubar is written from wrt_his.F by a call to subroutine nf_fwrite2d, the mask that is passed to that subroutine is umask_full. when rivers are implemented (using LuvSrc=T), this mask (created in set_masks.F) is identical to umask except that river inflow points are set to 1. so the umask_full array has more points flagged as water than does umask, which is the source of the error i encountered.

The solution:
I modified wpoints.F so that when subroutine wpoints_tile is called, the program passes umask_full and vmask_full to it (as opposed to umask and vmask as in the original code). in this way, the netcdf history file has the correct dimensions (for example, dimension xy_u now is equal to the number of wet u-points including those associated with river inflows.

As far as I can tell, the above solution results in netcdf history files with exactly the same values when reshaped onto the grid as do the history files from a run with WRITE_WATER disabled.

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

Re: History file output fails using WRITE_WATER and river in

#3 Post by arango »

Good catch, thank you :!: I updated the repository. Check the following :arrow: trac ticket for more information.

Post Reply