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

todorov@chapman.edu

 

 

&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,

/