C
C<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
C
C  INPUT DECK FOR TIGCM POST-MODEL PROCESSOR:
C  (mss executable = /FOSTER/binproc.ymp/tigcm.v1.o)
C
C NCASE = number of separate cases to be processed during this run
C         (full set of inputs needed for each case)
C LABEL = 32-character desciptive label for run (can leave blank)
C         (do NOT use colons (:))
C DYN = 1 -> dynamo case (for histories that have potential at 25 levels)
C HISTVOL(S) = history volumes (3 strings per volume)
C DDHHMM = Model day:hour:minute (md:mh:mm) where md,mh,mm are integers
C          and each time enclosed in quotes, e.g., '01:12:00'
C          (number of values usually multiple of 25)
C IYD = year-day for each day of run (e.g., 76080,76081,...)
C IDEN = flag for species densities:
C   0 -> species in number densities and total in mass density 
C   1 -> species and total in mass mixing ratio (as on history)
C   2 -> convert species and total to number densities (cm3)   
C   3 -> convert species and total to mass density (gm/cm3)
C LOG10PL = 1 -> plot log10 of fields where appropriate (densities)
C COLOR = 1 -> make color fill contours (otherwise, black and white)
C PLGLB = 1 -> Contour cylindrical equidistant global projection
C               at current ut and at selected pressures
C PLPOL = flag for stereographic (polar) projections at current
C          ut and at selected pressures as follows:
C   0 -> do not make polar plots
C   1 -> make polar plots for southern hemisphere only
C   2 -> make polar plots for northern hemisphere only
C   3 -> make polar plots for both hemispheres
C PLSATV = 1 -> Contour on satellite view projection
C               (use LATLON for projection orientation)
C PLLATZP = 1 -> Contour longitude slices (lat on x-axis, zp on y-axis)
C PLLATHT = 1 -> Contour longitude slices (lat on x-axis, ht on y-axis)
C PLUTLAT = 1 -> Contour ut on x-axis, latitude on y-axis
C PLUTZPM = 1 -> Contour global means with ut on x-axis, zp on y-axis
C                 (in this case, slon and spls are ignored)
C PLUTHTM = 1 -> Contour global means with ut on x-axis, height on y-axis
C                 (in this case, slon and spls are ignored)
C PLCOOL = 1 -> Include heating/cooling terms when iplglb=1,
C                iplpol=1, ipllon=1, or iplutlat=1
C SLON = selected longitudes (-180 to 180) (max of 5 lons)
C SSLT = selected solar local times (max of 5 times)
C SPLS = selected pressure levels (-7 to +5) (max of 5 pressures)
C SHTS = selected heights (will interpolate for global and polar plots)
C HTUB = upper boundary for linear height scale (must be multiple of 100.)
C        (if htub > htmax at grid point: zonal mean will extrapolate to
C         htub, but all others will leave plot blank above htmax)
C LATLON = lat,lon locations for center of satellite view projections
C          (lon = 999. means use local noon longitude)
C PLTIGCM -> choose fields to be contoured
C
C<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
C
 NCASE = 1
 ENDOFREAD
C
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C 3. EQUINOX WITH AURORA MEDIUM/LARGE
C    F107=150.,F107A=150.,POWER=16GW, PHI=70KV, DAY=91080
C    DATA=ROBLE/RGR87/IONSTS
C    SOURCE=ROBLE/RGR90/MARM03,3,0,0
C    HISTORY=ROBLE/RGR91/NICL02      FROM 0,0,0 TO 5,0,0  (24HR)
C    HISTORY=ROBLE/RGR91/NICP21,22   FROM 0,0,0 TO 1,0,0  (24HR)
C 6. SOLSTICE WITH AURORA MEDIUM/LARGE
C    F107=150.,F107A=150.,POWER=16GW, PHI=70KV, DAY=91355
C    DATA=ROBLE/RGR88/DVWSOL
C    SOURCE=ROBLE/RGR90/DECM03,3,0,0
C    HISTORY=ROBLE/RGR91/NICL05      FROM 0,0,0 TO 5,0,0  (24HR)
C    HISTORY=ROBLE/RGR91/NICP51,52   FROM 0,0,0 TO 1,0,0  (24HR)
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C
CLABEL = 'EQUINOX DAY 91080 (AUR MED/LAR) '
 LABEL = 'SOLSTICE DAY 91080 (AUR MED/LAR)'
 DYN = 0
