source: trunk/LMDZ.TITAN.old/libf/phytitan/tholin_cvd.F @ 3303

Last change on this file since 3303 was 1056, checked in by slebonnois, 11 years ago

SL: Titan runs ! see DOC/chantiers/commit_importants.log
File size: 12.2 KB
Line 
1      SUBROUTINE THOLIN_CVD(WLNMETERS,XNR,XNI)   ! original
2C THIS ROUTINE RETURNS THE OPTICAL CONSTANTS OF LABORATORY PRODUCED
3C FOR WAVELN ABOVE AND BELOW THE DATA, VALUES ARE EXTRAPOLATED AS
4CONSTANT
5      DIMENSION W(385),XN(385),XK(385)
6
7      data W/
8     &314.0367,303.2380,292.8110,282.7425,273.0199,263.6319,254.5665,
9     &245.8130,237.3603,229.3791,225.0075,220.7190,216.5123,212.3858,
10     &208.3380,204.3672,200.4722,196.6515,192.9034,189.2269,185.6204,
11     &182.0827,178.6123,175.2082,171.9691,169.6769,167.4153,165.1838,
12     &162.9822,160.8098,158.6664,156.5516,154.4649,152.4061,150.3746,
13     &148.3703,146.3927,144.4415,142.5163,140.6907,139.3035,137.9300,
14     &136.5701,135.2235,133.8902,132.5701,131.2629,129.9687,128.6872,
15     &127.4184,126.1621,124.9181,123.6865,122.4669,120.8540,117.6756,
16     &114.5806,111.5671,108.6329,105.7758,102.9938,100.2851, 97.6476,
17     & 95.0794, 92.5788, 90.1439, 87.7731, 85.4646, 83.2169, 81.0670,
18     & 79.0859, 77.1533, 75.2678, 73.4285, 71.6340, 69.8835, 68.1757,
19     & 66.5097, 64.8843, 63.2987, 61.7519, 60.2428, 58.7706, 57.3344,
20     & 55.8598, 54.2566, 52.6993, 51.1868, 49.7176, 48.2906, 46.9046,
21     & 45.5584, 44.2508, 42.9807, 41.7471, 40.5489, 39.3851, 38.2547,
22     & 37.1567, 36.3218, 35.9298, 35.5420, 35.1584, 34.7789, 34.4036,
23     & 34.0323, 33.6650, 33.3016, 32.9422, 32.5867, 32.2350, 31.8871,
24     & 31.5429, 31.2025, 30.7221, 30.0381, 29.3693, 28.7155, 28.0761,
25     & 27.4511, 26.8399, 26.2424, 25.6581, 25.0869, 24.5283, 23.9823,
26     & 23.4483, 22.9263, 22.4159, 22.0109, 21.7266, 21.4461, 21.1691,
27     & 20.8957, 20.6259, 20.3595, 20.0966, 19.8371, 19.5809, 19.3281,
28     & 19.0785, 18.8321, 18.5889, 18.3489, 18.2124, 18.1771, 18.1420,
29     & 18.1069, 18.0719, 18.0369, 18.0021, 17.9672, 17.9325, 17.8978,
30     & 17.8632, 17.8287, 17.7942, 17.7598, 17.7254, 17.5766, 17.3366,
31     & 17.0999, 16.8665, 16.6362, 16.4091, 16.1851, 15.9641, 15.7461,
32     & 15.5312, 15.3191, 15.1100, 14.9037, 14.7002, 14.4995, 14.3561,
33     & 14.2503, 14.1453, 14.0410, 13.9375, 13.8347, 13.7327, 13.6315,
34     & 13.5310, 13.4313, 13.3323, 13.2340, 13.1364, 13.0396, 12.9435,
35     & 12.8549, 12.7705, 12.6867, 12.6035, 12.5208, 12.4387, 12.3570,
36     & 12.2760, 12.1954, 12.1154, 12.0359, 11.9569, 11.8785, 11.8006,
37     & 11.7231, 11.6697, 11.6266, 11.5836, 11.5408, 11.4982, 11.4557,
38     & 11.4133, 11.3712, 11.3291, 11.2873, 11.2456, 11.2040, 11.1626,
39     & 11.1213, 11.0803, 11.0412, 11.0029, 10.9648, 10.9268, 10.8889,
40     & 10.8512, 10.8136, 10.7761, 10.7387, 10.7015, 10.6644, 10.6274,
41     & 10.5906, 10.5539, 10.5173, 10.4062, 10.2779, 10.1512, 10.0261,
42     &  9.9025,  9.7804,  9.6598,  9.5407,  9.4231,  9.3070,  9.1922,
43     &  9.0789,  8.9670,  8.8564,  8.7473,  8.6658,  8.5898,  8.5144,
44     &  8.4397,  8.3656,  8.2922,  8.2194,  8.1473,  8.0758,  8.0049,
45     &  7.9346,  7.8650,  7.7960,  7.7276,  7.6597,  7.6149,  7.5728,
46     &  7.5309,  7.4892,  7.4478,  7.4066,  7.3656,  7.3249,  7.2844,
47     &  7.2441,  7.2040,  7.1642,  7.1246,  7.0851,  7.0460,  7.0194,
48     &  6.9935,  6.9678,  6.9421,  6.9165,  6.8911,  6.8657,  6.8404,
49     &  6.8152,  6.7901,  6.7651,  6.7402,  6.7154,  6.6906,  6.6660,
50     &  5.7876,  5.4343,  5.1027,  4.7912,  4.4988,  4.2242,  3.9664,
51     &  3.7243,  3.4970,  3.2836,  3.0832,  2.8950,  2.7183,  2.5524,
52     &  2.3966,  2.2504,  2.1130,  1.9841,  1.8630,  1.7493,  1.6425,
53     &  1.5422,  1.4481,  1.3597,  1.2767,  1.1988,  1.1257,  1.0570,
54     &  0.9924,  0.9319,  0.8750,  0.8216,  0.7715,  0.7244,  0.6802,
55     &  0.6386,  0.5997,  0.5631,  0.5287,  0.4964,  0.4661,  0.4377,
56     &  0.4110,  0.3859,  0.3623,  0.3402,  0.3195,  0.3000,  0.2817,
57     &  0.2645,  0.2483,  0.2332,  0.2189,  0.2056,  0.1930,  0.1812,
58     &  0.1702,  0.1598,  0.1500,  0.1409,  0.1323,  0.1242,  0.1166,
59     &  0.1095,  0.1028,  0.0966,  0.0907,  0.0851,  0.0799,  0.0751,
60     &  0.0705,  0.0662,  0.0621,  0.0583,  0.0548,  0.0514,  0.0483,
61     &  0.0453,  0.0426,  0.0400,  0.0375,  0.0353,  0.0331,  0.0311,
62     &  0.0292,  0.0274,  0.0257,  0.0242,  0.0227,  0.0213,  0.0200/
63
64        data XN/
65     &  1.9168,  1.9212,  1.9258,  1.9306,  1.9353,  1.9397,  1.9428,
66     &  1.9421,  1.9430,  1.9444,  1.9454,  1.9469,  1.9476,  1.9474,
67     &  1.9463,  1.9439,  1.9387,  1.9312,  1.9260,  1.9210,  1.9160,
68     &  1.9101,  1.9037,  1.8982,  1.8931,  1.8894,  1.8856,  1.8818,
69     &  1.8779,  1.8738,  1.8693,  1.8646,  1.8593,  1.8530,  1.8475,
70     &  1.8423,  1.8369,  1.8314,  1.8247,  1.8189,  1.8149,  1.8112,
71     &  1.8076,  1.8040,  1.8003,  1.7964,  1.7930,  1.7901,  1.7874,
72     &  1.7849,  1.7826,  1.7805,  1.7785,  1.7768,  1.7745,  1.7715,
73     &  1.7688,  1.7652,  1.7613,  1.7574,  1.7535,  1.7490,  1.7442,
74     &  1.7399,  1.7363,  1.7337,  1.7317,  1.7308,  1.7283,  1.7248,
75     &  1.7198,  1.7164,  1.7132,  1.7066,  1.7027,  1.7001,  1.6959,
76     &  1.6923,  1.6897,  1.6877,  1.6859,  1.6847,  1.6840,  1.6835,
77     &  1.6831,  1.6839,  1.6875,  1.6914,  1.6930,  1.6934,  1.6951,
78     &  1.6949,  1.6925,  1.6932,  1.6943,  1.6962,  1.6976,  1.6983,
79     &  1.6997,  1.7005,  1.7008,  1.7012,  1.7018,  1.7024,  1.7028,
80     &  1.7032,  1.7037,  1.7044,  1.7051,  1.7055,  1.7058,  1.7061,
81     &  1.7066,  1.7064,  1.7044,  1.7014,  1.7018,  1.7014,  1.7019,
82     &  1.7027,  1.7029,  1.7060,  1.7076,  1.7074,  1.7083,  1.7091,
83     &  1.7106,  1.7116,  1.7096,  1.7104,  1.7114,  1.7112,  1.7114,
84     &  1.7120,  1.7116,  1.7121,  1.7144,  1.7165,  1.7159,  1.7126,
85     &  1.7098,  1.7112,  1.7108,  1.7105,  1.7115,  1.7117,  1.7119,
86     &  1.7121,  1.7121,  1.7121,  1.7121,  1.7120,  1.7120,  1.7119,
87     &  1.7118,  1.7117,  1.7117,  1.7116,  1.7116,  1.7116,  1.7118,
88     &  1.7121,  1.7127,  1.7135,  1.7152,  1.7170,  1.7164,  1.7131,
89     &  1.7127,  1.7140,  1.7158,  1.7193,  1.7226,  1.7255,  1.7276,
90     &  1.7285,  1.7293,  1.7297,  1.7298,  1.7295,  1.7285,  1.7270,
91     &  1.7243,  1.7193,  1.7132,  1.7094,  1.7064,  1.7033,  1.7016,
92     &  1.7010,  1.7015,  1.7032,  1.7037,  1.7040,  1.7040,  1.7039,
93     &  1.7037,  1.7035,  1.7032,  1.7029,  1.7026,  1.7023,  1.7019,
94     &  1.7015,  1.7013,  1.7011,  1.7009,  1.7008,  1.7006,  1.7003,
95     &  1.6998,  1.6992,  1.6986,  1.6977,  1.6964,  1.6955,  1.6949,
96     &  1.6944,  1.6940,  1.6936,  1.6933,  1.6930,  1.6927,  1.6925,
97     &  1.6923,  1.6922,  1.6921,  1.6920,  1.6919,  1.6918,  1.6918,
98     &  1.6917,  1.6917,  1.6917,  1.6918,  1.6921,  1.6926,  1.6933,
99     &  1.6942,  1.6953,  1.6966,  1.6981,  1.6999,  1.7020,  1.7048,
100     &  1.7082,  1.7107,  1.7131,  1.7153,  1.7171,  1.7187,  1.7204,
101     &  1.7221,  1.7240,  1.7259,  1.7278,  1.7312,  1.7359,  1.7411,
102     &  1.7464,  1.7494,  1.7514,  1.7526,  1.7522,  1.7499,  1.7421,
103     &  1.7382,  1.7424,  1.7539,  1.7647,  1.7797,  1.8034,  1.7889,
104     &  1.7552,  1.7328,  1.7341,  1.7439,  1.7669,  1.7864,  1.7968,
105     &  1.8076,  1.8169,  1.8204,  1.8160,  1.7935,  1.7718,  1.7534,
106     &  1.7302,  1.6941,  1.6947,  1.7192,  1.7730,  1.8033,  1.7930,
107     &  1.4623,  1.5506,  1.5865,  1.6135,  1.6100,  1.6293,  1.6393,
108     &  1.6479,  1.6500,  1.6500,  1.6229,  1.5865,  1.5951,  1.6102,
109     &  1.6200,  1.6279,  1.6300,  1.6300,  1.6320,  1.6400,  1.6400,
110     &  1.6400,  1.6400,  1.6400,  1.6447,  1.6496,  1.6500,  1.6500,
111     &  1.6515,  1.6557,  1.6599,  1.6651,  1.6704,  1.6758,  1.6813,
112     &  1.6875,  1.6938,  1.7001,  1.7053,  1.7105,  1.7157,  1.7149,
113     &  1.6874,  1.6587,  1.6376,  1.6349,  1.6442,  1.6568,  1.6684,
114     &  1.6796,  1.6800,  1.6777,  1.6711,  1.6645,  1.6590,  1.6556,
115     &  1.6523,  1.6557,  1.6731,  1.6906,  1.7102,  1.7323,  1.7500,
116     &  1.7189,  1.6778,  1.6211,  1.5510,  1.4599,  1.3682,  1.2730,
117     &  1.1930,  1.1130,  1.0330,  0.9586,  0.9234,  0.8881,  0.8528,
118     &  0.8176,  0.8053,  0.8083,  0.8189,  0.8295,  0.8401,  0.8507,
119     &  0.8614,  0.8721,  0.8829,  0.8936,  0.9044,  0.9151,  0.9200/
120
121        data XK/
122     &.139E-01,.173E-01,.214E-01,.266E-01,.330E-01,.409E-01,.507E-01,
123     &.601E-01,.676E-01,.759E-01,.810E-01,.869E-01,.942E-01,.102E+00,
124     &.111E+00,.120E+00,.130E+00,.135E+00,.140E+00,.145E+00,.150E+00,
125     &.155E+00,.159E+00,.162E+00,.165E+00,.167E+00,.170E+00,.172E+00,
126     &.174E+00,.177E+00,.180E+00,.182E+00,.185E+00,.186E+00,.187E+00,
127     &.188E+00,.188E+00,.189E+00,.190E+00,.188E+00,.187E+00,.186E+00,
128     &.185E+00,.184E+00,.183E+00,.181E+00,.179E+00,.177E+00,.175E+00,
129     &.173E+00,.171E+00,.169E+00,.167E+00,.165E+00,.163E+00,.158E+00,
130     &.156E+00,.154E+00,.152E+00,.150E+00,.148E+00,.146E+00,.143E+00,
131     &.140E+00,.137E+00,.133E+00,.130E+00,.128E+00,.128E+00,.127E+00,
132     &.126E+00,.123E+00,.123E+00,.120E+00,.115E+00,.112E+00,.109E+00,
133     &.104E+00,.989E-01,.939E-01,.890E-01,.841E-01,.792E-01,.747E-01,
134     &.695E-01,.627E-01,.571E-01,.552E-01,.546E-01,.528E-01,.513E-01,
135     &.523E-01,.488E-01,.452E-01,.423E-01,.401E-01,.391E-01,.376E-01,
136     &.363E-01,.356E-01,.351E-01,.347E-01,.342E-01,.341E-01,.340E-01,
137     &.337E-01,.334E-01,.332E-01,.335E-01,.338E-01,.342E-01,.344E-01,
138     &.351E-01,.362E-01,.379E-01,.356E-01,.333E-01,.322E-01,.300E-01,
139     &.294E-01,.270E-01,.258E-01,.271E-01,.273E-01,.267E-01,.268E-01,
140     &.265E-01,.289E-01,.288E-01,.274E-01,.280E-01,.285E-01,.282E-01,
141     &.286E-01,.287E-01,.274E-01,.273E-01,.296E-01,.334E-01,.351E-01,
142     &.327E-01,.311E-01,.319E-01,.301E-01,.300E-01,.300E-01,.302E-01,
143     &.304E-01,.307E-01,.309E-01,.310E-01,.311E-01,.311E-01,.312E-01,
144     &.312E-01,.312E-01,.312E-01,.311E-01,.311E-01,.309E-01,.305E-01,
145     &.302E-01,.299E-01,.296E-01,.292E-01,.316E-01,.345E-01,.344E-01,
146     &.322E-01,.304E-01,.293E-01,.286E-01,.305E-01,.326E-01,.355E-01,
147     &.378E-01,.403E-01,.430E-01,.461E-01,.493E-01,.528E-01,.565E-01,
148     &.604E-01,.646E-01,.636E-01,.621E-01,.607E-01,.587E-01,.556E-01,
149     &.528E-01,.503E-01,.494E-01,.495E-01,.496E-01,.496E-01,.497E-01,
150     &.498E-01,.499E-01,.499E-01,.500E-01,.501E-01,.501E-01,.502E-01,
151     &.503E-01,.503E-01,.504E-01,.504E-01,.506E-01,.509E-01,.512E-01,
152     &.515E-01,.518E-01,.521E-01,.525E-01,.522E-01,.516E-01,.510E-01,
153     &.504E-01,.498E-01,.493E-01,.488E-01,.483E-01,.477E-01,.472E-01,
154     &.467E-01,.462E-01,.457E-01,.453E-01,.448E-01,.443E-01,.438E-01,
155     &.434E-01,.429E-01,.424E-01,.411E-01,.395E-01,.380E-01,.365E-01,
156     &.351E-01,.338E-01,.324E-01,.312E-01,.300E-01,.288E-01,.277E-01,
157     &.281E-01,.289E-01,.297E-01,.305E-01,.311E-01,.317E-01,.323E-01,
158     &.329E-01,.336E-01,.342E-01,.342E-01,.337E-01,.348E-01,.365E-01,
159     &.415E-01,.475E-01,.529E-01,.589E-01,.656E-01,.705E-01,.744E-01,
160     &.645E-01,.558E-01,.515E-01,.532E-01,.550E-01,.806E-01,.121E+00,
161     &.139E+00,.112E+00,.921E-01,.777E-01,.696E-01,.819E-01,.914E-01,
162     &.102E+00,.121E+00,.145E+00,.173E+00,.199E+00,.204E+00,.208E+00,
163     &.213E+00,.190E+00,.147E+00,.114E+00,.111E+00,.157E+00,.220E+00,
164     &.428E-01,.109E-01,.807E-02,.141E-01,.147E-01,.100E-01,.128E-01,
165     &.190E-01,.321E-01,.252E-01,.273E-01,.124E-01,.719E-02,.719E-02,
166     &.719E-02,.719E-02,.719E-02,.719E-02,.719E-02,.719E-02,.719E-02,
167     &.718E-02,.709E-02,.701E-02,.720E-02,.699E-02,.759E-02,.723E-02,
168     &.694E-02,.723E-02,.777E-02,.864E-02,.102E-01,.126E-01,.160E-01,
169     &.204E-01,.264E-01,.346E-01,.435E-01,.551E-01,.701E-01,.889E-01,
170     &.109E+00,.129E+00,.147E+00,.168E+00,.188E+00,.205E+00,.226E+00,
171     &.250E+00,.275E+00,.294E+00,.299E+00,.303E+00,.309E+00,.318E+00,
172     &.328E+00,.339E+00,.353E+00,.368E+00,.407E+00,.484E+00,.584E+00,
173     &.698E+00,.788E+00,.873E+00,.950E+00,.101E+01,.107E+01,.104E+01,
174     &.988E+00,.942E+00,.899E+00,.850E+00,.758E+00,.676E+00,.603E+00,
175     &.538E+00,.461E+00,.389E+00,.326E+00,.274E+00,.230E+00,.193E+00,
176     &.164E+00,.140E+00,.119E+00,.101E+00,.861E-01,.733E-01,.681E-01/
177
178      XNR=XN(1)
179      XNI=XK(1)
180      WAVELN = WLNMETERS !* 1.E6  ! ARRAYS ARE IN MICRONS
181      IF (WAVELN .GE. W(1))  RETURN
182      XNR=XN(385)
183      XNI=XK(385)
184      IF (WAVELN .LE. W(385)) RETURN
185      DO 100 I=2,384
186      IF (WAVELN .GT. W(I) ) GOTO 101
187 100  CONTINUE
188 101  CONTINUE
189C ALL INTERPOLATION IS IN LOG LAMBDA
190      FACTOR= (alog(WAVELN) - alog(W(I)) ) / (alog(W(I-1)) - alog(W(I)))
191C REAL PART IS LINEARLY INTERPOLATED
192      XNR=XN(I) + FACTOR*(XN(I-1) - XN(I))
193C IMAGINARY PART IS LOG INTERPOLATED
194      XNI=alog(XK(I)) + FACTOR*(alog(XK(I-1)) - alog(XK(I)))
195      XNI=exp(XNI)
196
197      RETURN
198      END
199
200
Note: See TracBrowser for help on using the repository browser.