#
# QSUB -s /bin/csh
# QSUB -eo
# QSUB -lT 600
# QSUB -q prem
# QSUB -nr
# QSUB -lM 3MW
ja
cd $TMPDIR
set job = timesloc
batchname $job.$$.out
msread $job /FOSTER/bin/$job
chmod u+x $job
$job << 'EOF'
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 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
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 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 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 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 line plots)
C HT_SCALE = bottom,top,delta heights for y-axis, e.g., 100.,500.,10.
C            (used for PLUTVERT contours and UTLINE 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    If plotting e5577 at linear height scale (HT_SCALE) then doppler 
C    t,u,v and brightness will be calculated and line plots made 
C    (the doppler fields are vertically independent)
C   (e6300 is redline emissions, similiar to e5577 above)
C   (eo200 is nighttime O2 (0-0) atmospheric band emission, similiar to e5577)
C IONVEL = 1, 2, or 3 for calculation of ion drifts without efield,
C          with efield, or with electric potential, respectively
C  (note global or zonal means of ion drifts with ionvel=2 are not available)
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/RGR94/TS10DT2    ','/ROBLE/RGR94/TS10DT3    '
 DDHHMM = 25, 0, 0,25, 1, 0,25, 2, 0,25, 3, 0,25, 4, 0,25, 5, 0,
          25, 6, 0,25, 7, 0,25, 8, 0,25, 9, 0,25,10, 0,25,11, 0,
          25,12, 0,25,13, 0,25,14, 0,25,15, 0,25,16, 0,25,17, 0,
          25,18, 0,25,19, 0,25,20, 0,25,21, 0,25,22, 0,25,23, 0,
          26, 0, 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 = 0
 LOC =    42.6,-71.5
 LOCNAME = 'DURHAM N.H.     '
CLOCNAME = 'FT. COLLINS     '
CLOC =      42.5,-105.
CLOCNAME = 'HAUTE-PROVENCE  ','BISCARROSSE     '
CLOC =      44.,6.,           44.,-1.
CLOCNAME = 'FT. COLLINS+    ','FT. COLLINS+    ','FT. COLLINS     '
CLOC =      -57.5,-105.,      57.5,-105.,         42.5,-105.
CLOCNAME = 'HAUTE-PROVENCE+ ','BISCARROSSE+    ',
C          'HAUTE-PROVENCE+ ','BISCARROSSE+    '
CLOC =      59.,6.,           59.,-1.,
C           -59.,6.,          -59.,-1.
 BILIN = 1
 LOGLOC = 1
 ZP_RANGE = -17.,5.
CHT_SCALE = 80.,120.,1.
CHT_SCALE = 60.,150.,1.
CZPHTLINE = 97.
CUTLINE = 20,12,0
C
C utlines = 3->11 for Ft. Collins:
CUTLINES = 0, 3,0,  0, 4,0,  0, 5,0,  0, 6,0,  0, 7,0,  0, 8,0,  0, 9,0,
C          0,10,0,  0,11,0
C utlines = 1->4,20->0 for Haute-provence and Biscarrosse:
CUTLINES = 0, 1,0,  0, 2,0,  0, 3,0,  0, 4,0,  0,20,0,  0,21,0,  0,22,0,
C          0,23,0,  1, 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
 ZPHT_LON = -17.,-70., -4.,-70., 2.,-70.
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
 ZPHT_SLT = -17.,0., -4.,0., 2.,0.
 LOGUTLAT = 1
C
C          tn un vn o2 ox n4s noz co co2 h2o h2 hox o+ ch4 o21d
 FIELDS =   1, 1, 1, 1, 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 n2
            0,   0, 0, 1, 0,  0, 0, 0, 0,  0, 1, 0,  0, 0, 0,    1,
C          ui vi wi e5577 e6300 eo200
            0, 0, 0,  0,    0,    0
 IONVEL = 2
 WRASCII = 0
 DIRASCII = 'vishnu.hao.ucar.edu:/d/foster/timesloc/timesloc.dat     '
 DIRCGM =   'vishnu.hao.ucar.edu:/d/foster/timesloc/timesloc.cgm     '
CDIRCGM =   'vishnu.hao.ucar.edu:/d/foster/timesloc/t15diu.e5577.cgm '
CDIRCGM =   'vishnu.hao.ucar.edu:/d/foster/timesloc/tsfcexx.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 n2
C          1.,  1.,1., 1.,1.,1., 1.,1.,1.,1.,1.,1.,  1.,1.,1.
C          ui vi wi e5577 e6300 eo200
C          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 n2
C          0.,  0.,0., 0.,0.,0., 0.,0.,0.,0.,0.,0.,  0.,0.,0.
C          ui vi wi e5577 e6300 eo200
C          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 n2
C          0.,  0.,0., 0.,0.,0., 0.,0.,0.,0.,0.,0.,  0.,0.,0.
C          ui vi wi e5577 e6300 eo200
C          0.,0.,0., 0.,    0.,   0.
 ENDOFREAD
'EOF'
ja -st