Internal Error: Bad format node at get_date.f90?

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
rcxy
Posts: 31
Joined: Fri Jul 02, 2010 1:28 pm
Location: Guangdong Ocean University

Internal Error: Bad format node at get_date.f90?

#1 Post by rcxy » Sat Mar 05, 2011 11:23 am

Hi,

when i run the model after 3 days, it report an error:

Code: Select all

Process Information:

 Thread #  0 (pid=       0) is active.
 Thread #  1 (pid=       0) is active.

 Model Input Parameters:  ROMS/TOMS version 3.4  
                          Wednesday - March 2, 2011 -  3:30:11 PM
 -----------------------------------------------------------------------------

 zhubay Model

 Operating system : CYGWIN
 CPU/hardware     : i686
 Compiler system  : gfortran
 Compiler command : /usr/bin/gfortran-4
 Compiler flags   : -frepack-arrays -fopenmp -O3 -ffast-math -ffree-form -ffree-line-length-none

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

 Local Root    : /cygdrive/e/ROMS/CASE/NEW/trunk
 Header Dir    : /cygdrive/e/ROMS/CASE/NEW/trunk/ROMS/Include
 Header file   : zhubay_1225.h
 Analytical Dir: /cygdrive/e/ROMS/CASE/NEW/trunk/ROMS/Functionals

 Resolution, Grid 01: 0218x0238x020,  Parallel Threads:  2,  Tiling: 002x002


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

     259200  ntimes          Number of timesteps for 3-D equations.
      5.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.
        120  nRST            Number of timesteps between the writing of data
                               into restart fields.
          1  ninfo           Number of timesteps between print of information
                               to standard output.
        120  nSTA            Number of timesteps between the writing of data
                               the stations file.
          T  ldefout         Switch to create a new output NetCDF file(s).
        720  nHIS            Number of timesteps between the writing fields
                               into history file.
 2.0000E+01  nl_tnu2(01)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 01: temp
 2.0000E+01  nl_tnu2(02)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 02: salt
 1.0000E+01  nl_visc2        NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for momentum.
 1.0000E-06  Akt_bak(01)     Background vertical mixing coefficient (m2/s)
                               for tracer 01: temp
 1.0000E-06  Akt_bak(02)     Background vertical mixing coefficient (m2/s)
                               for tracer 02: salt
 1.0000E-05  Akv_bak         Background vertical mixing coefficient (m2/s)
                               for momentum.
 1.0000E-05  Akk_bak         Background vertical mixing coefficient (m2/s)
                               for turbulent energy.
 1.0000E-05  Akp_bak         Background vertical mixing coefficient (m2/s)
                               for turbulent generic statistical field.
 1.0000E-04  rdrg            Linear bottom drag coefficient (m/s).
 1.0000E-04  rdrg2           Quadratic bottom drag coefficient.
 2.0000E-02  Zob             Bottom roughness (m).
          1  Vtransform      S-coordinate transformation equation.
          1  Vstretching     S-coordinate stretching function.
 8.0000E-01  theta_s         S-coordinate surface control parameter.
 1.2000E+00  theta_b         S-coordinate bottom  control parameter.
      1.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).
