#
# QSUB -s /bin/csh
# QSUB -lT  2500
# QSUB -q prem
# QSUB -nr
# QSUB -lM 4MW
# QSUB -o flxproc.$$
# QSUB -eo
#
# 8/98: this must be run on paiute, as that is the only machine
#       that still has the ccm processor.
#
cd $TMPDIR
#
# RUNCCM2PROC = 1 -> execute ccm2 processor with below inputs to make
#   a "List Sorted Data" (lsd) save tape for use by flxproc below.
# (not necessary to run ccm2proc *if* ccmlsd (flxproc input) mss path 
#   pre-exists)
# (comment out or set to 0 to not run ccm2proc)
#
set RUNCCM2PROC = 1
set job = flxproc
#batchname $job.$$
#
# Run ccm2 processor and save lsd file:
#
if (! $?RUNCCM2PROC || $RUNCCM2PROC == 0) goto runflxproc
cat >! parms << 'END'
C
 TAPESA = '/LIUH/csm/cpl22mt/ccm3/hist/h0011'
 TITLEA = 'tgcm22mt/ccm3.6'
 DAYSA =  21.
C
C Get tgcm/t21 coupled fields from ccm2:
C
  FIELDA1 = 'T','U','V','Q','Z2','OMEGA'
C  FIELDA1 = 'V'
C
C 10 values of PRESSLE for T21: 
C (mb pressures corresponding to zp -20.5 to -16.0, every 0.5)
C
 PRESSLE =
  4.3292312622E+02, 2.6258111572E+02, 1.5926348877E+02, 
  9.6598197937E+01, 5.8589767456E+01, 3.5536491394E+01, 
  2.1553972244E+01, 1.3073143959E+01, 7.9292626381E+00, 
  4.8093414307E+00
C
C Run ccm2 proc in LSD mode:
C
 DRVRTYP = 'LSD1'
C
C RTLS1A = 2 -> ordering of horiz slices = field level, time, ensemble member
C
 RTLS1A = 2
C
C PKLS1A = 1 -> no packing
C
 PKLS1A = 1
C
C SAVLS1A = mss path at which to store lsd save tape
C (should be same as ccmlsd in processor inputs below)
C
 SAVLS1A = '/LIUH/csm/cpl22mt/ccm3/lsd/h0011'
ENDOFDATA
'END'
/ccm/proc/Processor || goto ccmproc_fail
echo ' '
echo '-------------------------------------------------------------------'
echo ' '
echo Completed ccm2 processor run at `date`
runflxproc:
ja
#
# Get flxproc executable from mss:
#
#msread $job /FOSTER/bin/$job
#chmod u+x $job
#
# Run flx processor:
#
cat >! $job.inp << 'EOF'
 &input_sut
