ROMS with monthly levitus climatology..

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
jvratnam
Posts: 10
Joined: Thu Oct 13, 2005 4:26 pm
Location: Frontier Research Center for Global Change

ROMS with monthly levitus climatology..

#1 Unread post by jvratnam »

Hi
I made a run of nearly 10 years with Levitus annual climatology. Now, when I am trying to use Levitus monthly climatology for further runs, the model is blowing up. Can someone please suggest how to get around this problem. I tried changing the time step but still the model blows up.

Thanks
Ratnam

User avatar
m.hadfield
Posts: 521
Joined: Tue Jul 01, 2003 4:12 am
Location: NIWA

#2 Unread post by m.hadfield »

I'm afraid that's nowhere near enough information for anyone else to daignose your problem.

User avatar
m.hadfield
Posts: 521
Joined: Tue Jul 01, 2003 4:12 am
Location: NIWA

#3 Unread post by m.hadfield »

Or even to diagnose it.

jvratnam
Posts: 10
Joined: Thu Oct 13, 2005 4:26 pm
Location: Frontier Research Center for Global Change

#4 Unread post by jvratnam »

Dear Dr. Hadfield

I am running ROMS for the domain 20s-40N, 40E-120E with a resolution of 60Km. When I change from annual levitus to monthly levitus I am getting the following output:

Process Information:

Thread # 0 (pid= 21443) is active.

Model Input Parameters: ROMS/TOMS version 2.2
Tuesday - June 19, 2007 - 9:30:34 AM
-----------------------------------------------------------------------------

ROMS/TOMS 2.0 - Indian Ocean , 60km Resolution

Operating system : Linux
CPU/hardware : i686
Compiler system : g95
Compiler command : g95
Compiler flags : -O3 -ffast-math

Resolution, Grid 01: 0136x0108x024, Parallel Threads: 1, Tiling: 004x004


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

1100000 ntimes Number of timesteps for 3-D equations.
1800.000 dt Timestep size (s) for 3-D equations.
30 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.
-1 nrrec Number of restart records to read from disk.
T LcycleRST Switch to recycle time-records in restart file.
48 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).
480 nHIS Number of timesteps between the writing fields
into history file.
480 ndefHIS Number of timesteps between creation of new
history files.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
480 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
480 ndefAVG Number of timesteps between creation of new
time-averaged file.
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.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
0.0000E+00 rdrg2 Quadratic bottom drag coefficient.
0.0000E+00 Zob Bottom roughness (m).
1 lmd_Jwt Jerlov water type.
5.0000E+00 theta_s S-coordinate surface control parameter.
4.0000E-01 theta_b S-coordinate bottom control parameter.
200.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
1.000 dstart Time-stamp assigned to model initialization (days).
-1.00 time_ref Reference time for units attribute (yyyymmdd.dd)
6.0000E+01 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
6.0000E+01 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(idTvar) Write out tracer 01: temp
T Hout(idTvar) Write out tracer 02: salt

Output/Input Files:

Output Restart File: IN_RST_60km.nc
Prefix for History Files: IN_HIS_60km
Prefix for Averages Files: IN_AVG_60km
Input Grid File: IND_OCN.nc
Input Nonlinear Initial File: IN_RST_s_60km.nc
Input Forcing File 01: IN_forc_60km.nc
IO Variable Information File: varinfo.dat

Tile partition information for Grid 01: 0136x0108x0024 tiling: 004x004

tile Istr Iend Jstr Jend Npts

0 1 34 1 27 22032
1 35 68 1 27 22032
2 69 102 1 27 22032
3 103 136 1 27 22032
4 1 34 28 54 22032
5 35 68 28 54 22032
6 69 102 28 54 22032
7 103 136 28 54 22032
8 1 34 55 81 22032
9 35 68 55 81 22032
10 69 102 55 81 22032
11 103 136 55 81 22032
12 1 34 82 108 22032
13 35 68 82 108 22032
14 69 102 82 108 22032
15 103 136 82 108 22032

