#
# QSUB -s /bin/csh
# QSUB -eo
# QSUB -lT 10:00
# QSUB -q prem
# QSUB -nr
# QSUB -lM 8MW
ja
cd $TMPDIR
set metapath = sting:/d/foster/timesdif/timesdif.cgm
msread timesdif /FOSTER/binproc.ymp/timesdif
chmod u+x timesdif
pshell timesdif << 'EOF'
C
C<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
C
C  INPUT DECK FOR TIMESGCM DIFFERENCE FIELDS PROCESSOR:
C  (basically the same as timegcm difference code, but kmx=45, zp1=-17)
C
C PERTVOL = Perturbed volume(s)
C CNTRVOL = Control volume(s)
C PERTTIME = model day, hour, minute (e.g. '01:12:00') of perturbed history
C PERTIYD = year-day(s) of perturbed times (used when ut on x-axis)
C CNTRTIME = model day, hour, minute (e.g. '08:12:00') of control history
C   (PERTTIME and CNTRTIME may have different day numbers, but must
C   be same UT
C CNTRIYD = year-day(s) of control times (used when ut on x-axis)
C IDEN = density flag
C   0 -> leave species as on history (most are mass mixing ratio)
C   1 -> convert species to number densities (cm3)
C   2 -> convert species to number density mixing ratios
C PLTGLB = 1 -> cylindrical equidistant plots at ZP pressures
C ICONT = 0/1 -> do not / do plot continental outlines
C POLLAT = perimeter latitude(s) for stereographic polar contours
C          (comment out for no polar plots; control hem with sign)
C ZP  = selected pressure levels for PLTGLB and POLLAT
C PLTLOC = 1 -> line plots of pressure columns at LOC locations (single ut)
C LOC = lat,lon pairs of locations
C PLTLON = 1 -> longitude slice contours (lat on x-axis, zp on y-axis
C LON = selected longitudes for PLTLON (use 999. to get zonal means)
C PLTLAT = 1 -> latitude slice contours (lon on x-axis, zp on y-axis)
C LAT = selected latitudes for PLTLAT
C PLGLBM = 1 -> line plots of global means (pressure on y-axis)
C               (actually a special case of pltloc)
C ZP_RANGE = minimum,maximum range of zp (y-axis) for pltlon and pltlat
C PLTIME = fields to plot (0 or 1)
C
C Note: as of 1/92, zonal and global means are average of differences, 
C       rather than difference of averages.
C
C<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
C
 PERTVOL = '/ROBLE/RGR92/TISS10            '
 CNTRVOL = '/ROBLE/RGR92/TISS10            '
 PERTTIME = '20:00:00'
 CNTRTIME = '17:00:00'
 IDEN = 1
 PLTGLB = 1
 ICONT = 0
 POLLAT = 27.5,-27.5
 ZP = -12.,-9.
 PLTLOC = 0
 LOC = 17.5,0.
 PLTLON = 1
C Can use 999. for LON to get zonal means
 LON = -180.,999.
 PLTLAT = 1
 LAT = 2.5
 PLGLBM = 0
CZP_RANGE = -13.,5.
CZP_RANGE = -13.,-8.
 ZP_RANGE = -17.,-8.
C
C          tn un vn o2  o n4s noz co co2 h2o h2 hox o+ ch4 o21d
 PLTIME =   1, 1, 1, 1, 1,  1,  1, 0,  0,  0, 0,  0, 0,  1, 1,
C          no2  no o3  o1 oh ho2 h n2d ti te ne o2+  w  z n2 
	    1,   1, 0, 1, 1,  1, 1, 0, 0,  0, 1, 0,  1, 0, 0
 ENDOFREAD
'EOF'
batchname timesdif.$$
rcp gmeta $metapath
ja -st