Page 1 of 1
ASCII to cdump file
Posted: October 19th, 2018, 4:45 am
by sreeds82
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
Re: ASCII to cdump file
Posted: October 19th, 2018, 10:48 am
by ariel.stein
In the exec directory there is a program called conmaxv that will calculate the maximum concentration at each concentration grid.
Re: ASCII to cdump file
Posted: October 19th, 2018, 10:50 am
by yaqiang
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)
Re: ASCII to cdump file
Posted: October 19th, 2018, 11:49 am
by sreeds82
Hi,
Many thanks. It's very useful info.
Sree