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
|
PROGRAM KAPROSSINTERP
DIMENSION AKAP1(5),AKAP2(5),AKAP3(5)
dimension potlo(42),xne(42),rho(42),xnatom(42)
read *,WT1,WT2
PRINT *,WT1,WT2
READ(1,5)
READ(1,5)
READ(2,5)
READ(2,5)
WRITE(3,22)
22 format(' ROSSELAND MASS ABSORPTION COEFFICENTS'/
1' logT logP 0km/s 1km/s 2km/s 4km/s 8km/s log Ne log Na l
2og rho IPloCM-1')
c 3.32-2.00 -5.200 -5.188 -5.169 -5.139 -5.104 5.29314 10.53958-13.12598 0.002
do 8 it=1,57
do 6 iP =1,42
READ(1,5)TABT,TABP,AKAP1
READ(2,5)TABT,TABP,AKAP2
potlo(ip)=1.
xne(ip)=0.
xnatom(ip)=0.
rho(ip)=0.
DO 4 IV=1,5
4 AKAP3(IV)=WT1*AKAP1(IV)+WT2*AKAP2(IV)
WRITE(3,5)TABT,TABP,AKAP3,xne(ip),xnatom(ip),rho(ip),potlo(ip)
5 format(f5.2,f5.2,5F7.3,f9.5,f9.5,f9.5,f8.3)
6 continue
8 continue
CALL EXIT
END
|
