GET_2DFLD - error while reading variable: sms_time

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
rainbowbluegz
Posts: 5
Joined: Sun Jan 06, 2008 4:54 pm
Location: SCSIO

GET_2DFLD - error while reading variable: sms_time

#1 Unread post by rainbowbluegz »

hi,everybady:
The following is the output information of ROMS2.2:
Process Information:

Thread # 0 (pid= 26018) is active.

Model Input Parameters: ROMS/TOMS version 2.2
Thursday - September 22, 2011 - 10:33:11 AM
-----------------------------------------------------------------------------

South China Sea, 1/10 degree resolution

Operating system : Linux
CPU/hardware : i686
Compiler system : ifort
Compiler command : ifort
Compiler flags : -ip -O3 -pc80 -xW

Resolution, Grid 01: 0155x0137x028, Parallel Threads: 1, Tiling: 001x001


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

2976 ntimes Number of timesteps for 3-D equations.
900.000 dt Timestep size (s) for 3-D equations.
60 ndtfast Number of timesteps for 2-D equations between
each 3D timestep.
1 ERstr Starting ensemble/perturbation run number.
2976 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.
2976 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
T ldefout Switch to create a new output NetCDF file(s).
960 nHIS Number of timesteps between the writing fields
into history file.
2976 ndefHIS Number of timesteps between creation of new
history files.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
192 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
5.0000E-01 tnu2(01) Horizontal, harmonic mixing coefficient (m2/s)
for tracer 01: temp
5.0000E-01 tnu2(02) Horizontal, harmonic mixing coefficient (m2/s)
for tracer 02: salt
0.0000E+00 visc4 Horizontal, biharmonic mixing coefficient (m4/s)
for momentum.
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-05 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
5.0000E-03 Zob Bottom roughness (m).
1 lmd_Jwt Jerlov water type.
6.0000E+00 theta_s S-coordinate surface control parameter.
0.0000E+00 theta_b S-coordinate bottom control parameter.
10.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
12538.000 dstart Time-stamp assigned to model initialization (days).
12538.000 tide_start Reference time origin for tidal forcing (days).
-2.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.
10.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(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(idVvis) Write out vertical viscosity coefficient.

Output/Input Files:

Output Restart File: roms_rst.nc
Prefix for History Files: roms_his
Output Averages File: roms_avg.nc
Input Grid File: roms_grd.nc
Input Nonlinear Initial File: roms_ini_ECCO_Y2010M1.nc
Input Forcing File 01: roms_frc_NCEP2_Y2010M1.nc
Biology Parameters File: /home/xie/roms-2.2/External/bioFasham_scs.in
IO Variable Information File: /home/xie/roms-2.2/varinfo.dat

Tile partition information for Grid 01: 0155x0137x0028 tiling: 001x001

tile Istr Iend Jstr Jend Npts

0 1 155 1 137 594580


Fasham Model Parameters, Grid: 01
=================================

1 BioIter Number of iterations for nonlinear convergence.
4.0000E-02 AttSW Light attenuation of seawater (m-1).
2.4860E-02 AttChl Light attenuation by chlorophyll (1/(mg_Chl m-2)).
4.3000E-01 PARfrac Fraction of shortwave radiation that is
photosynthetically active (nondimensional).
1.0000E+00 Vp0 Eppley temperature-limited growth parameter
(nondimensional).
9.5000E-03 I_thNH4 Radiation threshold for nitrification (W/m2).
1.0000E-01 D_p5NH4 Half-saturation radiation for nitrification (W/m2).
5.0000E-02 NitriR Nitrification rate (day-1).
2.0000E+00 K_NO3 Inverse half-saturation for phytoplankton NO3
uptake (1/(mmol_N m-3)).
2.0000E+00 K_NH4 Inverse half-saturation for phytoplankton NH4
uptake (1/(mmol_N m-3)).
2.0000E+00 K_Phy Zooplankton half-saturation constant for ingestion
(day-1).
5.3500E-02 Chl2C_m Maximum chlorophyll to carbon ratio (mg_Chl/mg_C).
1.0000E-03 ChlMin Chlorophyll minimum threshold (mg_Chl/m3).
6.6250E+00 PhyCN Phytoplankton Carbon:Nitrogen ratio (mol_C/mol_N).
1.5000E+00 PhyIP Phytoplankton NH4 inhibition parameter (1/mmol_N).
2.5000E-02 PhyIS Phytoplankton growth, initial slope of P-I curve
(mg_C/(mg_Chl Watts m-2 day)).
1.0000E-03 PhyMin Phytoplankton minimum threshold (mmol_N/m3).
1.5000E-01 PhyMR Phytoplankton mortality rate (day-1).
7.5000E-01 ZooAE_N Zooplankton nitrogen assimilation efficiency
(nondimensional).
1.0000E-01 ZooBM Rate for zooplankton basal metabolism (1/day).
6.6250E+00 ZooCN Zooplankton Carbon:Nitrogen ratio (mol_C/mol_N).
1.0000E-01 ZooER Zooplankton specific excretion rate (day-1).
6.0000E-01 ZooGR Zooplankton maximum growth rate (day-1).
1.0000E-03 ZooMin Zooplankton minimum threshold (mmol_N/m3).
1.0000E-01 ZooMR Zooplankton mortality rate (day-1).
1.0000E-02 LDeRRN Large detritus N re-mineralization rate (day-1).
1.0000E-02 LDeRRC Large detritus C re-mineralization rate (day-1).
5.0000E-03 CoagR Coagulation rate (day-1).
3.0000E-02 SDeRRN Remineralization rate for small detritus N (day-1).
3.0000E-02 SDeRRC Remineralization rate for small detritus C (day-1).
1.0000E-01 wPhy Phytoplankton sinking velocity (m/day).
1.0000E+00 wLDet Large detritus sinking velocity (m/day).
1.0000E-01 wSDet Small detritus sinking velocity (m/day).
0.0000E+00 tnu2(03) Horizontal, harmonic mixing coefficient (m2/s)
for tracer 03: NO3
0.0000E+00 tnu2(04) Horizontal, harmonic mixing coefficient (m2/s)
for tracer 04: NH4
0.0000E+00 tnu2(05) Horizontal, harmonic mixing coefficient (m2/s)
for tracer 05: chlorophyll
0.0000E+00 tnu2(06) Horizontal, harmonic mixing coefficient (m2/s)
for tracer 06: phytoplankton
0.0000E+00 tnu2(07) Horizontal, harmonic mixing coefficient (m2/s)
for tracer 07: zooplankton
0.0000E+00 tnu2(08) Horizontal, harmonic mixing coefficient (m2/s)
for tracer 08: LdetritusN
0.0000E+00 tnu2(09) Horizontal, harmonic mixing coefficient (m2/s)
for tracer 09: SdetritusN
1.0000E-06 Akt_bak(03) Background vertical mixing coefficient (m2/s)
for tracer 03: NO3
1.0000E-06 Akt_bak(04) Background vertical mixing coefficient (m2/s)
for tracer 04: NH4
1.0000E-06 Akt_bak(05) Background vertical mixing coefficient (m2/s)
for tracer 05: chlorophyll
1.0000E-06 Akt_bak(06) Background vertical mixing coefficient (m2/s)
for tracer 06: phytoplankton
1.0000E-06 Akt_bak(07) Background vertical mixing coefficient (m2/s)
for tracer 07: zooplankton
1.0000E-06 Akt_bak(08) Background vertical mixing coefficient (m2/s)
for tracer 08: LdetritusN
1.0000E-06 Akt_bak(09) Background vertical mixing coefficient (m2/s)
for tracer 09: SdetritusN
0.0000E+00 Tnudg(03) Nudging/relaxation time scale (days)
for tracer 03: NO3
0.0000E+00 Tnudg(04) Nudging/relaxation time scale (days)
for tracer 04: NH4
0.0000E+00 Tnudg(05) Nudging/relaxation time scale (days)
for tracer 05: chlorophyll
0.0000E+00 Tnudg(06) Nudging/relaxation time scale (days)
for tracer 06: phytoplankton
0.0000E+00 Tnudg(07) Nudging/relaxation time scale (days)
for tracer 07: zooplankton
0.0000E+00 Tnudg(08) Nudging/relaxation time scale (days)
for tracer 08: LdetritusN
0.0000E+00 Tnudg(09) Nudging/relaxation time scale (days)
for tracer 09: SdetritusN
T Hout(idTvar) Write out tracer 03: NO3
T Hout(idTvar) Write out tracer 04: NH4
T Hout(idTvar) Write out tracer 05: chlorophyll
T Hout(idTvar) Write out tracer 06: phytoplankton
T Hout(idTvar) Write out tracer 07: zooplankton
T Hout(idTvar) Write out tracer 08: LdetritusN
T Hout(idTvar) Write out tracer 09: SdetritusN

Activated C-preprocessing Options:

ANA_BIOLOGY Analytical biology initial conditions.
ANA_BPFLUX Analytical bottom passive tracers fluxes.
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_SPFLUX Analytical surface passive tracer fluxes.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged fields.
AVERAGES_AKS Writing out time-averaged vertical S-diffusion.
AVERAGES_AKT Writing out time-averaged vertical T-diffusion.
BIO_FASHAM Fasham type nitrogen-based model.
BIO_SEDIMENT Restore fallen particulate material to the nutrient pool.
DENITRIFICATION Add denitrification processes.
DIURNAL_SRFLUX Modulate shortwave radiation by the local diurnal cycle.
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_TRADIATION Eastern edge, tracers, radiation condition.
LMD_BKPP KPP bottom boundary layer mixing.
LMD_CONVEC LMD convective mixing due to shear instability.
LMD_MIXING Large/McWilliams/Doney interior mixing.
LMD_NONLOCAL LMD convective nonlocal transport.
LMD_RIMIX LMD diffusivity due to shear instability.
LMD_SKPP KPP surface boundary layer mixing.
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.
MIX_S_UV Mixing of momentum along constant S-surfaces.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
NORTH_FSGRADIENT Northern edge, free-surface, gradient condition.
NORTH_M2RADIATION Northern edge, 2D momentum, radiation condition.
NORTH_M3RADIATION Northern edge, 3D momentum, radiation condition.
NORTH_TRADIATION Northern edge, tracers, radiation condition.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLVE3D Solving 3D Primitive Equations.
SOUTH_FSGRADIENT Southern edge, free-surface, gradient condition.
SOUTH_M2RADIATION Southern edge, 2D momentum, radiation condition.
SOUTH_M3RADIATION Southern edge, 3D momentum, radiation condition.
SOUTH_TRADIATION Southern edge, tracers, radiation condition.
SPONGE Enhanced horizontal mixing in the sponge areas.
SSH_TIDES Add tidal elevation to SSH climatology.
TCLIMATOLOGY Processing tracer climatology data.
TCLM_NUDGING Nudging toward tracer climatology.
TS_U3HADVECTION Third-order upstream bias horizontal advection of tracers.
TS_A4VADVECTION Fouth-order Akima vertical advection of tracers.
TS_DIF2 Harmonic mixing of tracers.
UV_ADV Advection of momentum.
UV_COR Coriolis term.
UV_U3HADVECTION Third-order upstream bias advection of momentum.
UV_QDRAG Quadratic bottom stress.
UV_TIDES Add tidal currents to 2D momentum climatologies.
UV_VIS4 Biharmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
WESTERN_WALL Wall boundary at Western edge.
ZCLIMATOLOGY Processing sea surface height climatology data.
ZCLM_NUDGING Nudging toward sea surface height climatology.

INITIAL: Configurating and initializing forward nonlinear model ...


Vertical S-coordinate System:

level S-coord Cs-curve at_hmin over_slope at_hmax

28 0.0000000 0.0000000 0.000 0.000 0.000
27 -0.0357143 -0.0010705 -0.453 -3.076 -5.699
26 -0.0714286 -0.0021903 -0.911 -6.278 -11.644
25 -0.1071429 -0.0034111 -1.378 -9.736 -18.093
24 -0.1428571 -0.0047891 -1.860 -13.593 -25.326
23 -0.1785714 -0.0063879 -2.361 -18.011 -33.661
22 -0.2142857 -0.0082811 -2.888 -23.177 -43.465
21 -0.2500000 -0.0105560 -3.450 -29.312 -55.174
20 -0.2857143 -0.0133175 -4.056 -36.684 -69.311
19 -0.3214286 -0.0166928 -4.717 -45.614 -86.511
18 -0.3571429 -0.0208376 -5.447 -56.499 -107.551
17 -0.3928571 -0.0259429 -6.263 -69.824 -133.384
16 -0.4285714 -0.0322440 -7.188 -86.185 -165.183
15 -0.4642857 -0.0400314 -8.246 -106.323 -204.399
14 -0.5000000 -0.0496640 -9.470 -131.146 -252.823
13 -0.5357143 -0.0615858 -10.900 -161.785 -312.670
12 -0.5714286 -0.0763464 -12.585 -199.634 -386.683
11 -0.6071429 -0.0946261 -14.588 -246.422 -478.256
10 -0.6428571 -0.1172675 -16.983 -304.288 -591.593
9 -0.6785714 -0.1453143 -19.864 -375.884 -731.904
8 -0.7142857 -0.1800593 -23.348 -464.493 -905.639
7 -0.7500000 -0.2231040 -27.579 -574.184 -1120.789
6 -0.7857143 -0.2764325 -32.736 -709.996 -1387.255
5 -0.8214286 -0.3425030 -39.040 -878.172 -1717.304
4 -0.8571429 -0.4243610 -46.764 -1086.448 -2126.133
3 -0.8928571 -0.5257796 -56.249 -1344.409 -2632.569
2 -0.9285714 -0.6514336 -67.915 -1663.927 -3259.940
1 -0.9642857 -0.8071151 -82.283 -2059.715 -4037.147
0 -1.0000000 -1.0000000 -100.000 -2550.000 -5000.000

Time Splitting Weights: ndtfast = 60 nfast = 83

Primary Secondary Accumulated to Current Step

1-0.0003540175431531 0.0166666666666667-0.0003540175431531 0.0166666666666667
2-0.0006728864853186 0.0166725669590526-0.0010269040284717 0.0333392336257192
3-0.0009566069754562 0.0166837817338079-0.0019835110039279 0.0500230153595271
4-0.0012051795928536 0.0166997251833988-0.0031886905967815 0.0667227405429259
5-0.0014186058436584 0.0167198115099464-0.0046072964404400 0.0834425520528723
6-0.0015968888560239 0.0167434549406740-0.0062041852964638 0.1001860069935463
7-0.0017400342738662 0.0167700697549411-0.0079442195703301 0.1169560767484873
8-0.0018480513492362 0.0167990703261722-0.0097922709195663 0.1337551470746595
9-0.0019209542333025 0.0168298711819928-0.0117132251528688 0.1505850182566523
10-0.0019587634659484 0.0168618870858811-0.0136719886188171 0.1674469053425335
11-0.0019615076639817 0.0168945331436470-0.0156334962827989 0.1843414384861804
12-0.0019292254079571 0.0169272249380467-0.0175627216907560 0.2012686634242271
13-0.0018619673276114 0.0169593786948459-0.0194246890183674 0.2182280421190730
14-0.0017597983859117 0.0169904114836395-0.0211844874042792 0.2352184536027125
15-0.0016228003617172 0.0170197414567380-0.0228072877659964 0.2522381950594504
16-0.0014510745310524 0.0170467881294333-0.0242583622970488 0.2692849831888837
17-0.0012447445469948 0.0170709727049508-0.0255031068440436 0.2863559558938346
18-0.0010039595181742 0.0170917184474007-0.0265070663622178 0.3034476743412353
19-0.0007288972858858 0.0171084511060370-0.0272359636481036 0.3205561254472722
20-0.0004197678998158 0.0171205993941351-0.0276557315479194 0.3376767248414073
21-0.0000768172923796 0.0171275955257987-0.0277325488402990 0.3548043203672060
22 0.0002996688483262 0.0171288758140050-0.0274328799919727 0.3719331961812110
23 0.0007093610069602 0.0171238813331996-0.0267235189850126 0.3890570775144105
24 0.0011518815707581 0.0171120586497502-0.0255716374142545 0.4061691361641608
25 0.0016268003607435 0.0170928606235709-0.0239448370535109 0.4232619967877317
26 0.0021336299643244 0.0170657472842252-0.0218112070891865 0.4403277440719569
27 0.0026718208692771 0.0170301867848198-0.0191393862199094 0.4573579308567766
28 0.0032407563991181 0.0169856564369985-0.0158986298207914 0.4743435872937751
29 0.0038397474498617 0.0169316438303465-0.0120588823709297 0.4912752311241216
30 0.0044680270281660 0.0168676480395155-0.0075908553427637 0.5081428791636371
31 0.0051247445908649 0.0167931809223794-0.0024661107518988 0.5249360600860165
32 0.0058089601858883 0.0167077685125317 0.0033428494339895 0.5416438285985481
33 0.0065196383945684 0.0166109525094335 0.0098624878285579 0.5582547811079817
34 0.0072556420753336 0.0165022918695240 0.0171181299038915 0.5747570729775057
35 0.0080157259087899 0.0163813645016018 0.0251338558126814 0.5911384374791074
36 0.0087985297441882 0.0162477690697886 0.0339323855568696 0.6073862065488961
37 0.0096025717472802 0.0161011269073855 0.0435349573041498 0.6234873334562816
38 0.0104262413495604 0.0159410840449308 0.0539611986537102 0.6394284175012125
39 0.0112677919988951 0.0157673133557715 0.0652289906526052 0.6551957308569840
40 0.0121253337115396 0.0155795168224566 0.0773543243641448 0.6707752476794406
41 0.0129968254255413 0.0153774279272643 0.0903511497896861 0.6861526756067049
42 0.0138800671555305 0.0151608141701719 0.1042312169452167 0.7013134897768768
43 0.0147726919488984 0.0149294797175797 0.1190039088941151 0.7162429694944565
44 0.0156721576433617 0.0146832681850981 0.1346760665374768 0.7309262376795547
45 0.0165757384259147 0.0144220655577087 0.1512518049633915 0.7453483032372634
46 0.0174805161931688 0.0141458032506101 0.1687323211565603 0.7594941064878735
47 0.0183833717130782 0.0138544613140573 0.1871156928696385 0.7733485678019308
48 0.0192809755880536 0.0135480717855060 0.2063966684576921 0.7868966395874368
49 0.0201697790194623 0.0132267221923718 0.2265664474771544 0.8001233617798086
50 0.0210460043735164 0.0128905592087141 0.2476124518506708 0.8130139209885227
51 0.0219056355485466 0.0125397924691555 0.2695180873992175 0.8255537134576781
52 0.0227444081436646 0.0121746985433464 0.2922624955428821 0.8377284120010245
53 0.0235577994288118 0.0117956250742853 0.3158202949716939 0.8495240370753099
54 0.0243410181161952 0.0114029950838051 0.3401613130878892 0.8609270321591150
55 0.0250889939331105 0.0109973114485352 0.3652503070209996 0.8719243436076503
56 0.0257963669961523 0.0105791615496500 0.3910466740171519 0.8825035051573002
57 0.0264574769868114 0.0101492220997141 0.4175041510039634 0.8926527272570144
58 0.0270663521284595 0.0097082641499339 0.4445705031324229 0.9023609914069484
59 0.0276166979647200 0.0092571582811263 0.4721872010971429 0.9116181496880746
60 0.0281018859392274 0.0087968799817143 0.5002890870363703 0.9204150296697889
61 0.0285149417767724 0.0083285152160605 0.5288040288131427 0.9287435448858494
62 0.0288485336658352 0.0078532661864476 0.5576525624789779 0.9365968110722971
63 0.0290949602425051 0.0073724572920170 0.5867475227214829 0.9439692683643142
64 0.0292461383757877 0.0068875412879753 0.6159936610972706 0.9508568096522895
65 0.0292935907542990 0.0064001056483788 0.6452872518515697 0.9572569153006683
66 0.0292284332743468 0.0059118791358072 0.6745156851259165 0.9631687944364754
67 0.0290413622293988 0.0054247385812347 0.7035570473553153 0.9685935330177101
68 0.0287226413009382 0.0049407158774114 0.7322796886562536 0.9735342488951215
69 0.0282620883507064 0.0044620051890624 0.7605417770069599 0.9779962540841840
70 0.0276490620143323 0.0039909703832173 0.7881908390212922 0.9819872244674014
71 0.0268724480963495 0.0035301526829785 0.8150632871176416 0.9855173771503799
72 0.0259206457666001 0.0030822785480393 0.8409839328842418 0.9885996556984192
73 0.0247815535580256 0.0026502677852626 0.8657654864422674 0.9912499234836819
74 0.0234425551658448 0.0022372418926289 0.8892080416081122 0.9934871653763108
75 0.0218905050481195 0.0018465326398648 0.9110985466562317 0.9953336980161755
76 0.0201117138277064 0.0014816908890628 0.9312102604839381 0.9968153889052384
77 0.0180919334955969 0.0011464956586010 0.9493021939795350 0.9979618845638394
78 0.0158163424156431 0.0008449634336744 0.9651185363951780 0.9988068479975138
79 0.0132695301306719 0.0005813577267470 0.9783880665258499 0.9993882057242609
80 0.0104354819699858 0.0003601988912358 0.9888235484958356 0.9997484046154967
81 0.0072975634582501 0.0001862741917361 0.9961211119540858 0.9999346788072329
82 0.0038385045257687 0.0000646481340986 0.9999596164798545 0.9999993269413314
83 0.0000403835201456 0.0000006730586691 1.0000000000000002 1.0000000000000004

ndtfast, nfast = 60 83 nfast/ndtfast = 1.38333

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

1.000000000000 1.029032332991 0.514516166496 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.23667

Minimum barotropic Courant Number = 3.59609615E-02
Maximum barotropic Courant Number = 2.78404323E-01
Maximum Coriolis Courant Number = 5.75131924E-02


GET_STATE - Read state initial conditions, t = 0.0000
(Iter=0001, File: roms_ini_ECCO_Y2010M1.nc, Rec=0001, Index=1)
- free-surface
(Min = 3.42718750E-01 Max = 1.00076342E+00)
- vertically integrated u-momentum component
(Min = -6.13947819E-01 Max = 2.98077428E-01)
- vertically integrated v-momentum component
(Min = -3.51510533E-01 Max = 1.86907780E-01)
- u-momentum component
(Min = -1.25229418E+00 Max = 5.59096216E-01)
- v-momentum component
(Min = -5.26731347E-01 Max = 3.56355535E-01)
- potential temperature
(Min = 1.02650530E+00 Max = 3.00885906E+01)
- salinity
(Min = 3.10281887E+01 Max = 3.56634813E+01)
GET_NGFLD - tidal period
(Min = 4.30819210E+04 Max = 2.38071599E+06)
GET_2DFLD - tidal elevation amplitude
(Min = 6.42401830E-05 Max = 1.93443811E+01)
GET_2DFLD - tidal elevation phase angle
(Min = 8.20617395E-06 Max = 6.28275310E+00)
GET_2DFLD - tidal current inclination angle
(Min = 4.00291105E-06 Max = 6.28318364E+00)
GET_2DFLD - tidal current phase angle
(Min = 6.88083697E-06 Max = 6.28317974E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 2.20110333E-06 Max = 2.21910980E+00)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -4.86747945E-01 Max = 3.78358049E-01)

