what's wrong with Roms3.2?

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

what's wrong with Roms3.2?

#1 Unread post by rainbowbluegz »

hi, everyone:
The following is the output information:
Process Information:

Thread # 0 (pid= 13726) is active.

Model Input Parameters: ROMS/TOMS version 3.2
Wednesday - September 28, 2011 - 9:56:04 AM
-----------------------------------------------------------------------------

SCS, 1/6 degree resolution

Operating system : Linux
CPU/hardware : i686
Compiler system : ifort
Compiler command : /opt/intel/Compiler/11.1/073/bin/ia32/ifort
Compiler flags : -heap-arrays -fp-model precise -g -check uninit -ftrapuv -traceback -free

SVN Root URL :
SVN Revision :
Local Root : /home/xie/roms-3.3
Header Dir : /home/xie/roms3_3_project/scs
Header file : scs.h
Analytical Dir: /home/xie/roms3_3_project/scs

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.
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.
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).
180 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.
144 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
5.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
5.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
0.0000E+00 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.
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.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 Vtransform S-coordinate transformation equation.
1 Vstretching S-coordinate stretching function.
2.5000E+00 theta_s S-coordinate surface control parameter.
8.0000E-01 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)
0.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
0.0000E+00 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
0.0000E+00 Znudg Nudging/relaxation time scale (days)
for free-surface.
0.0000E+00 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
0.0000E+00 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
0.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.nc
Input Forcing File 01: roms_frc.nc
Input Climatology File: roms_clm.nc

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

tile Istr Iend Jstr Jend Npts

0 1 155 1 137 594580

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

tile Xmin Xmax Ymin Ymax grid

0 -0.50 156.50 -0.50 138.50 RHO-points

0 -0.50 156.50 -0.50 138.50 U-points

0 -0.50 156.50 -0.50 138.50 V-points

Activated C-preprocessing Options:

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.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged fields.
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.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
M2CLIMATOLOGY Processing 2D momentum climatology data.
M3CLIMATOLOGY Processing 3D momentum climatology data.
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.
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.
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.
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.
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.
SPLINES Conservative parabolic spline reconstruction.
SPONGE Enhanced horizontal mixing in the sponge areas.
SSH_TIDES Add tidal elevation to SSH climatology.
TCLIMATOLOGY Processing tracer climatology data.
TS_A4HADVECTION Fourth-order Akima horizontal advection of tracers.
TS_A4VADVECTION Fourth-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 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.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
WESTERN_WALL Wall boundary at Western edge.
ZCLIMATOLOGY Processing sea surface height climatology data.

