Source Term Geometry

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
Smaug
Posts: 55
Joined: July 1st, 2017, 8:28 am
Registered HYSPLIT User: No

Source Term Geometry

Post by Smaug » July 24th, 2017, 5:14 pm

Hi,

I am aware that HYSPLIT can do a stack emission release by doing multiple releases at different heights with the same longitude and latitude. Are they any more capabilities where HYSPLIT can do different source term geometry releases such as cuboid, ellipsiod, spherical or circle releases?

Thanks

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

Re: Source Term Geometry

Post by alicec » July 25th, 2017, 4:59 pm

You can create complex emission terms by specifying multiple sources in the control file or by using an
EMITIMES file.

https://ready.arl.noaa.gov/documents/Tu ... _emit.html

For example please find below part of an EMITIMES file which creates a cylindrical source term for a run which
simulates four particle sizes. The cylindrical source is approximated by specifying many line sources (form 300 to 14000 meters).

YYYY MM DD HH DURATION(hhhh) #RECORDS
YYYY MM DD HH MM DURATION(hhmm) LAT LON HGT(m) RATE(/h) AREA(m2) HEAT(w)
2008 08 22 00 9999 1096
2008 08 07 22 00 0100 52.169 -175.511 300 5.839e-05 0 -1
2008 08 07 22 00 0100 52.169 -175.511 300 4.964e-04 0 -1
2008 08 07 22 00 0100 52.169 -175.511 300 1.854e-03 0 -1
2008 08 07 22 00 0100 52.169 -175.511 300 4.891e-03 0 -1
2008 08 07 22 00 0100 52.169 -175.511 14000 0 0 -1
2008 08 07 22 00 0100 52.169 -175.511 14000 0 0 -1
2008 08 07 22 00 0100 52.169 -175.511 14000 0 0 -1
2008 08 07 22 00 0100 52.169 -175.511 14000 0 0 -1
2008 08 07 22 00 0100 52.169 -175.365 300 5.839e-05 0 -1
2008 08 07 22 00 0100 52.169 -175.365 300 4.964e-04 0 -1
2008 08 07 22 00 0100 52.169 -175.365 300 1.854e-03 0 -1
2008 08 07 22 00 0100 52.169 -175.365 300 4.891e-03 0 -1
2008 08 07 22 00 0100 52.169 -175.365 14000 0 0 -1
2008 08 07 22 00 0100 52.169 -175.365 14000 0 0 -1
2008 08 07 22 00 0100 52.169 -175.365 14000 0 0 -1
2008 08 07 22 00 0100 52.169 -175.365 14000 0 0 -1
2008 08 07 22 00 0100 52.227 -175.399 300 5.839e-05 0 -1
2008 08 07 22 00 0100 52.227 -175.399 300 4.964e-04 0 -1
2008 08 07 22 00 0100 52.227 -175.399 300 1.854e-03 0 -1
2008 08 07 22 00 0100 52.227 -175.399 300 4.891e-03 0 -1
2008 08 07 22 00 0100 52.227 -175.399 14000 0 0 -1
2008 08 07 22 00 0100 52.227 -175.399 14000 0 0 -1
2008 08 07 22 00 0100 52.227 -175.399 14000 0 0 -1
2008 08 07 22 00 0100 52.227 -175.399 14000 0 0 -1
2008 08 07 22 00 0100 52.258 -175.486 300 5.839e-05 0 -1
2008 08 07 22 00 0100 52.258 -175.486 300 4.964e-04 0 -1
2008 08 07 22 00 0100 52.258 -175.486 300 1.854e-03 0 -1
2008 08 07 22 00 0100 52.258 -175.486 300 4.891e-03 0 -1
2008 08 07 22 00 0100 52.258 -175.486 14000 0 0 -1
2008 08 07 22 00 0100 52.258 -175.486 14000 0 0 -1
2008 08 07 22 00 0100 52.258 -175.486 14000 0 0 -1
2008 08 07 22 00 0100 52.258 -175.486 14000 0 0 -1
2008 08 07 22 00 0100 52.247 -175.585 300 5.839e-05 0 -1
2008 08 07 22 00 0100 52.247 -175.585 300 4.964e-04 0 -1
2008 08 07 22 00 0100 52.247 -175.585 300 1.854e-03 0 -1
2008 08 07 22 00 0100 52.247 -175.585 300 4.891e-03 0 -1
2008 08 07 22 00 0100 52.247 -175.585 14000 0 0 -1
2008 08 07 22 00 0100 52.247 -175.585 14000 0 0 -1
2008 08 07 22 00 0100 52.247 -175.585 14000 0 0 -1
2008 08 07 22 00 0100 52.247 -175.585 14000 0 0 -1
2008 08 07 22 00 0100 52.200 -175.649 300 5.839e-05 0 -1
2008 08 07 22 00 0100 52.200 -175.649 300 4.964e-04 0 -1
2008 08 07 22 00 0100 52.200 -175.649 300 1.854e-03 0 -1
2008 08 07 22 00 0100 52.200 -175.649 300 4.891e-03 0 -1
2008 08 07 22 00 0100 52.200 -175.649 14000 0 0 -1
2008 08 07 22 00 0100 52.200 -175.649 14000 0 0 -1
2008 08 07 22 00 0100 52.200 -175.649 14000 0 0 -1
2008 08 07 22 00 0100 52.200 -175.649 14000 0 0 -1
2008 08 07 22 00 0100 52.139 -175.649 300 5.839e-05 0 -1
2008 08 07 22 00 0100 52.139 -175.649 300 4.964e-04 0 -1
2008 08 07 22 00 0100 52.139 -175.649 300 1.854e-03 0 -1
2008 08 07 22 00 0100 52.139 -175.649 300 4.891e-03 0 -1
2008 08 07 22 00 0100 52.139 -175.649 14000 0 0 -1
2008 08 07 22 00 0100 52.139 -175.649 14000 0 0 -1
2008 08 07 22 00 0100 52.139 -175.649 14000 0 0 -1
2008 08 07 22 00 0100 52.139 -175.649 14000 0 0 -1



It would be considered an advanced topic, but it is also possible to write your own particle dump (PARDUMP) file and
initialized your HYSPLIT run from that. The PARDUMP file is binary file which specifies the exact position and status of computational particles.
HYSPLIT can output this type of file at the end of the run and then you can use that file to restart your HYSPLIT run from the place you left off.
This is the most common use of the PARDUMP file. Information about the format of the file can be found in the HYSPLIT help in the file S442.htm. We do not provide any support for creating your own pardump file. You would have to write your own program to create a binary file in the PARDUMP format which specifies the status of each computational particle.

For most people, the emitimes file provides more than enough flexibility to create any type of source they need.
Here is another thread about the emitimes file.
viewtopic.php?t=1062

Smaug
Posts: 55
Joined: July 1st, 2017, 8:28 am
Registered HYSPLIT User: No

Re: Source Term Geometry

Post by Smaug » August 1st, 2019, 6:24 am

When creating my own PARDUMP file what variables in the CONTROL and SETUP.CFG file have to match the PARDUMP files? Is it the case that the total mass of the particles has to match the total mass calculated from the CONTROL file?

barbara.stunder
Posts: 330
Joined: November 9th, 2012, 4:23 pm
Registered HYSPLIT User: Yes

Re: Source Term Geometry

Post by barbara.stunder » August 5th, 2019, 10:32 am

No to your second question. See viewtopic.php?f=28&t=1496#p4914

Post Reply