GET_2DFLD - error while reading variable: sms_time at TIME index = 0

ROMS/TOMS - Input error ............. exit_flag: 2


Elapsed CPU time (seconds):

Thread # 0 CPU: 0.541
Total: 0.541



ncdump the forcing file is:
%% Dimensions:

nc('xi_u') = 156;
nc('eta_u') = 139;
nc('xi_v') = 157;
nc('eta_v') = 138;
nc('xi_rho') = 157;
nc('eta_rho') = 139;
nc('xi_psi') = 156;
nc('eta_psi') = 138;
nc('sms_time') = 124;
nc('shf_time') = 1;
nc('swf_time') = 1;
nc('sst_time') = 1;
nc('srf_time') = 1;
nc('sss_time') = 1;
nc('wind_time') = 124;
nc('tide_period') = 10;

%% Variables and attributes:

nc{'sms_time'} = ncdouble('sms_time'); %% 124 elements.
nc{'sms_time'}.long_name = ncchar(''surface momentum stress time'');
nc{'sms_time'}.units = ncchar(''days'');
nc{'sms_time'}.cycle_length = ncdouble(0);

nc{'shf_time'} = ncdouble('shf_time'); %% 1 element.
nc{'shf_time'}.long_name = ncchar(''surface heat flux time'');
nc{'shf_time'}.units = ncchar(''days'');
nc{'shf_time'}.cycle_length = ncdouble(0);

