Hello every RegCM3 users,

Now I am using RegCM3 to simulate regional climate. I chosen the 1990 year, in the simulation, ds = 10, clat = 48.5, clon = 117, iy = 90, jx = 100, DATTYP = EH5OM, SSTTYP = EH5RF. In this situations, I could run the RegCM3 successfully.

When I chosen the 2020 year, the simulation parameters as follow: ds = 10, clat = 48.5, clon = 117, iy = 90, jx = 100, DATTYP = EH5OM, SSTTYP = EHA1B. and other parameters were the same with the simulation in 1990. terrain.x and icbc.x could run successfully, but when the regcm.x run, it was found that “floating point exception”. Who can tell me the reasons? Best regards.

Part of information when the regcm.x run.

dtau = 7.500000000000000 15.00000000000000

NREC = 289276

IDATE1, IDATE2, dtmin, ktaur = 2080020100 2080021000

0.5000000000000000 0

READING HEADER FILE

DIMS 90 100 18

DOMAIN 10000.00 48.50000 117.0000 48.50000

117.0000 0.7155668

PROJLAMCON

SIGMA 0.000000 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 5.000000

OUTPUT 1 1

***** mdate = 2080020100

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 = 2 igcc = 2 ipptls =

1 iocnflx = 2 ipgf = 0 lakemod =

0 ichem = 0

model parameters

radfrq = 30.00000000000000 abatm = 90.00000000000000

abemh = 18.00000000000000 dt = 30.00000000000000

ncld = 1

HTSD

SATBRT

XLAT

XLONG

MSFX

MSFD

SNOWC

***************************************************

**** RegCM IS BEING RUN ON THE FOLLOWING GRID: ****

**** Map Projection: LAMCON ****

**** IX= 90 JX= 100 KX= 18

****

**** PTOP= 5.000000 DX= 10000.00 ****

**** CLAT= 48.50000 CLON= 117.0000 ****

***************************************************

Index of highest allowed pbl: kt = 10

AUTO-CONVERSION RATE: LAND= 5.0000002374872565E-004

OCEAN= 5.0000002374872565E-004

RELATIVE HUMIDITY THRESHOLDS: LAND= 0.8000000119209290

OCEAN= 0.8999999761581421

GULTEPE FACTORS: LAND= 0.4000000059604645 OCEAN=

0.4000000059604645

MAXIMUM CLOUD COVER FOR RADIATION: 0.8000000119209290

MAXIMUM RELATIVE HUMIDITY: 1.009999990463257

rh0 temperature threshold: 238.0000000000000

Raindrop Evaporation Rate 1.9999999494757503E-005

Raindrop Accretion Rate 6.000000000000000

*******************************************************

Grell Convection:

Max Shear: 0.5000000000000000

Min Shear: 0.2500000000000000

Max PPT eff: 1.000000000000000

Min PPT eff: 0.2500000000000000

Max PPT eff(o): 1.000000000000000

Min PPT eff(o): 0.000000000000000

Max PPT eff(x): 1.000000000000000

Min PPT eff(x): 0.2500000000000000

Max PBC (pbcmax): 150.0000000000000

Min Cloud Depth: 150.0000000000000

Max Cloud Base: 12 0.4000000059604645

Max Heating: 500.0000000000000

Min Heating: -250.0000000000000

Fritsch-Chappell (1980) Closure Assumption *****

ABE removal timescale: dtauc= 30.00000000000000

*******************************************************

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: 3.0000001424923539E-004

the surface energy budget is used to calculate the ground temperature. julday = 32 gmt = 0.0

the radiation is computed every 60 time steps.

relaxation boudnary conditions (exponential method) are used. fnudge = 0.16667E-02 gnudge = 0.66667E+05

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 = 9.000 days.

time step = 30.00 seconds

dx = 10000. meters

grid points (x,y) = (100, 90)

number of levels = 18

constant hor. diff. coef. = 0.50000E+04 m*m/s

maximum hor. diff. coef. = 0.52083E+05 m*m/s

READING INITAL CONDITIONS 2080020100

SO0

*** solar declination angle = -17.05 degrees.

ozone profiles

0.4963E-07

0.5239E-07

0.5267E-07

0.5596E-07

0.5775E-07

0.6359E-07

0.7041E-07

0.8196E-07

0.9969E-07

0.1265E-06

0.1744E-06

0.2778E-06

0.5098E-06

0.9142E-06

0.1556E-05

0.3323E-05

0.0000E+00

dectim = 1440.000000000000

BCs are ready from 2080020100 to 2080020106

dt, dtau = 30.00000000000000 7.500000000000000

15.00000000000000

0 linearization about standard atmosphere (lstand=.t.)

0sigmaf 0.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. 229. 239. 248. 255. 262. 268. 273.

277. 280. 283. 285. 286. 287. 287.

0ps mean 100.

0 vertical mode problem completed for kx= 18 0 errors detected (should be 0)

m, fac = 4 6.666666666666667

m, fac = 2 12.00000000000000

Writing output files in direct access format

******* OPENING NEW OUTPUT FILES: 2080020100

OPENING NEW OUT FILE: output/ATM.2080020100

OPENING NEW BAT FILE: output/SRF.2080020100

OPENING NEW RAD FILE: output/RAD.2080020100

OUT-history written date = 2080020100.000000

Floating point exception