source: trunk/LMDZ.TITAN/libf/phytitan/liqc2h6.F @ 1243

Last change on this file since 1243 was 175, checked in by slebonnois, 14 years ago

S.LEBONNOIS:

  • Revision majeure de la physique Titan => ajout des nuages version 10 bins (Jeremie Burgalat) Cette version reste a tester mais avec clouds=0, on reste sur l'ancienne.
  • Quelques ajouts dans la doc.
File size: 5.1 KB
Line 
1      SUBROUTINE LIQC2H6(WAVELN,XNR,XNI)
2
3*
4*
5* PERS. COMM E. QUIRICO ; eric.quirico@obs.ujf-grenoble.fr
6*
7*
8      DIMENSION W(154),XN(154),XK(154)
9        data W/
10     & 1000.  ,13.7832,12.5514,11.7883,11.1533,10.5841,10.0702,
11     & 9.6039, 9.1789, 8.7899, 8.4325, 8.1031, 7.7984, 7.5158,
12     & 7.2530, 7.0079, 6.7789, 6.5643, 6.3629, 6.1383, 5.8776,
13     & 5.6862, 5.3808, 4.4945, 3.8012, 3.6938, 3.6004, 3.5311,
14     & 3.4720, 3.4149, 3.3595, 3.3060, 3.2541, 3.2038, 3.1551,
15     & 3.1078, 3.0620, 3.0167, 2.9724, 2.9286, 2.8869, 2.8469,
16     & 2.8075, 2.7699, 2.7360, 2.7037, 2.6689, 2.6345, 2.5991,
17     & 2.5653, 2.5339, 2.5033, 2.4734, 2.4443, 2.4158, 2.3880,
18     & 2.3608, 2.3343, 2.3083, 2.2829, 2.2580, 2.2337, 2.2099,
19     & 2.1866, 2.1638, 2.1373, 2.1135, 2.0920, 2.0707, 2.0502,
20     & 2.0301, 2.0105, 1.9912, 1.9722, 1.9537, 1.9354, 1.9175,
21     & 1.9000, 1.8827, 1.8658, 1.8492, 1.8328, 1.8168, 1.8010,
22     & 1.7855, 1.7702, 1.7553, 1.7405, 1.7261, 1.7118, 1.6978,
23     & 1.6840, 1.6704, 1.6571, 1.6440, 1.6310, 1.6183, 1.6058,
24     & 1.5934, 1.5813, 1.5693, 1.5575, 1.5457, 1.5342, 1.5229,
25     & 1.5118, 1.5009, 1.4901, 1.4795, 1.4690, 1.4586, 1.4485,
26     & 1.4384, 1.4285, 1.4187, 1.4091, 1.3996, 1.3902, 1.3809,
27     & 1.3718, 1.3628, 1.3539, 1.3451, 1.3365, 1.3279, 1.3195,
28     & 1.3111, 1.3029, 1.2947, 1.2867, 1.2788, 1.2709, 1.2632,
29     & 1.2556, 1.2480, 1.2405, 1.2332, 1.2259, 1.2187, 1.2116,
30     & 1.2045, 1.1976, 1.1907, 1.1839, 1.1772, 1.1705, 1.1640,
31     & 1.1575, 1.1510, 1.1447, 1.1384, 1.1322, 1.1282, 0.0100
32     &/
33        data XN/
34     & 1.3267, 1.3272, 1.3294, 1.2944, 1.3174, 1.3197, 1.3211,
35     & 1.3217, 1.3225, 1.3229, 1.3235, 1.3242, 1.3255, 1.3286,
36     & 1.3266, 1.3408, 1.3029, 1.3101, 1.3149, 1.3171, 1.3188,
37     & 1.3196, 1.3203, 1.3226, 1.3273, 1.3298, 1.3334, 1.3432,
38     & 1.3569, 1.3637, 1.3311, 1.2866, 1.3017, 1.3067, 1.3099,
39     & 1.3120, 1.3133, 1.3139, 1.3148, 1.3153, 1.3159, 1.3162,
40     & 1.3167, 1.3170, 1.3371, 1.3600, 1.3600, 1.3600, 1.3600,
41     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600,
42     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600,
43     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600,
44     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600,
45     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600,
46     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600,
47     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600,
48     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600,
49     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600,
50     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600,
51     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600,
52     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600,
53     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600,
54     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600,
55     & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600
56     &/
57        data XK/
58     &.131E-02,.394E-03,.452E-03,.987E-03,.109E-02,.118E-02,.106E-02,
59     &.107E-02,.102E-02,.992E-03,.997E-03,.855E-03,.849E-03,.751E-03,
60     &.149E-02,.125E-02,.713E-02,.148E-02,.815E-03,.213E-03,.642E-04,
61     &.965E-04,.375E-04,.468E-04,.863E-04,.186E-05,.451E-04,.136E-02,
62     &.428E-02,.237E-01,.162E-01,.273E-02,.144E-02,.416E-03,.155E-03,
63     &.227E-03,.361E-03,.692E-04,.509E-04,.548E-04,.105E-03,.662E-04,
64     &.144E-03,.133E-03,.945E-04,.395E-04,.148E-03,.455E-05,.559E-04,
65     &.187E-04,.929E-04,.348E-04,.896E-04,.480E-03,.471E-03,.285E-03,
66     &.126E-03,.446E-03,.886E-03,.843E-03,.281E-03,.625E-04,.487E-04,
67     &.118E-04,.437E-04,.186E-05,.336E-04,.223E-05,.101E-04,.182E-04,
68     &.436E-04,.171E-04,.127E-04,.186E-04,.181E-04,.243E-04,.252E-04,
69     &.330E-04,.294E-04,.302E-04,.258E-04,.276E-04,.335E-04,.295E-04,
70     &.373E-04,.351E-04,.322E-04,.653E-04,.670E-04,.882E-04,.124E-03,
71     &.261E-04,.973E-05,.144E-04,.356E-05,.228E-05,.475E-05,.414E-05,
72     &.103E-05,.980E-06,.413E-06,.283E-06,.312E-06,.996E-06,.188E-05,
73     &.262E-05,.282E-05,.537E-05,.262E-05,.169E-05,.314E-05,.627E-05,
74     &.147E-04,.975E-05,.866E-05,.161E-04,.225E-04,.315E-04,.206E-04,
75     &.248E-04,.157E-04,.850E-05,.155E-05,.833E-06,.523E-06,.461E-06,
76     &.638E-06,.959E-06,.112E-05,.130E-05,.119E-05,.120E-05,.149E-05,
77     &.180E-05,.183E-05,.157E-05,.190E-05,.264E-05,.178E-05,.339E-05,
78     &.119E-04,.117E-04,.218E-04,.345E-04,.150E-04,.159E-04,.969E-05,
79     &.748E-05,.879E-05,.872E-05,.444E-05,.269E-05,.241E-05,.212E-05
80     &/
81*234567890123456789012345678901234567890123456789012345678901234567890
82*
83*
84      XNR=XN(1)
85      XNI=XK(1)
86      IF (WAVELN .GT. W(1))  RETURN
87      XNR=XN(154)
88      XNI=XK(154)
89      IF (WAVELN .LT. W(154)) RETURN
90      DO 100 I=2,154
91      IF (WAVELN .GT. W(I) ) GO TO 101
92 100  CONTINUE
93 101  CONTINUE
94C ALL INTERPOLATION IS IN LOG LAMBDA
95      FACTOR= (alog(WAVELN) - alog(W(I)) ) / (alog(W(I-1)) - alog(W(I)))
96C REAL PART IS LINEARLY INTERPOLATED
97      XNR=XN(I) + FACTOR*(XN(I-1) - XN(I))
98C IMAGINARY PART IS LOG INTERPOLATED
99      XNI=alog(XK(I)) + FACTOR*(alog(XK(I-1)) - alog(XK(I)))
100      XNI=exp(XNI)
101          iF(XNI.lt.1.e-7) XNI=1.e-7
102      RETURN
103      END
Note: See TracBrowser for help on using the repository browser.