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
C JULY 1982 DYNAMICS EXPLORER AND SME AND AIR FORCE SATELLITE STORM
C ANALYSIS
C 
C 1. Base steady state run before the storm,
C July 10,1982   day=82191 F107=176.,F107A=160. ,P=5 GW, PHI=30 KV
C DATA=ROBLE/RGR89/DVSSOL
C HISTORY=ROBLE/RGR90/JUST01,02 FROM 0,0,0 TO 1,0,0 (24 hour hist.)
C 				   1,0,0 to 2,0,0 ( 1 hour hist.)
C HISTORY=ROBLE/RGR90/JSTM01.......JSTM
C FORM 10,0,0 TO 19,0,0 WITH 1-HR TIME HISTORIES.
C
C As of Oct 2: RGR90/JSTM01 -> JSTM15
C
C   JSTM01: 10.00 -> 10:12 
C   JSTM02: 10.13 -> 11:01 
C   JSTM03: 11.02 -> 11:14 
C   JSTM04: 11.15 -> 12:03 
C   JSTM05: 12.04 -> 12:16 
C   JSTM06: 12.17 -> 13:05 
C   JSTM07: 13.06 -> 13:18 
C   JSTM08: 13.19 -> 14:07 
C   JSTM09: 14.08 -> 14:20 
C   JSTM10: 14.21 -> 15.09
C   JSTM11: 15.10 -> 15.22
C   JSTM12: 15.23 -> 16.11
C   JSTM13: 16.12 -> 17.00
C   JSTM14: 17.01 -> 17.13
C   JSTM15: 17.14 -> 18.00
C
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C
CLABEL = 'July 10-17 1982 (Hist JSTM01-15)'
 LABEL = 'July 13-14, 1982 (Storm)        '
 HISTVOL = 'ROBLE','RGR90','JSTM06','ROBLE','RGR90','JSTM07',
           'ROBLE','RGR90','JSTM08','ROBLE','RGR90','JSTM09',
           'ROBLE','RGR90','JSTM10'
CHISTVOL = 'ROBLE','RGR90','JSTM01','ROBLE','RGR90','JSTM02',
C          'ROBLE','RGR90','JSTM03','ROBLE','RGR90','JSTM04',
C          'ROBLE','RGR90','JSTM05','ROBLE','RGR90','JSTM06',
C          'ROBLE','RGR90','JSTM07','ROBLE','RGR90','JSTM08',
C          'ROBLE','RGR90','JSTM09','ROBLE','RGR90','JSTM10',
C          'ROBLE','RGR90','JSTM11','ROBLE','RGR90','JSTM12',
C          'ROBLE','RGR90','JSTM13','ROBLE','RGR90','JSTM14',
C          'ROBLE','RGR90','JSTM15'
C
CDDHHMM = '10:00:00','10:01:00','10:02:00','10:03:00','10:04:00',          
C         '10:05:00','10:06:00','10:07:00','10:08:00','10:09:00',      
C         '10:10:00','10:11:00','10:12:00','10:13:00','10:14:00',       
C         '10:15:00','10:16:00','10:17:00','10:18:00','10:19:00',       
C         '10:20:00','10:21:00','10:22:00','10:23:00',
C         '11:00:00','11:01:00','11:02:00','11:03:00','11:04:00',          
C         '11:05:00','11:06:00','11:07:00','11:08:00','11:09:00',      
C         '11:10:00','11:11:00','11:12:00','11:13:00','11:14:00',       
C         '11:15:00','11:16:00','11:17:00','11:18:00','11:19:00',       
C         '11:20:00','11:21:00','11:22:00','11:23:00',
C         '12:00:00','12:01:00','12:02:00','12:03:00','12:04:00',          
C         '12:05:00','12:06:00','12:07:00','12:08:00','12:09:00',      
C         '12:10:00','12:11:00','12:12:00','12:13:00','12:14:00',       
C         '12:15:00','12:16:00','12:17:00','12:18:00','12:19:00',       
C         '12:20:00','12:21:00','12:22:00','12:23:00',
C
 DDHHMM = '13:00:00','13:01:00','13:02:00','13:03:00','13:04:00',          
          '13:05:00','13:06:00','13:07:00','13:08:00','13:09:00',      
          '13:10:00','13:11:00','13:12:00','13:13:00','13:14:00',       
          '13:15:00','13:16:00','13:17:00','13:18:00','13:19:00',       
          '13:20:00','13:21:00','13:22:00','13:23:00',
          '14:00:00','14:01:00','14:02:00','14:03:00','14:04:00',          
          '14:05:00','14:06:00','14:07:00','14:08:00','14:09:00',      
          '14:10:00','14:11:00','14:12:00','14:13:00','14:14:00',       
          '14:15:00','14:16:00','14:17:00','14:18:00','14:19:00',       
          '14:20:00','14:21:00','14:22:00','14:23:00','15:00:00'
C
C         '15:00:00','15:01:00','15:02:00','15:03:00','15:04:00',          
C         '15:05:00','15:06:00','15:07:00','15:08:00','15:09:00',      
C         '15:10:00','15:11:00','15:12:00','15:13:00','15:14:00',       
C         '15:15:00','15:16:00','15:17:00','15:18:00','15:19:00',       
C         '15:20:00','15:21:00','15:22:00','15:23:00',
C         '16:00:00','16:01:00','16:02:00','16:03:00','16:04:00',          
C         '16:05:00','16:06:00','16:07:00','16:08:00','16:09:00',      
C         '16:10:00','16:11:00','16:12:00','16:13:00','16:14:00',       
C         '16:15:00','16:16:00','16:17:00','16:18:00','16:19:00',       
C         '16:20:00','16:21:00','16:22:00','16:23:00',
C         '17:00:00','17:01:00','17:02:00','17:03:00','17:04:00',          
C         '17:05:00','17:06:00','17:07:00','17:08:00','17:09:00',      
C         '17:10:00','17:11:00','17:12:00','17:13:00','17:14:00',       
C         '17:15:00','17:16:00','17:17:00','17:18:00','17:19:00',       
C         '17:20:00','17:21:00','17:22:00','17:23:00','18:00:00'
C
CIYD = 82191,82192,82193,82194,82195,82196,82197,82198
 IYD = 82194,82195
C
 IDEN = 2
 COLOR = 1
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 = 10.,15.,22.
 SPLS = -4.
 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,  1,  0,  0,  0,  0,  0,  0,  0,  0,  0,
C       FOF2 HMF2 RHO RAT 
          0,  0,  0,  0
C
C Fix contour intervals and/or min,max:
C At log10 z=-4: Ne=3.5->6.0,0.2; No=6.25->9.5,0.25
C
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.,6.85,1.,1.,  1., 1.,3.5, 1., 1., 1., 1.,
          1., 1., 1., 1.
 CONMAXS=0., 0., 0., 0., 0., 0.,9.6,0., 0.,  0., 0.,6.0, 0., 0., 0., 0.,
          0., 0., 0., 0.
 CONINTS=0., 0., 0., 0., 0., 0.,.25,0., 0.,  0., 0.,0.2, 0., 0., 0., 0.,
          0., 0., 0., 0.
 ENDOFREAD