[RegCNET] RegCM4 with dust

[Ismaila Diallo]

Dear Regcnet,

i plan to do a simulation of regcm4 with dust over west Africa. But i haven't 
idea about how kind of change must i do i my regcm.in. Hope in your help.
here is my regcm.in 

! This is a simple demo input file for regcm4.0 
! used in the RegCM Workshop in june 2010 
 &dimparam
 iy     = 120,
 jx     = 150,
 kz     = 18,
 nsg    = 1,
 /
 &geoparam
 iproj = 'NORMER',
 ds = 60.0,
 ptop = 5.0,
 clat = 10.39,
 clon = 9.48,
 plat = 10.39,
 plon = 9.48,
 truelatl = 30.0,
 truelath = 60,
 /
 &terrainparam
 domname = 'West_Africa',
 itype_in = 1,
 ntypec = 10,
 ntypec_s = 10,
 ifanal = .true.,
 smthbdy = .false.,
 lakadj = .false.,
 fudge_lnd   = .false.,
 fudge_lnd_s = .false.,
 fudge_tex   = .false.,
 fudge_tex_s = .false.,
 ntex = 17,
 h2opct = 75.,
 dirter = '/home/lpaosf/regcm/second_run/input',
 inpter = '/home/lpaosf/Data_Boundary',
 /
 &ioparam
 igrads = 1,
 ibigend = 1,
 ibyte = 4,
 iotyp = 1,
 /
 &debugparam
 debug_level = 1
 /
 &boundaryparam
 nspgx  = 12,
 nspgd  = 12,
 /
 &modesparam
 nsplit = 2,
 /
 &lakemodparam
 lkpts  = 10,
 /
 &globdatparam
 ssttyp = 'OISST',
 dattyp = 'EIN15',
 ehso4 = .false.,
 globidate1 = 1991010100,
 globidate2 = 1991123118,
 dirglob = '/home/lpaosf/regcm/second_run/input_icbc',
 inpglob = '/home/lpaosf/Data_Boundary',
 /
 &lsmparam
 lsmtyp = 'BATS',
 /
 &aerosolparam
 aertyp = 'AER00D0'
 ntr = 10,
 nbin = 4,
 /
 &restartparam
 ifrest  = .false. ,
 idate0  = 1991010100,
 idate1  = 1991010100,
 idate2  = 1991123118,
 nslice  =        120,
 /
 &timeparam
 radfrq  =    30.,
 abemh   =    18.,
 abatm   =   600.,
 dt      =   150.,
 ibdyfrq =          6,
 /
 &outparam
 ifsave  = .true. ,
   savfrq  =    48.,
 iftape  = .true. ,
   tapfrq  =     6.,
 ifrad   = .true. ,
   radisp  =     6.,
 ifbat   = .true. ,
 ifsub   = .true. ,
   batfrq  =     3.,
 ifprt   = .false.,
   prtfrq  =    12.,
   kxout   =          6,
   jxsex   =         40,
 ifchem  = .false.,
   chemfrq =     6.,
 dirout='/home/lpaosf/regcm/second_run/output',
 /
 &physicsparam
 iboudy  =          5,
 ibltyp  =          1,
 icup    =          4,
 igcc    =          1,
 ipptls  =          1,
 iocnflx =          2,
 ipgf    =          0,
 iemiss  =          0,
 lakemod =          0,
 ichem   =          0,
 /
 &subexparam
 qck1land=   .250E-03,
 qck1oce =   .250E-03,
 cevap   =   .100E-02,
 caccr   =      3.000,
 /
 &grellparam
 /
 &emanparam
 /
  &chemparam
 idirect   =    1,
 inpchtrname  = 'DUST','DUST','DUST','DUST', 
'SO2','SO4','BC_HB','BC_HL','OC_HB','OC_HL',
 inpchtrsol   =  0.10, 0.10, 0.10, 0.10, 0.20, 0.80, 0.05, 0.80, 0.05, 0.80
 inpchtrdpv(:,1) =0.,0.,0.,0., 0.00300,0.00200,0.00025,0.00025,0.00025,0.00025
 inpchtrdpv(:,2) = 0.,0.,0.,0., 0.00800,0.0020 
,0,0.00025,0.00250,0.00025,0.00250
 inpdustbsiz(:,1)  = 0.01, 1.00, 2.50, 5.00
 inpdustbsiz(:,2)= 1.00, 2.50, 5.00,20.00
 /
 &clmparam
 dirclm = '/set/this/to/where/your/input/clm/data/are',
 imask  = 1,
 clmfrq = 24.,
 /


thanks in advance

Cheers!!!!!
ismaila