#!/bin/csh
#QSUB -q prem -lt 600 -lT 600 -lm 6Mw -lM 6Mw -eo
ja
batchname mkcdf_mut.$$
cd $TMPDIR
msread mkcdf.x /FOSTER/bin/mkcdf_mut
chmod u+x mkcdf.x
mkcdf.x << 'EOF'
 &input
 ;
 ; histvol: history volumes (mss paths) from which to make new
 ;   netCDF file (may be tigcm, tiegcm, or timesgcm volumes)
 histvol = 
 '/ROBLE/RGR96/TE7EXT1','/ROBLE/RGR96/TE7EXT2',
 '/ROBLE/RGR96/TE7EXT3','/ROBLE/RGR96/TE7EXT4'
 mtimes =
 0,0,0,0,1,0,0,2,0,0,3,0,0,4,0,0,5,0,
 0,6,0,0,7,0,0,8,0,0,9,0,0,10,0,0,11,0,
 0,12,0,0,13,0,0,14,0,0,15,0,0,16,0,0,17,0,
 0,18,0,0,19,0,0,20,0,0,21,0,0,22,0,0,23,0,
 1,0,0
 ;
 ; selected longitudes (floats) (use 999. for zonal means):
 slon = 0.,-70.,180.,999.
 ;
 ; iden: species density conversion option
 ; iden=0 -> no conversion (most are mass mixing ratios)
 ; iden=1 -> convert to number densities (cm3)
 ; iden=2 -> convert to number density mixing ratios
 ; iden=3 -> convert to mass density (gm/cm3)
 iden =        1
 ;
 ; modelhts=0 -> calculate heights from tn, mean mass, and g0
 ; modelhts=1 -> return heights directly from model history
 modelhts =        1
 ;
 ; ionvel = 1 -> get ExB ion drifts
 ; ionvel = 2 -> get ExB+neutral ion drifts
 ; (ionvel ignored if UI,VI,WI not selected)
 ionvel =        1
 ;
 ifields = 
 ;   TN     UN     VN     O2     OX    N4S    NOZ     CO
      1      1      1      1      0      0      0      0
 ;  CO2    H2O     H2    HOX     O+    CH4   O21D    NO2
      0      0      0      0      0      0      0      0
 ;   NO     O3     O1     OH    HO2      H    N2D     TI
      0      0      0      0      0      0      0      0
 ;   TE     NE     O2      W      Z
      0      1      0      0      0
 ;   POTEN  UI     VI     WI      N2
      0      0      0      0      1
 unixdir = "vishnu.hao:/d/foster/tigcmdat/mut"
 mssdir = "/FOSTER/gencdf/mut"
 &end
'EOF'
ja -st
