Geoff:

I have added o/n2 to the new location processor (this was easier
to do than debug the older code). Below is a script you can use
for the new code. Note that it does the msread on 
/FOSTER/bin/tigcmloc.dev (this is the new file I've written
this morning). I've set this script up to use the same histories
you used on the bad job with the old code. Note also that
IDEN=2 (rather than 3) to get gm/cm3 with the new code.
Change the metapath in the script if you need to.

This new code also has capability to write data files
(see WRASCII, DIRASCII, etc). There are several other options
the old code did not have -- read the comments. It also makes
much better time axis labels, etc.

Let me know if you have problems with this...

--Ben


#
# QSUB -s /bin/csh
# QSUB -eo
# QSUB -lT 600
# QSUB -q prem
# QSUB -nr
# QSUB -lM 5MW
ja
cd $TMPDIR
set job = tigcmloc
set metapath = sting.hao.ucar.edu:/d/geoff/$job/$job.cgm
msread $job /FOSTER/bin/$job.dev
chmod u+x $job
pshell $job << 'EOF'
C
Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C
C INPUT FILE FOR TIGCMLOC CRAY PROCESSOR (zp = -7 to +5, kmx=25):
C (binary executable is on mss /FOSTER/bin/tigcmloc)
C (source is on sting:~foster/tigcmloc and mss /FOSTER/src/tigcmloc.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 = tigcm or tiegcm history volume(s) (24-char mss path per volume)
C           e.g., '/ROBLE/RGR92/LTCS01     ','/ROBLE/RGR92/LTCS02     '
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 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 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               (applies to contours of ut vs zp/ht)
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 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 LOGUTLAT = 1 -> plot log10 of density fields at locations, otherwise no
C               (applies to contours of ut vs lat)
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 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/tigcmloc                   '
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/RGR90/DECX23     ','/ROBLE/RGR90/DECX24     '
 DDHHMM = 3,0,0, 3,1,0, 3,2,0, 3,3,0, 3,4,0, 3,5,0, 3,6,0, 3,7,0,
          3,8,0, 3,9,0, 3,10,0, 3,11,0, 3,12,0, 3,13,0, 3,14,0, 
          3,15,0, 3,16,0, 3,17,0, 3,18,0, 3,19,0, 3,20,0, 3,21,0, 
          3,22,0, 3,23,0, 4,0,0
C
C IDEN=0 -> mass mixing ratios; IDEN=1 -> cm3; IDEN=2 -> gm/cm3
C
 IDEN = 2
 COLOR = 0
C
 PLUTVERT = 1
CLOC = 57.5,-70., 57.5,999., 60.,999., 999.,999.
 LOC = 74.7,-95.
 BILIN = 0
 LOGLOC = 0
CZP_RANGE = -7.,5.
 HT_SCALE = 90.,500.,10.
CZPHTLINE = 2.
CUTLINE = 2,18,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 = 2.,-70., -4.,999., 200.,-70., 200.,999. 
CZPHT_LON = 2.,-72.,2.,-70.
 LOGUTLAT = 1
C
C          tn  un  vn  o2  o1  n4s  no  o+  n2d  ti  te  ne  o2+  w  z
 FIELDS =   0,  0,  0,  0,  0,   0,  0,  0,   0,  0,  0,  0,   0, 0, 0,
C          pot  n2  ui  vi  wi fof2 hmf2 rho o/(o2+n2) o/n2
             0,  0,  0,  0,  0,   0,   0,  0,        1,   1
 WRASCII = 0
 DIRASCII = 'sting.hao.ucar.edu:/d/foster/tigcmloc                   '   
CDIRCGM =   'sting.hao.ucar.edu:/d/foster/tigcmloc                   '
C
C23456789012345678901234567890123456789012345678901234567890123456789012 
C        tn     un     vn      o2       o1      n4s      no       o+  
CCMIN =   1.,    1.,    1.,     1.,      1.,      1.,    1.,      1.,
C        n2d    ti     te      ne       o2+      w       z       pot  
C         1.,    1.,    1.,     1.,      1.,      1.,    1.,      1.,
C        n2     ui     vi      wi      fof2     hmf2     rho o/(o2+n2)
C         1.,    1.,    1.,     1.,      1.,      1.,    1.,      1.,
C        o/n2
C         1.
C
C23456789012345678901234567890123456789012345678901234567890123456789012 
C        tn     un     vn      o2       o1      n4s      no       o+  
CCMAX =   0.,    0.,    0.,     0.,      0.,      0.,    0.,      0.,
C        n2d    ti     te      ne       o2+      w       z       pot  
C         0.,    0.,    0.,     0.,      0.,      0.,    0.,      0.,
C        n2     ui     vi      wi      fof2     hmf2     rho o/(o2+n2)
C         0.,    0.,    0.,     0.,      0.,      0.,    0.,      0.,
C        o/n2
C         0.
C
C23456789012345678901234567890123456789012345678901234567890123456789012 
C        tn     un     vn      o2       o1      n4s      no       o+  
CCINT =   0.,    0.,    0.,     0.,      0.,      0.,    0.,      0.,
C        n2d    ti     te      ne       o2+      w       z       pot  
C         0.,    0.,    0.,     .2,      0.,      0.,    0.,      0.,
C        n2     ui     vi      wi      fof2     hmf2     rho o/(o2+n2)
C         0.,    0.,    0.,     0.,      0.,      0.,    0.,      0.,
C        o/n2
C         0.
 ENDOFREAD
'EOF'
batchname $job.$$
rcp gmeta $metapath
ja -st