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
ASCII to cdump file
-
- Posts: 660
- Joined: November 7th, 2012, 3:14 pm
- Registered HYSPLIT User: Yes
Re: ASCII to cdump file
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
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
Hi,
Many thanks. It's very useful info.
Sree
Many thanks. It's very useful info.
Sree