;FUNCTION GEOMF(I,J) FUNCTION geomf,j ; ; NOTE: The single parameter to the IDL function gives the ; satellite number (6..12). ; ; ADDED TO CODE 25 AUG 1992 ; ; I = CHANNELS 1 TO 40 ; J = ('F6', 'F7', 'F8', 'F9', 'F10', 'F11' AND 'F12') ; ; GEOM RETURNS THE GEOMETRIC FACTOR WITH UNITS (CM**2*STER*KEV) ; ; Converted from FORTRAN to IDL by DPS on 950519. ; ; SAVE COMMON geomfact,geof6,geof7,geof8,geof9,geof10,geof11,geof12 ; ; CHARACTER*3 J ; DIMENSION GEOF6(40),GEOF7(40),GEOF8(40),GEOF9(40),GEOF10(40) ; + ,GEOF11(40),GEOF12(40) ; geof6=[75.0,49.0,41.0,33.0,27.0,21.0,16.0,13.0,9.60, $ 7.60,1.40,1.00,0.80,0.51,0.34,0.20,0.10,0.05, $ 0.02,.0063,180.,120.,81.0,55.0,38.0,26.0,18.0, $ 12.0,8.10,5.60,230.,150.,140.,79.0,55.0,37.0, $ 25.0,18.0,12.0,8.00] geof7=[58.0,49.0,41.0,33.0,27.0,21.0,16.0,13.0,9.60, $ 7.60,3.30,2.50,2.00,1.30,0.91,0.58,0.34,0.20, $ 0.11,.057,240.,160.,110.,75.0,52.0,35.0,24.0, $ 17.0,11.0,7.60,210.,140.,98.0,66.0,46.0,31.0, $ 21.0,15.0,10.0,7.00] geof8=[32.6,27.5,23.0,18.5,15.2,11.8,8.99,7.30,5.39, $ 4.27,3.20,2.43,1.94,1.26,.884,.563,.330,.194, $ .107,.0553,115.,76.9,52.9,36.1,25.0,16.8,11.5, $ 8.17,5.29,3.65,111.,74.1,51.9,34.9,24.3,16.4, $ 11.1,7.94,5.29,3.70] geof9=[20.1,17.0,14.2,11.5,9.37,7.29,5.55,4.51,3.33, $ 2.64,3.93,2.98,2.38,1.55,1.08,0.69,.405,.238, $ .131,.0678,114.,76.2,52.4,35.7,24.8,16.7,11.4, $ 8.10,5.24,3.62,130.,86.4,60.5,40.7,28.4,19.1, $ 13.0,9.26,6.17,4.32] geof10=[46.2,38.9,30.2,24.3,19.4,15.8,12.6,10.2,8.24, $ 6.69,6.37,4.32,2.92,1.99,1.37,.930,.637,.438, $ .302,.206,56.7,38.7,26.0,17.5,11.9,8.18,5.54, $ 3.78,2.57,1.77,50.1,33.9,23.0,15.7,10.8,7.31, $ 5.01,3.45,2.38,1.62] geof11=[63.1,54.4,42.8,34.9,27.8,22.6,18.2,14.9,12.1, $ 9.50,10.1,6.85,4.71,3.10,2.09,1.43,.963,.663, $ .449,.300,54.9,38.1,26.0,17.7,12.1,8.24,5.64, $ 3.87,2.66,1.77,50.1,35.0,24.0,16.7,11.4,7.83, $ 5.27,3.71,2.56,1.77] geof12=[41.4,34.7,29.0,24.2,20.3,16.9,14.2,11.8,9.90, $ 8.28,7.30,4.91,3.30,2.22,1.49,1.00,.673,.453, $ .304,.204,70.6,48.2,32.9,22.5,15.3,10.5,7.15, $ 4.88,3.33,2.28,40.2,27.5,18.8,12.9,8.82,6.04, $ 4.13,2.83,1.94,1.32] geofac=REPLICATE(-1.,40) ; IF(J.EQ.' F6') THEN ; GEOMF=GEOF6(I) ; ELSE IF(J.EQ.' F7') THEN ; GEOMF=GEOF7(I) ; ELSE IF(J.EQ.' F8') THEN ; GEOMF=GEOF8(I) ; ELSE IF(J.EQ.' F9') THEN ; GEOMF=GEOF9(I) ; ELSE IF(J.EQ.'F10') THEN ; GEOMF=GEOF10(I) ; ELSE IF(J.EQ.'F11') THEN ; GEOMF=GEOF11(I) ; ELSE IF(J.EQ.'F12') THEN ; GEOMF=GEOF12(I) ; ENDIF CASE j OF 6: geofac=geof6 7: geofac=geof7 8: geofac=geof8 9: geofac=geof9 10: geofac=geof10 11: geofac=geof11 12: geofac=geof12 ELSE: MESSAGE,'Bad satellite: '+STRTRIM(j,2) ENDCASE ; IF geofac GT 0. THEN geofac=geofac/100000. ; RETURN,geofac END