Hello
I want to define a Vertical Emission Profile for a run in HYSPLIT. The dispersion run is a release from a nuclear power plant (Fukushima type accident). The release is spread in the vertical direction as follows: 0-50 m(15%); 50 - 300 m (15%) and 300 - 1000m (70%).
I'm looking for guidance on the best way to set up HYSPLIT for this type of scenario. I am currently thinking of using the same approach as for a volcanic (vertical line) release (https://ready.arl.noaa.gov/documents/Tu ... _emit.html). However, I just wondered if it is possible to do this with the 3 separate release ranges ( 0-50 m(15%); 50 - 300 m (15%) and 300 - 1000m (70%)). For example, I can easily amend the volcanic example above to cover the 0-50m but I just want advice on how I then input in the 50 - 300m and the 300-1000m line source into the EMITIMES file?
thanks
Kilian
EMITIMES file
Re: EMITIMES file
Here is an example of an EMITTIMES file.
This example emits 1 unit of mass over 5 minutes. 15% of the mass is emitted from 0-50m, 15% is emitted from 50-300m,
and 70% is emitted from 300-1000 m.
(DURATION is 0005= 5 minutes. To emit 1 unit mass total the rate is 12 unit mass/hour. 12*0.70 = 8.4 so rate for 300-1000m=8.5.
12*0.15=1.8 so rate for 0-50 and 50-300 is 1.8. The rate for the vertical line source is set by the bottom number.)
YYYY MM DD HH DURATION(hhhh) #RECORDS
YYYY MM DD HH MM DURATION(hhmm) LAT LON HGT(m) RATE(/h) AREA(m2) HEAT(w)
2015 01 09 03 0200 4
2015 01 09 03 00 0005 19.5 103 0 1.8 0 0
2015 01 09 03 00 0005 19.5 103 50 1.8 0 0
2015 01 09 03 00 0005 19.5 103 300 8.4 0 0
2015 01 09 03 00 0005 19.5 103 1000 0.0 0 0
This is the same as emitting 30% of mass from 0-300m and 70% mass from 300-1000m.
YYYY MM DD HH DURATION(hhhh) #RECORDS
YYYY MM DD HH MM DURATION(hhmm) LAT LON HGT(m) RATE(/h) AREA(m2) HEAT(w)
2015 01 09 03 0200 3
2015 01 09 03 00 0005 19.514 103.62 0 3.6 0 0
2015 01 09 03 00 0005 19.514 103.62 300 8.4 0 0
2015 01 09 03 00 0005 19.514 103.62 1000 0.0 0 0
These examples would be valid if one pollutant species is defined.
This example emits 1 unit of mass over 5 minutes. 15% of the mass is emitted from 0-50m, 15% is emitted from 50-300m,
and 70% is emitted from 300-1000 m.
(DURATION is 0005= 5 minutes. To emit 1 unit mass total the rate is 12 unit mass/hour. 12*0.70 = 8.4 so rate for 300-1000m=8.5.
12*0.15=1.8 so rate for 0-50 and 50-300 is 1.8. The rate for the vertical line source is set by the bottom number.)
YYYY MM DD HH DURATION(hhhh) #RECORDS
YYYY MM DD HH MM DURATION(hhmm) LAT LON HGT(m) RATE(/h) AREA(m2) HEAT(w)
2015 01 09 03 0200 4
2015 01 09 03 00 0005 19.5 103 0 1.8 0 0
2015 01 09 03 00 0005 19.5 103 50 1.8 0 0
2015 01 09 03 00 0005 19.5 103 300 8.4 0 0
2015 01 09 03 00 0005 19.5 103 1000 0.0 0 0
This is the same as emitting 30% of mass from 0-300m and 70% mass from 300-1000m.
YYYY MM DD HH DURATION(hhhh) #RECORDS
YYYY MM DD HH MM DURATION(hhmm) LAT LON HGT(m) RATE(/h) AREA(m2) HEAT(w)
2015 01 09 03 0200 3
2015 01 09 03 00 0005 19.514 103.62 0 3.6 0 0
2015 01 09 03 00 0005 19.514 103.62 300 8.4 0 0
2015 01 09 03 00 0005 19.514 103.62 1000 0.0 0 0
These examples would be valid if one pollutant species is defined.
Re: EMITIMES file
Thanks Alice for the response
That all makes perfect sense - I just have one (dumb!) question - can you explain where the 12 unit mass/hr comes from (To emit 1 unit mass total the rate is 12 unit mass/hour)?
thanks
That all makes perfect sense - I just have one (dumb!) question - can you explain where the 12 unit mass/hr comes from (To emit 1 unit mass total the rate is 12 unit mass/hour)?
thanks
Re: EMITIMES file
No problem - It's because I chose my emission duration as 5 minutes - 1/12 of an hour.
The duration of the emission is specified in the 6th column in format hhmm.
Then the rate of release per hour is specified in the 10th column.
This is a setup I often use for testing because it's easy to look in the MESSAGE file to verify that 1 unit of mass was released.
Cheers,
Alice
The duration of the emission is specified in the 6th column in format hhmm.
Then the rate of release per hour is specified in the 10th column.
This is a setup I often use for testing because it's easy to look in the MESSAGE file to verify that 1 unit of mass was released.
Cheers,
Alice
Re: EMITIMES file
thanks for the clarification - just one final question - is there anything I need to modify in the associated CONTROL file?
thanks
thanks
Re: EMITIMES file
To tell HYSPLIT to use the EMITIMES.txt file you need to put the following line in the SETUP.CFG file
efile='EMITIMES.txt'
So you can name your emitimes file anything you want; you just need to have this line in the SETUP.CFG file.
In the CONTROL file you will need to have as many sources lines as you have in the EMITIMES file.
So the first line should be the start time of the HYSPLIT run.
Then the next line should have the number of lines defining sources in the emitimes file.
Then there should be a dummy source line for each source in in the EMITIMES file - HYSPLIT will not use these so you can just put zeros if you like.
Below is the CONTROL file I used with the first EMITIMES example I posted earlier. It has 4 dummy sources in it.
If you don't have the correct number of dummy source lines in the CONTROL file, then HYSPLIT will stop running and you will get an error in the MESSAGE file which says
"ERROR emsmat: emission file mismatch with CONTROL"
Cheers,
Alice
15 01 09 03
4
0 0 0
0 0 0
0 0 0
0 0 0
3
0
32000.0
2
/pub/archives/gdas1/
gdas1.jan15.w1
/pub/archives/gdas1/
gdas1.jan15.w2
1
p006
1
01
00 00 00 00 00
1
55 , -145
0.5 0.50
90.0 360.0
./
out.bin
2
0 20000
00 00 00 00 00
00 00 00 00 00
00 01 00
1
0.6 2.5 1.0
0.0 0.0 0.0 0.0 0.0
0.0 4.0E+04 5.0E-06
0.0
0.0
~
efile='EMITIMES.txt'
So you can name your emitimes file anything you want; you just need to have this line in the SETUP.CFG file.
In the CONTROL file you will need to have as many sources lines as you have in the EMITIMES file.
So the first line should be the start time of the HYSPLIT run.
Then the next line should have the number of lines defining sources in the emitimes file.
Then there should be a dummy source line for each source in in the EMITIMES file - HYSPLIT will not use these so you can just put zeros if you like.
Below is the CONTROL file I used with the first EMITIMES example I posted earlier. It has 4 dummy sources in it.
If you don't have the correct number of dummy source lines in the CONTROL file, then HYSPLIT will stop running and you will get an error in the MESSAGE file which says
"ERROR emsmat: emission file mismatch with CONTROL"
Cheers,
Alice
15 01 09 03
4
0 0 0
0 0 0
0 0 0
0 0 0
3
0
32000.0
2
/pub/archives/gdas1/
gdas1.jan15.w1
/pub/archives/gdas1/
gdas1.jan15.w2
1
p006
1
01
00 00 00 00 00
1
55 , -145
0.5 0.50
90.0 360.0
./
out.bin
2
0 20000
00 00 00 00 00
00 00 00 00 00
00 01 00
1
0.6 2.5 1.0
0.0 0.0 0.0 0.0 0.0
0.0 4.0E+04 5.0E-06
0.0
0.0
~
Re: EMITIMES file
Hello again
would anyone have a copy of the various setup files (CONTROL, SETUP.CFG, CONC.CFG, EMITIMES) for a run using EMITIMES? I could then modify them for my own runs. I've just been trying to implement the changes suggested by alicec but still getting errors!!
I'm doing something silly but cant quite figure it out!!
thanks
Kilian
would anyone have a copy of the various setup files (CONTROL, SETUP.CFG, CONC.CFG, EMITIMES) for a run using EMITIMES? I could then modify them for my own runs. I've just been trying to implement the changes suggested by alicec but still getting errors!!
I'm doing something silly but cant quite figure it out!!
thanks
Kilian
-
- Posts: 660
- Joined: November 7th, 2012, 3:14 pm
- Registered HYSPLIT User: Yes
Re: EMITIMES file
Attached you will find the CONTROL, SETUP.CFG and EMITTIMES.txt files you will need to run 4 vertical line sources.
- Attachments
-
- CONTROL.txt
- (362 Bytes) Downloaded 290 times
-
- SETUP.txt
- (237 Bytes) Downloaded 259 times
-
- EMITTIMES.txt
- (571 Bytes) Downloaded 294 times