&input
 ;
 ; Read geophysical indices from mssfile and optionally plot and/or make
 ;   fortran data statements.
 ;
 ; f107a, f107d, and 3-hr Kp are obtained directly from the database.
 ; 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'
 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) )
 ; (note for dates, mo,id,yr must be 2-chars each)
 ;
 ;dates = '09/01/93','09/30/93'
 dates = '11/01/94','11/10/94'
 ;
 ;iyds = 85263,85263
 ;iyds = 86035,86055
 ;iyds = 92001,94120
 ;iyds = 93001,93031
 ;iyds = 92335,93061
 ;
 ; 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'
 ;
 ; iflex = 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,     1,     1,     1,         1,         1,    1
 sendlex = 'vishnu.hao:/d/foster/tgcminp/getndcs/rdgpi/rdgpi.dat'
 ;
 &end