20101028.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.
 2.0000E+01  obcfac          Factor between passive and active
                               open boundary conditions.
     12.000  T0              Background potential temperature (C) constant.
     35.000  S0              Background salinity (PSU) constant.
      1.000  gamma2          Slipperiness variable: free-slip (1.0) or 
                                                    no-slip (-1.0).
          T  Hout(idFsur)    Write out free-surface.
          T  Hout(idUbar)    Write out 2D U-momentum component.
          T  Hout(idVbar)    Write out 2D V-momentum component.
          T  Hout(idUvel)    Write out 3D U-momentum component.
          T  Hout(idVvel)    Write out 3D V-momentum component.
          T  Hout(idWvel)    Write out W-momentum component.
          T  Hout(idOvel)    Write out omega vertical velocity.
          T  Hout(idTvar)    Write out tracer 01: temp
          T  Hout(idTvar)    Write out tracer 02: salt

 Output/Input Files:

             Output Restart File:  ocean_rst.nc
             Output History File:  ocean_his.nc
            Output Stations File:  ocean_sta.nc
                 Input Grid File:  Data/ROMS/Grid/zhubay_grd.nc
    Input Nonlinear Initial File:  Data/ROMS/Initial/zhubay_ini.nc
          Station positions File:  ROMS/External/stations_zhubay0225.in

 Tile partition information for Grid 01:  0218x0238x0020  tiling: 002x002

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1      109        1      119   259420
        1      110      218        1      119   259420
        2        1      109      120      238   259420
        3      110      218      120      238   259420

 Tile minimum and maximum fractional grid coordinates:
   (interior points only)

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0    -0.50   220.50    -0.50   240.50  RHO-points
        1    -0.50   220.50    -0.50   240.50  RHO-points
        2    -0.50   220.50    -0.50   240.50  RHO-points
        3    -0.50   220.50    -0.50   240.50  RHO-points

        0    -0.50   220.50    -0.50   240.50    U-points
        1    -0.50   220.50    -0.50   240.50    U-points
        2    -0.50   220.50    -0.50   240.50    U-points
        3    -0.50   220.50    -0.50   240.50    U-points

        0    -0.50   220.50    -0.50   240.50    V-points
        1    -0.50   220.50    -0.50   240.50    V-points
        2    -0.50   220.50    -0.50   240.50    V-points
        3    -0.50   220.50    -0.50   240.50    V-points


 Stations Parameters, Grid: 01
 =============================

          7  Nstation        Number of stations to write out into stations file.
          T  Sout(idFsur)    Write out free-surface.
          T  Sout(idUbar)    Write out 2D U-momentum component.
          T  Sout(idVbar)    Write out 2D V-momentum component.
          T  Sout(idUvel)    Write out 3D U-momentum component.
          T  Sout(idVvel)    Write out 3D V-momentum component.
          T  Sout(idWvel)    Write out W-momentum component.
          T  Sout(idOvel)    Write out omega vertical velocity.
          T  Sout(idTvar)    Write out tracer 01: temp
          T  Sout(idTvar)    Write out tracer 02: salt

             Flag and positions for station 0001:  0    83.0000  144.0000
             Flag and positions for station 0002:  0    85.0000  131.0000
             Flag and positions for station 0003:  0    88.0000  143.0000
             Flag and positions for station 0004:  0   100.0000  132.0000
             Flag and positions for station 0005:  0   100.0000  144.0000
             Flag and positions for station 0006:  0   110.0000  136.0000
             Flag and positions for station 0007:  0   121.0000  140.0000

 Activated C-preprocessing Options:

 ZHUBAY              zhubay Model
 ADD_FSOBC           Adding tidal elevation to proccesed OBC data.
 ADD_M2OBC           Adding tidal currents to proccesed OBC data.
 ANA_BSFLUX          Analytical kinematic bottom salinity flux.
 ANA_BTFLUX          Analytical kinematic bottom temperature flux.
 ANA_FSOBC           Analytical free-surface boundary conditions.
 ANA_M2OBC           Analytical 2D momentum boundary conditions.
 ANA_SMFLUX          Analytical kinematic surface momentum flux.
 ANA_SSFLUX          Analytical kinematic surface salinity flux.
 ANA_STFLUX          Analytical kinematic surface temperature flux.
 ASSUMED_SHAPE       Using assumed-shape arrays.
 CURVGRID            Orthogonal curvilinear grid.
 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_M3RADIATION    Eastern edge, 3D momentum, radiation condition.
 EAST_TGRADIENT      Eastern edge, tracers, gradient condition.
 KANTHA_CLAYSON      Kantha and Clayson stability function formulation.
 MASKING             Land/Sea masking.
 MIX_GEO_TS          Mixing of tracers along geopotential surfaces.
 MIX_S_UV            Mixing of momentum along constant S-surfaces.
 MY25_MIXING         Mellor/Yamada Level-2.5 mixing closure.
 NONLINEAR           Nonlinear Model.
 NONLIN_EOS          Nonlinear Equation of State for seawater.
 NORTHERN_WALL       Wall boundary at Northern edge.
 N2S2_HORAVG         Horizontal smoothing of buoyancy and shear.
 OUT_DOUBLE          Double precision output fields in NetCDF files.
 _OPENMP             OpenMP parallel shared-memory directives.
 POWER_LAW           Power-law shape time-averaging barotropic filter.
 PROFILE             Time profiling activated .
 K_GSCHEME           Third-order upstream advection of TKE fields.
 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.
 SOUTHERN_WALL       Wall boundary at Southern edge.
 SPLINES             Conservative parabolic spline reconstruction.
 SPHERICAL           Spherical grid configuration.
 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_C4VADVECTION     Fourth-order centered vertical advection of momentum.
 UV_LDRAG            Linear bottom stress.
 UV_TIDES            Add tidal currents to 2D momentum climatologies.
 UV_VIS2             Harmonic mixing of momentum.
 VAR_RHO_2D          Variable density barotropic mode.
 WESTERN_WALL        Wall boundary at Western edge.

 INITIAL: Configuring and initializing forward nonlinear model ...


 Vertical S-coordinate System: 

 level   S-coord     Cs-curve          at_hmin  over_slope     at_hmax

    20   0.0000000  -0.0000000          -0.000      -0.000      -0.000
    19  -0.0500000  -0.0458417          -0.233      -0.637      -1.041
    18  -0.1000000  -0.0932041          -0.473      -1.294      -2.115
    17  -0.1500000  -0.1419507          -0.718      -1.968      -3.218
    16  -0.2000000  -0.1919290          -0.968      -2.658      -4.348
    15  -0.2500000  -0.2429724          -1.222      -3.362      -5.502
    14  -0.3000000  -0.2949008          -1.480      -4.077      -6.674
    13  -0.3500000  -0.3475226          -1.740      -4.801      -7.861
    12  -0.4000000  -0.4006364          -2.003      -5.531      -9.060
    11  -0.4500000  -0.4540334          -2.266      -6.265     -10.264
    10  -0.5000000  -0.5074993          -2.530      -7.000     -11.469
     9  -0.5500000  -0.5608171          -2.793      -7.733     -12.672
     8  -0.6000000  -0.6137698          -3.055      -8.461     -13.866
     7  -0.6500000  -0.6661423          -3.315      -9.181     -15.048
     6  -0.7000000  -0.7177247          -3.571      -9.892     -16.213
     5  -0.7500000  -0.7683138          -3.823     -10.590     -17.357
     4  -0.8000000  -0.8177159          -4.071     -11.273     -18.474
     3  -0.8500000  -0.8657483          -4.313     -11.938     -19.563
     2  -0.9000000  -0.9122412          -4.549     -12.583     -20.618
     1  -0.9500000  -0.9570387          -4.778     -13.207     -21.636
     0  -1.0000000  -1.0000000          -5.000     -13.807     -22.614

 Time Splitting Weights: ndtfast =  10    nfast =  15

    Primary            Secondary            Accumulated to Current Step

  1 0.0014850648765542 0.1000000000000000 0.0014850648765542 0.1000000000000000
  2 0.0059386031399401 0.0998514935123446 0.0074236680164944 0.1998514935123446
  3 0.0133457074936657 0.0992576331983506 0.0207693755101601 0.2991091267106951
  4 0.0236484051180389 0.0979230624489840 0.0444177806281990 0.3970321891596791
  5 0.0366959666818607 0.0955582219371801 0.0811137473100597 0.4925904110968592
  6 0.0521753389587947 0.0918886252689940 0.1332890862688544 0.5844790363658532
  7 0.0695217010484137 0.0866710913731146 0.2028107873172681 0.6711501277389678
  8 0.0878091442019233 0.0797189212682732 0.2906199315191914 0.7508690490072409
  9 0.1056214752525633 0.0709380068480809 0.3962414067717547 0.8218070558553218
 10 0.1209031436506851 0.0603758593228245 0.5171445504224398 0.8821829151781463
 11 0.1307902921035068 0.0482855449577560 0.6479348425259466 0.9304684601359023
 12 0.1314219308195445 0.0352065157474053 0.7793567733454911 0.9656749758833076
 13 0.1177312353577227 0.0220643226654509 0.8970880087032138 0.9877392985487585
 14 0.0832169680811588 0.0102911991296786 0.9803049767843727 0.9980304976784371
 15 0.0196950232156273 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

 Minimum X-grid spacing, DXmin =  2.02651619E-01 km
 Maximum X-grid spacing, DXmax =  5.12538984E-01 km
 Minimum Y-grid spacing, DYmin =  2.55004280E-01 km
 Maximum Y-grid spacing, DYmax =  1.54597709E+00 km
 Minimum Z-grid spacing, DZmin =  2.21845093E-01 m
 Maximum Z-grid spacing, DZmax =  1.20563425E+00 m

 Minimum barotropic Courant Number =  8.24929208E-03
 Maximum barotropic Courant Number =  2.29935423E-02
 Maximum Coriolis   Courant Number =  2.66051206E-04


 NLM: GET_STATE - Read state initial conditions,             t =     0 00:00:00
                   (File: zhubay_ini.nc, Rec=0001, Index=1)
                - free-surface
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - vertically integrated u-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - vertically integrated v-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - u-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - v-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - potential temperature
                   (Min =  2.24567986E+01 Max =  2.27305959E+01)
                - salinity
                   (Min =  3.41531511E+01 Max =  3.41673488E+01)
    GET_NGFLD   - tidal period
                   (Min =  4.30819210E+04 Max =  9.67262421E+04)
    GET_2DFLD   - tidal elevation amplitude
                   (Min =  5.73107538E-02 Max =  5.85056946E-01)
    GET_2DFLD   - tidal elevation phase angle
                   (Min =  1.36509702E+00 Max =  5.21011036E+00)
    GET_2DFLD   - tidal current inclination angle
                   (Min =  1.98049522E-06 Max =  6.28318432E+00)
    GET_2DFLD   - tidal current phase angle
                   (Min =  4.25854135E-03 Max =  6.28233562E+00)
    GET_2DFLD   - maximum tidal current, ellipse major axis
                   (Min =  9.60353901E-03 Max =  1.25420363E-01)
    GET_2DFLD   - minimum tidal current, ellipse minor axis
                   (Min = -1.56131073E-02 Max =  7.61347671E-03)

 Maximum grid stiffness ratios:  rx0 =   7.713563E-02 (Beckmann and Haidvogel)
                                 rx1 =   3.456655E+00 (Haney)


 Initial basin volumes: TotVolume =  3.2074008489E+10 m3
                        MinVolume =  1.8703482615E+04 m3
                        MaxVolume =  2.7515518927E+05 m3
                          Max/Min =  1.4711441443E+01

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00259200)

   STEP   Day HH:MM:SS  KINETIC_ENRG   POTEN_ENRG    TOTAL_ENRG    NET_VOLUME

      0     0 00:00:00  0.000000E+00  4.459277E+01  4.459277E+01  4.712004E+10
      DEF_HIS   - creating history file: ocean_his.nc
      WRT_HIS   - wrote history  fields (Index=1,1) into time record = 0000001
      DEF_STATION - creating stations file: ocean_sta.nc
      1     0 00:00:05  6.313774E-18  4.459277E+01  4.459277E+01  4.712004E+10
      2     0 00:00:10  1.065469E-11  4.459277E+01  4.459277E+01  4.712004E+10
      3     0 00:00:15  3.926059E-11  4.459277E+01  4.459277E+01  4.712004E+10
      4     0 00:00:20  8.142591E-11  4.459277E+01  4.459277E+01  4.712005E+10
      5     0 00:00:25  1.340193E-10  4.459278E+01  4.459278E+01  4.712005E+10
      6     0 00:00:30  1.954993E-10  4.459278E+01  4.459278E+01  4.712005E+10
      7     0 00:00:35  2.660043E-10  4.459279E+01  4.459279E+01  4.712005E+10
      8     0 00:00:40  3.472036E-10  4.459279E+01  4.459279E+01  4.712005E+10
      9     0 00:00:45  4.419565E-10  4.459280E+01  4.459280E+01  4.712006E+10
     10     0 00:00:50  5.538565E-10  4.459280E+01  4.459280E+01  4.712006E+10
