#
# QSUB -s /bin/csh
# QSUB -eo
# QSUB -lT 3:00
# QSUB -q prem
# QSUB -nr
# QSUB -lM 2MW
#
# Execute rdgpi to obtain and plot geophysical indices:
#
ja
cd $TMPDIR
setenv NCPUS 1
set job = rdgpi
batchname $job.$$
msread $job /FOSTER/bin/$job
chmod u+x $job
$job << 'EOF' || goto badx
 &input
E
 ;
 ; Read geophysical indices from database mssfile and optionally plot 
 ;   and/or make fortran data statements.
 ;
 ; 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-94080.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) )
 ; (note for dates, mo,id,yr must be 2-chars each)
 ;
 dates = '09/01/93','10/31/93'
 ;
 ;iyds = 85263,85263
 ;
 ; 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  3hrCtpoten aveCtpoten 3hrHP aveHP
 ifplt = 1,     1,     1,     1,     1,         1,         1,    1
 sendcgm = 'vishnu.hao:/d/foster/tgcminp/getndcs/rdgpi/rdgpi.cgm'
 ;
 ; ifdat = flags to make data statements of each field:
 ; senddat = machine:path to send file containing data statements
 ;
 ;       F107D  F107A  3hrKp  aveKp  3hrCtpoten aveCtpoten 3hrHP aveHP
 iflex = 1,     1,     0,     0,     1,         0,         1,    0
 sendlex = 'vishnu.hao:/d/foster/tgcminp/getndcs/rdgpi/rdgpi.lex'
 ;
 &end
'EOF'
ja -st
exit
badx:
echo 'Execution failed'
rcp $job.cgm vishnu.hao:/d/foster/tgcminp/getndcs/$job/$job.fail.cgm
exit