;
&proc_input
;
;histvols = '$TGCMDATA/TGCM.tiegcm.p001-2002-080.nc'
;histvols = 'BOUGHER.SWBV97.EQUMIN.nc'
;histvols = 'BOUGHER.SWBV97.EQUMED.nc'
;histvols = 'BOUGHER.SWBV97.EQUMAX.nc'
;mtimes = 5,12,0
;histvols = 'FOSTER.vtgcm.pequmin.nc'
;histvols = 'FOSTER.vtgcm.pequmin.nc'
;histvols = 'FOSTER.vtgcm.sequmin.nc'
 histvols = 'FOSTER.vtgcm.ptest001.nc'
;histvols = 'FOSTER.vtgcm.stest001.nc'
 mtimes = 7,0,0
;mtimes = 6,0,0,'to',7,0,0,'by',60
 vtgcm_fspath = '/fis/hao/home/bougher/VTGCM/VN2SNR'
;
maxdiskvols = 365
 cfields = 'TN','UN','VN','O2','O1','CO','N4S','N2D','ENO','W','Z'
 grid_levels = 'INTERFACES'
;
; addiag.F:
;cfields = 'TN','VN','VNa','diag_vc','N2a','barm1','barm2',
;          'TNa','W1a','W1b','ADIAG_Z'
;
; modelhts = 0 -> calculate geopotential height (non-constant gravity)
; modelhts = 1 -> use geopotential height on the history
;
modelhts = 0
icolor = 1
;
; iden=0 -> leave species as on history (most are mass mix ratios)
; iden=1 -> convert species to number densities (cm3)
; iden=2 -> convert species to number density mixing ratios
; iden=3 -> convert species to mass density (gm/cm3)
;
;iden = 1
 iden = 0
iplot = 1
outplt = 'cgm','ps'
;
; Make maps (global and/or polar):
;
ipltmaps = 1
  map_global = 1
  map_global_cenlon = 0
; map_global_censlt = 12
  map_polar = 0
  fmap_polar_perimlat = -40.,40.
; fmap_zpht = -10., -7., -4., 0., 2., 6.75, 7.25
  fmap_zpht = -10., -7., -4., 0., 2., 6.
  map_tn_unvn = 1
;
; Make longitude slices (lat vs zp/ht at selected lons):
;
ipltlon = 1
  flons = 'zm'
  fslts = 0.,12.
  flon_zprange = -99.,6.25
; flon_zprange = -99.,99.
; flon_htscale = 75.,150.,5.
;
; ilon_log10 = 0 -> do not plot log10 of any field
; ilon_log10 = 1 -> plot log10 of species densities only (default)
; ilon_log10 = 2 -> plot log10 of all fields
  ilon_log10=2		; plot log10 sp if ilat_log10 > 0
;
; Make latitude slices (lon vs zp/ht at selected lats):
;
ipltlat = 0
  flats = -40.,-60.,40.,60.
  flat_zprange = -99.,99.
; flat_htscale = 75.,150.,5.
;
; Make xy location line plots (vertical profiles):
;
ipltxyloc = 0
;
; xylocs = lat,lon locations
; Note 'zm' means plot zonal means
; Note 'zm','zm' means plot global means
;
  xylocs = 40.,'zm',60.,'slt12.',-60.,180.,'zm','zm'
  xyloc_locname = 'loc1','loc2','Global Mean'
  xyloc_zprange = -999.,999.
; xyloc_htscale = 50.,500.,10.
;
; Make ut vs zp/ht contours at selected locations:
;
ipltutvert = 0
  utvert_locs = 2.0,-157., 20.8,-156.5, 41.9,-111.4, 69.6,19.2
  utvert_locname = 'CHRISTMAS','MAUI','BEAR LAKE','EISCAT'
  utvert_zprange = -999.,999.
; utvert_htscale = 50.,100.,1.
;
; Make ut vs latitude contours at selected zp/ht and longitudes
;
ipltutlat = 0
  utlat_zphtlon = -7.,0., -7.,180., 2.,0., 2.,180.
;
; Make ut vs field line plots at selected locations:
;
ipltxyut = 0
  xyut_locs = 40.,-90., 40.,'zm', 40.,'slt12.', 'zm','zm'
; xyut_locname = 'LOC 1','LOC 2','LOC 3','LOC 4'
  xyut_zpht = -4.,120.,'hmf2'
;
; Example destinations for scp of cgm and ps output plot files:
; (If sendcgm and sendps are not provided, plots will be output
;   to local files tgcmproc.cgm and tgcmproc.ps)
;
;sendcgm = 'machine.hao:/e/foster/tgcmproc/tgcmproc.cgm'
;sendps  = 'machine.hao:/e/foster/tgcmproc/tgcmproc.ps'
;
 flnm_cgm = 'vtgcm.cgm'
 flnm_ps  = 'vtgcm.ps'
/