[3] | 1 | SUBROUTINE REGIS (FNU,ND1,NNU,T,ND2,NT,XN2N2,XCH4CH4,XN2CH4,XN2H2) |
---|
| 2 | C THIS SUBROUTINE RETURNS THE PRESSURE INDUCED ABSORBTION COEFFICIENTS |
---|
| 3 | C FOR N2-N2, CH4-CH4, N2-CH4 + CH4-N2, H2-N2 + N2-H2. |
---|
| 4 | C FNU IS THE WAVENUMBER ARRAY |
---|
| 5 | C T IS THE TEMPERATURE ARRAY |
---|
| 6 | C NNU IS THE NUMBER OF WAVENUMBERS |
---|
| 7 | C NT IS THE NUMBER OF TEMPERATURES |
---|
| 8 | C ND1 IS THE DIMENSION OF THE WAVENUMERS |
---|
| 9 | C ND2 IS THE DIMENSION OF THE TEMPERATURES |
---|
| 10 | C XN2N2(ND1,ND2) IS THE N2-N2 COEFFICIENT |
---|
| 11 | C XCH4CH4(ND1,ND2) IS THE CH4-CH4 COEFFICIENT |
---|
| 12 | C XN2CH4(ND1,ND2) IS THE N2-CH4 + CH4-N2 COEFFICIENT |
---|
| 13 | C XN2H2(ND1,ND2) IS THE N2-H2 + H2-N2 COEFFICIENT |
---|
| 14 | C FROM REGIS COURTIN ADAPTED BY CPM. |
---|
| 15 | IMPLICIT REAL (A-H,O-Z) |
---|
| 16 | DIMENSION FNU(NNU),T(NT), |
---|
| 17 | & XN2N2(ND1,ND2),XCH4CH4(ND1,ND2),XN2H2(ND1,ND2),XN2CH4(ND1,ND2) |
---|
| 18 | DIMENSION F8MN(100),F8NH(100),F8HN(100),F10(100), |
---|
| 19 | 2 XN2R(100),XN2T(100),XCH4R(100),XCH4T(100),XH2R(100),XH2T(100) |
---|
| 20 | DATA EPS1/71.4/,EPS2/148.6/,EPS3/36.8/ |
---|
| 21 | DO 1 IT=1,NT |
---|
| 22 | EPS4=SQRT(EPS1*EPS2) |
---|
| 23 | EPS5=SQRT(EPS1*EPS3) |
---|
| 24 | Y=4.*EPS4/T(IT) |
---|
| 25 | F8MN(IT)=FI8(Y) |
---|
| 26 | F10(IT)=FI10(Y) |
---|
| 27 | Y=4.*EPS5/T(IT) |
---|
| 28 | Z=FI8(Y) |
---|
| 29 | F8NH(IT)=Z |
---|
| 30 | F8HN(IT)=Z |
---|
| 31 | Y=4.*EPS1/T(IT) |
---|
| 32 | Z=FI8(Y) |
---|
| 33 | F8MN(IT)=F8MN(IT)/Z |
---|
| 34 | F8HN(IT)=F8HN(IT)/Z |
---|
| 35 | Y=4.*EPS2/T(IT) |
---|
| 36 | F10(IT)=F10(IT)/FI10(Y) |
---|
| 37 | Y=4.*EPS3/T(IT) |
---|
| 38 | F8NH(IT)=F8NH(IT)/FI8(Y) |
---|
| 39 | 1 CONTINUE |
---|
| 40 | CALL PIAN2(0.,NT,XN2R,XN2T,T) |
---|
| 41 | CALL PIACH4(0.,NT,XCH4R,XCH4T,T) |
---|
| 42 | CALL PIAH2(0.,NT,XH2R,XH2T,T) |
---|
| 43 | DO 2 INU=1,NNU |
---|
| 44 | CALL OPAN2(FNU(INU),NT,XN2R,XN2T,T) |
---|
| 45 | CALL OPACH4(FNU(INU),NT,XCH4R,XCH4T,T) |
---|
| 46 | CALL OPAH2(FNU(INU),NT,XH2R,XH2T,T) |
---|
| 47 | c!!!! CALL OPAN2(FNU(INU),NT,XN2R,XN2T) |
---|
| 48 | c!!!! CALL OPACH4(FNU(INU),NT,XCH4R,XCH4T) |
---|
| 49 | c!!!! CALL OPAH2(FNU(INU),NT,XH2R,XH2T) |
---|
| 50 | DO 11 IT=1,NT |
---|
| 51 | XN2N2(INU,IT)=XN2R(IT)+XN2T(IT) |
---|
| 52 | XCH4CH4(INU,IT)=XCH4R(IT)+XCH4T(IT) |
---|
| 53 | C |
---|
| 54 | C LINES 1 AND 3 ARE CORRECTION FACTORS INTRODUCED IN ORDER TO FIT |
---|
| 55 | C THE MEASUREMENTS OF DAGG ET AL (1986) BETWEEN 126 AND 212 K. |
---|
| 56 | C |
---|
| 57 | XN2CH4(INU,IT)=F8MN(IT)*(2.067*XN2R(IT)+2.865*XN2T(IT)) |
---|
| 58 | 1 * 2.48/(T(IT)**0.184) + |
---|
| 59 | 2 F10(IT)*(0.480*XCH4R(IT)+0.374*XCH4T(IT)) |
---|
| 60 | 3 * 1.54 |
---|
| 61 | C |
---|
| 62 | C LINE 3 IS A CORRECTION FACTOR INTRODUCED IN ORDER TO FIT |
---|
| 63 | C THE MEASUREMENTS OF DORE ET AL (1986) BETWEEN 91 AND 298 K. |
---|
| 64 | C OTHER CORRECTION FACTORS ARE IMBEDDED IN THE ROUTINE PIAH2. |
---|
| 65 | C |
---|
| 66 | XN2H2(INU,IT)= F8NH(IT)*(2.543*XH2R(IT)+1.445E1*XH2T(IT)) + |
---|
| 67 | 2 F8HN(IT)*(0.360*XN2R(IT)+0.246*XN2T(IT)) |
---|
| 68 | 3 * 0.30 |
---|
| 69 | 11 CONTINUE |
---|
| 70 | 2 CONTINUE |
---|
| 71 | RETURN |
---|
| 72 | END |
---|