Hi, Dear experts
I am expecting for your help.
The question is as follows.
I run the regcm-dust as follow steps,
1, datalinker.x
2, domain.param
I modified as follows,
      parameter(iproj='LAMCON') ! Map projection (LAMCON,NORMER,POLSTR,ROTMER)
      parameter(iy=96)          ! number of grid points in y direction (i)
      parameter(jx=96)          ! number of grid points in x direction (j)
      parameter(kz=18)          ! number of vertical levels (14, 18 or 23)
      parameter(nsg=1)          ! number of subgrids in 1 grid (one direction)
      parameter(ds=50.0)        ! grid point separation in km
      parameter(ptop= 5.0)      ! pressure of model top in cb.
      parameter(clat= 35.00)    ! Central latitude (North positive)
      parameter(clon= 105.00)    ! Central longitude (East positive)
      parameter(plat=clat)      ! Pole latitude  (used for ROTMER)
      parameter(plon=clon)      ! Pole longitude (used for ROTMER)
...
      parameter(AERTYP='AER11D1')    ! AEROSOL TYPE option,
      parameter(ntex=17)             ! Number of SOIL TEXTURE categories
      parameter(NPROC=16)             ! Number of CPU used for parallel job.
3, terrain.x
4, icbc.x
5, In the Main folder, I modified as follows.
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      integer ntr                      ! tracer number
      parameter( ntr = 10)
      integer nbin                     ! bins number ( dust)
      parameter(nbin = 4)
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
6, lrwxrwxrwx 1 yanlb yanlb 31 Feb 10 10:34 MAKECODE -> 0options/0_NODIAG_PARALLEL_CODE
Makefile
FC = mpif90
7, In my folder, regcm.in file
+++++++++
 $physicsparam
 iboudy  =          5,
 ibltyp  =          1,
 icup    =          2,
   igcc    =          2,
 ipptls  =          1,
 iocnflx =          2,
 ipgf    =          0,
 lakemod =          0,
 ichem   =          1,
 $
 $subexparam
 ncld    =  3,
 qck1land = 0.00025,
 qck1oce  = 0.00025,
 cevap    = 1.0e-3,
 caccr    = 3.0,
 $
 $grellparam
 $
 $emanparam
 $
 $chemparam
 idirect   =    2,
 chtrname  =    'SO2','SO4','BC_HB','BC_HL','OC_HB','OC_HL','DUST','DUST','DUST','DUST',
 chtrsol   =  0.20, 0.80, 0.05, 0.80, 0.05, 0.80, 0.50,  0.60,  0.70,  0.80,
 chtrdpv   = 0.00300,0.00200,0.00025,0.00025,0.00025,0.00025,0.00025,0.00050,0.00075,0.00100,0.00800,0.00200,0.00025,0.00250,0.00025,0.00250,0.00025,0.00050,0.00075,0.00100,
 dustbsiz  =  0.01, 1.00, 2.50, 5.00, 1.00, 2.50, 5.00,20.00,
 $
+++++++++
8, regcm.x
#!/bin/csh -f
set mydir=$PWD
cd ../Main
make clean
./MAKECODE
make
cd $mydir
mv ../Main/regcm .
/bin/ln -sf ../Input/DOMAIN.INFO fort.10
/bin/ln -sf ../Input/AERO.dat AERO.dat
/bin/ln -sf    ../Input/ICBC1994120100 fort.101
/bin/ln -sf    ../Input/ICBC1995010100 fort.102
...
9,mpirun -machinefile hostfile -np 16 -nolocal ./regcm<regcm.in >&log01
 
The follow is the output of my mpirun result in log01 file,
+++++++++++++++++++++++++
 process            0 of           16
 RESTARTPARAM READ IN
 TIMEPARAM READ IN
 OUTPARAM READ IN
 PHYSICSPARAM READ IN
 SUBEXPARAM READ IN
 GRELLPARAM READ IN
 CHEMPARAM READ IN
  dtau =     37.50000000000000         75.00000000000000
 NREC =        289276
 IDATE1, IDATE2, dtmin, ktaur =    1994120100   2000010100
    2.500000000000000                 0
 READING HEADER FILE
 DIMS           96           96           18
 DOMAIN    50000.00        35.00000        105.0000        35.00000
    105.0000       0.7155671
 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 =    1994120100
  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 =   T  chemfrq =    6.000000000000000
  physical parameterizations
  iboudy =             5  icup =             2  igcc =            2  ipptls =
            1  iocnflx =             2  ipgf =             0  lakemod =
            0  ichem =            1
  model parameters
  radfrq =     30.00000000000000       abatm =     450.0000000000000
  abemh =     18.00000000000000       dt =     150.0000000000000
  ncld =             3
 HT
 HTSD
 SATBRT
 XLAT
 XLONG
 MSFX
 MSFD
 F
 SNOWC
 ***************************************************
 ***************************************************
 **** RegCM IS BEING RUN ON THE FOLLOWING GRID: ****
 ****     Map Projection: LAMCON                ****
 ****     IX=           96  JX=           96  KX=           18
              ****
 ****     PTOP=    5.000000      DX=    50000.00            ****
 ****     CLAT=     35.00000      CLON=    105.0000         ****
 ***************************************************
  Index of highest allowed pbl:  kt =            10
 AUTO-CONVERSION RATE:  LAND=   2.5000000000000000E-004
                       OCEAN=   2.5000000000000000E-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.0000000000000000E-003
 Raindrop Accretion Rate    3.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 = 335   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 = 1857.000 days.
 time step =  150.00 seconds
 dx =  50000. meters
 grid points (x,y) = ( 96, 96)
 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   1994120100
 U
 V
 TA
 QV
 PS
 TS
           *** solar declination angle = -21.81 degrees.
 ozone profiles
   0.4963E-07
   0.4963E-07
   0.4963E-07
   0.4963E-07
   0.4963E-07
   0.5245E-07
   0.5282E-07
   0.5775E-07
   0.6289E-07
   0.7177E-07
   0.8547E-07
   0.1072E-06
   0.1498E-06
   0.2278E-06
   0.4169E-06
   0.8017E-06
   0.1419E-05
   0.3186E-05
   0.0000E+00
 dectim =     1440.000000000000
 BCs are ready from    1994120100   to    1994120106
 dt, dtau =     150.0000000000000         37.50000000000000
    75.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    33.33333333333334
 m, fac =             2    60.00000000000000
            1 SO2
            2 SO4
            3 BC_HB
            4 BC_HL
            5 OC_HB
            6 OC_HL
            7 DUST
            8 DUST
            9 DUST
           10 DUST
 Writing output files in direct access format
 ******* OPENING NEW OUTPUT FILES:   1994120100
 OPENING NEW OUT FILE: output/ATM.1994120100
 OPENING NEW BAT FILE: output/SRF.1994120100
 OPENING NEW RAD FILE: output/RAD.1994120100
 OPENING NEW CHEM FILE: output/CHE.1994120100
 OUT-history written date =     1994120100.000000
p0_26662:  p4_error: interrupt SIGFPE: 8
Killed by signal 2.
Killed by signal 2.
Killed by signal 2.
Killed by signal 2.
Killed by signal 2.
Killed by signal 2.
Killed by signal 2.
Killed by signal 2.
Killed by signal 2.
Killed by signal 2.
Killed by signal 2.
Killed by signal 2.
Killed by signal 2.
Killed by signal 2.
Killed by signal 2.
p0_26662: (18.421875) net_send: could not write to fd=4, errno = 32
+++++++++++++++++++++++++
 
How can i resolve it?
 
Yours,
Libin Yan
 
 
2009-02-18
Yanlb