| 1 | SUBROUTINE GAS2_NOHCN(J,K,T,P,U,TAUGAS) |
|---|
| 2 | #include "dimensions.h" |
|---|
| 3 | PARAMETER(NLAYER=llm,NLEVEL=NLAYER+1) |
|---|
| 4 | COMMON /STRATO/ C2H2(NLAYER),C2H6(NLAYER) |
|---|
| 5 | DIMENSION A(18),A1(18),A2(18),B(18),B1(18),B2(18) |
|---|
| 6 | DIMENSION C(18),C1(18),C2(18) |
|---|
| 7 | DATA P0/0.1/ |
|---|
| 8 | C P0 IS THE REFERENCE PRESSURE USED IN FORMULA |
|---|
| 9 | C PRESSURE IS IN BARS!!! ??FLAG? |
|---|
| 10 | C U IS IN MOLECULES PER CM2 |
|---|
| 11 | DATA a / |
|---|
| 12 | &-58.07850,-73.66951,-22.00411,-26.61158,-27.04048,-25.84024, |
|---|
| 13 | & 0.00000,-18.70291,-16.08370,-23.03309,-31.18288,-19.96538, |
|---|
| 14 | &-24.37518,-17.94444, 0.00000,-17.46335,-20.93769,-10.74428/ |
|---|
| 15 | DATA a1/ |
|---|
| 16 | & 31.24782, 45.09350, 3.53728, 7.34218, 7.09808, 5.35252, |
|---|
| 17 | & 0.00000, 1.58486, -2.62994, 2.88087, 11.38614, 0.33491, |
|---|
| 18 | & 3.36078, -3.94288, 0.00000, -1.57041, 0.28770,-10.54978/ |
|---|
| 19 | DATA a2/ |
|---|
| 20 | & -6.15763, -9.47812, -0.72611, -1.73805, -1.65804, -1.23760, |
|---|
| 21 | & 0.00000, -0.33713, 0.84310, -0.44196, -2.79135, 0.00277, |
|---|
| 22 | & -0.60197, 1.27842, 0.00000, 0.37686, 0.08828, 2.77701/ |
|---|
| 23 | DATA b/ |
|---|
| 24 | & 5.27820, 2.66675, 1.31741, -5.05287, 1.88008, 4.49298, |
|---|
| 25 | & 0.00000, -0.94667, 2.77477, 15.80634, 33.83065, -2.29512, |
|---|
| 26 | & 2.90790, 0.77265, 0.00000, -0.95471, 3.30490, 0.43886/ |
|---|
| 27 | DATA b1/ |
|---|
| 28 | & -4.00372, -1.49356, -0.80081, 5.86571, -0.95994, -3.38754, |
|---|
| 29 | & 0.00000, 1.24800, -1.75869,-13.07462,-29.97486, 2.67427, |
|---|
| 30 | & -1.73763, 0.87904, 0.00000, 1.35597, -2.09025, 1.40393/ |
|---|
| 31 | DATA b2/ |
|---|
| 32 | & 0.81995, 0.32331, 0.16170, -1.42280, 0.26043, 0.87377, |
|---|
| 33 | & 0.00000, -0.31593, 0.36935, 2.86686, 6.83669, -0.75690, |
|---|
| 34 | & 0.24220, -0.46079, 0.00000, -0.42681, 0.32905, -0.61443/ |
|---|
| 35 | DATA c/ |
|---|
| 36 | & -8.99733, 6.51736, -0.54907, -4.80457, 2.83010, -2.12019, |
|---|
| 37 | & 0.00000, -1.53211,-13.48286,-18.91273,-16.89597, -1.02362, |
|---|
| 38 | & 0.25710, -6.27612, 0.00000, -0.56746, 0.14651, -8.53711/ |
|---|
| 39 | DATA c1/ |
|---|
| 40 | & 8.43865, -6.62977, 0.48722, 4.43910, -3.21534, 1.92581, |
|---|
| 41 | & 0.00000, 1.33927, 12.54849, 17.74174, 14.03467, 0.84636, |
|---|
| 42 | & -0.30260, 5.98928, 0.00000, 0.40909, -0.17410, 8.13665/ |
|---|
| 43 | DATA c2/ |
|---|
| 44 | & -1.99052, 1.64435, -0.13150, -1.02683, 0.80781, -0.49989, |
|---|
| 45 | & 0.00000, -0.30569, -2.91842, -4.14473, -2.91521, -0.18050, |
|---|
| 46 | & 0.06981, -1.43493, 0.00000, -0.07578, 0.03729, -1.93800/ |
|---|
| 47 | XKV = A(K) + A1(K)*LOG10(T) + A2(K)*LOG10(T)**2 |
|---|
| 48 | &+ LOG10(P/P0)*(B(K) + B1(K)*LOG10(T) + B2(K)*LOG10(T)**2) |
|---|
| 49 | &+ (C(K) + C1(K)*LOG10(T) + C2(K)*LOG10(T)**2) |
|---|
| 50 | & *LOG10(P/P0)**2 |
|---|
| 51 | GASMIX=C2H2(J) |
|---|
| 52 | IF (K .GT. 11) GASMIX=C2H6(J) |
|---|
| 53 | XKV=10.**XKV |
|---|
| 54 | C CLIP OUT THE CLEAR INTERVALS. |
|---|
| 55 | IF (A(K) .EQ. 0.0) XKV=0.0 |
|---|
| 56 | TAUGAS=U*GASMIX*XKV |
|---|
| 57 | RETURN |
|---|
| 58 | END |
|---|
