Dear colleagues:
We are trying to run RegCM with the following input namelist (below) and are getting the following
SIGFPE: Floating-point exception – erroneous arithmetic operation. Screen shot is attached. Any idea what causes (and where) the SIGFPE? Any idea where the problem may be?
Thanks in advance,
Niko Grisel Todorov
Chapman University
(424) 262-5354
&dimparam
iy = 220, ! 250
jx = 200, ! 240
kz = 23,
/
&coreparam
idynamic = 1,
/
&geoparam
iproj = 'NORMER',
!iproj = 'ROTMER',
ds = 50,
ptop = 5.0,
clat = 00.00,
clon = 20.00, ! 15.00
plat = 00.00,
plon = 15.00,
/
&terrainparam
domname = 'had-rf', ! 'mp-rf'
smthbdy = .false.,
h2opct = 50.,
h2ohgt = .true.,
ismthlev = 1,
dirter = '/data/sharefolder/menaclimate/HadGEM2', ! '/sharefolder/menaclimate/shreef/input/'
inpter = '/data/sharefolder/menaclimate/RCM', ! '/sharefolder/menaclimate/RCM/'
/
&boundaryparam
nspgx = 20,
nspgd = 20,
high_nudge = 8.0D0,
medium_nudge = 6.0D0,
low_nudge = 4.0D0,
/
&globdatparam
ibdyfrq = 6,
ssttyp = 'HA_RF', ! 'MP_RF'
dattyp = 'HA_26', ! 'MP_RF'
gdate1 = 1964060100, ! 1949010100
gdate2 = 2005120100, ! 2005010100
dirglob = '/data/sharefolder/menaclimate/HadGEM2', ! '/sharefolder/menaclimate/shreef/input'
inpglob = '/data/sharefolder/menaclimate/RCM/AERGLOB', ! '/sharefolder/menaclimate/data'
calendar = '360_day', ! 'gregorian'
/
&restartparam
ifrest = .true., ! .false.
mdate0 = 1959120100, ! 1949010100
mdate1 = 1984060100, ! 1949010100
mdate2 = 2004121800, ! 2005010100
/
&timeparam
dt = 100.,
/
&outparam
ifsave = .true. ,
savfrq = 0.,
ifatm = .true. ,
atmfrq = 6.,
ifrad = .true. ,
radfrq = 6.,
ifsrf = .true. ,
ifsub = .true. ,
srffrq = 6.,
ifchem = .false.,
chemfrq = 6.,
dirout = '/data/sharefolder/menaclimate/HadGEM2', ! '/sharefolder/menaclimate/shreef/output'
/
&physicsparam
iboudy = 5,
isladvec = 0,
iqmsl = 1,
ibltyp = 1,
icup_lnd = 2,
icup_ocn = 6,
ipptls = 1,
iocnflx = 2,
iocnrough = 1,
iocnzoq = 1,
ipgf = 0,
iemiss = 0,
lakemod = 0,
ichem = 0,
!scenario = 'RPC4.5',
idcsst = 0,
iseaice = 0,
idesseas = 1,
iconvlwp = 1,
icldfrac = 0,
icldmstrat = 1,
icumcloud = 1,
irrtm = 0,
iclimao3 = 0,
isolconst = 1,
islab_ocean = 0,
itweak = 0,
/
&subexparam
qck1land = 0.1e-4, ! Autoconversion Rate for Land
qck1oce = 0.0005, ! Autoconversion Rate for Ocean
gulland = 0.80, ! Fract of Gultepe eqn (qcth) when prcp occurs (land)
guloce = 0.30, ! Fract of Gultepe eqn (qcth) for ocean
cevaplnd = 1.0e-5, ! Raindrop evap rate coef land [[(kg m-2 s-1)-1/2]/s]
cevapoce = 1.0e-5, ! Raindrop evap rate coef ocean [[(kg m-2 s-1)-1/2]/s]
caccrlnd = 6.0, ! Raindrop accretion rate land [m3/kg/s]
caccroce = 4.0, ! Raindrop accretion rate ocean [m3/kg/s]
conf = 1.00, ! Condensation efficiency
/
&grellparam
igcc = 2, ! Cumulus closure scheme
! 1 => Arakawa & Schubert (1974)
! 2 => Fritsch & Chappell (1980)
/
&tiedtkeparam
iconv = 4, ! Actual used scheme.
entrmax = 1.75e-3, ! Max entrainment iconv=[1,2,3]
entrdd = 3.0e-4, ! Entrainment rate for cumulus downdrafts
entrpen_lnd = 0.3e-2, ! Entrainment rate for penetrative convection
entrpen_ocn = 0.3e-2, ! Entrainment rate for penetrative convection
entrscv = 3.0e-4, ! Entrainment rate for shallow convection iconv=[1,2,3]
entrmid = 1.0e-4, ! Entrainment rate for midlevel convection iconv=[1,2,3]
cprcon = 1.0e-4, ! Coefficient for determining conversion iconv=[1,2,3]
detrpen_lnd = 0.75e-4, ! Detrainment rate for penetrative convection iconv=4
detrpen_ocn = 0.75e-4, ! Detrainment rate for penetrative convection iconv=
rcuc_lnd = 0.05, ! Convective cloud cover for rain evporation iconv=4
rcuc_ocn = 0.05, ! Convective cloud cover for rain evporation iconv=4
rcpec_lnd = 5.55e-5, ! Coefficient for rain evaporation below cloud iconv=4
rcpec_ocn = 5.55e-5, ! Coefficient for rain evaporation below cloud iconv=4
rhebc_lnd = 0.5, ! Critical rh below cloud for evaporation iconv=4
rhebc_ocn = 0.6, ! Critical rh below cloud for evaporation iconv=4
rprc_lnd = 6.0e-2, ! conversion coefficient from cloud water iconv=4
rprc_ocn = 3.0e-2, ! conversion coefficient from cloud water iconv=4
entshalp = 2.0, ! shallow entrainment factor for entrorg iconv=4
cmtcape = 3600.0, ! CAPE adjustment timescale iconv=4
lmfpen = .true., ! penetrative conv is switched on
lmfmid = .true., ! midlevel conv is switched on
lmfdd = .true., ! cumulus downdraft is switched on
lepcld = .true., ! prognostic cloud scheme is on
lmfdudv = .true., ! cumulus friction is switched on
lmfscv = .true., ! shallow convection is switched on
lmfuvdis = .true., ! use kinetic energy dissipation
lmftrac = .true., ! chemical tracer transport is on
lmfsmooth = .false., ! smoot of mass fluxes for tracers
lmfwstar = .false., ! Grant w* closure for shallow conv
/
&kfparam
kf_convrate = 0.03,
kf_entrate = 0.03, ! Entrainment rate
kf_min_pef = 0.2, ! Minimum precipitation efficiency
kf_max_pef = 0.9, ! Maximum precipitation efficiency
kf_dpp = 150.0, ! Start elevation for downdraft above cloud base (mb)
kf_tkemax = 5.0, ! Maximum turbolent kinetic energy in sub cloud layer
kf_min_dtcape = 1800.0, ! Consumption time of CAPE low limit
kf_max_dtcape = 3600.0, ! Consumption time of CAPE high limit
/
&holtslagparam
ricr_ocn = 0.25, ! Critical Richardson Number over Ocean
ricr_lnd = 0.25, ! Critical Richardson Number over Land
zhnew_fac = 0.25, ! Multiplicative factor for zzhnew in holtpbl
ifaholtth10 = 1, ! First approximation for obhukov length, th10 formula:
ifaholt = 1, ! th10 final adjustment:
/
&uwparam
iuwvadv = 0, ! 0=standard T/QV/QC advection, 1=GB01-style advection
atwo = 15.0, ! Efficiency of enhancement of entrainment by cloud evap.
rstbl = 1.5, ! Scaling parameter for stable boundary layer eddy length
czero = 5.869, ! Czero constant in UW PBL (eqn 44a and pgs 856-857)
nuk = 5.0D0, ! Multiplication factor for diffusion coefficients
/
&clm_inparm
fpftcon = 'pft-physiology.c130503.nc',
fsnowoptics = 'snicar_optics_5bnd_c090915.nc',
fsnowaging = 'snicar_drdt_bst_fit_60_c070416.nc',
/
&clm_soilhydrology_inparm
h2osfcflag = 1,
origflag = 0,
/
&clm_hydrology1_inparm
oldfflag = 0,
/