nc{'swf_time'} = ncdouble('swf_time'); %% 1 element.
nc{'swf_time'}.long_name = ncchar(''surface freshwater flux time'');
nc{'swf_time'}.units = ncchar(''days'');
nc{'swf_time'}.cycle_length = ncdouble(0);

nc{'sst_time'} = ncdouble('sst_time'); %% 1 element.
nc{'sst_time'}.long_name = ncchar(''sea surface temperature time'');
nc{'sst_time'}.units = ncchar(''days'');
nc{'sst_time'}.cycle_length = ncdouble(0);

nc{'sss_time'} = ncdouble('sss_time'); %% 1 element.
nc{'sss_time'}.long_name = ncchar(''sea surface salinity time'');
nc{'sss_time'}.units = ncchar(''days'');
nc{'sss_time'}.cycle_length = ncdouble(0);

nc{'srf_time'} = ncdouble('srf_time'); %% 1 element.
nc{'srf_time'}.long_name = ncchar(''solar shortwave radiation time'');
nc{'srf_time'}.units = ncchar(''days'');
nc{'srf_time'}.cycle_length = ncdouble(0);

nc{'wind_time'} = ncdouble('wind_time'); %% 124 elements.
nc{'wind_time'}.long_name = ncchar(''surface wind time'');
nc{'wind_time'}.units = ncchar(''days'');
nc{'wind_time'}.cycle_length = ncdouble(0);

