&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 = 'gpi_2002091-2002132.cnrc.nc'
 ; mssfile = 'gpi_1979001-2002059_2002091-2002132.nc'
 ; mssfile = '/TGCM/data/gpi_1979001-2002059.nc'
 ; mssfile = '/TGCM/data/gpi_1979001-2001031.nc'
 ; mssfile = '/FOSTER/tgcminp/79001-99111.gpi.dat'
 ; mssfile = '/FOSTER/tgcminp/79001-97212.gpi.dat'
   mssfile = '/TGCM/data/gpi_1979001-2002181.nc'
 ;
 ; 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 = 1993001,1993365
 ; iyds = 2001001,2001365
 ; iyds = 2001001,2002057
 ;
 ; iyds = 1999324,2000060
 ; iyds = 2000324,2001060
 ; iyds = 2002091,2002132
 ; iyds = 2002001,2002132 ; days 200260-90 are missing
 ; iyds = 2002001,2002181
   iyds = 2002106,2002115
 ;
 ; 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/rdgpi.1999324-2000060.cgm'
 ; sendcgm = 'vishnu.hao:/d/foster/rdgpi/rdgpi.2000324-2001060.cgm'
 ; sendcgm = 'vishnu.hao:/d/foster/rdgpi/rdgpi.2002091-2002132.cgm'
 ; sendcgm = 'vishnu.hao:/d/foster/rdgpi/rdgpi.2002001-2002132.cgm'
 ; sendcgm = 'vishnu.hao:/d/foster/rdgpi/rdgpi.2002001-2002181.cgm'
   sendcgm = 'vishnu.hao:/d/foster/rdgpi/rdgpi.2002106-2002115.cgm'
 ;
 ; sendcgm = 'vishnu.hao:/d/foster/rdgpi/rdgpi.2001001-2002057.cgm'
 ; sendcgm = 'vishnu.hao:/d/foster/rdgpi/rdgpi.yr1993.cgm'
 ; sendcgm = 'vishnu.hao:/d/foster/rdgpi/rdgpi.test.cgm'
 ; sendcgm = 'vishnu.hao:/d/foster/rdgpi/rdgpi.2000280-2000300.cgm'
 ; sendcgm = 'vishnu.hao:/d/foster/tiegcm/rdgpi.98075-80.cgm'
 ; sendcgm = 'vishnu.hao:tgcminp/getndcs/rdgpi/aug97.cgm'
 ; 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 = 0,     0,     0,     0,     0,    0,    0,    0
 ; sendlex = 'vishnu.hao:/d/foster/tiegcm/rdgpi.98075-80.lex'
 ; sendlex = 'vishnu.hao:tgcminp/getndcs/rdgpi/aug97.lex'
 ; 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