Activated C-preprocessing Options:

ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
AVERAGES Writing out time-averaged fields.
AVERAGES_AKS Writing out time-averaged vertical S-diffusion.
AVERAGES_AKT Writing out time-averaged vertical T-diffusion.
CURVGRID Orthogonal curvilinear grid.
DAMEE_4 North Atlantic DAMEE Application, Grid #4.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
EASTERN_WALL Wall boundary at Eastern edge.
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.
MASKING Land/Sea masking.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
NORTHERN_WALL Wall boundary at Northern edge.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
QCORRECTION Surface net heat flux correction.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SRELAXATION Surface salinity relaxation as surface flux.
SOLVE3D Solving 3D Primitive Equations.
SOUTHERN_WALL Wall boundary at Southern edge.
SPLINES Conservative parabolic spline reconstruction.
TCLIMATOLOGY Processing tracer climatology data.
TCLM_NUDGING Nudging toward tracer climatology.
TS_U3HADVECTION Third-order upstream bias horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection 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.
VAR_RHO_2D Variable density barotropic mode.
WESTERN_WALL Wall boundary at Western edge.

INITIAL: Configurating and initializing forward nonlinear model ...


Vertical S-coordinate System:

level S-coord Cs-curve at_hmin over_slope at_hmax

24 0.0000000 -0.0000000 0.000 -0.000 -0.000
23 -0.0416667 -0.0030852 -0.417 -8.094 -15.772
22 -0.0833333 -0.0069437 -0.833 -18.113 -35.393
21 -0.1250000 -0.0119914 -1.250 -31.091 -60.932
20 -0.1666667 -0.0187972 -1.667 -48.444 -95.221
19 -0.2083333 -0.0281303 -2.083 -72.086 -142.089
18 -0.2500000 -0.0409944 -2.500 -104.516 -206.531
17 -0.2916667 -0.0586177 -2.917 -148.788 -294.660
16 -0.3333333 -0.0823381 -3.333 -208.234 -413.134
15 -0.3750000 -0.1133185 -3.750 -285.746 -567.742
14 -0.4166667 -0.1520732 -4.167 -382.604 -761.042
13 -0.4583333 -0.1979504 -4.583 -497.188 -989.792
12 -0.5000000 -0.2489214 -5.000 -624.447 -1243.894
11 -0.5416667 -0.3020234 -5.417 -757.009 -1508.602
10 -0.5833333 -0.3543864 -5.833 -887.732 -1769.632
9 -0.6250000 -0.4042644 -6.250 -1012.272 -2018.293
8 -0.6666667 -0.4514899 -6.667 -1130.210 -2253.754
7 -0.7083333 -0.4972851 -7.083 -1244.589 -2482.095
6 -0.7500000 -0.5437741 -7.500 -1360.695 -2713.890
5 -0.7916667 -0.5935526 -7.917 -1484.987 -2962.056
4 -0.8333333 -0.6494565 -8.333 -1624.521 -3240.710
3 -0.8750000 -0.7145184 -8.750 -1786.846 -3564.942
2 -0.9166667 -0.7920448 -9.167 -1980.189 -3951.211
1 -0.9583333 -0.8857620 -9.583 -2213.823 -4418.063
0 -1.0000000 -1.0000000 -10.000 -2498.524 -4987.048

Time Splitting Weights: ndtfast = 30 nfast = 42

Primary Secondary Accumulated to Current Step

