segmentation fault trying to run realistic app

[das7105]

I got ROMS 3 up and running on a Mac OSX 10.5.6 using gfortran (i'm a new user). The upwelling case ran fine. Now I am trying to set up my own application, but I get segmentation faults that happen before the model even starts running (i.e I can 'make'). I have tried cutting down my grid size, messing with ulimit (although I couldn't get my stack size to umlimited- the max is 32768 kb), but to no avail. The log file looks like this:

Code: Select all

Process Information:

 Thread #  0 (pid=   61278) is active.

 Model Input Parameters:  ROMS/TOMS version 3.2  
                          Tuesday - January 27, 2009 -  2:39:36 PM
 -----------------------------------------------------------------------------

 Dave Test

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

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

 Local Root    : /Users/daves/ROMS
 Header Dir    : /Users/daves/ROMS/ROMS/Include
 Header file   : davetest.h
 Analytical Dir: /Users/daves/ROMS/ROMS/Functionals

 Resolution, Grid 01: 0009x0014x010,  Parallel Threads:  1,  Tiling: 001x001


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

       8640  ntimes          Number of timesteps for 3-D equations.
    300.000  dt              Timestep size (s) for 3-D equations.
         20  ndtfast         Number of timesteps for 2-D equations between
                               each 3D timestep.
          1  ERstr           Starting ensemble/perturbation run number.
          1  ERend           Ending ensemble/perturbation run number.
          0  nrrec           Number of restart records to read from disk.
          T  LcycleRST       Switch to recycle time-records in restart file.
       4320  nRST            Number of timesteps between the writing of data
                               into restart fields.
         10  ninfo           Number of timesteps between print of information
                               to standard output.
          T  ldefout         Switch to create a new output NetCDF file(s).
        864  nHIS            Number of timesteps between the writing fields
                               into history file.
        864  ndefHIS         Number of timesteps between creation of new
                               history files.
 5.0000E-01  tnu2(01)        Horizontal, harmonic mixing coefficient (m2/s)
                               for tracer 01: temp
 0.0000E+00  tnu2(02)        Horizontal, harmonic mixing coefficient (m2/s)
                               for 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.
 2.0000E-02  Zob             Bottom roughness (m).
 2.0000E-02  Zos             Surface roughness (m).
 1.0000E+01  blk_ZQ          Height (m) of surface air humidity measurement.
 1.0000E+01  blk_ZT          Height (m) of surface air temperature measurement.
 1.0000E+01  blk_ZW          Height (m) of surface winds measurement.
          5  lmd_Jwt         Jerlov water type.
 5.0000E+00  theta_s         S-coordinate surface control parameter.
 8.0000E-01  theta_b         S-coordinate bottom  control parameter.
      0.000  Tcline          S-coordinate surface/bottom layer width (m) used
                               in vertical coordinate stretching.
   1023.700  rho0            Mean density (kg/m3) for Boussinesq approximation.
    120.000  dstart          Time-stamp assigned to model initialization (days).
  -5114.000  tide_start      Reference time origin for tidal forcing (days).