nc{'sustr'} = ncdouble('sms_time', 'eta_u', 'xi_u'); %% 2688816 elements.
nc{'sustr'}.long_name = ncchar(''surface u-momentum stress'');
nc{'sustr'}.units = ncchar(''Newton meter-2'');

nc{'svstr'} = ncdouble('sms_time', 'eta_v', 'xi_v'); %% 2686584 elements.
nc{'svstr'}.long_name = ncchar(''surface v-momentum stress'');
nc{'svstr'}.units = ncchar(''Newton meter-2'');

nc{'Uwind'} = ncdouble('wind_time', 'eta_u', 'xi_u'); %% 2688816 elements.
nc{'Uwind'}.long_name = ncchar(''surface u component'');
nc{'Uwind'}.units = ncchar(''Meter day-1'');

nc{'Vwind'} = ncdouble('wind_time', 'eta_v', 'xi_v'); %% 2686584 elements.
nc{'Vwind'}.long_name = ncchar(''surface v component'');
nc{'Vwind'}.units = ncchar(''Meter day-1'');

nc{'shflux'} = ncdouble('shf_time', 'eta_rho', 'xi_rho'); %% 21823 elements.
nc{'shflux'}.long_name = ncchar(''surface net heat flux'');
nc{'shflux'}.units = ncchar(''Watts meter-2'');

