[RegCNET] How to resolve it?
Yanlb
yanlibin at ieecas.cn
Wed Feb 18 06:07:47 CET 2009
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20090218/55527c35/attachment.html>
More information about the RegCNET
mailing list