20060101.00  time_ref        Reference time for units attribute (yyyymmdd.dd)
 3.0000E+00  Tnudg(01)       Nudging/relaxation time scale (days)
                               for tracer 01: temp
 3.0000E+00  Tnudg(02)       Nudging/relaxation time scale (days)
                               for tracer 02: salt
 3.0000E+00  Znudg           Nudging/relaxation time scale (days)
                               for free-surface.
 3.0000E+00  M2nudg          Nudging/relaxation time scale (days)
                               for 2D momentum.
 3.0000E+00  M3nudg          Nudging/relaxation time scale (days)
                               for 3D momentum.
 4.0000E+00  obcfac          Factor between passive and active
                               open boundary conditions.
      4.000  T0              Background potential temperature (C) constant.
     32.000  S0              Background salinity (PSU) constant.
      1.000  gamma2          Slipperiness variable: free-slip (1.0) or 
                                                    no-slip (-1.0).
          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(idTvar)    Write out tracer 01: temp
          T  Hout(idTvar)    Write out tracer 02: salt
          T  Hout(idUsms)    Write out surface U-momentum stress.
          T  Hout(idVsms)    Write out surface V-momentum stress.
          T  Hout(idUbms)    Write out bottom U-momentum stress.
          T  Hout(idVbms)    Write out bottom V-momentum stress.
          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(idDano)    Write out density anomaly.
          T  Hout(idVvis)    Write out vertical viscosity coefficient.
          T  Hout(idSdif)    Write out vertical S-diffusion coefficient.
          T  Hout(idMtke)    Write out turbulent kinetic energy.
          T  Hout(idMtls)    Write out turbulent generic length-scale.

 Output/Input Files:

             Output Restart File:  OUT/ocean_rst.nc
        Prefix for History Files:  OUT/ocean_his
                 Input Grid File:  grid.nc
    Input Nonlinear Initial File:  Ocn/ocean_ini.nc
           Input Forcing File 01:  Ocn/tides.nc
           Input Forcing File 02:  Atm/Uwind.nc
           Input Forcing File 03:  Atm/lwrad_down.nc
           Input Forcing File 04:  Atm/Pair.nc
           Input Forcing File 05:  Atm/rain.nc
           Input Forcing File 06:  Atm/swrad.nc
           Input Forcing File 07:  Atm/Tair.nc
           Input Forcing File 08:  Atm/Qair.nc
           Input Forcing File 09:  Atm/Vwind.nc
          Input Climatology File:  Ocn/ocean_clm.nc
             Input Boundary File:  Ocn/ocean_bry.nc

 Tile partition information for Grid 01:  0009x0014x0010  tiling: 001x001

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1        9        1       14     1260

 Tile minimum and maximum fractional grid coordinates:

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     1.00     9.00     1.00    14.00  RHO-points

        0     1.50     8.50     1.00    14.00    U-points

        0     1.00     9.00     1.50    13.50    V-points

 Activated C-preprocessing Options:

  DAVETEST           Dave Test
  ANA_BSFLUX         Analytical kinematic bottom salinity flux.
  ANA_BTFLUX         Analytical kinematic bottom temperature flux.
  ANA_RAIN           Analytical rain fall rate.
  ANA_SSFLUX         Analytical kinematic surface salinity flux.
  ASSUMED_SHAPE      Using assumed-shape arrays.
  BULK_FLUXES        Surface bulk fluxes parametererization.
  CANUTO_A           Canuto A-stability function formulation.
  DJ_GRADPS          Parabolic Splines density Jacobian (Shchepetkin, 2002).
  DOUBLE_PRECISION   Double precision arithmetic.
  EAST_FSCHAPMAN     Eastern edge, free-surface, Chapman condition.
  EAST_M2FLATHER     Eastern edge, 2D momentum, Flather condition.
  EAST_M3NUDGING     Eastern edge, 3D momentum, passive/active outflow/inflow.
  EAST_M3RADIATION   Eastern edge, 3D momentum, radiation condition.
  EAST_TNUDGING      Eastern edge, tracers, passive/active outflow/inflow.
  EAST_TRADIATION    Eastern edge, tracers, radiation condition.
  GLS_MIXING         Generic Length-Scale turbulence closure.
  LONGWAVE_OUT       Compute outgoing longwave radiation internally.
  M2CLIMATOLOGY      Processing 2D momentum climatology data.
  M2CLM_NUDGING      Nudging toward 2D momentum climatology.
  M3CLIMATOLOGY      Processing 3D momentum climatology data.
  M3CLM_NUDGING      Nudging toward 3D momentum climatology.
  MASKING            Land/Sea masking.
  MIX_GEO_TS         Mixing of tracers along geopotential surfaces.
  NONLINEAR          Nonlinear Model.
  NONLIN_EOS         Nonlinear Equation of State for seawater.
  NORTH_FSCHAPMAN    Northern edge, free-surface, Chapman condition.
  NORTH_M2FLATHER    Northern edge, 2D momentum, Flather condition.
  NORTH_M3NUDGING    Northern edge, 3D momentum, passive/active outflow/inflow.
  NORTH_M3RADIATION  Northern edge, 3D momentum, radiation condition.
  NORTH_TNUDGING     Northern edge, tracers, passive/active outflow/inflow.
  NORTH_TRADIATION   Northern edge, tracers, radiation condition.
  N2S2_HORAVG        Horizontal smoothing of buoyancy and shear.
  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.
  !RST_SINGLE        Double precision fields in restart NetCDF file.
  SALINITY           Using salinity.
  SOLAR_SOURCE       Solar Radiation Source Term.
  SOLVE3D            Solving 3D Primitive Equations.
  SOUTH_FSCHAPMAN    Southern edge, free-surface, Chapman condition.
  SOUTH_M2FLATHER    Southern edge, 2D momentum, Flather condition.
  SOUTH_M3NUDGING    Southern edge, 3D momentum, passive/active outflow/inflow.
  SOUTH_M3RADIATION  Southern edge, 3D momentum, radiation condition.
  SOUTH_TNUDGING     Southern edge, tracers, passive/active outflow/inflow.
  SOUTH_TRADIATION   Southern edge, tracers, radiation condition.
  SPLINES            Conservative parabolic spline reconstruction.
  SPHERICAL          Spherical grid configuration.
  SSH_TIDES          Add tidal elevation to SSH climatology.
  TCLIMATOLOGY       Processing tracer climatology data.
  TCLM_NUDGING       Nudging toward tracer climatology.
  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_C4VADVECTION    Fourth-order centered vertical advection of momentum.
  UV_QDRAG           Quadratic bottom stress.
  UV_TIDES           Add tidal currents to 2D momentum climatologies.
  VAR_RHO_2D         Variable density barotropic mode.
  WEST_FSCHAPMAN     Western edge, free-surface, Chapman condition.
  WEST_M2FLATHER     Western edge, 2D momentum, Flather condition.
  WEST_M3NUDGING     Western edge, 3D momentum, passive/active outflow/inflow.
  WEST_M3RADIATION   Western edge, 3D momentum, radiation condition.
  WEST_TNUDGING      Western edge, tracers, passive/active outflow/inflow.
  WEST_TRADIATION    Western edge, tracers, radiation condition.
  ZCLIMATOLOGY       Processing sea surface height climatology data.
  ZCLM_NUDGING       Nudging toward sea surface height climatology.

 INITIAL: Configuring and initializing forward nonlinear model ...


 Vertical S-coordinate System: 

 level   S-coord     Cs-curve          at_hmin  over_slope     at_hmax

    10   0.0000000  -0.0000000          -0.000      -0.000      -0.000
     9  -0.1000000  -0.0105618          -0.032     -15.408     -30.785
     8  -0.2000000  -0.0361960          -0.109     -52.806    -105.504
     7  -0.3000000  -0.0969683          -0.291    -141.466    -282.641
     6  -0.4000000  -0.2224207          -0.667    -324.488    -648.308
     5  -0.5000000  -0.4163071          -1.249    -607.347   -1213.445
     4  -0.6000000  -0.6143559          -1.843    -896.278   -1790.714
     3  -0.7000000  -0.7533581          -2.260   -1099.067   -2195.875
     2  -0.8000000  -0.8405260          -2.522   -1226.236   -2449.950
     1  -0.9000000  -0.9121394          -2.736   -1330.712   -2658.688
     0  -1.0000000  -1.0000000          -3.000   -1458.891   -2914.782

 Time Splitting Weights: ndtfast =  20    nfast =  29

    Primary            Secondary            Accumulated to Current Step

  1-0.0009651193358779 0.0500000000000000-0.0009651193358779 0.0500000000000000
  2-0.0013488780126037 0.0500482559667939-0.0023139973484816 0.1000482559667939
  3-0.0011514592651645 0.0501156998674241-0.0034654566136460 0.1501639558342179
  4-0.0003735756740661 0.0501732728306823-0.0038390322877122 0.2003372286649002
  5 0.0009829200513762 0.0501919516143856-0.0028561122363360 0.2505291802792858
  6 0.0029141799764308 0.0501428056118168 0.0000580677400948 0.3006719858911026
  7 0.0054132615310267 0.0499970966129952 0.0054713292711215 0.3506690825040978
  8 0.0084687837865132 0.0497264335364439 0.0139401130576347 0.4003955160405417
  9 0.0120633394191050 0.0493029943471183 0.0260034524767397 0.4496985103876600
 10 0.0161716623600090 0.0486998273761630 0.0421751148367486 0.4983983377638230
 11 0.0207585511322367 0.0478912442581626 0.0629336659689853 0.5462895820219856
 12 0.0257765478740990 0.0468533167015507 0.0887102138430843 0.5931428987235363
 13 0.0311633730493853 0.0455644893078458 0.1198735868924696 0.6387073880313821
 14 0.0368391158442262 0.0440063206553765 0.1567127027366958 0.6827137086867585
 15 0.0427031802506397 0.0421643648631652 0.1994158829873354 0.7248780735499237
 16 0.0486309868367616 0.0400292058506332 0.2480468698240970 0.7649072794005569
 17 0.0544704302037591 0.0375976565087951 0.3025173000278562 0.8025049359093520
 18 0.0600380921294285 0.0348741349986072 0.3625553921572847 0.8373790709079592
 19 0.0651152103984763 0.0318722303921358 0.4276706025557610 0.8692513013000949
 20 0.0694434033194839 0.0286164698722119 0.4971140058752449 0.8978677711723068
 21 0.0727201499285569 0.0251442997062377 0.5698341558038018 0.9230120708785445
 22 0.0745940258796570 0.0215082922098099 0.6444281816834588 0.9445203630883544
 23 0.0746596950216180 0.0177785909158270 0.7190878767050768 0.9622989540041814
 24 0.0724526566618460 0.0140456061647461 0.7915405333669228 0.9763445601689276
 25 0.0674437485167025 0.0104229733316538 0.8589842818836253 0.9867675335005814
 26 0.0590334053485720 0.0070507859058187 0.9180176872321973 0.9938183194064002
 27 0.0465456732896125 0.0040991156383901 0.9645633605218099 0.9979174350447904
 28 0.0292219798521905 0.0017718319739095 0.9937853403740005 0.9996892670186999
 29 0.0062146596259994 0.0003107329813000 0.9999999999999999 0.9999999999999999

 ndtfast, nfast =   20  29   nfast/ndtfast = 1.45000

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

    1.000000000000 1.060707743385 0.530353871693 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.28400   0.14200

 Minimum X-grid spacing, DXmin =  2.97356675E+01 km
 Maximum X-grid spacing, DXmax =  2.97356675E+01 km
 Minimum Y-grid spacing, DYmin =  3.07439608E+01 km
 Maximum Y-grid spacing, DYmax =  3.07439608E+01 km
 Minimum Z-grid spacing, DZmin =  3.16853172E-02 m
 Maximum Z-grid spacing, DZmax =  5.77269100E+02 m

 Minimum barotropic Courant Number =  3.80717940E-03
 Maximum barotropic Courant Number =  1.18671316E-01
 Maximum Coriolis   Courant Number =  3.28944953E-02


 NLM: GET_STATE - Read state initial conditions,             t =   120 00:00:00
                   (File: ocean_ini.nc, Rec=0001, Index=1)
                - free-surface
                   (Min =  4.29925293E-03 Max =  1.39852402E-01)
                - vertically integrated u-momentum component
                   (Min = -1.58153090E-01 Max =  1.67690075E-01)
                - vertically integrated v-momentum component
                   (Min = -8.34430335E-02 Max =  1.30607023E-01)
                - u-momentum component
                   (Min = -1.98971007E-01 Max =  2.69165120E-01)
                - v-momentum component
                   (Min = -1.22494426E-01 Max =  2.07887629E-01)
                - potential temperature
                   (Min =  1.74229923E+00 Max =  1.50000000E+01)
                - salinity
                   (Min =  0.00000000E+00 Max =  3.46305472E+01)
    GET_NGFLD   - tidal period
                   (Min = -7.20002575E+03 Max = -7.20002575E+03)
    GET_2DFLD   - tidal elevation amplitude
                   (Min =  1.00000000E+37 Max = -1.00000000E+37)
    GET_2DFLD   - tidal elevation phase angle
                   (Min =  1.00000000E+37 Max = -1.00000000E+37)
    GET_2DFLD   - tidal current inclination angle
                   (Min =  1.00000000E+37 Max = -1.00000000E+37)
    GET_2DFLD   - tidal current phase angle
                   (Min =  1.00000000E+37 Max = -1.00000000E+37)
    GET_2DFLD   - maximum tidal current, ellipse major axis
                   (Min =  1.00000000E+37 Max = -1.00000000E+37)
    GET_2DFLD   - minimum tidal current, ellipse minor axis
                   (Min =  1.00000000E+37 Max = -1.00000000E+37)
    GET_2DFLD   - surface u-wind component,                  t =   120 00:00:00
                   (File: Uwind.nc, Rec=0714, Index=1)
                   (Tmin=         90.2917 Tmax=        273.2500)
                   (Min =  6.85470033E-01 Max =  8.74236333E+00)
    GET_2DFLD   - surface v-wind component,                  t =   120 00:00:00
                   (File: Vwind.nc, Rec=0714, Index=1)
                   (Tmin=         90.2917 Tmax=        273.2500)
                   (Min = -1.06013147E+01 Max =  6.27439000E-01)
    GET_2DFLD   - solar shortwave radiation flux,            t =   120 00:00:00
                   (File: swrad.nc, Rec=0714, Index=1)
                   (Tmin=         90.2917 Tmax=        273.2500)
                   (Min =  5.66614492E-05 Max =  1.40407721E-04)
    GET_2DFLD   - downwelling longwave radiation flux,       t =   120 00:00:00
                   (File: lwrad_down.nc, Rec=0714, Index=1)
                   (Tmin=         90.2917 Tmax=        273.2500)
                   (Min =  5.23297392E-05 Max =  7.78363329E-05)
    GET_2DFLD   - surface air pressure,                      t =   120 00:00:00
                   (File: Pair.nc, Rec=0714, Index=1)
                   (Tmin=         90.2917 Tmax=        273.2500)
                   (Min =  9.41421267E+02 Max =  1.02793233E+03)
    GET_2DFLD   - surface air temperature,                   t =   120 00:00:00
                   (File: Tair.nc, Rec=0714, Index=1)
                   (Tmin=         90.2917 Tmax=        273.2500)
                   (Min =  8.46280000E+00 Max =  1.24724000E+01)
    GET_2DFLD   - surface air relative humidity,             t =   120 00:00:00
                   (File: Qair.nc, Rec=0714, Index=1)
                   (Tmin=         90.2917 Tmax=        273.2500)
                   (Min =  4.82636058E-01 Max =  8.69312700E-01)
    GET_NGFLD   - free-surface western boundary condition,   t =   120 00:00:00
                   (File: ocean_bry.nc, Rec=0001, Index=2)
                   (Tmin=        120.0000 Tmax=        150.0000)
                   (Min =  4.29925293E-03 Max =  1.14858652E-01)
    GET_NGFLD   - free-surface eastern boundary condition,   t =   120 00:00:00
                   (File: ocean_bry.nc, Rec=0001, Index=2)
                   (Tmin=        120.0000 Tmax=        150.0000)
                   (Min =  6.69948349E-02 Max =  1.39385963E-01)
    GET_NGFLD   - free-surface southern boundary condition,  t =   120 00:00:00
                   (File: ocean_bry.nc, Rec=0001, Index=2)
                   (Tmin=        120.0000 Tmax=        150.0000)
                   (Min =  4.29925293E-03 Max =  7.11807813E-02)
    GET_NGFLD   - free-surface northern boundary condition,  t =   120 00:00:00
                   (File: ocean_bry.nc, Rec=0001, Index=2)
                   (Tmin=        120.0000 Tmax=        150.0000)
                   (Min =  1.11451183E-01 Max =  1.24161614E-01)
    GET_NGFLD   - 2D u-momentum western boundary condition,  t =   120 00:00:00
                   (File: ocean_bry.nc, Rec=0001, Index=2)
                   (Tmin=        120.0000 Tmax=        150.0000)
                   (Min = -2.23290309E-02 Max = -8.31978531E-03)
    GET_NGFLD   - 2D v-momentum western boundary condition,  t =   120 00:00:00
                   (File: ocean_bry.nc, Rec=0001, Index=2)
                   (Tmin=        120.0000 Tmax=        150.0000)
                   (Min =  5.34520762E-03 Max =  1.98139420E-02)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,  t =   120 00:00:00
                   (File: ocean_bry.nc, Rec=0001, Index=2)
                   (Tmin=        120.0000 Tmax=        150.0000)
                   (Min = -1.50792976E-01 Max =  1.67690075E-01)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,  t =   120 00:00:00
                   (File: ocean_bry.nc, Rec=0001, Index=2)
                   (Tmin=        120.0000 Tmax=        150.0000)
                   (Min = -8.34430335E-02 Max =  1.16547539E-01)
    GET_NGFLD   - 2D u-momentum southern boundary condition, t =   120 00:00:00
                   (File: ocean_bry.nc, Rec=0001, Index=2)
                   (Tmin=        120.0000 Tmax=        150.0000)
                   (Min = -1.01661586E-02 Max =  4.89103950E-03)
    GET_NGFLD   - 2D v-momentum southern boundary condition, t =   120 00:00:00
                   (File: ocean_bry.nc, Rec=0001, Index=2)
                   (Tmin=        120.0000 Tmax=        150.0000)
                   (Min =  1.35854326E-02 Max =  1.17949375E-01)
    GET_NGFLD   - 2D u-momentum northern boundary condition, t =   120 00:00:00
                   (File: ocean_bry.nc, Rec=0001, Index=2)
                   (Tmin=        120.0000 Tmax=        150.0000)
                   (Min = -1.63991366E-02 Max =  1.67690075E-01)
    GET_NGFLD   - 2D v-momentum northern boundary condition, t =   120 00:00:00
                   (File: ocean_bry.nc, Rec=0001, Index=2)
                   (Tmin=        120.0000 Tmax=        150.0000)
                   (Min = -8.34430335E-02 Max =  3.28114659E-02)
    GET_NGFLD   - 3D u-momentum western boundary condition,  t =   120 00:00:00
                   (File: ocean_bry.nc, Rec=0001, Index=2)
                   (Tmin=        120.0000 Tmax=        150.0000)
                   (Min = -1.62868211E-01 Max =  1.27311642E-02)
    GET_NGFLD   - 3D v-momentum western boundary condition,  t =   120 00:00:00
                   (File: ocean_bry.nc, Rec=0001, Index=2)
                   (Tmin=        120.0000 Tmax=        150.0000)
                   (Min = -1.06481557E-02 Max = 

any ideas? the only thing that pops out at me is why my tidal variables are 1e37! I've checked the input netcdf file and they're correct there. Thanks!

Dave Sutherland

[kate]

Can you run it in gdb? Recompile with USE_DEBUG first. I don't think the 1.e37 values should cause a seg fault, though they may cause other troubles. You may have two completely different problems, though both are happening so early in the run that a debugger should be useful for finding them. Or lots of print statements of the "made it past xyz" sort.

[das7105]

Hi Kate,

Thanks for the reply and suggestions. I recompiled with the use_debug flag 'on', in the same configuration as before. When I ran the code this time I got the error:

"At line 151 of file get_cycle.f90
Fortran runtime error: Array reference out of bounds for array 'tval', lower bound of dimension 1 exceeded (0 < 1)"

while looking at the ocean_clm.nc file for the ssh climatology. I assume this means my time referencing is off somehow. I'll dig through some more and see if I can find the error. It did get farther along than before though (i.e. the log file made it through more variables). Thanks again,

Dave

[das7105]

I finally got this run to make it- thanks for the help Kate. I think the problem was in the defining of the time variables I had in my input netcdf files.

One more question though- I can't seem to get rid of those 1.e37 values as min/max for my tidal values, and the weird thing is that for tidal period, it is giving me a min/max of -7.2e3 in the log! I've checked the input netcdf file, varinfo.dat, etc. but can't seem to find the problem. Plus my case ran, so is it a typo in the log, or is tide_period reading something else? here's the snippet from the log file with the info I'm talking about:
========
NLM: GET_STATE - Read state initial conditions, t = 122 00:00:00
(File: ocean_ini.nc, Rec=0001, Index=1)
- free-surface
(Min = -7.73503553E-03 Max = 1.04913927E-01)
- vertically integrated u-momentum component
(Min = -1.22036214E-01 Max = 1.95703370E-01)
- vertically integrated v-momentum component
(Min = -1.03021563E-01 Max = 1.21610556E-01)
- u-momentum component
(Min = -2.42570592E-01 Max = 2.32679226E-01)
- v-momentum component
(Min = -2.08331570E-01 Max = 2.22899850E-01)
- potential temperature
(Min = 1.74042434E+00 Max = 1.50000000E+01)
- salinity
(Min = 0.00000000E+00 Max = 3.46323478E+01)
GET_NGFLD - tidal period
(Min = -7.20002575E+03 Max = -7.20002575E+03)
GET_2DFLD - tidal elevation amplitude
(Min = 1.00000000E+37 Max = -1.00000000E+37)
GET_2DFLD - tidal elevation phase angle
(Min = 1.00000000E+37 Max = -1.00000000E+37)
GET_2DFLD - tidal current inclination angle
(Min = 1.00000000E+37 Max = -1.00000000E+37)
GET_2DFLD - tidal current phase angle
(Min = 1.00000000E+37 Max = -1.00000000E+37)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 1.00000000E+37 Max = -1.00000000E+37)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = 1.00000000E+37 Max = -1.00000000E+37)
GET_2DFLD - surface u-wind component, t = 122 00:00:00
=====================
any suggestions? I assume the 1.e37 are related to the FillValue for the netcdf files, but I can't figure out the tide period thing. Thanks,