1-0.0008094437383769 0.0333333333333333-0.0008094437383769 0.0333333333333333
2-0.0014053566728197 0.0333603147912792-0.0022148004111966 0.0666936481246126
3-0.0017877524645903 0.0334071600137065-0.0040025528757869 0.1001008081383191
4-0.0019566842408176 0.0334667517625262-0.0059592371166046 0.1335675599008453
5-0.0019122901320372 0.0335319745705535-0.0078715272486418 0.1670995344713988
6-0.0016548570247459 0.0335957175749547-0.0095263842733877 0.2006952520463535
7-0.0011849025289723 0.0336508794757796-0.0107112868023601 0.2343461315221331
8-0.0005032751608632 0.0336903762267453-0.0112145619632232 0.2680365077488784
9 0.0003887272597151 0.0337071520654408-0.0108258347035082 0.3017436598143192
10 0.0014892209965583 0.0336941944901169-0.0093366137069499 0.3354378543044362
11 0.0027955815694920 0.0336445537902317-0.0065410321374578 0.3690824080946679
12 0.0043042707117221 0.0335513677379153-0.0022367614257357 0.4026337758325831
13 0.0060106451121704 0.0334078920475245 0.0037738836864347 0.4360416678801076
14 0.0079087469427945 0.0332075372104522 0.0116826306292292 0.4692492050905598
15 0.0099910761708919 0.0329439123123590 0.0216737068001212 0.5021931174029188
16 0.0122483446563884 0.0326108764399960 0.0339220514565096 0.5348039938429148
17 0.0146692120341107 0.0322025982847830 0.0485912634906203 0.5670065921276978
18 0.0172400033810439 0.0317136245503127 0.0658312668716642 0.5987202166780105
19 0.0199444086685725 0.0311389577709445 0.0857756755402366 0.6298591744489550
20 0.0227631639997064 0.0304741441486588 0.1085388395399430 0.6603333185976138
21 0.0256737146312910 0.0297153720153352 0.1342125541712341 0.6900486906129490
22 0.0286498597812016 0.0288595815276255 0.1628624139524357 0.7189082721405746
23 0.0316613792205219 0.0279045862015855 0.1945237931729576 0.7468128583421600
24 0.0346736416507075 0.0268492068942348 0.2291974348236651 0.7736620652363948
25 0.0376471948657328 0.0256934188392112 0.2668446296893979 0.7993554840756060
26 0.0405373376992232 0.0244385123436867 0.3073819673886211 0.8237939964192927
27 0.0432936737565710 0.0230872677537126 0.3506756411451921 0.8468812641730054
28 0.0458596469320356 0.0216441452951603 0.3965352880772277 0.8685254094681656
29 0.0481720587108284 0.0201154903974257 0.4447073467880561 0.8886408998655914
30 0.0501605672561820 0.0185097551070648 0.4948679140442381 0.9071506549726561
31 0.0517471682814030 0.0168377361985254 0.5466150823256410 0.9239883911711815
32 0.0528456577069106 0.0151128305891453 0.5994607400325517 0.9391012217603267
33 0.0533610761022577 0.0133513086655816 0.6528218161348094 0.9524525304259084
34 0.0531891349131380 0.0115726061288397 0.7060109510479473 0.9640251365547481
35 0.0522156244733761 0.0097996349650684 0.7582265755213233 0.9738247715198165
36 0.0503158038019031 0.0080591141492892 0.8085423793232264 0.9818838856691057
37 0.0473537721847154 0.0063819206892258 0.8558961515079417 0.9882658063583315
38 0.0431818225418188 0.0048034616164019 0.8990779740497605 0.9930692679747334
39 0.0376397765791564 0.0033640675316746 0.9367177506289169 0.9964333355064080
40 0.0305543017255207 0.0021094083123694 0.9672720523544376 0.9985427438187774
41 0.0217382098544504 0.0010909315881854 0.9890102622088880 0.9996336754069628
42 0.0109897377911118 0.0003663245930371 0.9999999999999998 0.9999999999999999

ndtfast, nfast = 30 42 nfast/ndtfast = 1.40000

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

1.000000000000 1.047601458608 0.523800729304 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.18933

Minimum barotropic Courant Number = 1.57296216E-02
Maximum barotropic Courant Number = 3.04070270E-01
Maximum Coriolis Courant Number = 1.69537754E-01


