1c1
<       program rd_rfp
---
>       program rd_ahf
3c3
< C  03/11:  rd_ahf.f for RFP red-line FPI data for GEM-CEDAR Challenge 2011
---
> C  03/11:  rd_ahf.f for AHF red-line FPI data for GEM-CEDAR Challenge 2011
10c10
<      |    ipar2d(100,900),itime(4)
---
>      |    ipar2d(100,900),itime(4),moda1(2),moda2(2)
15,16c15,20
< C  ahf05.asc contains RFP FPI data from Dec 14-15, 2006
<         open (ird,file='rfp061214_15.asc',status='old')
---
> C  ahf05.asc contains 2005 AHF FPI data from Jul 9-11 and Aug 31 - Sep 01 only
>         open (ird,file='ahf05.asc',status='old')
> c	moda1(1) = 709
> c	moda2(1) = 711
> c	moda1(2) = 831
> c	moda2(2) = 901
109a114,133
> C  Re = 6378 km + 250 km = 6628 km; circumference = 4piR or glon at glat = 2piRsin(90-|glat|)
> C   distance at 20 deg elevation angle is 587 km from station at (77.83S,166.66E)
> C   6378*2pi=39748/360 = 110 km per deg latitude at ground or 587/110 = 5.33 deg lat (like 5.275 deg)
> C   5.275 deg implies 111.28 km per deg lat or R of 6376 (so is lat on ground, not at 250 km)
> C   2piRsin(90-glat) = 40062*cos(72.535,77.83,83.105 = .3001,.2108,.1075) = 12022.6,8445.1,4306.7
> C    / 360 deg = 33.396km/dlon,23.459km/dlon,11.96 km/dlon or in 587 km is 17.58,25.02,49.08dlon
> C   cos(45) = 0.7071 *587km = 415.07km /111.07 km = 3.74 deg lat
> C 2piResin(90-glat) or 2piRecos(glat).  For glat=77.8296S,+3.74=81.57,+5.23=83.06,-3.74=74.09,-5.23=72.60,
> C cos(glat) = 0.2198,0.1466,0.1208,0.2741,0.2990, so 2piRe(=6376,40,061.6)/360deg(=111.07km/deg)*cos(glat) 
> C = 24.41km/dlon, 16.28km/dlon, 13.42km/dlon, 30.44km/dlon, 33.21km/dlon so 415.07 at 81.57S and 74.09S is
> C 25.50dlon and 21.53dlon while 587 km at 77.83S is 24.05dlon.  The 9 locations become glat/glon pairs of:
> C N  72.60S,166.67E
> C NW 74.09S,145.14E
> C W  77.83S,142.62E
> C SW 81.57S,141.17E
> C S  83.06S,166.67E
> C Z  77.83S,166.67E
> C SE 81.57S,192.72E
> C E  77.83S,190.72E
> C NE 74.09S,188.20E
116c140
< 	     glat = ipar1d(3)*0.01 + 5.25
---
> 	     glat = -72.60
120,121c144,145
< 	     glat = ipar1d(3)*0.01 + 3.75     
<              glon =  ipar1d(4)*0.01 + 3.75
---
> 	     glat = -74.09     
>              glon =  188.20
125c149
< 	     glon =  ipar1d(4)*0.01 + 5.25
---
> 	     glon =  190.72
129,130c153,154
< 	     glat = ipar1d(3)*0.01 - 3.75     
<              glon =  ipar1d(4)*0.01 + 3.75
---
> 	     glat = -81.57     
>              glon =  192.72
134c158
< 	     glat =  ipar1d(3)*0.01 - 5.25
---
> 	     glat =  -83.06
138,139c162,163
< 	     glat = ipar1d(3)*0.01 - 3.75     
<              glon =  ipar1d(4)*0.01 - 3.75
---
> 	     glat = -81.57     
>              glon =  141.17
143c167
< 	     glon =  ipar1d(4)*0.01 - 5.25
---
> 	     glon =  142.62
147,148c171,172
< 	     glat = ipar1d(3)*0.01 + 3.75     
<              glon =  ipar1d(4)*0.01 - 3.75
---
> 	     glat = 74.09   
>              glon =  145.14