C
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C
C INPUT FILE FOR TIMESLOC CRAY PROCESSOR (zp = -7 to +5, kmx=25):
C (binary executable is on mss /FOSTER/bin/timesloc)
C (source is on sting:~foster/timesloc and mss /FOSTER/src/timesloc.tar)
C
C This is a multi-ut processor, i.e., most plots are made with ut on x-axis
C
C This processor will accept tigcm or tiegcm histories, however:
C If electric potential or ion drift plots are requested (pot,ui,vi,wi),
C   and the history is tigcm (not tiegcm), a notice will be printed 
C   to stdout (unit 6), and potential and/or ion drifts will not be plotted.
C
C All values below DDHHMM except FIELDS may be commented out to get the
C   default value, or to suppress the plotting option. Always comment out 
C   in column one with an upper case C.
C All uncommented keyword = value lexical read lines must start in 
C   column two and be <= 72 columns per line. (you can use the ruler at 
C   CMIN,CMAX,CINT for reference). New lines may be used if a comma 
C   appears at the end of the previous line.
C Do not use tabs as white space in lexical reads.
C
C HISTVOL = timesgcm history volume(s) (24-char mss path per volume)
C           e.g., '/ROBLE/RGR92/TISED1     ','/ROBLE/RGR92/TISED2     '
C DDHHMM = model day,hour,min of each desired history 
C   (3 comma-separated integers per history, e.g., 3,0,0, 3,1,0)
C   (can use 99 as one of the integers to determine histories on volume;
C    99 will not be found, and times found will be printed to stdout)
C   (HISTVOL and DDHHMM are required inputs)
C RMVOLS = 1 -> remove history volumes from disk after reads (default = 1)
C IDEN = 0 -> leave species as on history (most are mass mixing ratios) 
C IDEN = 1 -> convert species to number densities (cm3)   
C IDEN = 2 -> convert species to mass density (gm/cm3) 
C             (except ne, o+, and o2+ which remain number density)
C   (default IDEN = 1) 
C   Note o+, ne, and o2+ are number densities on the history, and will
C     be left as such (cm-3) even when IDEN = 0
C COLOR = 1 -> make color fill contours, monochrome otherwise
C   (default COLOR = 0)
C LABELS = 1,2 for simple, normal labeling (default is 2)
C PLUTVERT = 1 -> contour with ut on x-axis, zp/ht on y-axis at LOC or LOCSLT
C   locations. ZP_RANGE and/or HT_SCALE are used to determine y-axis
C   (set PLUTVERT=0 for no such contours)
C LOC = lat,lon pairs of selected locations. May use lat,999. to get
C       zonal means at selected lat, or 999.,999. to get global means
C       (for PLUTVERT, UTLINE, UTLINES, ZPHTLINES)
C LOCNAME = 16-char location name. Optional, but if given, there must be
C           one name per lat,lon pair in LOC (will appear on plot labels)
C           (blank or dummy names may be used for zonal or global mean locs)
C LOCSLT = lat,slt pairs at which to contour zp/ht vs ut. Like LOC, except
C          at fixed local time rather than longitude (LOCNAME not used)
C          (for PLUTVERT, UTLINE, UTLINES, ZPHTLINES)
C NOTE: YOU CANNOT SPECIFY BOTH LOC and LOCSLT IN THE SAME RUN
C BILIN = 1 -> do bilinear interpolation to desired location; if BILIN=0,
C       use nearest grid point. (No interp is done for global means.
C       If doing zonal means and BILIN=1, interp is done to desired lat)
C       (if LOCSLT is set, bilin interps to the slt longitudes)
C LOGLOC = 1 -> plot log10 of density fields at locations, otherwise no
C ZP_RANGE = bottom,top range of pressures for y-axis (-7 -> +5)
C            (e.g.: ZP_RANGE= -4.,2.)
C            (used for PLUTVERT contours and UTLINE and UTLINES line plots)
C HT_SCALE = bottom,top,delta heights for y-axis, e.g., 100.,500.,10.
C            (used for PLUTVERT contours and UTLINE and UTLINES line plots)
C ZPHTLINE = selected zp(s) and/or ht(s) at which to make line plots at
C            LOC locations with ut on x-axis, field on y-axis
C            (is assumed a height if > 5., otherwise is assumed zp)
C UTLINE = selected times (integer triplets) at which to make line plots
C          at LOC locations with field on x-axis, zp/ht on y-axis
C          (must be subset of DDHHMM) (ZP_RANGE and/or HT_SCALE will be
C          used for y-axis)
C UTLINES = like UTLINE except plot lines for each specified ut on the
C           same frame for each field
C ZPHT_LON = pairs of zp/ht,lon for contours with ut on x-axis, lat on y-axis
C (first of each pair may be pressure or height, 2nd of each pair may be
C  a longitude, or 999. for zonal means)
C ZPHT_SLT = like ZPHT_LON, except second of each pair is local time instead
C  of longitude (i.e., latitude vs ut at selected zp/ht and fixed local 
C  times)
C
C          --------------- FIELDS AND CONTOUR VALUES ----------------
C
C FIELDS = fields to plot (integers 0 or 1) 
C   (this is required, and must provide a value for each field)
C   (rho is total density o2 + o1 + n2)
C   (fof2 and hmf2 will not be plotted where zp/ht is on y-axis)
C   (e5577 is volume emission rate for greenline (5577A))
C   (e6300 is redline emissions, similiar to e5577 above)
C   (eo200 is nighttime O2 (0-0) atmospheric band emission, similiar to e5577)
C   If plotting emissions at linear height scale (HT_SCALE) then 
C     doppler t,u,v and brightness will be calculated and line plots made 
C     (the doppler fields are vertically independent)
C IONVEL = flag for calculation of ion drift vel (for ui,vi,wi only)
C   = 1 -> calculate ExB drifts 
C   = 2 -> calculate magnetic ExB + effects of neutral atmos
C WRASCII = 0 -> do not make ascii data file(s)
C WRASCII = 1 -> make an ascii data file for each frame
C   (file names will be tigcmproc.fr001.dat, tigcmproc.fr002.dat, ...)
C WRASCII = 2 -> make single ascii file containing data from each frame
C   (file name will be tigcmproc.dat)
C DIRASCII = 56-char remote directory to which ascii data file(s) may be
C            remote copied, e.g.:
C e.g., DIRASCII = 'sting.hao.ucar.edu:/d/foster/timesloc                   '
C DIRCGM = 56-char remote directory to which metacode file (gmeta) may be
C          remote copied.
C
C CMIN = contour minimums (if > CMAX, will use min of plotted array)
C CMAX = contour maximums (if < CMIN, will use max of plotted array)
C CINT = contour intervals (if = 0., will choose interval in software)
C   (CMIN,CMAX,CINT are floating point values)
C   (defaults for CMIN,CMAX,CINT = 1.,0.,0. for all fields)
C   
C   Note CMIN must be < CMAX, and CINT > 0. for CMIN,CMAX to take effect
C        however, CINT > 0. will take effect regardless of CMIN,CMAX
C   CMIN,CMAX,CINT apply to all plots. Remember that these will change
C   radically on mapped projections depending on LOG10. Note also that a
C   CINT that is appropriate for mapped projections may not be appropriate
C   for vertical slices.
C
C Please direct comments, bugs, and suggestions for additional options to:
C Ben Foster (hao) (foster@ncar.ucar.edu) 303-497-1595
C
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C
 HISTVOL = '/ROBLE/RGR95/TSCLX16    '
 DDHHMM = 110,0,0, 110,12,0, 111,0,0, 111,12,0
 RMVOLS = 0
