p = WRFEnviron()
meteogram = 0
plot2D = 1
nens = 9
itime = (2015,2,7,0,0,0)
ftime = (2015,2,8,1,0,0)
hourly = 1
times = utils.generate_times(itime,ftime,hourly*60*60)
ncfiles = [os.path.join(ncroot,'SGp{0:02d}'.format(n)) for n in range(1,11)]
ncfiles.append(os.path.join(ncroot,'Gm4km'))
vrbls = ['T2','wind10']
loclist = [{'Stockton':(54.57,-1.32)},
{'Norwich':(52.63,1.30)},
{'FortWill':(56.82,-5.11)},
{'Wick':(58.45,-3.09)},
]
if meteogram:
for loc in loclist:
for vrbl in vrbls:
p.meteogram(vrbl,loc,ncfiles,outdir=outroot)
if plot2D:
for t in times:
p.plot2D('T2',t,ncdir=ncfiles[nens],outdir=outroot,cb=True,clvs=N.arange(263,284,1),f_suffix='{0}'.format(nens))