nc{'swflux'} = ncdouble('swf_time', 'eta_rho', 'xi_rho'); %% 21823 elements.
nc{'swflux'}.long_name = ncchar(''surface freshwater flux (E-P)'');
nc{'swflux'}.units = ncchar(''centimeter day-1'');
nc{'swflux'}.positive = ncchar(''net evaporation'');
nc{'swflux'}.negative = ncchar(''net precipitation'');

nc{'SST'} = ncdouble('sst_time', 'eta_rho', 'xi_rho'); %% 21823 elements.
nc{'SST'}.long_name = ncchar(''sea surface temperature'');
nc{'SST'}.units = ncchar(''Celsius'');

nc{'SSS'} = ncdouble('sss_time', 'eta_rho', 'xi_rho'); %% 21823 elements.
nc{'SSS'}.long_name = ncchar(''sea surface salinity'');
nc{'SSS'}.units = ncchar(''PSU'');

nc{'dQdSST'} = ncdouble('sst_time', 'eta_rho', 'xi_rho'); %% 21823 elements.
nc{'dQdSST'}.long_name = ncchar(''surface net heat flux sensitivity to SST'');
nc{'dQdSST'}.units = ncchar(''Watts meter-2 Celsius-1'');

nc{'swrad'} = ncdouble('srf_time', 'eta_rho', 'xi_rho'); %% 21823 elements.
nc{'swrad'}.long_name = ncchar(''solar shortwave radiation'');
nc{'swrad'}.units = ncchar(''Watts meter-2'');
nc{'swrad'}.positive = ncchar(''downward flux, heating'');
nc{'swrad'}.negative = ncchar(''upward flux, cooling'');

