Time varying source

General questions and postings pertaining to the use of HYSPLIT regardless of the platform. For platform specific questions, use the HYSPLIT Platform forums.
Post Reply
Fabio.Dioguardi
Posts: 4
Joined: August 17th, 2018, 12:36 pm
Registered HYSPLIT User: Yes

Time varying source

Post by Fabio.Dioguardi » November 22nd, 2018, 7:30 am

Hi,

is it possible to set up a time-varying source with a time variation resolution < 1 hour, both in terms of emission rate and location (e.g. height)?
I am aware of the options explained in the manual at page 193 ("Time Variation of the Emission Rate") but:
- with the EFILE option, apparently I can act on the emission rate, not on the locations (hence height)
- with the restart option, I can again change the emission rate up to a time resolution of minutes, but this does not apply to the locations, since these are controlled from the Concentration Setup run (Total run time). I can run for 1 hour, then change the location heights and then restart, but apparently it does not work if I specify a fraction of hour in the field "Total run time". Unless I am completely wrong...

Thank you
Fabio

alicec
Posts: 101
Joined: February 8th, 2016, 12:56 pm
Registered HYSPLIT User: Yes

Re: Time varying source

Post by alicec » November 26th, 2018, 11:23 am

An EMITIMES file can handle almost any emission scenario.
https://www.ready.noaa.gov/hysplitusersguide/S417.htm

Constructing the more complicated ones can be confusing though.
Here is a short example of what you want to do (If I understand your
description correctly).


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 01 6
2015 01 09 03 00 0005 19.5 103.6 50 5 0 0
2015 01 09 03 35 0005 19.5 103.6 1000 1 0 0
2015 01 09 03 40 0010 19.5 103.6 300 8 0 0
2015 01 09 03 30 0005 19.0 103.0 50 5 0 0
2015 01 09 03 35 0005 19.0 103.0 1000 1 0 0
2015 01 09 03 40 0005 19.0 103.0 300 8 0 0

EMIT 5 units/hr from 19.5 103.6 at height of 50 m for 5 minutes starting at 09 03 00.
EMIT 1 units/hr from 19.5 103.6 at height of 1000 m for 5 minutes starting at 09 03 35.
EMIT 8 units/hr from 19.5 103.6 at height of 300 m for 10 minutes starting at 09 03 40.
Then three lines which are emit at a different location (19.0, 103.0)

A few words of caution

you may need to set the DELT in the SETUP.CFG file when using
time-varying emissions with resolution < 1 hour.
Otherwise when you allow HYSPLIT to set the time step, it may be set such that it misses one of your emissions.
For instance, for the above example I could set DELT=5 or DELT=1.

If are using multiple particle sizes, then you need to repeat lines for each particle size.
For instance, the first line would be emissions for your first particle size.
If you want same emission for second particle size, you would need to repeat the line.
HYSPLIT will expect that each particle size has it's own line to specify emissions for it.

If you have two lines that ONLY differ by the height (emission time and location are the same), then HYSPLIT will assume you want
a vertical line source there.

Fabio.Dioguardi
Posts: 4
Joined: August 17th, 2018, 12:36 pm
Registered HYSPLIT User: Yes

Re: Time varying source

Post by Fabio.Dioguardi » November 26th, 2018, 12:37 pm

Hi,

thank you very much.
In the meantime I found out, and I was able to run with the following EMITIMES file:

YYYY MM DD HH DURATION(hhhh) #RECORDS
YYYY MM DD HH MM DURATION(hhmm) LAT LON HGT(m) RATE(/h) AREA(m2) HEAT(w)
2018 11 22 14 0001 12
2018 11 22 14 00 0100 37.5 15.0 3300 74270634.42 0.0 0.0
2018 11 22 14 00 0100 37.5 15.0 3300 54803370.70 0.0 0.0
2018 11 22 14 00 0100 37.5 15.0 3300 24480582.40 0.0 0.0
2018 11 22 14 00 0100 37.5 15.0 3300 12274808.17 0.0 0.0
2018 11 22 14 00 0100 37.5 15.0 3300 3538068.61 0.0 0.0
2018 11 22 14 00 0100 37.5 15.0 3300 357057.07 0.0 0.0
2018 11 22 14 00 0100 37.5 15.0 7745 74270634.42 0.0 0.0
2018 11 22 14 00 0100 37.5 15.0 7745 54803370.70 0.0 0.0
2018 11 22 14 00 0100 37.5 15.0 7745 24480582.40 0.0 0.0
2018 11 22 14 00 0100 37.5 15.0 7745 12274808.17 0.0 0.0
2018 11 22 14 00 0100 37.5 15.0 7745 3538068.61 0.0 0.0
2018 11 22 14 00 0100 37.5 15.0 7745 357057.07 0.0 0.0
2018 11 22 15 0001 12
2018 11 22 15 00 0100 37.5 15.0 3300 36097370.60 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 3300 26635797.55 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 3300 11898170.28 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 3300 5965861.25 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 3300 1719589.10 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 3300 173538.59 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 6829 36097370.60 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 6829 26635797.55 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 6829 11898170.28 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 6829 5965861.25 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 6829 1719589.10 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 6829 173538.59 0.0 0.0