E
 ;
 ;----------------------------------------------------------------------
 ;
 ; tgcm/t21:
 ;
 ; histvols = time-gcm history (-16 to +5 coupled run)
 ;
   histvols = '/LIUH/cpl22mt/prim14', mtimes= 21,0,0
 ;
 ; ccmlsd = mss path of ccm2 lsd save tape written by ccm2 processor
 ; (should be same as SAVLS1A above if ccm2 processor was executed)
 ;
   ccmlsd = '/LIUH/csm/cpl22mt/ccm3/lsd/h0011'
 ;
 ;----------------------------------------------------------------------
 ;
 ; Request fields: cfields is string array with the below valid values
 ;   (spelling is important)
 ; Coupled fields are bracketed, e.g., [TN] (CH4 coupled only with maccm)
 ; Quotes around cfield field names are not necessary. 
 ; Do not include brackets or parens in cfields.
 ;
 ; [TN] [UN] [VN] O2   OX   N4S  NOZ  CO   CO2  [H2O] H2   HOX O+ (CH4) 
 ;  AR   HE  O21D NO2  NO   O3   O1   OH   HO2  H    N2D  TI   TE   NE   
 ;  O2+ [W]  [Z] POTEN UI   VI   WI   N2   RHO  UNVN UIVI FOF2 HMF2
 ;
 ; cfields = TN,UN,VN,H2O,W,Z
  cfields = VN
 ; 
 ;fmnmxint = 'TN',1.,0.,5.
 ;
 ; iden = 0 -> mass mix ratios (for coupled H2O)
 ;
 iden = 0
 ionvel = 0
 modelhts = 1
 ;
 ;----------------------------------------------------------------------
 ;
 icolor = 0             ; make color fill contours if icolor=1
 ibox_clabs = 0		; box contour line labels
 ;outplt = 'cgm','ps'	; plot output type(s) ('cgm' and/or 'ps')
 ;outplt = 'ps'		; plot output type(s) ('cgm' and/or 'ps')
 ;outplt = 'cgm'		; plot output type(s) ('cgm' and/or 'ps')
 outplt = 'cgm'
 ;psmode = 'port'	; for ps: 'port' for portrait or 'land' for landscape
 psmode = 'land'	; for ps: 'port' for portrait or 'land' for landscape
 ;
 ; multiplt = 1 -> put multiple plots on each frame
 ; ipltrowcol = nrows,ncols of plots on each frame (multiplt > 0 only)
 ; multiadvfr = 0 -> advance frame *only* when page is full according to
 ;                   ipltrowcol.
 ; multiadvfr = 1 -> advance frame with full page, *and* between plot
 ;                   types (maps, slices, xy plots) and ut's
 ;
 multiplt = 0           ; multiple plots per frame (0 or 1 toggle)
   ipltrowcol = 2,2	; # rows,cols of plots per frame (multiplt > 0)
   multiadvfr = 1	; frame adv flag for multiplt
 ;
 ; iyaxright: single integer flag for extra right-hand y-axes. 
 ;   (applies to ipltlat and ipltlon vertical slices only)
 ; For now, assume left hand y-axis is always linear Zp, then:
 ; iyaxright = 0 -> add no extra right-hand y-axes 
 ; iyaxright = 1 -> add extra (non-linear) right-hand y-axis in height (km)
 ; iyaxright = 2 -> add extra right-hand y-axes in height *and* pressure (mb)
 ;
 iyaxright = 2
 ;
 ;----------------------------------------------------------------------
 ; Maps:
 ; (contour over map projections if ipltmaps=1)
 ; (if ipltmaps=0, remaining map parameters are ignored)
 ; (as of 11/3/95, interpolation to height surface not yet available --
 ;  if fmap_zpht > 5., it will be reset to 5.)
 ;
 ipltmaps = 1
 ;   fmap_zpht = -20.,-18.,-16.,-14.,-12.,-10.,-8.  ; 1 int: make map projection
   fmap_zpht = -20.,-19.5,-19.,-18.5,-18.,-17.5,-17.,-16.5,-16.,-15.5,-15.,-14.5,-14.
   map_continents = 0                           ; 1 int: continent outlines
   map_tn_unvn = 1                              ; 1 int: add unvn to tn
   map_ht_unvn = 1                              ; 1 int: add unvn to ht (z)
   map_ep_uivi = 0                              ; 1 int: add uivi to epot
   ivec_label = 0                               ; 1 int: label vec mags
 ; map_top_anno = 'ANNOTATION FOR MAPS'
 ;
 ; Map projections and related options:
 ;
   map_global = 1                               ; 1 int: CE projection
     map_global_cenlon = 0                      ; 1 int: center longitude
 ;   map_global_censlt = 12                     ; 1 int: center local time
   map_polar = 0                                ; 1 int: ST projection
     fmap_polar_perimlat = -30., 30.            ; n floats: perimeter lat(s)
   map_satview = 0                              ; 1 int: SATV projection
     fmap_satview_latlon = -55.,0.              ; 2 floats: center of satv
 ;   fmap_satview_eradii = 1.5                  ; 1 float: earth radii dist 
 ;
 ;----------------------------------------------------------------------
 ;
 ; Latitude slices (contour longitude vs zp or ht at selected latitudes)
 ; (more than one pair of flat_zprange is allowed for vertical scales)
 ; (if flat_zprange(1) is < bottom of model(s), or flat_zprange(2) is above
 ;  model(s), it will be adjusted to match bottom and top of model(s),
 ;  e.g., flat_zprange = -99.,99. will be adjusted to full vertical range)
 ; flat_htscale = linear height scale(s) for interpolation on y-axis
 ;   (triplets of bot,top,delta heights for y-axis)
 ;   (more than one triplet for flat_htscale is allowed for height scales)
 ;
 ipltlat = 0
   flats = -72.5,-67.5,-42.5,-17.5,2.5,17.5,   ; n floats: selected latitudes
	    42.5,67.5,72.5
 ;  flat_zprange = -99.,5., -99.,-8.
   flat_htscale = 5.,500.,10., 5.,100.,1., 5.,50.,.5
   ilat_log10 = 1                        ; 1 int: plot log10 of species
 ; lat_top_anno = 'TOP ANNOTATION FOR LAT SLICES'
 ;
 ;----------------------------------------------------------------------
 ;
 ; Longitude slices (contour latitude vs zp or ht at selected longitudes)
 ; (more than one pair of flon_zprange is allowed for vertical scales)
 ; (if flon_zprange(1) is < bottom of model(s), or flon_zprange(2) is above
 ;  model(s), it will be adjusted to match bottom and top of model(s),
 ;  e.g., flon_zprange = -99.,99. will be adjusted to full vertical range)
 ; flon_htscale = linear height scale(s) for interpolation on y-axis
 ;   (triplets of bot,top,delta heights for y-axis)
 ;   (more than one triplet for flon_htscale is allowed for height scales)
 ; 
 ; (flons = selected longitudes. Use flons = 'zm' for zonal means,
 ;  e.g. flons = -70.,'zm' will make lon slices at 70W, and zonal means)
 ;                                         1/flxproc.day168
 ipltlon = 0
   flons = -180.,0.,'zm'       ; n floats: selected longitudes,'zm'
 ; flon_zprange = -99.,5., -99.,-8.
   flon_htscale = 5.,500.,10., 5.,100.,1., 5.,50.,.5
   ilon_log10 = 1                        ; 1 int: plot log10 of species
 ; lon_top_anno = 'TOP ANNOTATION FOR LON SLICES'
 ;
 ;----------------------------------------------------------------------
 ;
 ; ipltutvert = 1 -> contour ut vs zp/ht at utvert_locs
 ;                  ('ZM' not allowed, but 'LTxx.x' is allowed)
 ipltutvert = 0
 ;
 ; utvert_locs = lat,lon pairs of locations (lon can be 'LTxx.x')
 ; E.g.: utvert_locs = 0.,'LT12.' -> lat=0., local noon
 ;       utvert_locs = 60.,-70.   -> lat=60.,lon = -70.
 ;
 ;utvert_locs =     2.0,-157.,  20.8,-156.5,  41.9,-111.4,  74.4,-94.5,
 ;                  40.,-105.5, -44.,170.4, -36.,139.
 ;utvert_locname = 'CHRISTMAS ISL', 'MAUI',  'BEAR LAKE',  'RESOLUTE BAY',
 ;                 'FORT COLLINS','MTJOHNS','ADELAIDE'
 ;utvert_locs =     2.0,-157.,  20.8,-156.5,  41.9,-156.5,  41.9,-111.4
 ;utvert_locname = 'CHRISTMAS', 'MAUI',       'UWOMI',  'BEAR LAKE'
 utvert_locs =     2.0,-157.
 utvert_locname = 'CHRISTMAS'
 ;utvert_zprange = -99.,99.,-99.,-8.		; bot,top for y-axis zp
 utvert_zprange = -99.,-8.
 ;utvert_htscale = 30.,100.,1.		; bot,top,delta for y-axis ht
 iutvert_log10 = 1                      ; log10 of species
 ;
 ;----------------------------------------------------------------------
 ; Send Files:
 ;
 ;sendxdr = 'worsel.hao:/local/e/liuh/tgcmproc/f90/xdr/flxproc.021.V.xdr'
 ;sendcgm = 'worsel.hao:/local/e/liuh/tgcmproc/f90/flxproc.d013.cgm'
 ;sendps = 'ouray.hao:/d/roble/flxcplr/t21/flxproc.ps'
 ;
 sendcgm = 'vishnu.hao:/d/foster/tgcmvis/flxproc/flxproc.cgm'
 ;sendps = 'vishnu.hao:/d/foster/tgcmvis/flxproc/flxproc.ps'
 sendxdr = 'vishnu.hao:/d/foster/tgcmvis/flxproc/xdr/flxproc.xdr'
 ;
 ;----------------------------------------------------------------------
 &end_sut
'EOF'
echo Executing $job..
#$job < $job.inp || goto xfail
flxproc_j90 < $job.inp || goto xfail
ja -st
exit
#
ccmproc_fail:
echo \>\>\> ccm2 processor execution failed
exit
#
mkproc_fail:
echo \>\>\> make of flx processor failed
exit
#
xfail:
echo \>\>\> Execution of flxproc failed
exit