nc{'tide_period'} = ncdouble('tide_period'); %% 10 elements.
nc{'tide_period'}.long_name = ncchar(''Tide angular period'');
nc{'tide_period'}.units = ncchar(''Hours'');

nc{'tide_Ephase'} = ncdouble('tide_period', 'eta_rho', 'xi_rho'); %% 218230 elements.
nc{'tide_Ephase'}.long_name = ncchar(''Tidal elevation phase angle'');
nc{'tide_Ephase'}.units = ncchar(''Degrees'');

nc{'tide_Eamp'} = ncdouble('tide_period', 'eta_rho', 'xi_rho'); %% 218230 elements.
nc{'tide_Eamp'}.long_name = ncchar(''Tidal elevation amplitude'');
nc{'tide_Eamp'}.units = ncchar(''Meter'');

nc{'tide_Cmin'} = ncdouble('tide_period', 'eta_rho', 'xi_rho'); %% 218230 elements.
nc{'tide_Cmin'}.long_name = ncchar(''Tidal current ellipse semi-minor axis'');
nc{'tide_Cmin'}.units = ncchar(''Meter second-1'');

nc{'tide_Cmax'} = ncdouble('tide_period', 'eta_rho', 'xi_rho'); %% 218230 elements.
nc{'tide_Cmax'}.long_name = ncchar(''Tidal current, ellipse semi-major axis'');
nc{'tide_Cmax'}.units = ncchar(''Meter second-1'');

nc{'tide_Cangle'} = ncdouble('tide_period', 'eta_rho', 'xi_rho'); %% 218230 elements.
nc{'tide_Cangle'}.long_name = ncchar(''Tidal current inclination angle'');
nc{'tide_Cangle'}.units = ncchar(''Degrees between semi-major axis and East'');

nc{'tide_Cphase'} = ncdouble('tide_period', 'eta_rho', 'xi_rho'); %% 218230 elements.
nc{'tide_Cphase'}.long_name = ncchar(''Tidal current phase angle'');
nc{'tide_Cphase'}.units = ncchar(''Degrees'');

endef(nc)
close(nc)


what's the problem? thanks for any help!!

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

Re: GET_2DFLD - error while reading variable: sms_time

#2 Unread post by arango »

Do you really want us to tell you what is the problem? Really, you are using a version of the code that was released on Jun 30, 2003. That's a code that it is 8 years old. ROMS have changed and evolved substantially since them.

:idea: I highly recommend to use at much newer version of the code. I assume that you will be able to update but I will assume that you will have a lot of conflicts if you customize the code. It is very unlikely that anyone can help you. The logic is get_2dfld is very complex and has changes a lot through the years...

Post Reply