CHISTVOL = 'ROBLE','RGR91','NICP22'
 HISTVOL = 'ROBLE','RGR91','NICP52'
C
 DDHHMM = '02:00:00'
CIYD = 91080
 IYD = 91355
 IDEN = 0
 LOG10PL = 0
 COLOR = 0
C
C Plots at single ut:
 PLGLB = 0
 PLPOL = 2
 PLSATV = 0
 PLLATZP = 0
 PLLATHT = 0
C
C Plots with ut on x-axis:
C (End with M refers to global means)
 PLUTLAT = 0
 PLUTZPM = 0
 PLUTHTM = 0
C
C Cooling/heating terms:
 PLCOOL = 0
C
C Selected grid:
 SLON = -105.
CSSLT = 0.,6.,12.,18.
CSSLT = 0.,12.
CSPLS = 2.
 SHTS = 120.
 HTUB = 600.
 LATLON = -70.0,999.,+70.,999.
C
C         T   U   V  O2   O N4S  NO  O+ N2D  TI  TE  NE O2+   W   Z  N2
 PLTIGCM= 1,  0,  0,  1,  1,  0,  1,  0,  0,  0,  0,  0,  0,  0,  0,  0,
C       FOF2 HMF2 RHO RAT 
          0,  0,  0,  0
C
C           POTEN UI VI WI SUMPOT HPOTEN HUI HVI HWI
 PLTIEGCM =   1,  1, 1, 1,   1,     1,    1,  1,  1
C
C23456789012345678901234567890123456789012345678901234567890123456789012 
C        T   U   V  O2   O N4S  NO  O+ N2D  TI  TE  NE O2+   W   Z  N2
C      FOF2 HMF2 RHO RAT 
 CONMINS=1., 1., 1., 1., 1., 1., 1.,1., 1., 1.,  1., 1., 1., 1., 1., 1.,
          1., 1., 1., 1.
 CONMAXS=0., 0., 0., 0., 0., 0., 0.,0., 0., 0.,  0., 0., 0., 0., 0., 0.,
          0., 0., 0., 0.
 CONINTS= 10., 0., 0.,0.1,0.1, 0.,0.1,0., 0., 0., 0., 0., 0., 0., 0.,0.,
          0., 0., 0., 0.
 ENDOFREAD
C
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C 1. BASE CASE , EQUINOX (SOLAR ONLY)
C   F107=150.,F107A=150.,POWER=0GW, PHI=0KV, DAY=91080
C   DATA=ROBLE/RGR87/IONSTS
C   SOURCE=ROBLE/RGR90/MARM03,3,0,0
C   HISTORY=ROBLE/RGR91/NICL00      FROM 0,0,0 TO 5,0,0  (24HR)
C   HISTORY=ROBLE/RGR91/NICP01,02   FROM 0,0,0 TO 1,0,0  (24HR)
C 4. BASE CASE , SOLSTICE (SOLAR ONLY)
C    F107=150.,F107A=150.,POWER=0GW, PHI=0KV, DAY=91355
C    DATA=ROBLE/RGR88/DVWSOL
C    SOURCE=ROBLE/RGR90/DECM03,3,0,0
C    HISTORY=ROBLE/RGR91/NICL03      FROM 0,0,0 TO 5,0,0  (24HR)
C    HISTORY=ROBLE/RGR91/NICP31,32   FROM 0,0,0 TO 1,0,0  (24HR)
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C
CLABEL = 'EQUINOX DAY 91080 (SOLAR ONLY)  '
 LABEL = 'SOLSTICE DAY 91355 (SOLAR ONLY) '
 DYN = 0
CHISTVOL = 'ROBLE','RGR91','NICP02'
 HISTVOL = 'ROBLE','RGR91','NICP32'
C
 DDHHMM = '02:00:00'
CIYD = 91080
 IYD = 91355
 IDEN = 0
 LOG10PL = 0
 COLOR = 0
C
C Plots at single ut:
 PLGLB = 0
 PLPOL = 2
 PLSATV = 0
 PLLATZP = 0
 PLLATHT = 0
C
C Plots with ut on x-axis:
C (End with M refers to global means)
 PLUTLAT = 0
 PLUTZPM = 0
 PLUTHTM = 0
C
C Cooling/heating terms:
 PLCOOL = 0
