diff options
| author | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-02-16 12:40:45 -0500 |
|---|---|---|
| committer | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-02-16 12:40:45 -0500 |
| commit | 0373ffdfaaaa3845306ca71243d535fdffd941d4 (patch) | |
| tree | 194c3c278d7e352e39d555d31aae93c0be2dfc03 /kappa9/kaprossinterp.for | |
| parent | 01b51f73bd06b2d6eabb776ba6cc69e4abfaa0b3 (diff) | |
| download | kasym-0373ffdfaaaa3845306ca71243d535fdffd941d4.tar.gz | |
Initial commit
Diffstat (limited to 'kappa9/kaprossinterp.for')
| -rw-r--r-- | kappa9/kaprossinterp.for | 30 |
1 files changed, 30 insertions, 0 deletions
diff --git a/kappa9/kaprossinterp.for b/kappa9/kaprossinterp.for new file mode 100644 index 0000000..3bbc161 --- /dev/null +++ b/kappa9/kaprossinterp.for @@ -0,0 +1,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 |