ASCII to cdump file

Topics about the HYSPLIT dispersion model.
Post Reply
sreeds82
Posts: 13
Joined: December 2nd, 2015, 5:57 am
Registered HYSPLIT User: Yes

ASCII to cdump file

Post by sreeds82 » October 19th, 2018, 4:45 am

Dear HYSPLIT community,

I have a query, is it possible to display the maximum concentration (single display) value of whole domain considering all time steps from cdump file?

Is it possible to prepare an ascii text to cdump file, in order to display using hysplit display utility? Any help in this regard is appreciated.

Thanks in advacnce
Sree

ariel.stein
Posts: 635
Joined: November 7th, 2012, 3:14 pm
Registered HYSPLIT User: Yes
Contact:

Re: ASCII to cdump file

Post by ariel.stein » October 19th, 2018, 10:48 am

In the exec directory there is a program called conmaxv that will calculate the maximum concentration at each concentration grid.

yaqiang
Posts: 44
Joined: October 5th, 2015, 10:20 pm
Registered HYSPLIT User: Yes

Re: ASCII to cdump file

Post by yaqiang » October 19th, 2018, 10:50 am

You may try MeteoInfoLab script to get the max value and location from cdump file, and plot it on map.

Code: Select all

#Read data array from HYSPLIT concentration output data
fn = r'D:\Working\MIScript\JapanNuclear\cdump15'
f = addfile_hyconc(fn)
data = f['C137'][:,0,:,:]  # 3 dimension array - time, lat, lon

#Get the maximum value and location
v_max = data.max()
idx = data.argmax()
idx = unravel_index(idx, data.shape)
t = f.gettime(idx[0])
lat = data.dimvalue(1)[idx[1]]
lon = data.dimvalue(2)[idx[2]]

#Plot the maximum value
axesm()
geoshow('country', edgecolor='k')
scatterm(lon, lat, color='r', size=10)
text(lon + 2, lat - 2, '%s\n%s' % (str(v_max), t.strftime('%Y-%m-%d %H')))
xlim(110, 200)
ylim(15, 70)
hyconc_max.png

sreeds82
Posts: 13
Joined: December 2nd, 2015, 5:57 am
Registered HYSPLIT User: Yes

Re: ASCII to cdump file

Post by sreeds82 » October 19th, 2018, 11:49 am

Hi,
Many thanks. It's very useful info.
Sree

Post Reply