I used to running HYSPLIT Model to simulate Nuclear Power Plant accident.
I used kg/hr as a unit of emission rate.
I attached my default_conc below and i emphasize emission rate with Bold font and under line.
Code: Select all
19 07 23 01
1
35.41 126.42 50.0
96
0
10000.0
1
D:/G-CIEM_test/
0723_5DAY
1
C137
[b][u]4.65E+15[/u][/b]
96
19 07 23 01
1
35.88 127.5
0.03 0.05
5.6 5.0
D:/G-CIEM_test/
0723_5DAY_
1
0
19 07 23 01 00
19 07 27 01 00
00 12 00
1
0.4 1.9 1.0
0.002 136.96 0.0 0.0 0.0
0.0 3.2E+05 5.0E-05
10960.0
1.0E-06
But I found that total deposited mass is larger than total amount of emission mass.
I attached calculation history below.
Code: Select all
Total deposited mass = Σ(Concentration of grid[kg/m2])*(grid area[m2])
Total emission mass = (Emission rate[kg/hr])*(emission duration[hr])
1. Are there any error or mistake in my mass calculation?
2. Why the situation that Total deposited mass is larger than Total emission mass, occur?
3. To calculate MASS BALANCE with HYSPLIT Result, what equation should i use? or Which section in the manual should i read?
And i also attached CONC.CFG and SETUP.CFG below.
CONC.CFG
Code: Select all
&SETUP
tratio = 0.75,
initd = 0,
kpuff = 0,
khmax = 9999,
kmixd = 0,
kmix0 = 250,
kzmix = 0,
kdef = 0,
kbls = 1,
kblt = 0,
conage = 48,
numpar = -2500,
qcycle = 0.0,
efile = '',
tkerd = 0.18,
tkern = 0.18,
ninit = 1,
ndump = 0,
ncycl = 0,
pinpf = 'PARINIT',
poutf = 'PARDUMP',
mgmin = 10,
kmsl = 0,
maxpar = 10000,
cpack = 1,
cmass = 0,
dxf = 1.0,
dyf = 1.0,
dzf = 0.01,
ichem = 0,
maxdim = 1,
kspl = 1,
krnd = 6,
frhs = 1.0,
frvs = 0.01,
frts = 0.10,
frhmax = 3.0,
splitf = 1.0,
/
Code: Select all
&SETUP
tratio = 0.75,
initd = 0,
kpuff = 0,
khmax = 9999,
kmixd = 0,
kmix0 = 250,
kzmix = 0,
kdef = 0,
kbls = 1,
kblt = 0,
conage = 48,
numpar = -2500,
qcycle = 0.0,
efile = '',
tkerd = 0.18,
tkern = 0.18,
ninit = 1,
ndump = 0,
ncycl = 0,
pinpf = 'PARINIT',
poutf = 'PARDUMP',
mgmin = 10,
kmsl = 0,
maxpar = 10000,
cpack = 1,
cmass = 0,
dxf = 1.0,
dyf = 1.0,
dzf = 0.01,
ichem = 0,
maxdim = 1,
kspl = 1,
krnd = 6,
frhs = 1.0,
frvs = 0.01,
frts = 0.10,
frhmax = 3.0,
splitf = 1.0,
/
Thank you in advanced.