dave

[m.hadfield]

The faulty tidal data sounds like this

https://www.myroms.org/projects/src/ticket/269

The recent IO updates appear to have broken a tidally-forced application of mine. The tidal period and parameters (elevation amplitude, etc) are not being read from the file, as illustrated by the following output:

GET_NGFLD - tidal period
(Min = 3.11285267-303 Max = 3.11285267-303)
GET_2DFLD - tidal elevation amplitude
(Min = 1.00000000E+37 Max = -1.00000000E+37)
GET_2DFLD - tidal elevation phase angle
(Min = 1.00000000E+37 Max = -1.00000000E+37)

The problem seems to be that in ROMS/Nonlinear/get_idata.F at line 101, NTC(ng) (i.e. the number of tidal consttituents) is set to 0

It's been fixed in the current version.

[m.hadfield]

Or this

https://www.myroms.org/projects/src/ticket/274

[das7105]

thanks- does this mean I need to fix get_idata.F myself, or just checkout the latest version? I'm running roms 3.2, svn revision 294M, which I got in early Jan.

[m.hadfield]

You should get the latest version, as there are fixes to files other than get_idata.F. You should be able to do this by running "svn update".

[antoinemarie]

Hi Dave,

I've seen you manage to find your first error:

"At line 151 of file get_cycle.f90
Fortran runtime error: Array reference out of bounds for array 'tval', lower bound of dimension 1 exceeded (0 < 1)"

I have the same problem but I don't manage to find my error. Can you tell me how did you do in your case please. I think it can help me.

Thanks,

Antoine MARIE