MeteoInfoLab script - plot HYSPLIT concentration data

Post Reply
Posts: 32
Joined: October 5th, 2015, 10:20 pm
Registered HYSPLIT User: Yes

MeteoInfoLab script - plot HYSPLIT concentration data

Post by yaqiang » March 9th, 2018, 1:26 am

Get a HYSPLIT output concentration data file object using addfile_hyocnc function, and then read 2D array for filled contour plot. Following is an example of cesium-137 simulation during 2011 Fukushima nuclear disaster.

Code: Select all

fn = r'D:\Working\MIScript\JapanNuclear\cdump15'
f = addfile_hyconc(fn)
tidx = 5
data = f['C137'][tidx,0,:,:]

geoshow('cn_province', edgecolor='b')
geoshow('country', edgecolor='b')
levs = [1e-14,1e-13,1e-12,1e-11]
cols = ['w',(0,255,255),'g','y','r']
layer = contourfm(data, levs, colors=cols, smooth=False)
scatterm(141.0, 37.5, size=10, marker='S', color='r')
text(142.5, 39, 'Fukushima', bold=True)
xlim(110, 200)
ylim(15, 70)
yticks(arange(20, 61, 20))
colorbar(layer, label='C137', labelloc='top', shrink=0.8)
t = f.gettime(tidx)
title('Simulation of cesium-137 emitted from Fukushima (%s)' % t.strftime('%Y-%m-%d'))

Post Reply