GET_STATE - Read state initial conditions, t = 3103.0000
(File: IN_RST_s_60km.nc, Rec=0002, Index=1)
- free-surface
(Min = -3.55579064E-01 Max = 2.65218695E-01)
- vertically integrated u-momentum component
(Min = -3.74503450E-01 Max = 5.53708394E-01)
- vertically integrated v-momentum component
(Min = -2.61824144E-01 Max = 9.56240395E-01)
- u-momentum component
(Min = -1.15985571E+00 Max = 1.92306246E+00)
- v-momentum component
(Min = -1.73017521E+00 Max = 1.73017521E+00)
- potential temperature
(Min = -2.77201824E-03 Max = 3.63656878E+01)
- salinity
(Min = -1.14554432E-04 Max = 3.85475131E+01)
GET_2DFLD - surface u-momentum stress, t = 195.0000
(File: IN_forc_60km.nc, Rec=0007, Index=1)
(Min = -1.22533484E-04 Max = 2.10668168E-04)
GET_2DFLD - surface v-momentum stress, t = 195.0000
(File: IN_forc_60km.nc, Rec=0007, Index=1)
(Min = -3.00937855E-05 Max = 1.99176829E-04)
GET_2DFLD - solar shortwave radiation flux, t = 195.0000
(File: IN_forc_60km.nc, Rec=0007, Index=1)
(Min = 3.31832774E-05 Max = 9.48760414E-05)
GET_2DFLD - surface net heat flux, t = 195.0000
(File: IN_forc_60km.nc, Rec=0007, Index=1)
(Min = -3.80253898E-05 Max = 6.19588329E-05)
GET_2DFLD - sea surface temperature climatology, t = 195.0000
(File: IN_forc_60km.nc, Rec=0007, Index=1)
(Min = 2.16377163E+01 Max = 3.27661476E+01)
GET_2DFLD - surface net heat flux sensitivity to SST, t = 195.0000
(File: IN_forc_60km.nc, Rec=0007, Index=1)
(Min = -2.30225741E-05 Max = -6.84274295E-06)
GET_2DFLD - sea surface salinity climatology, t = 195.0000
(File: IN_forc_60km.nc, Rec=0007, Index=1)
(Min = 1.78007412E+01 Max = 4.01702118E+01)
GET_3DFLD - potential temperature climatology, t = 195.0000
(File: IN_Lclm_mon_60km.nc, Rec=0007, Index=1)
(Min = 0.00000000E+00 Max = 3.23528691E+01)
GET_3DFLD - salinity climatology, t = 195.0000
(File: IN_Lclm_mon_60km.nc, Rec=0007, Index=1)
(Min = 0.00000000E+00 Max = 4.15420839E+01)

Maximum grid stiffness ratios: rx0 = 1.832900E-01 (Beckmann and Haidvogel)
rx1 = 4.231252E+00 (Haney)


Initial basin volumes: TotVolume = 1.04222658837077E+17
MinVolume = 1.06710667228717E+09
MaxVolume = 2.36735455036476E+12
Max/Min = 2.21847975637779E+03

NL ROMS/TOMS: started time-stepping:( TimeSteps: 00148897 - 01100000)

