latbounds = [90, 90]
nyrs = 0 # length of run in years
xx = -12 * nyrs # (until end in months)
l = False
levs = np.arange(0, 25, 1)
outdir = '/home/ih280/Analysis/fixedOH/plots/zonal_mean/individual_plots/'
if not os.path.exists(outdir):
os.makedirs(outdir)
# labels
ylab = r'Height / km'
ylim = [0, 20]
# labels
label = plots.labels()
# Colours
c = ccol.custom_colors('grads')
cmap = cm.bwr
cnorm = clrs.Normalize(cmap, clip=False)
cmap.set_under(color=cmap(0.0), alpha=1.0)
cmap.set_over(color=cmap(1.0), alpha=1.0)
#************************************************************
# TROPOPAUSE
trophgt = ncdf.Dataset(
'/' + disk + '/ih280/um/' + jobid + '/' + jobid + '_trophgt.nc', 'r')
tropheight = trophgt.variables['tropopause_height'] #time,lat,lon
trop = np.mean(np.mean(\
tropheight[xx:,:,:]\
unit = 'ppbv'
nyrs = 0 # entire run length
plottitle = jobs2.metas[jobs2.jobid.index(jobid)]
var = jobs2.variable[jobs2.stash.index(stash)]
# range for contour limits
#levs = np.arange(1.65,1.95,0.025)
plotname = jobid+stash+'_'+localtime+'.png'
outdir = 'xx'
if not os.path.exists(outdir):
os.makedirs(outdir)
#************************************************************
# label format
fmtlab = plots.labels()
# Colours
c=ccol.custom_colors('grads')
cmap = cm.bwr
cnorm=clrs.Normalize(cmap,clip=False)
cmap.set_under(color=cmap(0.0),alpha=1.0)
cmap.set_over(color=cmap(1.0),alpha=1.0)
#************************************************************
# Read in files
# time,model_level_number,latitude,longitude
file1 = '/'+disk+'/ih280/um/'+jobid+'/'+jobid+stash+'.nc'
ncbase = ncdf.Dataset(file1,'r')
data = ncbase.variables[var]
nyrs = 0 # entire run length
plottitle = jobid + '\nTropospheric Average'
var = jobs2.variable[jobs2.stash.index(stash)]
# range for contour limits
levs = np.arange(500, 1200, 100)
#levs = np.arange(0,700,100)
plotname = jobid + stash + '_' + localtime + '.png'
outdir = '/home/ih280/Analysis/fixedOH/plots/global_distribution/individual_plots/'
if not os.path.exists(outdir):
os.makedirs(outdir)
#************************************************************
# label format
fmtlab = plots.labels()
# Colours
c = ccol.custom_colors('grads')
cmap = cm.bwr
cnorm = clrs.Normalize(cmap, clip=False)
cmap.set_under(color=cmap(0.0), alpha=1.0)
cmap.set_over(color=cmap(1.0), alpha=1.0)
#***********************************************************
# Tropospheric mask: 1 for troposphere, 0 other
#INPUT:
#ncfile = file containing tropopause height
#LTROP = 1 for troposphere (i.e. mask out stratosphere); 0 otherwise
#tmean = 0 (no time meaning),
# 1 (climatological monthly mean),
latbounds = [90,90]
nyrs = 0 # length of run in years
xx = -12*nyrs # (until end in months)
l = False
levs = np.arange(0,25,1)
outdir ='/home/ih280/Analysis/fixedOH/plots/zonal_mean/individual_plots/'
if not os.path.exists(outdir):
os.makedirs(outdir)
# labels
ylab = r'Height / km'
ylim = [0,20]
# labels
label=plots.labels()
# Colours
c=ccol.custom_colors('grads')
cmap = cm.bwr
cnorm=clrs.Normalize(cmap,clip=False)
cmap.set_under(color=cmap(0.0),alpha=1.0)
cmap.set_over(color=cmap(1.0),alpha=1.0)
#************************************************************
# TROPOPAUSE
trophgt = ncdf.Dataset('/'+disk+'/ih280/um/'+jobid+'/'+jobid+'_trophgt.nc','r')
tropheight = trophgt.variables['tropopause_height'] #time,lat,lon
trop = np.mean(np.mean(\
tropheight[xx:,:,:]\
, axis=2, dtype=np.float64)\