#
# QSUB -s /bin/csh
# QSUB -eo
# QSUB -lT 5:00
# QSUB -q ded_1
# QSUB -nr
# QSUB -lM 2MW
#
# Compile and execute rdgpi to obtain and plot geophysical indices.
# Submit to an SGI machine (e.g., dataproc)
# (source is at hao: /home/tgcm/rdgpi)
#
cd $TMPDIR
set job = rdgpi
echo Begin $job job at `date`
#
# Get source files, Makefile, and mkdep script:
#
INCLUDE -h /home/tgcm/rdgpi/*.f /home/tgcm/rdgpi/*.F 
INCLUDE -h /home/tgcm/rdgpi/Makefile
INCLUDE -h /home/tgcm/bld/mkdep
#
# Build executable (exec must be same as EXEC in Makefile)
#
perl mkdep || echo "WARNING: error from mkdep."
gmake || echo "WARNING: gmake failed."
set exec = rdgpi.irix
if (! -e $exec) then
  echo "WARNING: cannot find exec file $exec (ls follows):"
  /usr/bin/ls -l
endif
#
$exec << '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 = '/TGCM/data/gpi_1979001-2002059.nc' ; new file made 3/20/02
 ;mssfile = '/TGCM/data/gpi_1979001-2001031.nc'
  mssfile = '/TGCM/data/gpi_1979001-2002090_2002091-2002132.nc' ; with cnrc/isgi data
 ;
 ; 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 = '08/01/00','08/15/00'
 ; iyds = 1999360,2000005
 ; iyds = 2000280,2000300
 ; iyds = 2000101,2001365
 ;
 ; iyds = 1999324,2000060
 ; iyds = 2000324,2001060
   iyds = 2002075,2002132
 ;
 ; 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:/d/foster/rdgpi/gpi.1999324-2000060.cgm'
 ; sendcgm = 'vishnu.hao:/d/foster/rdgpi/gpi.2000324-2001060.cgm'
 ; sendcgm = 'vishnu.hao:/d/foster/rdgpi/rdgpi.yr2001.cgm'
 ; sendcgm = 'vishnu.hao:/d/foster/rdgpi/rdgpi.test.cgm'
   sendcgm = 'vishnu.hao:/d/foster/rdgpi/gpi.2002075-2002132.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 = 0,     0,     0,     0,     0,    0,    0,    0
 ; sendlex = 'vishnu.hao:/d/foster/tiegcm/rdgpi.98075-80.lex'
 ;
 &end
'EOF'
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