コード例 #1
0
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:,:,:]\
コード例 #2
0
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]
コード例 #3
0
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),
コード例 #4
0
ファイル: zonal_mean_plot_84.py プロジェクト: jws52/iris_pyle
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)\