[3] | 1 | SUBROUTINE XMIE(RAD,RREAL,XIMG,QEXT,QSCT,QABS,CBAR,WAVEN) |
---|
| 2 | C CHANGED TO WORK WITH THE VAX MIE VERSION |
---|
| 3 | COMPLEX ACAP(7000) |
---|
| 4 | c COMPLEX ACAP(700000) |
---|
| 5 | c COMPLEX ACAP(2000000) |
---|
| 6 | REAL ELTRMX(4,1,2),PIE(3,1),TAU(3,1),CSTHT(1),SI2THT(1) |
---|
| 7 | REAL RO,REALO,XIMGO,QE,QS,CB,RD2,RN,CN,TPILAM |
---|
| 8 | c print*,RAD,RREAL,XIMG,QEXT,QSCT,QABS,CBAR,WAVEN |
---|
| 9 | RO=RAD |
---|
| 10 | REALO=RREAL |
---|
| 11 | XIMGO=XIMG |
---|
| 12 | RD2=0. |
---|
| 13 | RN=1.0 |
---|
| 14 | CN=0. |
---|
| 15 | WAVEL=1.E4/WAVEN |
---|
| 16 | TPILAM=2.0*3.14159/WAVEL |
---|
| 17 | CALL DMIESS (RO, REALO, XIMGO, 0.0, 1, |
---|
| 18 | * QE, QS,CB, |
---|
| 19 | * ELTRMX,PIE,TAU,CSTHT,SI2THT,ACAP,1, |
---|
| 20 | c * 2000000, RD2, RN, CN, TPILAM) |
---|
| 21 | * 7000, RD2, RN, CN, TPILAM) |
---|
| 22 | CBAR=CB/QS |
---|
| 23 | QEXT=QE*3.14159*RAD*RAD*1.E-8 |
---|
| 24 | QSCT=QS*3.14159*RAD*RAD*1.E-8 |
---|
| 25 | QABS=QEXT-QSCT |
---|
| 26 | SIZE=3.14159*RAD*RAD*1.E-8 |
---|
| 27 | RETURN |
---|
| 28 | END |
---|