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
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
At line 496 of file get_date.f90
Internal Error: Bad format node
Thank you !