GET_2DFLD - surface u-momentum stress, t = 225.0000
(File: IN_forc_60km.nc, Rec=0008, Index=2)
(Min = -1.27057826E-04 Max = 1.69477477E-04)
GET_2DFLD - surface v-momentum stress, t = 225.0000
(File: IN_forc_60km.nc, Rec=0008, Index=2)
(Min = -3.75618186E-05 Max = 1.89059624E-04)
GET_2DFLD - solar shortwave radiation flux, t = 225.0000
(File: IN_forc_60km.nc, Rec=0008, Index=2)
(Min = 3.95470418E-05 Max = 8.92582593E-05)
GET_2DFLD - surface net heat flux, t = 225.0000
(File: IN_forc_60km.nc, Rec=0008, Index=2)
(Min = -3.01855186E-05 Max = 5.09043850E-05)
GET_2DFLD - sea surface temperature climatology, t = 225.0000
(File: IN_forc_60km.nc, Rec=0008, Index=2)
(Min = 2.21497269E+01 Max = 3.36901932E+01)
GET_2DFLD - surface net heat flux sensitivity to SST, t = 225.0000
(File: IN_forc_60km.nc, Rec=0008, Index=2)
(Min = -2.09121310E-05 Max = -8.16135130E-06)
GET_2DFLD - sea surface salinity climatology, t = 225.0000
(File: IN_forc_60km.nc, Rec=0008, Index=2)
(Min = 1.79777184E+01 Max = 3.99104347E+01)
GET_3DFLD - potential temperature climatology, t = 225.0000
(File: IN_Lclm_mon_60km.nc, Rec=0008, Index=2)
(Min = 0.00000000E+00 Max = 3.29817480E+01)
GET_3DFLD - salinity climatology, t = 225.0000
(File: IN_Lclm_mon_60km.nc, Rec=0008, Index=2)
(Min = 0.00000000E+00 Max = 4.09810503E+01)

STEP time[DAYS] KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME trd

148896 3103.000000 2.835926E-03 1.783567E+04 1.783567E+04 1.090724E+17 0
DEF_HIS - creating history file: IN_HIS_60km_0000000311.nc
148897 3103.020833 2.836060E-03 1.783560E+04 1.783560E+04 1.090724E+17 0
148898 3103.041667 2.838052E-03 1.783553E+04 1.783553E+04 1.090724E+17 0
148899 3103.062500 2.850864E-03 1.783545E+04 1.783546E+04 1.090724E+17 0
148900 3103.083333 2.887540E-03 1.783538E+04 1.783539E+04 1.090724E+17 0
148901 3103.104167 2.966744E-03 1.783532E+04 1.783532E+04 1.090724E+17 0
148902 3103.125000 3.111943E-03 1.783525E+04 1.783526E+04 1.090724E+17 0
148903 3103.145833 3.351165E-03 1.783519E+04 1.783519E+04 1.090724E+17 0
148904 3103.166667 3.716828E-03 1.783513E+04 1.783513E+04 1.090724E+17 0
148905 3103.187500 4.245130E-03 1.783507E+04 1.783507E+04 1.090724E+17 0
148906 3103.208333 4.975982E-03 1.783501E+04 1.783502E+04 1.090724E+17 0
148907 3103.229167 5.952659E-03 1.783495E+04 1.783496E+04 1.090724E+17 0
148908 3103.250000 7.220932E-03 1.783490E+04 1.783491E+04 1.090724E+17 0
148909 3103.270833 8.829156E-03 1.783485E+04 1.783486E+04 1.090724E+17 0
148910 3103.291667 1.082868E-02 1.783481E+04 1.783482E+04 1.090724E+17 0
148911 3103.312500 1.327266E-02 1.783476E+04 1.783477E+04 1.090724E+17 0
148912 3103.333333 1.621659E-02 1.783471E+04 1.783473E+04 1.090724E+17 0
148913 3103.354167 1.971809E-02 1.783467E+04 1.783469E+04 1.090724E+17 0
148914 3103.375000 2.383728E-02 1.783462E+04 1.783464E+04 1.090724E+17 0
148915 3103.395833 2.863725E-02 1.783456E+04 1.783459E+04 1.090724E+17 0
148916 3103.416667 3.418468E-02 1.783451E+04 1.783454E+04 1.090724E+17 0
148917 3103.437500 4.055132E-02 1.783444E+04 1.783448E+04 1.090724E+17 0
148918 3103.458333 4.781575E-02 1.783438E+04 1.783442E+04 1.090724E+17 0
148919 3103.479167 5.606596E-02 1.783431E+04 1.783436E+04 1.090724E+17 0
148920 3103.500000 6.540042E-02 1.783423E+04 1.783430E+04 1.090724E+17 0
148921 3103.520833 7.593576E-02 1.783416E+04 1.783424E+04 1.090724E+17 0
148922 3103.541667 8.781377E-02 1.783409E+04 1.783417E+04 1.090724E+17 0
148923 3103.562500 1.012036E-01 1.783401E+04 1.783411E+04 1.090724E+17 0
148924 3103.583333 1.163092E-01 1.783393E+04 1.783404E+04 1.090724E+17 0
148925 3103.604167 1.333730E-01 1.783385E+04 1.783398E+04 1.090724E+17 0
148926 3103.625000 1.526997E-01 1.783377E+04 1.783392E+04 1.090724E+17 0
148927 3103.645833 1.746556E-01 1.783368E+04 1.783386E+04 1.090724E+17 0
148928 3103.666667 1.997377E-01 1.783360E+04 1.783380E+04 1.090724E+17 0
148929 3103.687500 2.291772E-01 1.783352E+04 1.783374E+04 1.090724E+17 0
148930 3103.708333 2.707614E-01 1.783343E+04 1.783370E+04 1.090724E+17 0
148931 3103.729167 4.051067E-01 1.783334E+04 1.783375E+04 1.090724E+17 0
148932 3103.750000 2.167770E+00 1.783328E+04 1.783545E+04 1.090724E+17 0
148933 3103.770833 2.165873E+02 9.884504E+04 9.906163E+04 1.090724E+17 0
148934 3103.791667 NaN NaN NaN NaN 0