INITIAL: Configuring 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.0156922 -1.769 -40.215 -78.661
26 -0.0714286 -0.0334181 -3.722 -85.596 -167.470
25 -0.1071429 -0.0534230 -5.880 -136.766 -267.652
24 -0.1428571 -0.0759436 -8.263 -194.325 -380.387
23 -0.1785714 -0.1011932 -10.893 -258.816 -506.740
22 -0.2142857 -0.1293455 -13.784 -330.680 -647.577
21 -0.2500000 -0.1605159 -16.946 -410.210 -803.474
20 -0.2857143 -0.1947408 -20.384 -497.499 -974.614
19 -0.3214286 -0.2319599 -24.091 -592.393 -1160.694
18 -0.3571429 -0.2720013 -28.052 -694.455 -1360.858
17 -0.3928571 -0.3145746 -32.240 -802.948 -1573.656
16 -0.4285714 -0.3592742 -36.620 -916.842 -1797.064
15 -0.4642857 -0.4055940 -41.146 -1034.852 -2028.557
14 -0.5000000 -0.4529542 -45.766 -1155.504 -2265.242
13 -0.5357143 -0.5007368 -50.423 -1277.229 -2504.034
12 -0.5714286 -0.5483273 -55.064 -1398.466 -2741.868
11 -0.6071429 -0.5951560 -59.635 -1517.768 -2975.900
10 -0.6428571 -0.6407336 -64.095 -1633.892 -3203.689
9 -0.6785714 -0.6846786 -68.407 -1745.869 -3423.332
8 -0.7142857 -0.7267318 -72.549 -1853.042 -3633.535
7 -0.7500000 -0.7667600 -76.508 -1955.070 -3833.632
6 -0.7857143 -0.8047488 -80.285 -2051.919 -4023.553
5 -0.8214286 -0.8407893 -83.885 -2143.819 -4203.753
4 -0.8571429 -0.8750596 -87.327 -2231.223 -4375.119
3 -0.8928571 -0.9078056 -90.631 -2314.755 -4538.878
2 -0.9285714 -0.9393221 -93.825 -2395.164 -4696.503
1 -0.9642857 -0.9699374 -96.937 -2473.284 -4849.631
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.0010269040284716 0.0333392336257192
3-0.0009566069754562 0.0166837817338079-0.0019835110039279 0.0500230153595271
4-0.0012051795928536 0.0166997251833988-0.0031886905967814 0.0667227405429259
5-0.0014186058436584 0.0167198115099464-0.0046072964404399 0.0834425520528722
6-0.0015968888560238 0.0167434549406740-0.0062041852964637 0.1001860069935462
7-0.0017400342738662 0.0167700697549411-0.0079442195703299 0.1169560767484873
8-0.0018480513492362 0.0167990703261722-0.0097922709195661 0.1337551470746595
9-0.0019209542333024 0.0168298711819928-0.0117132251528686 0.1505850182566522
10-0.0019587634659484 0.0168618870858811-0.0136719886188169 0.1674469053425334
11-0.0019615076639817 0.0168945331436469-0.0156334962827987 0.1843414384861803
12-0.0019292254079571 0.0169272249380466-0.0175627216907558 0.2012686634242269
13-0.0018619673276114 0.0169593786948459-0.0194246890183672 0.2182280421190729
14-0.0017597983859117 0.0169904114836395-0.0211844874042789 0.2352184536027123
15-0.0016228003617172 0.0170197414567380-0.0228072877659961 0.2522381950594503
16-0.0014510745310524 0.0170467881294333-0.0242583622970485 0.2692849831888836
17-0.0012447445469948 0.0170709727049508-0.0255031068440433 0.2863559558938343
18-0.0010039595181742 0.0170917184474007-0.0265070663622175 0.3034476743412350
19-0.0007288972858858 0.0171084511060370-0.0272359636481033 0.3205561254472720
20-0.0004197678998158 0.0171205993941351-0.0276557315479191 0.3376767248414071
21-0.0000768172923796 0.0171275955257987-0.0277325488402987 0.3548043203672058
22 0.0002996688483262 0.0171288758140050-0.0274328799919725 0.3719331961812107
23 0.0007093610069602 0.0171238813331995-0.0267235189850123 0.3890570775144103
24 0.0011518815707581 0.0171120586497502-0.0255716374142542 0.4061691361641604
25 0.0016268003607435 0.0170928606235709-0.0239448370535106 0.4232619967877314
26 0.0021336299643244 0.0170657472842252-0.0218112070891863 0.4403277440719565
27 0.0026718208692771 0.0170301867848198-0.0191393862199091 0.4573579308567763
28 0.0032407563991181 0.0169856564369985-0.0158986298207911 0.4743435872937748
29 0.0038397474498617 0.0169316438303465-0.0120588823709294 0.4912752311241213
30 0.0044680270281659 0.0168676480395155-0.0075908553427635 0.5081428791636368
31 0.0051247445908649 0.0167931809223794-0.0024661107518986 0.5249360600860161
32 0.0058089601858883 0.0167077685125316 0.0033428494339897 0.5416438285985478
33 0.0065196383945684 0.0166109525094335 0.0098624878285581 0.5582547811079812
34 0.0072556420753336 0.0165022918695240 0.0171181299038917 0.5747570729775052
35 0.0080157259087899 0.0163813645016018 0.0251338558126816 0.5911384374791070
36 0.0087985297441882 0.0162477690697886 0.0339323855568698 0.6073862065488956
37 0.0096025717472802 0.0161011269073855 0.0435349573041500 0.6234873334562812
38 0.0104262413495604 0.0159410840449308 0.0539611986537104 0.6394284175012119
39 0.0112677919988951 0.0157673133557715 0.0652289906526054 0.6551957308569835
40 0.0121253337115396 0.0155795168224566 0.0773543243641450 0.6707752476794401
41 0.0129968254255413 0.0153774279272642 0.0903511497896863 0.6861526756067043
42 0.0138800671555305 0.0151608141701719 0.1042312169452169 0.7013134897768762
43 0.0147726919488984 0.0149294797175797 0.1190039088941153 0.7162429694944560
44 0.0156721576433617 0.0146832681850981 0.1346760665374769 0.7309262376795540
45 0.0165757384259147 0.0144220655577087 0.1512518049633916 0.7453483032372628
46 0.0174805161931688 0.0141458032506101 0.1687323211565604 0.7594941064878729
47 0.0183833717130782 0.0138544613140573 0.1871156928696386 0.7733485678019302
48 0.0192809755880535 0.0135480717855060 0.2063966684576921 0.7868966395874362
49 0.0201697790194623 0.0132267221923718 0.2265664474771545 0.8001233617798079
50 0.0210460043735164 0.0128905592087141 0.2476124518506709 0.8130139209885220
51 0.0219056355485466 0.0125397924691555 0.2695180873992175 0.8255537134576775
52 0.0227444081436646 0.0121746985433464 0.2922624955428820 0.8377284120010239
53 0.0235577994288118 0.0117956250742853 0.3158202949716938 0.8495240370753092
54 0.0243410181161952 0.0114029950838051 0.3401613130878890 0.8609270321591143
55 0.0250889939331104 0.0109973114485352 0.3652503070209994 0.8719243436076495
56 0.0257963669961522 0.0105791615496500 0.3910466740171517 0.8825035051572995
57 0.0264574769868114 0.0101492220997141 0.4175041510039630 0.8926527272570136
58 0.0270663521284594 0.0097082641499339 0.4445705031324225 0.9023609914069476
59 0.0276166979647200 0.0092571582811263 0.4721872010971425 0.9116181496880739
60 0.0281018859392274 0.0087968799817143 0.5002890870363698 0.9204150296697882
61 0.0285149417767724 0.0083285152160605 0.5288040288131421 0.9287435448858486
62 0.0288485336658351 0.0078532661864476 0.5576525624789772 0.9365968110722963
63 0.0290949602425050 0.0073724572920170 0.5867475227214822 0.9439692683643134
64 0.0292461383757876 0.0068875412879753 0.6159936610972698 0.9508568096522887
65 0.0292935907542990 0.0064001056483788 0.6452872518515688 0.9572569153006675
66 0.0292284332743468 0.0059118791358072 0.6745156851259155 0.9631687944364746
67 0.0290413622293988 0.0054247385812347 0.7035570473553143 0.9685935330177093
68 0.0287226413009382 0.0049407158774114 0.7322796886562525 0.9735342488951207
69 0.0282620883507063 0.0044620051890625 0.7605417770069588 0.9779962540841832
70 0.0276490620143322 0.0039909703832173 0.7881908390212909 0.9819872244674006
71 0.0268724480963495 0.0035301526829785 0.8150632871176404 0.9855173771503791
72 0.0259206457666001 0.0030822785480393 0.8409839328842404 0.9885996556984185
73 0.0247815535580255 0.0026502677852627 0.8657654864422659 0.9912499234836811
74 0.0234425551658447 0.0022372418926289 0.8892080416081106 0.9934871653763100
75 0.0218905050481194 0.0018465326398648 0.9110985466562300 0.9953336980161748
76 0.0201117138277064 0.0014816908890628 0.9312102604839364 0.9968153889052376
77 0.0180919334955968 0.0011464956586011 0.9493021939795332 0.9979618845638386
78 0.0158163424156430 0.0008449634336744 0.9651185363951762 0.9988068479975131
79 0.0132695301306719 0.0005813577267471 0.9783880665258481 0.9993882057242601
80 0.0104354819699857 0.0003601988912359 0.9888235484958338 0.9997484046154960
81 0.0072975634582500 0.0001862741917361 0.9961211119540838 0.9999346788072321
82 0.0038385045257686 0.0000646481340986 0.9999596164798524 0.9999993269413306
83 0.0000403835201471 0.0000006730586691 0.9999999999999994 0.9999999999999998

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 X-grid spacing, DXmin = 1.66385662E+01 km
Maximum X-grid spacing, DXmax = 1.84770211E+01 km
Minimum Y-grid spacing, DYmin = 1.66703769E+01 km
Maximum Y-grid spacing, DYmax = 1.84751077E+01 km
Minimum Z-grid spacing, DZmin = 1.76944522E+00 m
Maximum Z-grid spacing, DZmax = 2.38792282E+02 m

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


NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(File: roms_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 = 1.03284264E+00 Max = 2.87353796E+01)
- salinity
(Min = 3.13632231E+01 Max = 3.50847621E+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 - surface u-momentum stress, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.69066110E-04 Max = 2.88435117E-06)
GET_2DFLD - surface v-momentum stress, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.31120433E-04 Max = -7.53048503E-06)
GET_2DFLD - surface net heat flux, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -8.74117223E-05 Max = 9.91936039E-06)
GET_2DFLD - surface net freswater flux, (E-P), t = 345 00:00:00
(Rec=0012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -9.18419336E-08 Max = 9.79046117E-08)
GET_2DFLD - sea surface height climatology, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.36335665E-01 Max = 5.20850485E-01)
GET_2DFLD - vertically integrated u-momentum climatologt = 345 00:00:00
(Rec=0012, Index=1, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -5.68134895E+00 Max = 2.74522526E+00)
GET_2DFLD - vertically integrated v-momentum climatologt = 345 00:00:00
(Rec=0012, Index=1, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.54390195E+00 Max = 7.90076624E+00)
GET_3DFLD - potential temperature climatology, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.03284264E+00 Max = 2.92403694E+01)
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
oceanG 081A3E4D get_cycle_ 169 get_cycle.f90
oceanG 0819F1FE get_3dfld_ 207 get_3dfld.f90
oceanG 08138F1E get_data_ 125 get_data.f90
oceanG 0804C287 initial_ 218 initial.f90
oceanG 0804A662 ocean_control_mod 94 ocean_control.f90
oceanG 0804A370 MAIN__ 75 master.f90
oceanG 0804A2E1 Unknown Unknown Unknown
libc.so.6 009AAE9C Unknown Unknown Unknown
oceanG 0804A1D1 Unknown Unknown Unknown

