ASCII to cdump file

[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

[ariel.stein]

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

[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)

[sreeds82]

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