source: trunk/libf/phytitan/gas2.F @ 104

Last change on this file since 104 was 3, checked in by slebonnois, 15 years ago

Creation de repertoires:

  • chantiers : pour communiquer sur nos projets de modifs
  • documentation : pour stocker les docs

Ajout de:

  • libf/phytitan : physique de Titan
  • libf/chimtitan: chimie de Titan
  • libf/phyvenus : physique de Venus
File size: 4.9 KB
RevLine 
[3]1      SUBROUTINE GAS2(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      COMMON /STRAT2/ HCN(NLAYER)
6      DIMENSION A(18),A1(18),A2(18),B(18),B1(18),B2(18)
7      DIMENSION C(18),C1(18),C2(18)
8      DIMENSION AHCN(11),A1HCN(11),A2HCN(11),BHCN(11),B1HCN(11)
9      DIMENSION B2HCN(11),CHCN(11),C1HCN(11),C2HCN(11)
10      DATA P0/0.1/
11
12C P0 IS THE REFERENCE PRESSURE USED IN FORMULA
13C PRESSURE IS IN BARS!!! ??FLAG?
14C U IS IN MOLECULES PER CM2
15
16      DATA a /
17     &-58.07850,-73.66951,-22.00411,-26.61158,-27.04048,-25.84024,
18     &  0.00000,-18.70291,-16.08370,-23.03309,-31.18288,-19.96538,
19     &-24.37518,-17.94444,  0.00000,-17.46335,-20.93769,-10.74428/
20      DATA a1/
21     & 31.24782, 45.09350,  3.53728,  7.34218,  7.09808,  5.35252,
22     &  0.00000,  1.58486, -2.62994,  2.88087, 11.38614,  0.33491,
23     &  3.36078, -3.94288,  0.00000, -1.57041,  0.28770,-10.54978/
24      DATA a2/
25     & -6.15763, -9.47812, -0.72611, -1.73805, -1.65804, -1.23760,
26     &  0.00000, -0.33713,  0.84310, -0.44196, -2.79135,  0.00277,
27     & -0.60197,  1.27842,  0.00000,  0.37686,  0.08828,  2.77701/
28      DATA b/
29     &  5.27820,  2.66675,  1.31741, -5.05287,  1.88008,  4.49298,
30     &  0.00000, -0.94667,  2.77477, 15.80634, 33.83065, -2.29512,
31     &  2.90790,  0.77265,  0.00000, -0.95471,  3.30490,  0.43886/
32      DATA b1/
33     & -4.00372, -1.49356, -0.80081,  5.86571, -0.95994, -3.38754,
34     &  0.00000,  1.24800, -1.75869,-13.07462,-29.97486,  2.67427,
35     & -1.73763,  0.87904,  0.00000,  1.35597, -2.09025,  1.40393/
36      DATA b2/
37     &  0.81995,  0.32331,  0.16170, -1.42280,  0.26043,  0.87377,
38     &  0.00000, -0.31593,  0.36935,  2.86686,  6.83669, -0.75690,
39     &  0.24220, -0.46079,  0.00000, -0.42681,  0.32905, -0.61443/
40      DATA c/
41     & -8.99733,  6.51736, -0.54907, -4.80457,  2.83010, -2.12019,
42     &  0.00000, -1.53211,-13.48286,-18.91273,-16.89597, -1.02362,
43     &  0.25710, -6.27612,  0.00000, -0.56746,  0.14651, -8.53711/
44      DATA c1/
45     &  8.43865, -6.62977,  0.48722,  4.43910, -3.21534,  1.92581,
46     &  0.00000,  1.33927, 12.54849, 17.74174, 14.03467,  0.84636,
47     & -0.30260,  5.98928,  0.00000,  0.40909, -0.17410,  8.13665/
48      DATA c2/
49     & -1.99052,  1.64435, -0.13150, -1.02683,  0.80781, -0.49989,
50     &  0.00000, -0.30569, -2.91842, -4.14473, -2.91521, -0.18050,
51     &  0.06981, -1.43493,  0.00000, -0.07578,  0.03729, -1.93800/
52
53      DATA ahcn/
54     & -33.41732,-35.36144,-19.76105,-16.45201,-13.09335,-18.84868,
55     &   0.00000,-18.605,  -13.479,  -20.717,  -17.126/
56      DATA a1hcn/
57     &  11.42543, 12.36695,  2.43240, -2.78746, -6.36678, -0.49211,
58     &   0.00000,  1.0009,  -4.9550,   2.0424,  -2.0131/
59      DATA a2hcn/
60     &  -2.20538, -2.46746, -0.66486,  0.78681,  1.54548, -0.05627,
61     &   0.00000, -0.33381,  1.0814,  -0.73052,  0.24167/
62      DATA bhcn/
63     &  -1.14762,  0.90664,  1.02184, -0.33594, -3.11059, -2.64450,
64     &   0.00000, -0.66681, -5.9630,  -1.7662,  -1.6016/
65      DATA b1hcn/
66     &   2.24061,  0.49483, -0.80223,  1.42836,  4.46001,  4.38418,
67     &   0.00000,  0.88177,  6.7881,   2.8273,   2.6971/
68      DATA b2hcn/
69     &  -0.69153, -0.23858,  0.21351, -0.45182, -1.19719, -1.21663,
70     &   0.00000, -0.20593, -1.6747,  -0.70656, -0.62399/
71      DATA chcn/
72     &  -0.41463,  5.42790,  1.78780, -0.84225, -3.68653,  4.96968,
73     &   0.00000,  0.46294, -5.4664,   3.1475,  -2.2176/
74      DATA c1hcn/
75     &   0.35005, -5.78523, -2.04583,  0.78874,  2.98507, -6.43269,
76     &   0.00000, -0.69540,  5.1038,  -4.0865,   1.5738/
77      DATA c2hcn/
78     &  -0.08749,  1.48015,  0.54974, -0.19827, -0.60826,  1.84860,
79     &   0.00000,  0.19795, -1.1936,   1.1768,  -0.30869/
80
81       IF (K .GT. 11) THEN
82          XKV =         A(K) + A1(K)*LOG10(T) + A2(K)*LOG10(T)**2
83     &   + LOG10(P/P0)*(B(K) + B1(K)*LOG10(T) + B2(K)*LOG10(T)**2)
84     &   +             (C(K) + C1(K)*LOG10(T) + C2(K)*LOG10(T)**2)
85     &    *LOG10(P/P0)**2
86          XKV=10.**XKV
87C CLIP OUT THE CLEAR INTERVALS.
88          IF (A(K) .EQ. 0.0) XKV=0.0
89          TAUGAS=U*C2H6(J)*XKV
90       ELSE
91          XKV =         A(K) + A1(K)*LOG10(T) + A2(K)*LOG10(T)**2
92     &   + LOG10(P/P0)*(B(K) + B1(K)*LOG10(T) + B2(K)*LOG10(T)**2)
93     &   +             (C(K) + C1(K)*LOG10(T) + C2(K)*LOG10(T)**2)
94     &    *LOG10(P/P0)**2
95          XKV=10.**XKV
96          IF (A(K) .EQ. 0.0) XKV=0.0
97          TAUGAS=U*C2H2(J)*XKV
98          XKVhcn =   Ahcn(K) + A1hcn(K)*LOG10(T) + A2hcn(K)*LOG10(T)**2
99     &+ LOG10(P/P0)*(Bhcn(K) + B1hcn(K)*LOG10(T) + B2hcn(K)*LOG10(T)**2)
100     &+             (Chcn(K) + C1hcn(K)*LOG10(T) + C2hcn(K)*LOG10(T)**2)
101     & *LOG10(P/P0)**2
102          XKVhcn=10.**XKVhcn
103          IF (Ahcn(K) .EQ. 0.0) XKVhcn=0.0
104          TAUGAShcn=U*HCN(J)*XKVhcn
105
106c          if(j.eq.25) then
107c           print*,'J=',j,' P=',P,' HCN=',HCN(J),taugashcn,
108c     &' C2H2=',C2H2(J),taugas,' K=',K
109c          endif
110
111          taugas = taugas+taugashcn
112
113       ENDIF
114       RETURN
115       END
Note: See TracBrowser for help on using the repository browser.