Python read_region_mask_popdiag Examples

Example #1

File: vert_zonal_avgs.py Project: NCAR/CESM_postprocessing

    myrend   = yrend - 181 + 1948

year_string = '%.0f-%.0f'%(myrstart,myrend)

#---------------------------------------------------------------------------

# set some plot resources
pl.rc('contour',negative_linestyle='solid') # make negative contour lines solid
pl.rc('font',**{'size':14})

#---------------------------------------------------------------------------

varlist = sorted(variables)

if POPDIAGPY == 'TRUE':
  reg_mask, nreg, reg_lname, reg_sname = read_region_mask_popdiag(mgrid)
else:
  reg_mask, nreg, reg_lname, reg_sname = read_region_mask('gx1v6')

# add global to region names
reg_sname.insert(0,'glo')
reg_lname.insert(0,'Global')
nreg = nreg + 1

if POPDIAGPY == 'TRUE':
    indir = WORKDIRPY
    outdir = WORKDIRPY
else:
    indir  = '/bonaire/data2/ivan/ccsm_output'
    outdir = '/fiji/home/ivan/www/ccsm_log/cases/%s/plots'%case

Example #2

    myrstart = yrstart - 181 + 1948
    myrend = yrend - 181 + 1948
    year_string_exp = '[%.0f-%.0f]' % (myrstart, myrend)

#---------------------------------------------------------------------------

# set some plot resources
pl.rc('contour', negative_linestyle='solid')  # solid negative contour lines
pl.rc('font', **{'size': 14})

#-----------------------------------------------------------------------------

varlist = sorted(variables)

if POPDIAGPY == 'TRUE':
    reg_mask, nreg, reg_lname, reg_sname = read_region_mask_popdiag(mgridexp)
else:
    reg_mask, nreg, reg_lname, reg_sname = read_region_mask('gx1v6')

# add global to region names
reg_sname.insert(0, 'glo')
reg_lname.insert(0, 'Global')
nreg = nreg + 1

if POPDIAGPY == 'TRUE':
    expdir = WORKDIRPY
    cntrldir = WORKDIRPY
    outdir = WORKDIRPY
else:
    indir = '/bonaire/data2/ivan/ccsm_output'
    outdir = '/fiji/home/ivan/www/ccsm_log/cases/%s/plots' % case