..................
  66775     3 20:44:35  4.771685E-02  4.386175E+01  4.390947E+01  4.547124E+10
  66776     3 20:44:40  4.777835E-02  4.386318E+01  4.391096E+01  4.547301E+10
  66777     3 20:44:45  4.783993E-02  4.386460E+01  4.391244E+01  4.547478E+10
  66778     3 20:44:50  4.790159E-02  4.386603E+01  4.391393E+01  4.547655E+10
  66779     3 20:44:55At line 496 of file get_date.f90
Internal Error: Bad format node
  4.796333E-02  4.386746E+01  4.391542E+01  4.547832E+10
  66780     3 20:45:00  4.802515E-02  4.386888E+01  4.391691E+01  4.548009E+10
  66781     3 20:45:05  4.808705E-02  4.387031E+01  4.391840E+01  4.548187E+10
  66782     3 20:45:10  4.814903E-02  4.387174E+01  4.391989E+01  4.548365E+10
  66783     3 20:45:15  4.821110E-02  4.387317E+01  4.392138E+01  4.548542E+10
  66784     3 20:45:20  4.827324E-02  4.387461E+01  4.392288E+01  4.548720E+10
  66785     3 20:45:25  4.833546E-02  4.387604E+01  4.392437E+01  4.548899E+10
  66786     3 20:45:30  4.839776E-02  4.387747E+01  4.392587E+01  4.549077E+10
  66787     3 20:45:35  4.846014E-02  4.387891E+01  4.392737E+01  4.549255E+10
  66788     3 20:45:40  4.852261E-02  4.388034E+01  4.392887E+01  4.549434E+10
i open the get_date.f90, the code is :

Code: Select all

490      iday=INT(time/86400)
491      isec=INT(time-iday*86400.0_r8)
492      ihour=isec/3600
493      minute=MOD(isec,3600)/60
494      isec=MOD(isec,60)
495!
496      WRITE (text,10) iday, ihour, minute, isec
497 10   FORMAT (i5,1x,i2.2,':',i2.2,':',i2.2)
498      time_code=text
Why the error :

At line 496 of file get_date.f90
Internal Error: Bad format node


Thank you !

Post Reply