C
C IDEN=0 -> mass mixing ratios; IDEN=1 -> cm3; IDEN=2 -> gm/cm3
C
 IDEN = 1
 MODELHTS = 1
 COLOR = 0
 LABELS = 2
C
 PLUTVERT = 1
C
C Note: cannot do both LOC and LOCSLT in same run
C       (LOCNAME not used if LOCSLT is set)
C       (if LOC is set, number of LOCNAME and LOC pairs must be the same)
C
CLOCSLT =    42.6,0., 42.6,6., 42.6,12., 42.6,18.
 LOC = -35.,138., 18.3,-66.8
 LOCNAME = 'ADELAIDE        ','ARECIBO         '
 BILIN = 1
 LOGLOC = 1
 ZP_RANGE = -15.5,5.
CHT_SCALE = 30.,200.,5.
CHT_SCALE = 60.,150.,1.
C
C ZPHTLINE = selected zp or hts for line plots of field vs ut:
CZPHTLINE = 120.
C
C UTLINE(S) = selected ut's for line plots of zp/ht vs LOC or LOCSLT
CUTLINE = 20,12,0
CUTLINES = 25,0,0, 25,6,0, 25,12,0, 25,18,0, 26,0,0 
C
C ZPHT_LON = pairs of zp/ht,lon for contours with ut on x-axis, lat on y-axis
C (first of each pair may be pressure or height, 2nd of each pair may be
C  a longitude or 999. for zonal means)
C
CZPHT_LON = 120.,-70., 120.,0.
C
C ZPHT_SLT = like ZPHT_LON above, except 2nd of each pair is local time
C rather than longitude (i.e., lat vs ut at selected zp/ht and local times)
C
CZPHT_SLT = 120.,0., 120.,6., 120.,12., 120.,18.
 LOGUTLAT = 1