I donot know what cause the problem and how to solve.
Thanks any help!

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: what's wrong with Roms3.2?

#2 Unread post by kate »

rainbowbluegz wrote: GET_3DFLD - potential temperature climatology, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.03284264E+00 Max = 2.92403694E+01)
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
oceanG 081A3E4D get_cycle_ 169 get_cycle.f90
oceanG 0819F1FE get_3dfld_ 207 get_3dfld.f90
oceanG 08138F1E get_data_ 125 get_data.f90
oceanG 0804C287 initial_ 218 initial.f90
oceanG 0804A662 ocean_control_mod 94 ocean_control.f90
oceanG 0804A370 MAIN__ 75 master.f90
oceanG 0804A2E1 Unknown Unknown Unknown
libc.so.6 009AAE9C Unknown Unknown Unknown
oceanG 0804A1D1 Unknown Unknown Unknown

I donot know what cause the problem and how to solve.
Thanks any help!
It has given you the stack trace and even the line numbers (I think). Which field is get_data.f90 reading on line 125? What does it read after potential temperature climatology? Salt climatology perhaps? It's having trouble in get_cycle, which has to do with the time axis for said variable.

Have you considered trying again with a newer ROMS? Some of the I/O has changed since 3.2 so the bugs have changed too. In any case, a seg fault could be an array out of bounds, so I would recompile with bounds checking - just temporarily for figuring this out.

ymamoutos
Posts: 71
Joined: Fri Nov 19, 2010 2:33 pm
Location: University of Aegean

