1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
|
PROGRAM FLUXAVERAGE1A
IMPLICIT REAL*8 (A-H,O-Z)
DIMENSION FLAM1(100000),FLAM2(100000),FLAMC(100000),
1 AVERHSURF(100000)
DIMENSION TITLE(74),WLEDGE(200),XMU(20),AVERCONTIN(100000)
DATA NPTI/100000/
C
C FIX THE WAVELENGHT RANGE IN WHICH WE WANT AVERAGED FLUXES
C
READ(5,56)WAVE1,WAVE2
56 FORMAT(3F10.4)
TYPE*,WAVE1,WAVE2
C
READ(5,56)FLAM1(1),FLAM2(1),FLAMC(1)
TYPE*,FLAM1(1),FLAM2(1),FLAMC(1)
C FLAM1(1)=362.0
C FLAM2(1)=362.1
C FLAMC(1)=362.05
DO 550 I=2,NPTI
FLAM1(I)=FLAM1(I-1)+0.100
FLAM2(I)=FLAM2(I-1)+0.100
FLAMC(I)=FLAM1(I)+0.0500
C TYPE*,N,FLAM1(I),FLAM2(I),FLAMC(I)
c if(flam2(i).gt.wave2)go to 501
if(flam2(i).ge.wave2)go to 501
550 CONTINUE
C
501 npti=i
c type*,npti,flam1(npti),flam2(npti),flamc(npti)
do 1001 k=1,npti
if(flam1(k).lt.wave1)go to 1001
nlam1=k
go to 101
1001 continue
101 CONTINUE
c type*,nlam1,wave1,flam1(nlam1)
OPEN(UNIT=1,FORM='UNFORMATTED',READONLY,STATUS='OLD')
READ(1)TEFF,GLOG,TITLE,WBEGIN,RESOLU,NWL,IFSURF,NMU,XMU,NEDGE,
1WLEDGE
2 FORMAT(F10.1,F10.3,/6HTITLE ,74A1/F10.3,F10.1,I10,I5,I5/
1 10F8.4/10F8.4/I10/(5F16.5))
WRITE(2,3)TEFF,GLOG
3 FORMAT(5X,F7.0,9X,F8.5,5X)
WRITE(2,334)TITLE
c334 FORMAT(5HTITLE,74A1)
334 FORMAT(74A1/)
write(2,335)
335 format(4x,'Wavelength',10x,'Flux'/5x,
1 'Angstrom', 2x,'erg/sec/cm2/Angstrom')
c type*,nmu
K=nlam1
NAVER=0
SUMHSURF=0.
SUMCONTIN=0.
NIWL=0
DO 100 I=1,1000000000
READ(1)HSURF,CONTIN
NIWL=NIWL+1
WAVE=WBEGIN*(1.+1./RESOLU)**(NIWL-1)
5 FORMAT(1X,F10.4,2E10.4)
IF(WAVE.LT.FLAM1(K))GO TO 100
IF(WAVE.GE.FLAM1(K).AND.WAVE.LE.FLAM2(K)) THEN
NAVER=NAVER+1
SUMHSURF=SUMHSURF+HSURF
SUMCONTIN=SUMCONTIN+CONTIN
c TYPE*,naver,WAVE,hsurf,sumhsurf
END IF
IF(WAVE.GE.FLAM2(K)-0.001) THEN
AVERHSURF(K)=SUMHSURF/NAVER
AVERCONTIN(K)=SUMCONTIN/NAVER
c IF(K.GT.100000)type*,i,k,flam2(k),averhsurf(k)
K=K+1
NAVER=0
SUMHSURF=0.
SUMCONTIN=0.
END IF
IF(WAVE.GE.wave2-0.001) GO TO 98
100 CONTINUE
98 DO 200 I=1,NPTI
FLAMC(I)=10.*FLAMC(I)
freq=2.997925e18/flamc(i)
FLUXLAM=4.*FREQ/FLAMC(I)*AVERHSURF(I)
WRITE(2,67)FLAMC(I),FLUXLAM
67 FORMAT(3X,F9.1,3X,1PE18.12)
200 CONTINUE
STOP
END
|