C
C Selected grid:
 SLON = -105.
CSSLT = 0.,6.,12.,18.
CSSLT = 0.,12.
CSPLS = 2.
 SHTS = 120.
 HTUB = 600.
 LATLON = -70.0,999.,+70.,999.
C
C         T   U   V  O2   O N4S  NO  O+ N2D  TI  TE  NE O2+   W   Z  N2
 PLTIGCM= 1,  0,  0,  1,  1,  0,  1,  0,  0,  0,  0,  0,  0,  0,  0,  0,
C       FOF2 HMF2 RHO RAT 
          0,  0,  0,  0
C
C           POTEN UI VI WI SUMPOT HPOTEN HUI HVI HWI
 PLTIEGCM =   1,  1, 1, 1,   1,     1,    1,  1,  1
C
C23456789012345678901234567890123456789012345678901234567890123456789012 
C        T   U   V  O2   O N4S  NO  O+ N2D  TI  TE  NE O2+   W   Z  N2
C      FOF2 HMF2 RHO RAT 
 CONMINS=1., 1., 1., 1., 1., 1., 1.,1., 1., 1.,  1., 1., 1., 1., 1., 1.,
          1., 1., 1., 1.
 CONMAXS=0., 0., 0., 0., 0., 0., 0.,0., 0., 0.,  0., 0., 0., 0., 0., 0.,
          0., 0., 0., 0.
 CONINTS= 10., 0., 0., .1,.1,0.,.05,0., 0., 0., 0., 0., 0., 0., 0.,0.,
          0., 0., 0., 0.
 ENDOFREAD
C
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C 2. EQUINOX WITH AURORA SMALL
C    F107=150.,F107A=150.,POWER=5GW, PHI=30KV, DAY=91080
C    DATA=ROBLE/RGR87/IONSTS
C    SOURCE=ROBLE/RGR90/MARM03,3,0,0
C    HISTORY=ROBLE/RGR91/NICL01      FROM 0,0,0 TO 5,0,0  (24HR)
C    HISTORY=ROBLE/RGR91/NICP11,12   FROM 0,0,0 TO 1,0,0  (24HR)
C  5. SOLSTICE WITH AURORA SMALL
C   F107=150.,F107A=150.,POWER=5GW, PHI=30KV, DAY=91355
C   DATA=ROBLE/RGR88/DVWSOL
C   SOURCE=ROBLE/RGR90/DECM03,3,0,0
C   HISTORY=ROBLE/RGR91/NICL04      FROM 0,0,0 TO 5,0,0  (24HR)
C   HISTORY=ROBLE/RGR91/NICP41,42   FROM 0,0,0 TO 1,0,0  (24HR)
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C
CLABEL = 'EQUINOX DAY 91080 (AURORA SMALL)'
 LABEL = 'SOLSTICE DAY 91355 AURORA SMALL'
 DYN = 0
CHISTVOL = 'ROBLE','RGR91','NICP12'
 HISTVOL = 'ROBLE','RGR91','NICP42'
C
 DDHHMM = '02:00:00'
CIYD = 91080
 IYD = 91355
 IDEN = 0
 LOG10PL = 0
 COLOR = 0
C
C Plots at single ut:
 PLGLB = 0
 PLPOL = 2
 PLSATV = 0
 PLLATZP = 0
 PLLATHT = 0
C
C Plots with ut on x-axis:
C (End with M refers to global means)
 PLUTLAT = 0
 PLUTZPM = 0
 PLUTHTM = 0
C
C Cooling/heating terms:
 PLCOOL = 0
C
C Selected grid:
 SLON = -105.
CSSLT = 0.,6.,12.,18.
CSSLT = 0.,12.
CSPLS = 2.
 SHTS = 120.
 HTUB = 600.
 LATLON = -70.0,999.,+70.,999.
C
C         T   U   V  O2   O N4S  NO  O+ N2D  TI  TE  NE O2+   W   Z  N2
 PLTIGCM= 1,  0,  0,  1,  1,  0,  1,  0,  0,  0,  0,  0,  0,  0,  0,  0,
C       FOF2 HMF2 RHO RAT 
          0,  0,  0,  0
C
C           POTEN UI VI WI SUMPOT HPOTEN HUI HVI HWI
 PLTIEGCM =   1,  1, 1, 1,   1,     1,    1,  1,  1
 ENDOFREAD