Basically I release from a line source 6 species with different emission rates. The height of release changes after 1 hour.
But when I try to increase the time resolution to 30 minutes, I have some problems. The calculation simply does not begin and no error messages are displayed. Here is the version of EMITIMES I am trying to use now:
YYYY MM DD HH DURATION(hhhh) #RECORDS
YYYY MM DD HH MM DURATION(hhmm) LAT LON HGT(m) RATE(/h) AREA(m2) HEAT(w)
2018 11 22 14 0001 24
2018 11 22 14 00 0030 37.5 15.0 3300 74270634.42 0.0 0.0
2018 11 22 14 00 0030 37.5 15.0 3300 54803370.70 0.0 0.0
2018 11 22 14 00 0030 37.5 15.0 3300 24480582.40 0.0 0.0
2018 11 22 14 00 0030 37.5 15.0 3300 12274808.17 0.0 0.0
2018 11 22 14 00 0030 37.5 15.0 3300 3538068.61 0.0 0.0
2018 11 22 14 00 0030 37.5 15.0 3300 357057.07 0.0 0.0
2018 11 22 14 00 0030 37.5 15.0 7745 74270634.42 0.0 0.0
2018 11 22 14 00 0030 37.5 15.0 7745 54803370.70 0.0 0.0
2018 11 22 14 00 0030 37.5 15.0 7745 24480582.40 0.0 0.0
2018 11 22 14 00 0030 37.5 15.0 7745 12274808.17 0.0 0.0
2018 11 22 14 00 0030 37.5 15.0 7745 3538068.61 0.0 0.0
2018 11 22 14 00 0030 37.5 15.0 7745 357057.07 0.0 0.0
2018 11 22 14 30 0030 37.5 15.0 3300 742706.42 0.0 0.0
2018 11 22 14 30 0030 37.5 15.0 3300 548070.70 0.0 0.0
2018 11 22 14 30 0030 37.5 15.0 3300 244805.40 0.0 0.0
2018 11 22 14 30 0030 37.5 15.0 3300 122748.17 0.0 0.0
2018 11 22 14 30 0030 37.5 15.0 3300 35380.61 0.0 0.0
2018 11 22 14 30 0030 37.5 15.0 3300 3570.07 0.0 0.0
2018 11 22 14 30 0030 37.5 15.0 4500 742706.42 0.0 0.0
2018 11 22 14 30 0030 37.5 15.0 4500 548033.70 0.0 0.0
2018 11 22 14 30 0030 37.5 15.0 4500 244805.40 0.0 0.0
2018 11 22 14 30 0030 37.5 15.0 4500 122748.17 0.0 0.0
2018 11 22 14 30 0030 37.5 15.0 4500 35380.61 0.0 0.0
2018 11 22 14 30 0030 37.5 15.0 4500 3570.07 0.0 0.0
2018 11 22 15 0001 12
2018 11 22 15 00 0100 37.5 15.0 3300 36097370.60 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 3300 26635797.55 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 3300 11898170.28 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 3300 5965861.25 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 3300 1719589.10 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 3300 173538.59 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 6829 36097370.60 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 6829 26635797.55 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 6829 11898170.28 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 6829 5965861.25 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 6829 1719589.10 0.0 0.0
2018 11 22 15 00 0100 37.5 15.0 6829 173538.59 0.0 0.0

Thank you
Fabio

alicec
Posts: 101
Joined: February 8th, 2016, 12:56 pm
Registered HYSPLIT User: Yes

Re: Time varying source

Post by alicec » November 26th, 2018, 12:57 pm

Each Emission Cycle must have the same number of lines. So if you have 24 lines in your
first emission cycle, you must also have 24 lines in your second emission cycle.
If you don't need 24 lines in your second emission cycle, you can add "dummy" sources
by setting the rate and duration to 0 for some lines.
This is necessary due to the way the arrays are allocated.

the other thing to check is that now you should have 4 source location lines in your CONTROL file instead
of 2.

Are you running from command line or GUI?
Check to see if a MESSAGE file is generated and see what that says.

Fabio.Dioguardi
Posts: 4
Joined: August 17th, 2018, 12:36 pm
Registered HYSPLIT User: Yes

Re: Time varying source

Post by Fabio.Dioguardi » November 27th, 2018, 7:03 am

OK now it works, thanks. I am running with the GUI.

However, when reverting back to the previous case (emissions changing with a 1 hr cycle), it does not work anymore. I am really getting mad this time, I don't see why. I am using exactly the same EMITIMES file that worked previously, and everything in CONTROL and SETUP.CFG seem ok. I attach them if you want to take a look...

Thanks
Fabio
Attachments
SETUP.CFG.txt
(568 Bytes) Downloaded 1 time
CONTROL.txt
(814 Bytes) Downloaded 3 times
EMITIMES.txt
(1.5 KiB) Downloaded 3 times

alicec
Posts: 101
Joined: February 8th, 2016, 12:56 pm
Registered HYSPLIT User: Yes

Re: Time varying source

Post by alicec » November 27th, 2018, 9:54 am

You should only have two source locations in your control file.
The number of source locations in the control file is important because array allocation
occurs before the EMITIMES file is read.
So the number of source locations in your control file multiplied by the number of
particle sizes that you have should be equal to the number of emission lines in each emission
cycle in your EMITIMES file.

Post Reply