ROMS/TOMS - Blows up ................ exit_flag: 1


Blowing-up: Saving latest model state into RESTART file

WRT_RST - wrote re-start fields (Index=1,1) into time record = 0000001

Elapsed CPU time (seconds):

Thread # 0 CPU: 156.993
Total: 156.993

Nonlinear model elapsed time profile:

Initialization ................................... 0.318 ( 0.2025 %)
Reading of input data ............................ 0.093 ( 0.0592 %)
Processing of input data ......................... 1.033 ( 0.6579 %)
Processing of output time averaged data .......... 3.215 ( 2.0476 %)
Computation of vertical boundary conditions ...... 0.064 ( 0.0408 %)
Computation of global information integrals ...... 2.754 ( 1.7539 %)
Writing of output data ........................... 0.033 ( 0.0210 %)
Model 2D kernel .................................. 63.306 (40.3243 %)
2D/3D coupling, vertical metrics ................. 1.709 ( 1.0884 %)
Omega vertical velocity .......................... 2.241 ( 1.4272 %)
Equation of state for seawater ................... 9.744 ( 6.2063 %)
KPP vertical mixing parameterization ............. 18.354 (11.6911 %)
3D equations right-side terms .................... 4.220 ( 2.6882 %)
3D equations predictor step ...................... 8.347 ( 5.3166 %)
Pressure gradient ................................ 3.952 ( 2.5176 %)
Corrector time-step for 3D momentum .............. 15.365 ( 9.7868 %)
Corrector time-step for tracers .................. 19.720 (12.5611 %)
Total: 154.467 98.3906

ROMS/TOMS - number of time records written in history file, Grid 01, 00000000
number of time records written in restart file, Grid 01, 00000001
number of time records written in averages file, Grid 01, 00000000

MAIN: Abnormal termination: BLOWUP.

ROMS/TOMS: DONE... Tuesday - June 19, 2007 - 9:33:23 AM

Thanks
Ratnam

jvratnam
Posts: 10
Joined: Thu Oct 13, 2005 4:26 pm
Location: Frontier Research Center for Global Change

#5 Unread post by jvratnam »

Hi

I found the problem for the blowing of the model to be in the monthly climatology I created. Actually the problem was in the .hdat files read by the OA program. The .hdat files have levels as 0 9 99.....149999 instead of 0 10 100....150000. I made the change in the .hdat files then the OA program was taking levels properly. Now roms is running properly.

Thanks
Ratnam

Post Reply