C
C          tn un vn o2 ox n4s noz co co2 h2o h2 hox o+ ch4 o21d
 FIELDS =   1, 1, 1, 0, 0,  0,  0, 0,  0,  0, 0,  0, 0,  0, 0,
C          no2  no o3  o1 oh ho2 h n2d ti te ne o2+  w  z  poten
            0,   0, 0, 1, 0,  0, 0, 0, 0,  0, 1, 0,  0, 0, 0,
C          ui vi wi  n2 e5577 e6300 eo200
            0, 0, 0,  1,  1,    1,    1
 IONVEL = 2
 WRASCII = 0
 DIRASCII = 'vishnu.hao:/d/foster/timesloc/timeslocz/timeslocz.dat   '
 DIRCGM =   'vishnu.hao:/d/foster/timesloc/timeslocz/timeslocz.cgm   ' 
C
C23456789012345678901234567890123456789012345678901234567890123456789012 
C          tn un vn o2 ox n4s noz co co2 h2o h2 hox o+ ch4 o21d
CCMIN =    1.,1.,1.,1.,1.,1., 1., 1.,1., 1., 1.,1., 1.,1., 1., 
C          no2  no o3  o1 oh ho2 h n2d ti te ne o2+  w  z 
C          1.,  1.,1., 1.,1.,1., 1.,1.,1.,1.,1.,1.,  1.,1.,
C          ui vi wi n2 e5577 e6300 eo200
C          1.,1.,1., 1., 1.,   1.,   1.
C          tn un vn o2 ox n4s noz co co2 h2o h2 hox o+ ch4 o21d
CCMAX =    0.,0.,0.,0.,0.,0., 0., 0.,0., 0., 0.,0., 0.,0., 0., 
C          no2  no o3  o1 oh ho2 h n2d ti te ne o2+  w  z
C          0.,  0.,0., 0.,0.,0., 0.,0.,0.,0.,0.,0.,  0.,0.,
C          ui vi wi n2 e5577 e6300 eo200
C          0.,0.,0., 0., 0.,    0.,   0.
C          tn un vn o2 ox n4s noz co co2 h2o h2 hox o+ ch4 o21d
CCINT =    0.,0.,0.,0.,0.,0., 0., 0.,0., 0., 0.,0., 0.,0., 0., 
C          no2  no o3  o1 oh ho2 h n2d ti te ne o2+  w  z
C          0.,  0.,0., 0.,0.,0., 0.,0.,0.,0.,0.,0.,  0.,0.,
C          ui vi wi n2 e5577 e6300 eo200
C          0.,0.,0., 0., 0.,    0.,   0.
 ENDOFREAD