;
&proc_input
;
; Days 81-106 (daily):
;
;histvols = 'TGCM.tiegcm1.95.p_whi2008_heelis_gpi_001.nc','to', 
;           'TGCM.tiegcm1.95.p_whi2008_heelis_gpi_003.nc','by','1' 
;mtimes = 81,0,0,'to',106,0,0,'by',1440
;flnm_ps  = 's_whi2008_heelis_gpi_daily.ps'
;
 histvols = 'TGCM.tiegcm1.95.p_whi2008_weimer_imf_001.nc','to', 
            'TGCM.tiegcm1.95.p_whi2008_weimer_imf_003.nc','by','1' 
 mtimes = 81,0,0,'to',106,0,0,'by',1440
 flnm_ps  = 's_whi2008_weimer_imf_daily.ps'
;
;
;cfields = 'TN','UN','VN','O2','O1','N2','N4S','NO','OP','N2D',
;          'TI','TE','NE','O2P','OMEGA','POTEN','Z'
 cfields = 'TN','NE'

 fmnmxint = 'TN',0.,0.,50., 'NE',0.,0.,0.10
;
; 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
iplot = 1
outplt = 'ps'
;
; Local ascii output file:
;flnm_dat = 'tgcmproc.dat'
;
; Make maps (global and/or polar):
;
ipltmaps = 0
  map_global = 1
  map_polar = 0
  fmap_polar_perimlat = -40.,40.
  fmap_zpht = -4., 2.
;
; Make longitude slices (lat vs zp/ht at selected lons):
;
ipltlon = 0
  flons = 'zm'
  fslts = 12.
  flon_zprange = -99.,99.
; flon_htscale = 75.,150.,5.
;
; Make latitude slices (lon vs zp/ht at selected lats):
;
ipltlat = 0
  flats = -60.,-40.,0.,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 = 'zm','zm'
  xyloc_locname = 'Global Mean'
  xyloc_zprange = -999.,999.
; xyloc_htscale = 50.,500.,10.
;
; Make lat vs zp/ht plots of zm at selected slt
; (must have 1 day time series):
;
ipltzmslt = 0
; zmslt = 0.,6.,12.,18.
zmslt = 0.,12.
zmslt_zprange = -999.,999.
zmslt_htscale = 90.,400.,5.
;
; Make ut vs zp/ht contours at selected locations:
;
ipltutvert = 1
  utvert_locs = 0.,180., 60.,180.
; 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 = 1
  utlat_zphtlon = -4.,180., 2.,180.
;
; Make ut vs field line plots at selected locations:
;
ipltxyut = 1
  xyut_locs = 'zm','zm'
; xyut_locname = 'LOC 1','LOC 2','LOC 3','LOC 4'
  xyut_zpht = -4.,2.
/