[RegCNET] Floating point error in RegCM
Jonathan Winter
jwinter at MIT.EDU
Tue Sep 20 17:50:53 CEST 2005
Jagadish,
I thought I'd reply because I use redhat enterprise linux v3.0 with
regcm3 and it works fine. It's not exactly the same as 9.0, but close.
a few things you might try:
1. boot to the regular kernel, anything "smp" is an acronym for
symmetric multiprocessing. the smp kernel is the default because it
should have no problems with one processor (I'm assuming there's only 1
proc on this machine), but maybe some piece of hardware is just not
playing with the software correctly. I doubt this is the problem, but
it wouldn't hurt to try.
2. for more info on your crash, use the gnu debugger in redhat, it
helps with finding the problem (what subroutine, what matrix, why,
etc.). the commands go:
> gdb ../Main/regcm
> run < regcm.in
you can also add flags in your makefile to make the gdb more sensitive,
return more information. I can send you a makefile with all the flags I
use if you'd like. When it crashes, typing simple questions like
"where", etc. gives you the info. google gnu debugger, there should be
some docs somewhere.
hope this helps,
jonathan
On Tue, 2005-09-20 at 07:18, jagadish karmacharya wrote:
>
> Dear Jeremy and all,
>
> Recently we have upgraded Redhat 9.0 (kernel 2.4.20-6smp) on
> out office computers.
>
> I have tried to run the simualtion with different version
> of RegCM on those mechine but they generate floating point
> exception error all the time.
>
> I have tried all the options i can think of... like decreaing
> the time step, changing the domain size, different the driving
> data, different convection scheme etc but with no success.
>
> Model crashes imidately or after few time step.
>
> I am running the model in different computer with Linux
> Mandrike 10 with same model configuration.
>
> Is this problem related with linux version or something else?
>
> I have attched two messages of simulation for trouble
> shooting.
>
> Any help/ suggestion will be appriciated
>
> Best Reagards
>
> Jagadish
>
>
> __________________________________________________
> Do You Yahoo!?
> Tired of spam? Yahoo! Mail has the best spam protection around
> http://mail.yahoo.com
>
> ______________________________________________________________________
> [root at cc1 exp_blhnew]# ./regcm<./regcm.in
> RESTARTPARAM READ IN
> TIMEPARAM READ IN
> OUTPARAM READ IN
> PHYSICSPARAM READ IN
> SUBEXPARAM READ IN
> GRELLPARAM READ IN
> EMANPARAM READ IN
> CHEMPARAM READ IN
> dtau = 37.50000000000000 75.00000000000000
> NREC = 289276
> IDATE1, IDATE2, dtmin, ktaur = 2000052400 2000080100
> 2.500000000000000 0
> READING HEADER FILE
> DIMS 111 120 18
> DOMAIN 50000.00 21.00000 80.00000 21.00000
> 80.00000 0.0000000E+00
> PROJROTMER
> SIGMA 0.0000000E+00 5.0000001E-02 0.1000000 0.1600000
> 0.2300000 0.3100000 0.3900000 0.4700000 0.5500000
> 0.6300000 0.7100000 0.7800000 0.8400000 0.8900000
> 0.9300000 0.9600000 0.9800000 0.9900000 1.000000
> PTOP 1.000000
> OUTPUT 1 1
> ***** mdate = 2000052400
> input/output parameters
> ifsave = T savfrq = 48.00000000000000 iftape = F tapfrq =
> 6.000000000000000 ifprt = F prtfrq = 12.00000000000000
> kxout = 6 jxsex = 40 radisp =
> 6.000000000000000 batfrq = 3.000000000000000 nslice =
> 120 ifchem = F chemfrq = 6.000000000000000
>
> physical parameterizations
> iboudy = 5 icup = 4 igcc = 2 ipptls =
> 1 iocnflx = 2 ipgf = 0 lakemod = 0
> ichem = 0
>
> model parameters
> radfrq = 30.00000000000000 abatm = 450.0000000000000 abemh =
> 18.00000000000000 dt = 150.0000000000000
>
> ncld = 1
>
> HT
> HTSD
> SATBRT
> XLAT
> XLONG
> MSFX
> MSFD
> F
> SNOWC
>
> ***************************************************
> ***************************************************
> **** RegCM IS BEING RUN ON THE FOLLOWING GRID: ****
> **** Map Projection: ROTMER ****
> **** IX= 111 JX= 120 KX= 18 ****
> **** PTOP= 1.000000 DX= 50000.00 ****
> **** CLAT= 21.00000 CLON= 80.00000 ****
> ***************************************************
>
> Index of highest allowed pbl: kt = 10
>
> AUTO-CONVERSION RATE: LAND= 5.0000000000000000E-004
> OCEAN= 5.0000000000000000E-004
> RELATIVE HUMIDITY THRESHOLDS: LAND= 0.9900000000000000
> OCEAN= 0.9900000000000000
> GULTEPE FACTORS: LAND= 1.000000000000000 OCEAN=
> 1.000000000000000
> MAXIMUM CLOUD COVER FOR RADIATION: 0.8000000000000000
> MAXIMUM RELATIVE HUMIDITY: 1.010000000000000
> rh0 temperature threshold: 238.0000000000000
> Raindrop Evaporation Rate 2.0000000000000000E-005
> Raindrop Accretion Rate 6.000000000000000
>
>
> EMANUEL (1991) CONVECTION V4.3C (20 May, 2002)
> MIN CONVECTION ORIGIN (minsig/orig): 0.9500000000000000 3
> AUTOCONVERSION THERSHOLD (elcrit): 1.0000000000000000E-003
> AUTOCONVERSION THRESHOLD TO ZERO (tlcrit): -55.00000000000000
> ENTRAINMENT COEFFICIENT (entp): 1.500000000000000
> FRACTIONAL AREA OF UNSATURATED DNDRAFT (sigd): 5.0000000000000000E-002
> PRECIP FRACTION OUTSIDE OF CLOUD (sigs): 0.1200000000000000
> FALL SPEED OF RAIN (omtrain): 50.00000000000000
> FALL SPEED OF SNOW (omtsnow): 5.500000000000000
> RAIN EVAPORATION COEFFICIENT (coeffr): 1.000000000000000
> SNOW EVAPORATION COEFFICIENT (coeffs): 0.8000000000000000
> CONVECTIVE MOMENTUM TRANSPORT COEFFICIENT (cu): 0.7000000000000000
> DOWNDRAFT VELOCITY SCALE (betae): 10.00000000000000
> MAX NEGATIVE PERTURBATION BELOW LFC (dtmax): 0.9000000000000000
> QUASI-EQUILIBRIUM APPROACH RATE (alphae): 1.0000000000000000E-002
> QUASI-EQUILIBRIUM APPROACH RATE (damp): 5.0000000000000000E-002
>
>
>
> CONVECTIVE CLOUD FRACTION/WATER
> Maximum Convective Cloud Cover
> before resolution scaling: 0.2500000000000000
> Maximum Convective Cloud Cover
> after resolution scaling: 0.2500000000000000
> Convective Cloud Water: 4.9999998736893758E-005
> the surface energy budget is used to calculate the ground temperature. julday = 145 gmt = 0.0
> the radiation is computed every 12 time steps.
>
> relaxation boudnary conditions (exponential method) are used. fnudge = 0.33333E-03 gnudge = 0.33333E+06
> 0 k sigma(k) a(k) dsigma(k) twt(k,1) twt(k,2) qcon(k)
>
> 1 0.0000 0.0250 0.0500 0.0000 0.0000 0.0000
> 2 0.0500 0.0750 0.0500 0.5000 0.5000 0.5000
> 3 0.1000 0.1300 0.0600 0.4545 0.5455 0.5455
> 4 0.1600 0.1950 0.0700 0.4615 0.5385 0.5385
> 5 0.2300 0.2700 0.0800 0.4667 0.5333 0.5333
> 6 0.3100 0.3500 0.0800 0.5000 0.5000 0.5000
> 7 0.3900 0.4300 0.0800 0.5000 0.5000 0.5000
> 8 0.4700 0.5100 0.0800 0.5000 0.5000 0.5000
> 9 0.5500 0.5900 0.0800 0.5000 0.5000 0.5000
> 10 0.6300 0.6700 0.0800 0.5000 0.5000 0.5000
> 11 0.7100 0.7450 0.0700 0.5333 0.4667 0.4667
> 12 0.7800 0.8100 0.0600 0.5385 0.4615 0.4615
> 13 0.8400 0.8650 0.0500 0.5455 0.4545 0.4545
> 14 0.8900 0.9100 0.0400 0.5556 0.4444 0.4444
> 15 0.9300 0.9450 0.0300 0.5714 0.4286 0.4286
> 16 0.9600 0.9700 0.0200 0.6000 0.4000 0.4000
> 17 0.9800 0.9850 0.0100 0.6667 0.3333 0.3333
> 18 0.9900 0.9950 0.0100 0.5000 0.5000 0.5000
> 19 1.0000
>
>
> maximum time = 69.000 days.
> time step = 150.00 seconds
> dx = 50000. meters
> grid points (x,y) = (120,111)
> number of levels = 18
> constant hor. diff. coef. = 0.25000E+05 m*m/s
> maximum hor. diff. coef. = 0.26042E+06 m*m/s
> READING INITAL CONDITIONS 2000052400
> U
> V
> TA
> QV
> PS
> TS
> *** initial total air = 0.30585E+18 kg, total water = 0.10863E+16 kg in large domain.
> *** solar declination angle = 20.99 degrees.
> ozone profiles
> 0.4963E-07 0.4963E-07 0.4963E-07
> 0.4963E-07 0.4963E-07 0.4963E-07
> 0.4963E-07 0.4963E-07 0.4963E-07
> 0.4963E-07 0.4963E-07 0.4963E-07
> 0.4981E-07 0.4985E-07 0.4985E-07
> 0.5267E-07 0.5267E-07 0.5267E-07
> 0.5367E-07 0.5370E-07 0.5371E-07
> 0.5775E-07 0.5775E-07 0.5775E-07
> 0.6413E-07 0.6417E-07 0.6417E-07
> 0.7402E-07 0.7406E-07 0.7407E-07
> 0.9089E-07 0.9101E-07 0.9104E-07
> 0.1196E-06 0.1197E-06 0.1197E-06
> 0.1658E-06 0.1661E-06 0.1661E-06
> 0.2862E-06 0.2866E-06 0.2867E-06
> 0.5890E-06 0.5896E-06 0.5897E-06
> 0.1156E-05 0.1157E-05 0.1157E-05
> 0.2574E-05 0.2576E-05 0.2576E-05
> 0.7294E-05 0.7296E-05 0.7297E-05
> 0.0000E+00 0.0000E+00 0.0000E+00
> dectim = 1440.000000000000
> BCs are ready from 2000052400 to 2000052406
> dt, dtau = 150.0000000000000 37.50000000000000
> 75.00000000000000
> 0 linearization about standard atmosphere (lstand=.t.)
> 0sigmaf 0.000E+00 5.000E-02 0.100 0.160 0.230 0.310 0.390 0.470 0.550 0.630 0.710
> 0.780 0.840 0.890 0.930 0.960 0.980 0.990 1.00
> 0t mean 218. 218. 218. 218. 225. 236. 245. 253. 260. 267. 272.
> 276. 280. 282. 284. 286. 287. 287.
> 0ps mean 100.
> 0 vertical mode problem completed for kx= 18 0 errors detected (should be 0)
> m, fac = 4 33.33333333333334
> m, fac = 2 60.00000000000000
> Writing output files in direct access format
>
> ******* OPENING NEW OUTPUT FILES: 2000052400
> Note that rotated Mercartor (ROTMER) projections are not supported by GrADS.
> Although not exact, the eta.u projection in GrADS is somewhat similar.
> FERRET, however, does support this projection.
> OPENING NEW BAT FILE: output/SRF.2000052400
> Note that rotated Mercartor (ROTMER) projections are not supported by GrADS.
> Although not exact, the eta.u projection in GrADS is somewhat similar.
> FERRET, however, does support this projection.
> BATS variables written at 2000052400 2.500000000000000
> at day = 0.0851, ktau = 50 : 1st, 2nd time deriv of ps = 0.21323E-03 0.24245E-06, no. of points w/convection = 3816
> BATS variables written at 2000052403 180.0000000000000
> at day = 0.1719, ktau = 100 : 1st, 2nd time deriv of ps = 0.30562E-03 0.91678E-06, no. of points w/convection = 3388
> Floating point exception
> [root at cc1 exp_blhnew]#
>
> ______________________________________________________________________
> [root at cc4 test_final]# ./regcm <./regcm.in
> RESTARTPARAM READ IN
> TIMEPARAM READ IN
> OUTPARAM READ IN
> PHYSICSPARAM READ IN
> SUBEXPARAM READ IN
> GRELLPARAM READ IN
> EMANPARAM READ IN
> CHEMPARAM READ IN
> dtau = 37.50000000000000 75.00000000000000
> NREC = 289276
> IDATE1, IDATE2, dtmin, ktaur = 2000052400 2000060100
> 2.500000000000000 0
> READING HEADER FILE
> DIMS 111 120 18
> DOMAIN 50000.00 21.00000 80.00000 21.00000
> 80.00000 0.0000000E+00
> PROJROTMER
> SIGMA 0.0000000E+00 5.0000001E-02 0.1000000 0.1600000
> 0.2300000 0.3100000 0.3900000 0.4700000 0.5500000
> 0.6300000 0.7100000 0.7800000 0.8400000 0.8900000
> 0.9300000 0.9600000 0.9800000 0.9900000 1.000000
> PTOP 1.000000
> OUTPUT 1 1
> ***** mdate = 2000052400
> input/output parameters
> ifsave = T savfrq = 48.00000000000000 iftape = T tapfrq =
> 6.000000000000000 ifprt = F prtfrq = 12.00000000000000
> kxout = 6 jxsex = 40 radisp =
> 6.000000000000000 batfrq = 3.000000000000000 nslice =
> 120 ifchem = F chemfrq = 6.000000000000000
>
> physical parameterizations
> iboudy = 5 icup = 4 igcc = 2 ipptls =
> 1 iocnflx = 2 ipgf = 0 lakemod = 0
> ichem = 0
>
> model parameters
> radfrq = 30.00000000000000 abatm = 450.0000000000000 abemh =
> 18.00000000000000 dt = 150.0000000000000
>
> ncld = 1
>
> HT
> HTSD
> SATBRT
> XLAT
> XLONG
> MSFX
> MSFD
> F
> SNOWC
>
> ***************************************************
> ***************************************************
> **** RegCM IS BEING RUN ON THE FOLLOWING GRID: ****
> **** Map Projection: ROTMER ****
> **** IX= 111 JX= 120 KX= 18 ****
> **** PTOP= 1.000000 DX= 50000.00 ****
> **** CLAT= 21.00000 CLON= 80.00000 ****
> ***************************************************
>
> Index of highest allowed pbl: kt = 10
>
> AUTO-CONVERSION RATE: LAND= 5.0000000000000000E-004
> OCEAN= 5.0000000000000000E-004
> RELATIVE HUMIDITY THRESHOLDS: LAND= 0.8000000000000000
> OCEAN= 0.9000000000000000
> GULTEPE FACTORS: LAND= 0.4000000000000000 OCEAN=
> 0.4000000000000000
> MAXIMUM CLOUD COVER FOR RADIATION: 0.8000000000000000
> MAXIMUM RELATIVE HUMIDITY: 1.010000000000000
> rh0 temperature threshold: 238.0000000000000
> Raindrop Evaporation Rate 2.0000000000000000E-005
> Raindrop Accretion Rate 6.000000000000000
>
>
> EMANUEL (1991) CONVECTION V4.3C (20 May, 2002)
> MIN CONVECTION ORIGIN (minsig/orig): 0.9500000000000000 3
> AUTOCONVERSION THERSHOLD (elcrit): 1.0000000000000000E-003
> AUTOCONVERSION THRESHOLD TO ZERO (tlcrit): -55.00000000000000
> ENTRAINMENT COEFFICIENT (entp): 1.500000000000000
> FRACTIONAL AREA OF UNSATURATED DNDRAFT (sigd): 5.0000000000000000E-002
> PRECIP FRACTION OUTSIDE OF CLOUD (sigs): 0.1200000000000000
> FALL SPEED OF RAIN (omtrain): 50.00000000000000
> FALL SPEED OF SNOW (omtsnow): 5.500000000000000
> RAIN EVAPORATION COEFFICIENT (coeffr): 1.000000000000000
> SNOW EVAPORATION COEFFICIENT (coeffs): 0.8000000000000000
> CONVECTIVE MOMENTUM TRANSPORT COEFFICIENT (cu): 0.7000000000000000
> DOWNDRAFT VELOCITY SCALE (betae): 10.00000000000000
> MAX NEGATIVE PERTURBATION BELOW LFC (dtmax): 0.9000000000000000
> QUASI-EQUILIBRIUM APPROACH RATE (alphae): 1.0000000000000000E-002
> QUASI-EQUILIBRIUM APPROACH RATE (damp): 5.0000000000000000E-002
>
>
>
> CONVECTIVE CLOUD FRACTION/WATER
> Maximum Convective Cloud Cover
> before resolution scaling: 0.2500000000000000
> Maximum Convective Cloud Cover
> after resolution scaling: 0.2500000000000000
> Convective Cloud Water: 4.9999998736893758E-005
> the surface energy budget is used to calculate the ground temperature. julday = 145 gmt = 0.0
> the radiation is computed every 12 time steps.
>
> relaxation boudnary conditions (exponential method) are used. fnudge = 0.33333E-03 gnudge = 0.33333E+06
> 0 k sigma(k) a(k) dsigma(k) twt(k,1) twt(k,2) qcon(k)
>
> 1 0.0000 0.0250 0.0500 0.0000 0.0000 0.0000
> 2 0.0500 0.0750 0.0500 0.5000 0.5000 0.5000
> 3 0.1000 0.1300 0.0600 0.4545 0.5455 0.5455
> 4 0.1600 0.1950 0.0700 0.4615 0.5385 0.5385
> 5 0.2300 0.2700 0.0800 0.4667 0.5333 0.5333
> 6 0.3100 0.3500 0.0800 0.5000 0.5000 0.5000
> 7 0.3900 0.4300 0.0800 0.5000 0.5000 0.5000
> 8 0.4700 0.5100 0.0800 0.5000 0.5000 0.5000
> 9 0.5500 0.5900 0.0800 0.5000 0.5000 0.5000
> 10 0.6300 0.6700 0.0800 0.5000 0.5000 0.5000
> 11 0.7100 0.7450 0.0700 0.5333 0.4667 0.4667
> 12 0.7800 0.8100 0.0600 0.5385 0.4615 0.4615
> 13 0.8400 0.8650 0.0500 0.5455 0.4545 0.4545
> 14 0.8900 0.9100 0.0400 0.5556 0.4444 0.4444
> 15 0.9300 0.9450 0.0300 0.5714 0.4286 0.4286
> 16 0.9600 0.9700 0.0200 0.6000 0.4000 0.4000
> 17 0.9800 0.9850 0.0100 0.6667 0.3333 0.3333
> 18 0.9900 0.9950 0.0100 0.5000 0.5000 0.5000
> 19 1.0000
>
>
> maximum time = 8.000 days.
> time step = 150.00 seconds
> dx = 50000. meters
> grid points (x,y) = (120,111)
> number of levels = 18
> constant hor. diff. coef. = 0.25000E+05 m*m/s
> maximum hor. diff. coef. = 0.26042E+06 m*m/s
> READING INITAL CONDITIONS 2000052400
> U
> V
> TA
> QV
> PS
> TS
> *** initial total air = 0.30585E+18 kg, total water = 0.10863E+16 kg in large domain.
> *** solar declination angle = 20.99 degrees.
> ozone profiles
> 0.4963E-07 0.4963E-07 0.4963E-07
> 0.4963E-07 0.4963E-07 0.4963E-07
> 0.4963E-07 0.4963E-07 0.4963E-07
> 0.4963E-07 0.4963E-07 0.4963E-07
> 0.4981E-07 0.4985E-07 0.4985E-07
> 0.5267E-07 0.5267E-07 0.5267E-07
> 0.5367E-07 0.5370E-07 0.5371E-07
> 0.5775E-07 0.5775E-07 0.5775E-07
> 0.6413E-07 0.6417E-07 0.6417E-07
> 0.7402E-07 0.7406E-07 0.7407E-07
> 0.9089E-07 0.9101E-07 0.9104E-07
> 0.1196E-06 0.1197E-06 0.1197E-06
> 0.1658E-06 0.1661E-06 0.1661E-06
> 0.2862E-06 0.2866E-06 0.2867E-06
> 0.5890E-06 0.5896E-06 0.5897E-06
> 0.1156E-05 0.1157E-05 0.1157E-05
> 0.2574E-05 0.2576E-05 0.2576E-05
> 0.7294E-05 0.7296E-05 0.7297E-05
> 0.0000E+00 0.0000E+00 0.0000E+00
> dectim = 1440.000000000000
> BCs are ready from 2000052400 to 2000052406
> dt, dtau = 150.0000000000000 37.50000000000000
> 75.00000000000000
> 0 linearization about standard atmosphere (lstand=.t.)
> 0sigmaf 0.000E+00 5.000E-02 0.100 0.160 0.230 0.310 0.390 0.470 0.550 0.630 0.710
> 0.780 0.840 0.890 0.930 0.960 0.980 0.990 1.00
> 0t mean 218. 218. 218. 218. 225. 236. 245. 253. 260. 267. 272.
> 276. 280. 282. 284. 286. 287. 287.
> 0ps mean 100.
> 0 vertical mode problem completed for kx= 18 0 errors detected (should be 0)
> m, fac = 4 33.33333333333334
> m, fac = 2 60.00000000000000
> Writing output files in direct access format
>
> ******* OPENING NEW OUTPUT FILES: 2000052400
> Note that rotated Mercartor (ROTMER) projections are not supported by GrADS.
> Although not exact, the eta.u projection in GrADS is somewhat similar.
> FERRET, however, does support this projection.
> OPENING NEW OUT FILE: output/ATM.2000052400
> Note that rotated Mercartor (ROTMER) projections are not supported by GrADS.
> Although not exact, the eta.u projection in GrADS is somewhat similar.
> FERRET, however, does support this projection.
> OPENING NEW BAT FILE: output/SRF.2000052400
> Note that rotated Mercartor (ROTMER) projections are not supported by GrADS.
> Although not exact, the eta.u projection in GrADS is somewhat similar.
> FERRET, however, does support this projection.
> OPENING NEW RAD FILE: output/RAD.2000052400
> Note that rotated Mercartor (ROTMER) projections are not supported by GrADS.
> Although not exact, the eta.u projection in GrADS is somewhat similar.
> FERRET, however, does support this projection.
> OUT-history written date = 2000052400.000000
> BATS variables written at 2000052400 2.500000000000000
> Writing rad fields at ktau = 1 2000052400
> at day = 0.0851, ktau = 50 : 1st, 2nd time deriv of ps = 0.59293E-04 0.15078E-06, no. of points w/convection = 4470
> BATS variables written at 2000052403 180.0000000000000
> at day = 0.1719, ktau = 100 : 1st, 2nd time deriv of ps = 0.41453E-04 0.91110E-07, no. of points w/convection = 5652
> OUT-history written date = 2000052406.000000
> BATS variables written at 2000052406 0.0000000000000000E+000
> Writing rad fields at ktau = 144 2000052406
> BCs are ready from 2000052406 to 2000052412
> at day = 0.2587, ktau = 150 : 1st, 2nd time deriv of ps = 0.38419E-04 0.58503E-07, no. of points w/convection = 6507
> at day = 0.3455, ktau = 200 : 1st, 2nd time deriv of ps = 0.27137E-04 0.43160E-07, no. of points w/convection = 6954
> BATS variables written at 2000052409 180.0000000000000
> at day = 0.4323, ktau = 250 : 1st, 2nd time deriv of ps = 0.22644E-04 0.31854E-07, no. of points w/convection = 6982
> OUT-history written date = 2000052412.000000
> BATS variables written at 2000052412 0.0000000000000000E+000
> Writing rad fields at ktau = 288 2000052412
> BCs are ready from 2000052412 to 2000052418
> at day = 0.5191, ktau = 300 : 1st, 2nd time deriv of ps = 0.17762E-04 0.26758E-07, no. of points w/convection = 6238
> at day = 0.6059, ktau = 350 : 1st, 2nd time deriv of ps = 0.16011E-04 0.21594E-07, no. of points w/convection = 5079
> BATS variables written at 2000052415 180.0000000000000
> at day = 0.6927, ktau = 400 : 1st, 2nd time deriv of ps = 0.17237E-04 0.17195E-07, no. of points w/convection = 4570
> OUT-history written date = 2000052418.000000
> BATS variables written at 2000052418 0.0000000000000000E+000
> Writing rad fields at ktau = 432 2000052418
> BCs are ready from 2000052418 to 2000052500
> at day = 0.7795, ktau = 450 : 1st, 2nd time deriv of ps = 0.11075E-04 0.15838E-07, no. of points w/convection = 4624
> at day = 0.8663, ktau = 500 : 1st, 2nd time deriv of ps = 0.13804E-04 0.12813E-07, no. of points w/convection = 4897
> BATS variables written at 2000052421 180.0000000000000
> at day = 0.9531, ktau = 550 : 1st, 2nd time deriv of ps = 0.10921E-04 0.11440E-07, no. of points w/convection = 5027
> *** solar declination angle = 21.16 degrees.
> dectim = 2880.000000000000
> OUT-history written date = 2000052500.000000
> BATS variables written at 2000052500 0.0000000000000000E+000
> Writing rad fields at ktau = 576 2000052500
> BCs are ready from 2000052500 to 2000052506
> at day = 1.0399, ktau = 600 : 1st, 2nd time deriv of ps = 0.10631E-04 0.12301E-07, no. of points w/convection = 5154
> BATS variables written at 2000052503 180.0000000000000
> at day = 1.1267, ktau = 650 : 1st, 2nd time deriv of ps = 0.11342E-04 0.12299E-07, no. of points w/convection = 5889
> at day = 1.2135, ktau = 700 : 1st, 2nd time deriv of ps = 0.73936E-05 0.11815E-07, no. of points w/convection = 6707
> OUT-history written date = 2000052506.000000
> BATS variables written at 2000052506 0.0000000000000000E+000
> Writing rad fields at ktau = 720 2000052506
> BCs are ready from 2000052506 to 2000052512
> at day = 1.3003, ktau = 750 : 1st, 2nd time deriv of ps = 0.12899E-04 0.10888E-07, no. of points w/convection = 7218
> BATS variables written at 2000052509 180.0000000000000
> at day = 1.3872, ktau = 800 : 1st, 2nd time deriv of ps = 0.12873E-04 0.10021E-07, no. of points w/convection = 7451
> at day = 1.4740, ktau = 850 : 1st, 2nd time deriv of ps = 0.64977E-05 0.96710E-08, no. of points w/convection = 7176
> OUT-history written date = 2000052512.000000
> BATS variables written at 2000052512 0.0000000000000000E+000
> Writing rad fields at ktau = 864 2000052512
> BCs are ready from 2000052512 to 2000052518
> at day = 1.5608, ktau = 900 : 1st, 2nd time deriv of ps = 0.11829E-04 0.10404E-07, no. of points w/convection = 6055
> BATS variables written at 2000052515 180.0000000000000
> at day = 1.6476, ktau = 950 : 1st, 2nd time deriv of ps = 0.10732E-04 0.10311E-07, no. of points w/convection = 5195
> at day = 1.7344, ktau = 1000 : 1st, 2nd time deriv of ps = 0.69576E-05 0.94423E-08, no. of points w/convection = 4978
> OUT-history written date = 2000052518.000000
> BATS variables written at 2000052518 0.0000000000000000E+000
> Writing rad fields at ktau = 1008 2000052518
> BCs are ready from 2000052518 to 2000052600
> at day = 1.8212, ktau = 1050 : 1st, 2nd time deriv of ps = 0.81357E-05 0.10495E-07, no. of points w/convection = 4803
> BATS variables written at 2000052521 180.0000000000000
> at day = 1.9080, ktau = 1100 : 1st, 2nd time deriv of ps = 0.86746E-05 0.93758E-08, no. of points w/convection = 4885
> at day = 1.9948, ktau = 1150 : 1st, 2nd time deriv of ps = 0.48543E-05 0.93718E-08, no. of points w/convection = 4930
> *** solar declination angle = 21.33 degrees.
> dectim = 4320.000000000000
> OUT-history written date = 2000052600.000000
> BATS variables written at 2000052600 0.0000000000000000E+000
> Writing rad fields at ktau = 1152 2000052600
> SAVTMP RESTART WRITTEN: idatex= 2000052600ktau= 1152
> BCs are ready from 2000052600 to 2000052606
> at day = 2.0816, ktau = 1200 : 1st, 2nd time deriv of ps = 0.10597E-04 0.10068E-07, no. of points w/convection = 5089
> BATS variables written at 2000052603 180.0000000000000
> at day = 2.1684, ktau = 1250 : 1st, 2nd time deriv of ps = 0.70683E-05 0.10852E-07, no. of points w/convection = 5959
> OUT-history written date = 2000052606.000000
> BATS variables written at 2000052606 0.0000000000000000E+000
> Writing rad fields at ktau = 1296 2000052606
> BCs are ready from 2000052606 to 2000052612
> at day = 2.2552, ktau = 1300 : 1st, 2nd time deriv of ps = 0.10318E-04 0.12261E-07, no. of points w/convection = 6717
> at day = 2.3420, ktau = 1350 : 1st, 2nd time deriv of ps = 0.12736E-04 0.90837E-08, no. of points w/convection = 6981
> BATS variables written at 2000052609 180.0000000000000
> at day = 2.4288, ktau = 1400 : 1st, 2nd time deriv of ps = 0.76357E-05 0.92424E-08, no. of points w/convection = 7027
> OUT-history written date = 2000052612.000000
> BATS variables written at 2000052612 0.0000000000000000E+000
> Writing rad fields at ktau = 1440 2000052612
> BCs are ready from 2000052612 to 2000052618
> at day = 2.5156, ktau = 1450 : 1st, 2nd time deriv of ps = 0.66111E-05 0.11154E-07, no. of points w/convection = 6249
> at day = 2.6024, ktau = 1500 : 1st, 2nd time deriv of ps = 0.11121E-04 0.89689E-08, no. of points w/convection = 5414
> BATS variables written at 2000052615 180.0000000000000
> at day = 2.6892, ktau = 1550 : 1st, 2nd time deriv of ps = 0.80934E-05 0.90144E-08, no. of points w/convection = 4903
> OUT-history written date = 2000052618.000000
> BATS variables written at 2000052618 0.0000000000000000E+000
> Writing rad fields at ktau = 1584 2000052618
> BCs are ready from 2000052618 to 2000052700
> at day = 2.7760, ktau = 1600 : 1st, 2nd time deriv of ps = 0.75041E-05 0.10191E-07, no. of points w/convection = 4787
> at day = 2.8628, ktau = 1650 : 1st, 2nd time deriv of ps = 0.96583E-05 0.86815E-08, no. of points w/convection = 4957
> BATS variables written at 2000052621 180.0000000000000
> at day = 2.9497, ktau = 1700 : 1st, 2nd time deriv of ps = 0.62292E-05 0.87005E-08, no. of points w/convection = 5127
> *** solar declination angle = 21.50 degrees.
> dectim = 5760.000000000000
> OUT-history written date = 2000052700.000000
> BATS variables written at 2000052700 0.0000000000000000E+000
> Writing rad fields at ktau = 1728 2000052700
> BCs are ready from 2000052700 to 2000052706
> at day = 3.0365, ktau = 1750 : 1st, 2nd time deriv of ps = 0.90043E-05 0.96675E-08, no. of points w/convection = 5304
> at day = 3.1233, ktau = 1800 : 1st, 2nd time deriv of ps = 0.10972E-04 0.95897E-08, no. of points w/convection = 5860
> BATS variables written at 2000052703 180.0000000000000
> at day = 3.2101, ktau = 1850 : 1st, 2nd time deriv of ps = 0.54809E-05 0.99352E-08, no. of points w/convection = 6735
> OUT-history written date = 2000052706.000000
> BATS variables written at 2000052706 0.0000000000000000E+000
> Writing rad fields at ktau = 1872 2000052706
> BCs are ready from 2000052706 to 2000052712
> Floating point exception
> [root at cc4 test_final]#
>
>
>
> ______________________________________________________________________
> _______________________________________________
> RegCNET mailing list
> RegCNET at lists.ictp.it
> https://lists.ictp.it/mailman/listinfo/regcnet
More information about the RegCNET
mailing list