Re: what's wrong with Roms3.2?

#3 Unread post by ymamoutos »

Hi to everyone,

I've got "similar" problem.My ROMS version is 3.6. Kate you suggest to recompile with bounds checking. How can i do that?

Thanks for any help and for your time.
Attachments
log.txt
(20.5 KiB) Downloaded 223 times

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: what's wrong with Roms3.2?

#4 Unread post by kate »

It will depend on your compiler so you should ideally check its man pages. Mine says:
Debugging Options
-C Add array bounds checking; the same as -Mbounds.
It can give you info such as this:

Code: Select all

0: Subscript out of range for array a (tri_exchange_2d.f90: 501)
    subscript=42, lower bound=78, upper bound=123, dimension=1
0: Subscript out of range for array a (tri_exchange_2d.f90: 501)
    subscript=42, lower bound=118, upper bound=164, dimension=1

ymamoutos
Posts: 71
Joined: Fri Nov 19, 2010 2:33 pm
Location: University of Aegean

Re: what's wrong with Roms3.2?

#5 Unread post by ymamoutos »

Thanks for your advise Kate. My compiler is gfortran 4.4.3 and i used -fbounds-check flag.
The result for my current cpp definitions file (myapp.h) is:

Code: Select all

At line 74 of file get_data.f90
Fortran runtime error: Array reference out of bounds for array 'frc', upper bound of dimension 1 exceeded (1 > 0)
If i undefine the ANA_BSFLUX, ANA_BTFLUX, ANA_SMFLUX and ANA_STFLUX
because i don't use an external forcing file i get again runtime error for frc array:

Code: Select all

At line 2156 of file read_phypar.f90
Fortran runtime error: Array reference out of bounds for array 'frc', upper bound of dimension 1 exceeded (1 > 0)
Any suggestion/help would be appreciated. Thanks again for your time and your help.

Yiannis
Attachments
ocean.in
(97.07 KiB) Downloaded 202 times
myapp.h
(921 Bytes) Downloaded 210 times

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: what's wrong with Roms3.2?

#6 Unread post by kate »

ymamoutos wrote:

Code: Select all

At line 74 of file get_data.f90
Fortran runtime error: Array reference out of bounds for array 'frc', upper bound of dimension 1 exceeded (1 > 0)
So this is with ANA_BSFLUX and friends? What is it trying to read on line 74 of get_data.f90?
If i undefine the ANA_BSFLUX, ANA_BTFLUX, ANA_SMFLUX and ANA_STFLUX
because i don't use an external forcing file i get again runtime error for frc array:
You want these all #defined, I believe. You need to provide a bottom salinity flux one way or another - most of us analytically set it to zero. Same for the rest, unless we're explicitly providing the surface fluxes.

ymamoutos
Posts: 71
Joined: Fri Nov 19, 2010 2:33 pm
Location: University of Aegean

Re: what's wrong with Roms3.2?

#7 Unread post by ymamoutos »

Kate hi,

on line 74 get_data.f90 calls subroutine get_2dfld for surface wind stress components.
You need to provide a bottom salinity flux one way or another - most of us analytically set it to zero.
So if am not providing an external file i have to modify ana_btflux.h and the other files if am understanding correct.

Thanks again for your patience and help.

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: what's wrong with Roms3.2?

#8 Unread post by kate »

Yes, you have to set ANA_SMFLUX to tell it not to read surface momentum flux from a file. For these common ones, you can probably just leave it be the default of zero. If you aren't sure, you can read the analytical.f90 file to see what got picked up during the C preprocessor phase.

ymamoutos
Posts: 71
Joined: Fri Nov 19, 2010 2:33 pm
Location: University of Aegean

Re: what's wrong with Roms3.2?

#9 Unread post by ymamoutos »

Hi Kate and greetings from Greece,

My problem is solved.I modified the ana_*.h files
that i wanted for my test case following your suggestion
(i set the variables that i needed equal to zero) and
my model configuration is running smoothly until now.

Thanks again for your help and have a pleasant evening.

Yiannis

Post Reply