C
C<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
C
C  INPUT DECK FOR TIGCM POST-MODEL PROCESSOR:
C
C NCASE = number of separate cases to be processed during this run
C         (full set of inputs needed for each case)
C
C LABEL = 32-character desciptive label for run (can leave blank)
C         (do NOT use colons (:))
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 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 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 PLTIGCM -> choose fields to be contoured
C
C<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
C
 NCASE = 1
 ENDOFREAD
C
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C Case 1:
C August 21,1979 solar maximum, day=79233, F107=243.,F107A=215.,P=5GW,
C                                  PHI=30KV
C DATA=ROBLE/RGR89/DVSSOL
C HISTORY=ROBLE/RGR90/AUGX03,04 FROM 2,0,0 to 3,0,0 (24 hour hist.)
C                                    3,0,0 to 4,0,0 (1  hour hist.)
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C
 LABEL = 'August 21, 1979 Solar Max       ' 
 HISTVOL = 'ROBLE','RGR90','AUGX03','ROBLE','RGR90','AUGX04'
CHISTVOL = 'ROBLE','RGR90','AUGX04'
CDDHHMM = '04:00:00'
 DDHHMM = '03:00:00','03:01:00','03:02:00','03:03:00','03:04:00',    
          '03:05:00','03:06:00','03:07:00','03:08:00','03:09:00'
C         '03:10:00','03:11:00','03:12:00','03:13:00','03:14:00',       
C         '03:15:00','03:16:00','03:17:00','03:18:00','03:19:00',       
C         '03:20:00','03:21:00','03:22:00','03:23:00','04:00:00'
 IYD = 79233
C
 IDEN = 2
 COLOR = 0
C
C Plots at single ut:
 PLGLB = 0
 PLPOL = 0
 PLSATV = 1
 PLLATZP = 0
 PLLATHT = 0
C
C Plots with ut on x-axis:
C (End with M means global means)
 PLUTLAT = 0
 PLUTZPM = 0
 PLUTHTM = 0
C
C Cooling/heating terms:
 PLCOOL = 0
C
C Selected grid:
CSLON = -70.
CSSLT = 3.,15.
 SPLS = -4.
CSHTS = 200.
 LATLON = 67.5,999.
C
C         T   U   V  O2   O N4S  NO  O+ N2D  TI  TE  NE O2+   W   Z  N2
 PLTIGCM= 0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  1,  0,  0,  0,  0,
C       FOF2 HMF2 RHO RAT 
          0,  0,  0,  0
 CONINTS= 50.,19:0.
 ENDOFREAD
C
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C Case 2:
C 10.October 21,1979 solar maximum, day=79294, F107=243.,F107A=215.,P=5GW,
C                                  PHI=30KV
C DATA=ROBLE/RGR89/DVSSOL
C HISTORY=ROBLE/RGR90/OCTX03,04 FROM 2,0,0 to 3,0,0 (24 hour hist.)
C                                    3,0,0 to 4,0,0 (1  hour hist.)
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C
 LABEL = 'October 21, 1979 Solar Max      '
 HISTVOL = 'ROBLE','RGR90','OCTX04'
C
 DDHHMM = '04:00:00'
 IYD = 79294
C
 IDEN = 2
 COLOR = 0
C
C Plots at single ut:
 PLGLB = 1
 PLPOL = 3
 PLLATZP = 0
 PLLATHT = 0
C
C Plots with ut on x-axis:
C (End with M means global means)
 PLUTLAT = 0
 PLUTZPM = 0
 PLUTHTM = 0
C
C Cooling/heating terms:
 PLCOOL = 0
C
C Selected grid:
CSLON = -70.
 SSLT = 3.,15.
 SPLS = 1.,+2.
C
C         T   U   V  O2   O N4S  NO  O+ N2D  TI  TE  NE O2+   W   Z  N2
 PLTIGCM= 1,  1,  1,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,
C       FOF2 HMF2 RHO RAT 
          0,  0,  0,  0
 ENDOFREAD
C
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C Case 3:
C 7.March 21,1979 solar maximum, day=79080, F107=243.,F107A=215.,P=5GW,
C                                  PHI=30KV
C DATA=ROBLE/RGR87/IONSTS
C HISTORY=ROBLE/RGR90/MARX03,04 FROM 2,0,0 to 3,0,0 (24 hour hist.)
C                                    3,0,0 to 4,0,0 (1  hour hist.)
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C
 LABEL = 'March 21, 1979 Solar Max        '
 HISTVOL = 'ROBLE','RGR90','MARX04'
C
 DDHHMM = '04:00:00'
 IYD = 79080
C
 IDEN = 2
 COLOR = 0
C
C Plots at single ut:
 PLGLB = 1
 PLPOL = 3
 PLLATZP = 0
 PLLATHT = 0
C
C Plots with ut on x-axis:
C (End with M means global means)
 PLUTLAT = 0
 PLUTZPM = 0
 PLUTHTM = 0
C
C Cooling/heating terms:
 PLCOOL = 0
C
C Selected grid:
CSLON = -70.
 SSLT = 3.,15.
 SPLS = 1.,+2.
C
C         T   U   V  O2   O N4S  NO  O+ N2D  TI  TE  NE O2+   W   Z  N2
 PLTIGCM= 1,  1,  1,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,
C       FOF2 HMF2 RHO RAT 
          0,  0,  0,  0
 ENDOFREAD