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