Context Navigation

lell.F @ 537

Last change on this file since 537 was 106, checked in by slebonnois, 14 years ago
SL : mise a jour de phytitan pour etre conforme aux sources actuelles utilisees sur gnome.
File size: 14.1 KB

Line
1	SUBROUTINE LELL(NLEVEL,Z,RHCH4L,FH2L,FARGON,TEMP,PRESS,DEN,XMU,
2	& CH4,H2,XN2,AR,IPRINT)
3	C THIS SUBROUTINE SETS UP THE INTITAL ATMOSPHERIC PROFILE FOR TITAN
4	C BASED ON THE LELLOUCH AT AL DATA. THE ROUTINE STARTS WITH INPUTS
5	C INPUTS:
6	C NLEVEL: NUMBER OF ALTITUDE LEVELS, J=1 IS AT THE TOP
7	C P PRESS GRID IN BARS
8	C
9	C ASSUMES: FARGON = 0
10	C : FH2 = 0.03
11	C
12	C OUTPUTS: RHCH4 AT SURFACE
13	C AND AT EACH LEVEL (NOT LAYER AVERAGES)
14	C TEMP (K), PRESS(BARS), DEN(CM-3), XMU = MEAN MOLECUALR WEIGHT
15	C CH4, H2, XN2, AR ARE THE NUMBER MIXING RATIOS OF THE GASES
16	C
17	C DATA: IS THE LELLOUCH ET AL VALUES DENOTED BY __LE(J)
18	C
19	C NEW VERSION FOR GCM : Z E [0,1200]
20	INTEGER ITIME
21	DIMENSION Z(NLEVEL),TEMP(NLEVEL),PRESS(NLEVEL),
22	$ DEN(NLEVEL),XMU(NLEVEL)
23	DIMENSION CH4(NLEVEL),H2(NLEVEL),XN2(NLEVEL),AR(NLEVEL)
24	C
25	DIMENSION ZLE(148), XN2LE(148), CH4LE(148), TLE(148), PLE(148),
26	& DLE(148), XMULE(148), DENMLE(148)
27
28	save ZLE,XN2LE,CH4LE,TLE,PLE,DLE,XMULE,DENMLE
29	save itime
30	data itime/0/
31	C
32	C
33	DATA ZLE/
34	& 1265., 1215., 1165., 1116., 1050., 1000., 950.0, 900.0, 880.0,
35	& 820.0, 800.0, 760.0, 700.0, 675.0, 650.0, 625.0, 600.0, 575.0,
36	& 550.0, 525.0, 500.0, 475.0, 450.0, 435.0, 420.0, 410.0, 400.0,
37	& 390.0, 380.0, 350.0, 340.0, 330.0, 320.0, 310.0, 300.0, 280.0,
38	& 260.0, 250.0, 240.0, 230.0, 220.0, 210.0,
39	& 200.0, 198.0, 196.0, 194.0, 192.0, 190.0, 188.0, 186.0, 184.0,
40	& 182.0, 180.0, 178.0, 176.0, 174.0, 172.0, 170.0, 168.0, 166.0,
41	& 164.0, 162.0, 160.0, 158.0, 156.0, 154.0, 152.0, 150.0, 148.0,
42	& 146.0, 144.0, 142.0, 140.0, 138.0, 136.0, 134.0, 132.0, 130.0,
43	& 128.0, 126.0, 124.0, 122.0, 120.0, 118.0, 116.0, 114.0, 112.0,
44	& 110.0, 108.0, 106.0, 104.0, 102.0, 100.0, 98.0, 96.0, 94.0,
45	& 92.0, 90.0, 88.0, 86.0, 84.0, 82.0, 80.0, 78.0, 76.0,
46	& 74.0, 72.0, 70.0, 68.0, 66.0, 64.0, 62.0, 60.0, 58.0,
47	& 56.0, 54.0, 52.0, 50.0, 48.0, 46.0, 44.0, 42.0, 40.0,
48	& 38.0, 36.0, 34.0, 32.0, 30.0, 28.0, 26.0, 24.0, 22.0,
49	& 20.0, 18.0, 16.0, 14.0, 12.0, 10.0, 8.0, 6.0, 5.0,
50	& 4.0, 3.0, 2.0, 1.5, 1.0, 0.5, 0.0/
51	! N2 MIXING RATIO
52	DATA XN2LE/
53	& 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985,
54	& 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985,
55	& 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985,
56	& 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985,
57	& 0.985, 0.985, 0.985, 0.985, 0.985, 0.985,
58	& 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985,
59	& 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985,
60	& 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985,
61	& 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985,
62	& 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985,
63	& 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985,
64	& 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985,
65	& 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985,
66	& 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985, 0.985,
67	& 0.985, 0.985, 0.985, 0.984, 0.983, 0.983, 0.982, 0.980, 0.979,
68	& 0.977, 0.974, 0.971, 0.966, 0.960, 0.955, 0.945, 0.935, 0.926,
69	& 0.920, 0.920, 0.920, 0.920, 0.920, 0.920, 0.920/
70	! CH4 MIXING RATIO
71	DATA CH4LE/
72	& 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015,
73	& 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015,
74	& 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015,
75	& 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015,
76	& 0.015, 0.015, 0.015, 0.015, 0.015, 0.015,
77	& 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015,
78	& 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015,
79	& 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015,
80	& 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015,
81	& 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015,
82	& 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015,
83	& 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015,
84	& 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015,
85	& 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015, 0.015,
86	& 0.015, 0.015, 0.015, 0.016, 0.017, 0.017, 0.018, 0.020, 0.021,
87	& 0.023, 0.026, 0.029, 0.034, 0.040, 0.045, 0.055, 0.065, 0.074,
88	& 0.080, 0.080, 0.080, 0.080, 0.080, 0.080, 0.080/
89	! TEMPERATURE IN k
90	DATA TLE/
91	& 183.0, 181.0, 178.0, 174.7, 169.0, 164.0, 158.0, 150.0, 147.0,
92	& 136.0, 135.0, 137.0, 143.0, 146.5, 150.0, 154.5, 159.0, 163.5,
93	& 166.5, 169.0, 171.1, 172.8, 174.3, 175.1, 175.6, 175.9, 176.0,
94	& 176.0, 175.9, 175.8, 175.7, 175.6, 175.5, 175.4, 175.3, 175.2,
95	& 175.1, 175.0, 174.9, 174.8, 174.6, 174.4,
96	& 174.0, 173.9, 173.8, 173.6, 173.4, 173.2, 173.0, 172.9, 172.7,
97	& 172.5, 172.4, 172.2, 172.0, 171.8, 171.6, 171.4, 171.2, 171.0,
98	& 170.8, 170.4, 169.9, 169.5, 169.0, 168.5, 167.9, 167.4, 166.7,
99	& 166.0, 165.3, 164.5, 163.7, 163.0, 162.4, 161.8, 161.2, 160.5,
100	& 159.6, 158.6, 157.5, 156.3, 155.3, 154.3, 153.4, 152.6, 151.9,
101	& 151.0, 150.1, 149.0, 147.8, 146.5, 145.1, 143.7, 142.1, 140.6,
102	& 139.1, 137.6, 135.8, 133.8, 131.5, 128.8, 125.8, 122.7, 119.7,
103	& 116.6, 111.9, 106.1, 100.5, 92.6, 85.9, 80.9, 77.6, 75.4,
104	& 73.8, 72.8, 72.2, 71.7, 71.5, 71.4, 71.2, 71.1, 71.1,
105	& 71.2, 71.4, 71.5, 71.8, 72.3, 73.1, 73.9, 74.7, 75.7,
106	& 76.7, 78.0, 79.2, 80.6, 82.1, 83.6, 85.5, 87.3, 88.5,
107	& 89.5, 90.5, 91.5, 92.1, 92.9, 93.7, 94.5/
108	! & 89.5, 90.5, 91.5, 92.1, 92.7, 93.3, 93.9/
109	! PRESSURE IN MILLIBARS
110	DATA PLE/
111	& 8.08e-09, 1.39e-08, 2.45e-08, 4.41e-08, 1.03e-07, 2.04e-07,
112	& 4.27e-07, 9.44e-07, 1.32e-06, 3.92e-06, 5.77e-06, 1.26e-05,
113	& 4.11e-05, 6.68e-05, 1.08e-04, 1.75e-04, 2.80e-04, 4.46e-04,
114	& 7.08e-04, 1.12e-03, 1.78e-03, 2.84e-03, 4.54e-03, 6.03e-03,
115	& 8.01e-03, 9.70e-03, 1.18e-02, 1.43e-02, 1.73e-02, 3.13e-02,
116	& 3.83e-02, 4.69e-02, 5.74e-02, 7.05e-02, 8.66e-02, 1.32e-01,
117	& 2.01e-01, 2.49e-01, 3.09e-01, 3.84e-01, 4.78e-01, 5.96e-01, !
118	& 7.59E-01, 7.98E-01, 8.37E-01, 8.77E-01, 9.20E-01, 9.64E-01,
119	& 1.01E+00, 1.06E+00, 1.11E+00, 1.16E+00, 1.22E+00, 1.28E+00,
120	& 1.35E+00, 1.41E+00, 1.48E+00, 1.55E+00, 1.63E+00, 1.71E+00,
121	& 1.79E+00, 1.88E+00, 1.97E+00, 2.07E+00, 2.17E+00, 2.29E+00,
122	& 2.40E+00, 2.52E+00, 2.65E+00, 2.78E+00, 2.93E+00, 3.08E+00,
123	& 3.24E+00, 3.41E+00, 3.59E+00, 3.78E+00, 3.90E+00, 4.19E+00,
124	& 4.42E+00, 4.66E+00, 4.91E+00, 5.19E+00, 5.48E+00, 5.78E+00,
125	& 6.10E+00, 6.45E+00, 6.82E+00, 7.22E+00, 7.63E+00, 8.08E+00,
126	& 8.56E+00, 9.06E+00, 9.61E+00, 1.02E+01, 1.08E+01, 1.15E+01,
127	& 1.22E+01, 1.30E+01, 1.38E+01, 1.47E+01, 1.57E+01, 1.68E+01,
128	& 1.80E+01, 1.93E+01, 2.07E+01, 2.23E+01, 2.40E+01, 2.60E+01,
129	& 2.83E+01, 3.10E+01, 3.42E+01, 3.79E+01, 4.23E+01, 4.75E+01,
130	& 5.34E+01, 6.01E+01, 6.79E+01, 7.67E+01, 8.67E+01, 9.81E+01,
131	& 1.11E+02, 1.26E+02, 1.42E+02, 1.61E+02, 1.83E+02, 2.07E+02,
132	& 2.35E+02, 2.65E+02, 3.00E+02, 3.40E+02, 3.83E+02, 4.32E+02,
133	& 4.87E+02, 5.47E+02, 6.14E+02, 6.88E+02, 7.70E+02, 8.59E+02,
134	& 9.57E+02, 1.06E+03, 1.12E+03, 1.18E+03, 1.24E+03, 1.30E+03,
135	& 1.34E+03, 1.38E+03, 1.43E+03, 1.48E+03/
136	! & 1.34E+03, 1.37E+03, 1.40E+03, 1.44E+03/
137	! NUMBER DENSITY
138	DATA DLE/
139	& 3.20e+08, 5.55e+08, 9.96e+08, 1.83e+09, 4.40e+09, 9.03e+09,
140	& 1.96e+10, 4.56e+10, 6.51e+10, 2.09e+11, 3.09e+11, 6.68e+11,
141	& 2.08e+12, 3.30e+12, 5.23e+12, 8.19e+12, 1.28e+13, 1.98e+13,
142	& 3.08e+13, 4.81e+13, 7.55e+13, 1.19e+14, 1.89e+14, 2.49e+14,
143	& 3.31e+14, 3.99e+14, 4.84e+14, 5.87e+14, 7.14e+14, 1.29e+15,
144	& 1.58e+15, 1.93e+15, 2.37e+15, 2.91e+15, 3.58e+15, 5.44e+15,
145	& 8.31e+15, 1.03e+16, 1.28e+16, 1.59e+16, 1.98e+16, 2.47e+16, !
146	& 3.11E+16, 3.30E+16, 3.47E+16, 3.70E+16, 3.90E+16, 4.09E+16,
147	& 4.25E+16, 4.48E+16, 4.67E+16, 4.93E+16, 5.16E+16, 5.45E+16,
148	& 5.68E+16, 5.97E+16, 6.27E+16, 6.56E+16, 6.91E+16, 7.20E+16,
149	& 7.55E+16, 7.98E+16, 8.41E+16, 8.86E+16, 9.34E+16, 9.83E+16,
150	& 1.04E+17, 1.09E+17, 1.15E+17, 1.22E+17, 1.28E+17, 1.36E+17,
151	& 1.44E+17, 1.52E+17, 1.61E+17, 1.69E+17, 1.79E+17, 1.89E+17,
152	& 2.00E+17, 2.13E+17, 2.26E+17, 2.41E+17, 2.56E+17, 2.72E+17,
153	& 2.88E+17, 3.06E+17, 3.26E+17, 3.46E+17, 3.68E+17, 3.92E+17,
154	& 4.19E+17, 4.48E+17, 4.80E+17, 5.14E+17, 5.51E+17, 5.92E+17,
155	& 6.38E+17, 6.86E+17, 7.36E+17, 7.95E+17, 8.65E+17, 9.45E+17,
156	& 1.04E+18, 1.14E+18, 1.26E+18, 1.39E+18, 1.56E+18, 1.78E+18,
157	& 2.04E+18, 2.43E+18, 2.89E+18, 3.40E+18, 3.96E+18, 4.57E+18,
158	& 5.25E+18, 6.00E+18, 6.83E+18, 7.77E+18, 8.83E+18, 1.00E+19,
159	& 1.14E+19, 1.29E+19, 1.46E+19, 1.66E+19, 1.87E+19, 2.12E+19,
160	& 2.39E+19, 2.69E+19, 3.02E+19, 3.38E+19, 3.78E+19, 4.21E+19,
161	& 4.68E+19, 5.19E+19, 5.75E+19, 6.34E+19, 6.97E+19, 7.66E+19,
162	& 8.36E+19, 9.11E+19, 9.48E+19, 9.88E+19, 1.03E+20, 1.07E+20,
163	& 1.09E+20, 1.11E+20, 1.14E+20, 1.17E+20/
164	! & 1.09E+20, 1.11E+20, 1.13E+20, 1.15E+20/
165	! MEAN MOLECULAR WEIGHT
166	DATA XMULE/
167	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8,
168	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8,
169	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8,
170	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8,
171	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8,
172	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8,
173	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8,
174	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8,
175	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8,
176	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8,
177	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8,
178	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8,
179	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8,
180	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8,
181	& 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.8, 27.7,
182	& 27.7, 27.7, 27.6, 27.6, 27.5, 27.5, 27.3, 27.2, 27.1,
183	& 27.0, 27.0, 27.0, 27.0, 27.0, 27.0, 27.0/
184	! DENSITY IN GRAMS/CUBIT CENTIMETER
185	DATA DENMLE/
186	& 1.40e-14, 2.46e-14, 4.46e-14, 8.25e-14, 2.01e-13, 4.13e-13,
187	& 8.99e-13, 2.10e-12, 3.01e-12, 9.63e-12, 1.43e-11, 3.08e-11,
188	& 9.60e-11, 1.53e-10, 2.42e-10, 3.78e-10, 5.89e-10, 9.12e-10,
189	& 1.42e-09, 2.22e-09, 3.49e-09, 5.50e-09, 8.71e-09, 1.15e-08,
190	& 1.53e-08, 1.84e-08, 2.23e-08, 2.71e-08, 3.29e-08, 5.96e-08,
191	& 7.29e-08, 8.92e-08, 1.09e-07, 1.34e-07, 1.65e-07, 2.51e-07,
192	& 3.84e-07, 4.75e-07, 5.90e-07, 7.34e-07, 9.15e-07, 1.14e-06, !
193	& 1.43E-06, 1.52E-06, 1.60E-06, 1.71E-06, 1.80E-06, 1.89E-06,
194	& 1.96E-06, 2.07E-06, 2.16E-06, 2.28E-06, 2.38E-06, 2.52E-06,
195	& 2.63E-06, 2.76E-06, 2.89E-06, 3.03E-06, 3.19E-06, 3.33E-06,
196	& 3.49E-06, 3.69E-06, 3.88E-06, 4.09E-06, 4.32E-06, 4.54E-06,
197	& 4.79E-06, 5.04E-06, 5.31E-06, 5.62E-06, 5.91E-06, 6.27E-06,
198	& 6.63E-06, 7.01E-06, 7.42E-06, 7.83E-06, 8.26E-06, 8.71E-06,
199	& 9.26E-06, 9.84E-06, 1.04E-05, 1.11E-05, 1.18E-05, 1.26E-05,
200	& 1.33E-05, 1.41E-05, 1.51E-05, 1.60E-05, 1.70E-05, 1.81E-05,
201	& 1.94E-05, 2.07E-05, 2.22E-05, 2.38E-05, 2.55E-05, 2.73E-05,
202	& 2.95E-05, 3.17E-05, 3.40E-05, 3.67E-05, 4.00E-05, 4.37E-05,
203	& 4.78E-05, 5.26E-05, 5.80E-05, 6.40E-05, 7.20E-05, 8.22E-05,
204	& 9.45E-05, 1.12E-04, 1.33E-04, 1.57E-04, 1.83E-04, 2.11E-04,
205	& 2.43E-04, 2.77E-04, 3.16E-04, 3.59E-04, 4.08E-04, 4.62E-04,
206	& 5.26E-04, 5.96E-04, 6.75E-04, 7.65E-04, 8.65E-04, 9.79E-04,
207	& 1.10E-03, 1.24E-03, 1.39E-03, 1.56E-03, 1.74E-03, 1.94E-03,
208	& 2.16E-03, 2.39E-03, 2.64E-03, 2.91E-03, 3.19E-03, 3.49E-03,
209	& 3.79E-03, 4.12E-03, 4.27E-03, 4.43E-03, 4.62E-03, 4.80E-03,
210	& 4.89E-03, 4.98E-03, 5.08E-03, 5.17E-03/
211	C
212	print*,'press ',press
213	c write(77,*) press
214	C RETURNS PRESSURE IN BARS
215	IF (ITIME.EQ.0) THEN
216	DO 201 I=1,148
217	PLE(I)=PLE(I)*0.001
218	ITIME=1
219	201 CONTINUE
220	ENDIF
221	C AND SET ARGON AND HYDROGEN
222	FARGON=0.
223	FH2=0.003
224	C
225	DO 202 J=1,NLEVEL
226	H2(J)=FH2
227	AR(J)=0.0
228	202 CONTINUE
229	C
230	DO 100 J=1,NLEVEL
231	C EXTRAPOLATE WITH ISOTHERMAL ATM ABOVE DATA POINTS
232	ISTART=1
233	IF (PRESS(J) .LT. PLE(1) ) THEN
234	TEMP(J)=TLE(1)
235	XMU(J)=XMULE(1)
236	CH4(J)=CH4LE(1)
237	XN2(J)=XN2LE(1)
238	Z(J)=ZLE(1)+ALOG(PLE(1)/PRESS(J))*40.43
239	DEN(J)=DLE(1)*PRESS(J)/PLE(1)
240	ELSE
241	DO 101 I=ISTART,147
242	C INTERPOLATE LINEAR IN LOP P
243	IF (PRESS(J) .GT. PLE(I+1) ) GO TO 101
244	FACTOR= ALOG(PRESS(J)/PLE(I) )/ALOG(PLE(I+1)/PLE(I))
245	TEMP(J)=TLE(I) + FACTOR*(TLE(I+1) - TLE(I))
246	XMU(J)=XMULE(I) + FACTOR*(XMULE(I+1) - XMULE(I))
247	CH4(J)=CH4LE(I) + FACTOR*(CH4LE(I+1) - CH4LE(I))
248	XN2(J)=XN2LE(I) + FACTOR*(XN2LE(I+1) - XN2LE(I))
249	Z(J)=ZLE(I) + FACTOR*(ZLE(I+1) - ZLE(I))
250	DEN(J)=EXP(ALOG(DLE(I))+FACTOR*(ALOG(DLE(I+1))-ALOG(DLE(I))))
251	ISTART=I
252	GO TO 100
253	101 CONTINUE
254	ENDIF
255	100 CONTINUE
256
257	C SET RHCH4 AT SURFACE...
258	RHCH4 =CH4(NLEVEL)*PRESS(NLEVEL) / PCH4(TEMP(NLEVEL))
259	C
260	IF (IPRINT .LT. 0) RETURN
261	WRITE (6,139)RHCH4,FH2,FARGON
262	DO 135 J=1,NLEVEL-1
263	WRITE(6,140)J,Z(J),PRESS(J),DEN(J),TEMP(J),
264	& CH4(J)*PRESS(J)/PCH4(TEMP(J))
265	& ,CH4(J)100.,XN2(J)100.,H2(J)100.,AR(J)100.,XMU(J)
266	& ,(TEMP(J+1)-TEMP(J))/(Z(J+1)-Z(J))
267	135 CONTINUE
268	J=NLEVEL
269	WRITE(6,140)J,Z(J),PRESS(J),DEN(J),TEMP(J),
270	& CH4(J)*PRESS(J)/PCH4(TEMP(J))
271	& ,CH4(J)100.,XN2(J)100.,H2(J)100.,AR(J)100.,XMU(J)
272	139 FORMAT(///' BACKGROUND ATMOSPHERE AT LEVELS (LELLOUCH ET AL)'/
273	& ' SURFACE HUMIDITY OF CH4:',F5.3,' H2 MIXING RATIO:',F6.4,
274	& ' ARGON SETTING:',F8.4/' LELLOUCH ET AL RESULTS '/
275	&' LVL ALTITUDE P(BARS) DEN(CM-3) TEMP RH-CH4'
276	& , ' %CH4 %N2 %H2 %AR MU DT/DZ' )
277	140 FORMAT(1X,I3,F8.3,1P2E10.3,0PF7.2,F5.2,2F6.2,2F5.2,4F6.2)
278	C
279	RETURN
280	END

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format