#
# QSUB -s /bin/csh
# QSUB -eo
# QSUB -lT 9:00
# QSUB -q prem
# QSUB -nr
# QSUB -lM 2MW
#
# Execute rdgpi to obtain and plot geophysical indices:
#
ja
cd $TMPDIR
set job = rdgpi
batchname $job.$$
msread $job /FOSTER/bin/$job
chmod u+x $job
$job << 'EOF' || goto badx
 &input
 ;
 ; Read geophysical indices from database mssfile and optionally plot 
 ;   and/or make lexical read statements for model input.
 ;
 ; f107a, f107d, and 3-hr Kp are obtained directly from the database file.
 ; Hemispheric Power is calculated from Kp as follows:
 ;   hpower = amax1(0., -2.78+9.33*fkp)
 ; Cross-cap potential is calculated from Kp as follows:
 ;   ctpoten = 29.+11.*fkp
 ;
 mssfile = '/FOSTER/tgcminp/79001-97212.gpi.dat'
 ; mssfile = '/FOSTER/tgcminp/79001-97212.gpi.dat'
 ;
 ; Provide beginning and end of time period of desired indices:
 ; Can give either dates (two strings, each mo/id/yr), *or*
 ;  year-days (two 5-digit integers yyddd) )
 ; (for dates, mo,id,yr must be 2-chars each)
 ; (min and max dates are given in the mssfile name)
 ;
 ;dates = '11/ 1/94','11/10/94'
 ;dates = '11/ 1/94','11/10/94'
 ;dates = ' 8/1/98',' 8/15/98'
   iyds = 94305,94314
 ; iyds = 97214,97212
 ;
 ; ifplt = flags to make line plots of each field:
 ; sendcgm = machine:path to send plots
 ; (fields with "ave" prefix are daily averaged from the 3-hr values)
 ;
 ;       F107D  F107A  3hrKp  aveKp  3hrCP aveCP 3hrHP aveHP
 ifplt = 1,     1,     1,     1,     1,    1,    1,    1
   sendcgm = 'vishnu.hao:tgcminp/getndcs/rdgpi/mahrsi.94305-94314.cgm'
 ; sendcgm = 'ouray.hao:timegcm/tgcm22/mahrsi.94305-94314.cgm'
 ; sendcgm = 'ouray.hao:timegcm/tgcm22/mahrsi.110194-111094.cgm'
 ; sendcgm = 'ouray.hao:timegcm/tgcm22/mahrsi.080197-081598.cgm'
 ; sendcgm = 'ouray.hao:timegcm/kibo/strtwarm.93.indicies'
 ; sendcgm = 'ouray.hao:timegcm/kibo/sissi3.aug.90.cgm'
 ;
 ; ifdat = flags to make data statements of each field:
 ; senddat = machine:path to send file containing data statements
 ;
 ;       F107D  F107A  3hrKp  aveKp  3hrCP aveCP 3hrHP aveHP
 iflex = 1,     1,     0,     0,     1,    0,    1,    0
   sendlex = 'vishnu.hao:tgcminp/getndcs/rdgpi/mahrsi1.lex'
 ; sendlex = 'ouray.hao:timegcm/tgcm22/mahrsi1.lex'
 ; sendlex = 'ouray.hao:/d/roble/timegcm/mahrsi.lex'
 ; sendlex = 'ouray.hao:tgcm/tgcm12/may.96.lex'
 ; sendlex = 'ouray.hao:timegcm/kibo/sissi3.aug.90.lex'
 ;
 &end
'EOF'
ja -st
exit
badx:
echo 'Execution failed'
#rcp gmeta ouray.hao:/d/roble/timegcm/$job.fail.cgm
#rcp rdgpi.lex ouray.hao:/d/roble/timegcm/$job.fail.lex
exit