indexes,weights = TL.select(req)
VARCONV=CONVERSION_DICT[var]
# setting up filelist for requested period -----------------
filelist=[]
for k in indexes:
t = TL.Timelist[k]
filename = INPUTDIR + "ave." + t.strftime("%Y%m%d-%H:%M:%S") + ".nc"
print filename
filelist.append(filename)
# ----------------------------------------------------------
M3d = TimeAverager3D(filelist, weights, var, TheMask)
for il,layer in enumerate(LAYERLIST):
De = DataExtractor(TheMask,rawdata=M3d)
if args.optype=='integral':
integrated = MapBuilder.get_layer_integral(De, layer)
else:
integrated = MapBuilder.get_layer_average(De, layer)
integrated=integrated * VARCONV
# mask200=TheMask.mask_at_level(200)
mask200=TheMask.mask_at_level(LIMIT_PER_MASK[il])
# clim = [M3d[TheMask.mask].min(), M3d[TheMask.mask].max()]
clim=CLIM_DICT[var]
integrated200=integrated*mask200 # taglio il costiero
integrated200[integrated200==0]=np.nan # sostituisco gli 0 con i NAN
#pl.set_cmap('gray_r') #changes the colormap
fig,ax = mapplot({'varname':var, 'clim':clim, 'layer':layer, 'data':integrated200, 'date':''},fig=None,ax=None,mask=TheMask,coastline_lon=clon,coastline_lat=clat)
ax.set_xlim([-5,36])
