Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

init_be.F90 @ 1663

Last change on this file since 1663 was 1279, checked in by Laurent Fairhead, 15 years ago
Merged LMDZ4-dev branch changes r1241:1278 into the trunk Running trunk and LMDZ4-dev in LMDZOR configuration on local machine (sequential) and SX8 (4-proc) yields identical results (restart and restartphy are identical binarily) Log history from r1241 to r1278 is available by switching to source:LMDZ4/branches/LMDZ4-dev-20091210
File size: 24.8 KB

Line
1	!$Id $
2
3	SUBROUTINE init_be(pctsrf,masktr,tautr,vdeptr,scavtr,srcbe)
4
5	USE dimphy
6	USE comgeomphy
7	USE infotrac, ONLY : nbtr
8
9	IMPLICIT NONE
10	!=====================================================================
11	! Objet : prescription d'une source de Beryllium 7
12	! pour 19 niveaux verticaux
13	! (d'apres le diagramme de Lal and Peters, 1967)
14	!
15	!
16	! written by : O. Coindreau (CEA/LDG) 05/2005
17	! last modified by : A. Jamelot (LMD/CEA) 04/03/2009
18	!=====================================================================
19
20	INCLUDE "YOMCST.h"
21	INCLUDE "YOECUMF.h"
22	INCLUDE "indicesol.h"
23
24	!
25	! Input Arguments
26	!
27	REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: pctsrf !Pourcentage de sol (f(nature du sol))
28	!
29	! Output Arguments
30	!
31	REAL,DIMENSION(klon),INTENT(OUT) :: masktr ! Masque de l'echange avec la surface (possible => 1 )
32	REAL,INTENT(OUT) :: tautr ! Constante de decroissance radioactive
33	REAL,INTENT(OUT) :: vdeptr ! Vitesse de depot sec dans la couche Brownienne
34	REAL,INTENT(OUT) :: scavtr ! Coefficient de lessivage
35	REAL,DIMENSION(klon,klev),INTENT(OUT) :: srcbe ! source volumique de 7Be
36	!
37	! Local Variables
38	!
39	REAL,DIMENSION(klon) :: rlatgeo ! latitudes geomagnetiques de la grille
40	REAL :: glt ! latitude du pole geomagnetique
41	REAL :: glg ! longitude du pole geomagnetique
42	REAL :: latgeo,qcos
43	INTEGER :: k,i
44
45	WRITE(,)'PASSAGE init_be ...'
46
47	! Source actuellement definie pour klev = 19 et klev >= 39
48	IF (klev /= 19 .AND. klev<39) CALL abort_gcm("init_be","Source du be7 necessite klev=19 ou klev>=39",1)
49	!
50	! Definition des constantes
51	! -------------------------
52	tautr = 6645000.
53	vdeptr = 1.E-3
54	scavtr = 0.5
55
56	WRITE(,) '-------------- SOURCE DE BERYLLIUM ------------------- '
57	WRITE(,)'Decroissance (s): ', tautr
58	WRITE(,)'Vitesse de depot sec: ',vdeptr
59	WRITE(,)'Facteur de lessivage: ',scavtr
60
61	DO i = 1,klon
62	masktr(i) = 0.
63	IF ( NINT(pctsrf(i,1)) .EQ. 1 ) masktr(i) = 1.
64	END DO
65
66	! Premiers niveaux: source nulle
67	! ------------------------------
68	DO k = 1,6
69	DO i = 1,klon
70	srcbe(i,k) = 0.
71	END DO
72	END DO
73	!
74	! Pour les autres niveaux:
75	! 1-passer des coordonnees geographiques a la latitude geomagnetique
76	! 2-prescrire la source de Be (en 10exp5 at/g/s) dans ce repere
77	! 3-mettre la source de Be ds la bonne unite (en at/kgA/s)
78	!
79	glt=78.5*rpi/180.
80	glg=69.0*rpi/180.
81
82	DO i = 1,klon
83	qcos=sin(glt)*sin(rlatd(i))
84	qcos=qcos+cos(glt)cos(rlatd(i))cos(rlond(i)+glg)
85	IF ( qcos .LT. -1.) qcos = -1.
86	IF ( qcos .GT. 1.) qcos = 1.
87	rlatgeo(i)=rpi/2.-acos(qcos)
88	ENDDO
89
90	!===========================
91	! Cas 19 niveaux verticaux
92	!===========================
93	IF (klev.eq.19) then
94	DO k = 1,klev
95	DO i = 1,klon
96	latgeo=(180./rpi)*abs(rlatgeo(i))
97	IF ( k .EQ. 1 ) THEN
98	IF (latgeo.GE.50.0) srcbe(i,k)=0.1
99	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.09
100	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.07
101	END IF
102	IF ( k .EQ. 2 ) THEN
103	IF (latgeo.GE.50.0) srcbe(i,k)=0.12
104	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.1
105	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.09
106	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.07
107	END IF
108	IF ( k .EQ. 3 ) THEN
109	IF (latgeo.GE.50.0) srcbe(i,k)=0.14
110	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.12
111	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.1
112	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.09
113	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.07
114	END IF
115	IF ( k .EQ. 4 ) THEN
116	IF (latgeo.GE.50.0) srcbe(i,k)=0.175
117	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.16
118	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.14
119	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.12
120	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.1
121	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.09
122	END IF
123	IF ( k .EQ. 5 ) THEN
124	IF (latgeo.GE.50.0) srcbe(i,k)=0.28
125	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.26
126	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.23
127	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.175
128	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.14
129	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.12
130	END IF
131	IF ( k .EQ. 6 ) THEN
132	IF (latgeo.GE.50.0) srcbe(i,k)=0.56
133	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.49
134	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.42
135	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.28
136	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.26
137	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.245
138	END IF
139	IF ( k .EQ. 7 ) THEN
140	IF (latgeo.GE.50.0) srcbe(i,k)=1.05
141	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.875
142	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.7
143	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.52
144	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.44
145	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.385
146	END IF
147	IF ( k .EQ. 8 ) THEN
148	IF (latgeo.GE.50.0) srcbe(i,k)=2.
149	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=1.8
150	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=1.5
151	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=1.
152	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.8
153	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.75
154	END IF
155	IF ( k .EQ. 9 ) THEN
156	IF (latgeo.GE.50.0) srcbe(i,k)=4.
157	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=3.5
158	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=3.
159	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=2.5
160	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=1.8
161	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=1.4
162	END IF
163	IF ( k .EQ. 10 ) THEN
164	IF (latgeo.GE.50.0) srcbe(i,k)=8.5
165	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=8.
166	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=7.
167	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=4.5
168	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=3.5
169	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=3.
170	END IF
171	IF ( k .EQ. 11 ) THEN
172	IF (latgeo.GE.50.0) srcbe(i,k)=17.
173	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=15.
174	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=11.
175	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=8.
176	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=5.
177	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=4.
178	END IF
179	IF ( k .EQ. 12 ) THEN
180	IF (latgeo.GE.50.0) srcbe(i,k)=25.
181	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=22.
182	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=17.
183	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=11.
184	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=7.5
185	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=7.
186	END IF
187	IF ( k .EQ. 13 ) THEN
188	IF (latgeo.GE.60.0) srcbe(i,k)=33.
189	IF (latgeo.GE.50.0 .and. latgeo.LT.60.0) srcbe(i,k)=32.
190	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=30.
191	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=22.
192	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=15.
193	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=11.
194	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=8.
195	END IF
196	IF ( k .EQ. 14 ) THEN
197	IF (latgeo.GE.60.0) srcbe(i,k)=48.
198	IF (latgeo.GE.50.0 .and. latgeo.LT.60.0) srcbe(i,k)=45.
199	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=36.
200	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=26.
201	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=17.5
202	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=12.5
203	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=10.
204	END IF
205	IF ( k .EQ. 15 ) THEN
206	IF (latgeo.GE.70.0) srcbe(i,k)=58.
207	IF (latgeo.GE.60.0 .and. latgeo.LT.70.0) srcbe(i,k)=57.
208	IF (latgeo.GE.50.0 .and. latgeo.LT.60.0) srcbe(i,k)=50.
209	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=38.
210	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=25.
211	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=15.
212	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=12.5
213	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=10.
214	END IF
215	IF ( k .EQ. 16 ) THEN
216	IF (latgeo.GE.70.0) srcbe(i,k)=70.
217	IF (latgeo.GE.60.0 .and. latgeo.LT.70.0) srcbe(i,k)=65.
218	IF (latgeo.GE.50.0 .and. latgeo.LT.60.0) srcbe(i,k)=50.
219	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=32.
220	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=20.
221	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=13.
222	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=9.
223	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=7.5
224	END IF
225	IF ( k .GE. 17 ) THEN
226	IF (latgeo.GE.70.0) srcbe(i,k)=80.
227	IF (latgeo.GE.60.0 .and. latgeo.LT.70.0) srcbe(i,k)=70.
228	IF (latgeo.GE.50.0 .and. latgeo.LT.60.0) srcbe(i,k)=45.
229	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=27.
230	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=17.5
231	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=12.
232	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=8.
233	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=7.
234	END IF
235	END DO
236	END DO
237	END IF ! fin de 19 niveaux verticaux
238
239	!================================
240	! Cas 39 niveaux verticaux
241	!================================
242	IF (klev .ge. 39) then
243	DO k = 1,klev
244	DO i = 1,klon
245	latgeo=(180./rpi)*abs(rlatgeo(i))
246	IF ( k .LE. 4 ) THEN
247	IF (latgeo.GE.50.0) srcbe(i,k)=0.07
248	END IF
249	IF ( k .EQ. 5 ) THEN
250	IF (latgeo.GE.50.0) srcbe(i,k)=0.1
251	IF (latgeo.GE.20.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.09
252	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.07
253	END IF
254	IF ( k .EQ. 6 ) THEN
255	IF (latgeo.GE.50.0) srcbe(i,k)=0.14
256	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.12
257	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.1
258	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.09
259	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.07
260	END IF
261	IF ( k .EQ. 7 ) THEN
262	IF (latgeo.GE.50.0) srcbe(i,k)=0.16
263	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.16
264	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.14
265	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.12
266	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.1
267	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.09
268	END IF
269	IF ( k .EQ. 8 ) THEN
270	IF (latgeo.GE.50.0) srcbe(i,k)=0.175
271	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.16
272	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.14
273	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.12
274	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.1
275	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.1
276	END IF
277	IF ( k .EQ. 9 ) THEN
278	IF (latgeo.GE.50.0) srcbe(i,k)=0.245
279	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.21
280	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.175
281	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.14
282	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.12
283	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.12
284	END IF
285	IF ( k .EQ. 10 ) THEN
286	IF (latgeo.GE.50.0) srcbe(i,k)=0.31
287	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.28
288	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.245
289	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.21
290	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.16
291	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.14
292	END IF
293	IF ( k .EQ. 11 ) THEN
294	IF (latgeo.GE.50.0) srcbe(i,k)=0.35
295	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.3
296	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.3
297	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.2
298	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.18
299	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.16
300	END IF
301	IF ( k .EQ. 12 ) THEN
302	IF (latgeo.GE.40.0) srcbe(i,k)=0.5
303	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.4
304	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.35
305	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.3
306	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.25
307	END IF
308	IF ( k .EQ. 13 ) THEN
309	IF (latgeo.GE.50.0) srcbe(i,k)=0.8
310	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=0.7
311	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.6
312	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.5
313	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.4
314	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.35
315	END IF
316	IF ( k .EQ. 14 ) THEN
317	IF (latgeo.GE.50.0) srcbe(i,k)=1.2
318	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=1.
319	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=0.75
320	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.6
321	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.5
322	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.4
323	END IF
324	IF ( k .EQ. 15 ) THEN
325	IF (latgeo.GE.60.0) srcbe(i,k)=1.75
326	IF (latgeo.GE.50.0 .and. latgeo.LT.60.0) srcbe(i,k)=1.8
327	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=1.6
328	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=1.4
329	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=0.9
330	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=0.75
331	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.65
332	END IF
333	IF ( k .EQ. 16 ) THEN
334	IF (latgeo.GE.50.0) srcbe(i,k)=3.
335	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=2.5
336	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=1.8
337	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=1.5
338	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=1.2
339	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=0.9
340	END IF
341	IF ( k .EQ. 17 ) THEN
342	IF (latgeo.GE.50.0) srcbe(i,k)=4.
343	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=3.
344	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=2.5
345	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=2.
346	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=1.6
347	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=1.4
348	END IF
349	IF ( k .EQ. 18 ) THEN
350	IF (latgeo.GE.50.0) srcbe(i,k)=7.
351	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=6.
352	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=4.5
353	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=3.5
354	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=3.
355	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=2.
356	END IF
357	IF ( k .EQ. 19 ) THEN
358	IF (latgeo.GE.50.0) srcbe(i,k)=8.5
359	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=8.
360	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=7.
361	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=4.
362	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=3.5
363	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=3.
364	END IF
365	IF ( k .EQ. 20 ) THEN
366	IF (latgeo.GE.50.0) srcbe(i,k)=12.5
367	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=12.
368	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=8.
369	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=6.
370	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=4.
371	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=3.5
372	END IF
373	IF ( k .EQ. 21 ) THEN
374	IF (latgeo.GE.50.0) srcbe(i,k)=16.
375	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=13.
376	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=10.
377	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=7.5
378	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=4.5
379	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=4.
380	END IF
381	IF ( k .EQ. 22 ) THEN
382	IF (latgeo.GE.50.0) srcbe(i,k)=20.
383	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=17.5
384	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=12.5
385	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=9.
386	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=6.
387	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=4.5
388	END IF
389	IF ( k .EQ. 23 ) THEN
390	IF (latgeo.GE.50.0) srcbe(i,k)=25.
391	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=22.
392	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=15.
393	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=10.
394	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=7.5
395	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=6.
396	END IF
397	IF ( k .EQ. 24 ) THEN
398	IF (latgeo.GE.50.0) srcbe(i,k)=28.
399	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=26.
400	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=18.
401	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=12.
402	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=8.5
403	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=7.
404	END IF
405	IF ( k .EQ. 25 ) THEN
406	IF (latgeo.GE.50.0) srcbe(i,k)=33.
407	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=28.
408	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=20.
409	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=14.
410	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=10.
411	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=8.5
412	END IF
413	IF ( k .EQ. 26 ) THEN
414	IF (latgeo.GE.60.0) srcbe(i,k)=38.
415	IF (latgeo.GE.50.0 .and. latgeo.LT.60.0) srcbe(i,k)=36.
416	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=32.
417	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=24.
418	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=15.
419	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=11.5
420	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=10.
421	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=6.
422	END IF
423	IF ( k .EQ. 27 ) THEN
424	IF (latgeo.GE.60.0) srcbe(i,k)=46.
425	IF (latgeo.GE.50.0 .and. latgeo.LT.60.0) srcbe(i,k)=44.
426	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=35.
427	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=25.
428	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=16.
429	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=12.5
430	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=10.
431	END IF
432	IF ( k .EQ. 28 ) THEN
433	IF (latgeo.GE.60.0) srcbe(i,k)=53.
434	IF (latgeo.GE.50.0 .and. latgeo.LT.60.0) srcbe(i,k)=48.
435	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=37.
436	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=25.
437	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=16.
438	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=12.5
439	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=10.
440	END IF
441	IF ( k .EQ. 29 ) THEN
442	IF (latgeo.GE.70.0) srcbe(i,k)=58.
443	IF (latgeo.GE.60.0 .and. latgeo.LT.70.0) srcbe(i,k)=56.
444	IF (latgeo.GE.50.0 .and. latgeo.LT.60.0) srcbe(i,k)=50.
445	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=36.
446	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=24.
447	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=15.
448	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=11.5
449	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=10.
450	END IF
451	IF ( k .EQ. 30 ) THEN
452	IF (latgeo.GE.70.0) srcbe(i,k)=65.
453	IF (latgeo.GE.60.0 .and. latgeo.LT.70.0) srcbe(i,k)=60.
454	IF (latgeo.GE.50.0 .and. latgeo.LT.60.0) srcbe(i,k)=50.
455	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=35.
456	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=22.
457	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=14.
458	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=10.
459	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=9.
460	END IF
461	IF ( k .EQ. 31 ) THEN
462	IF (latgeo.GE.70.0) srcbe(i,k)=70.
463	IF (latgeo.GE.60.0 .and. latgeo.LT.70.0) srcbe(i,k)=62.
464	IF (latgeo.GE.50.0 .and. latgeo.LT.60.0) srcbe(i,k)=48.
465	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=32.
466	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=21.
467	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=13.
468	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=9.
469	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=7.6
470	END IF
471	IF ( k .EQ. 32 ) THEN
472	IF (latgeo.GE.70.0) srcbe(i,k)=80.
473	IF (latgeo.GE.60.0 .and. latgeo.LT.70.0) srcbe(i,k)=60.
474	IF (latgeo.GE.50.0 .and. latgeo.LT.60.0) srcbe(i,k)=46.
475	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=30.
476	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=17.5
477	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=11.
478	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=8.
479	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=7.4
480	END IF
481	IF ( k .GE. 33 ) THEN
482	IF (latgeo.GE.70.0) srcbe(i,k)=80.
483	IF (latgeo.GE.60.0 .and. latgeo.LT.70.0) srcbe(i,k)=70.
484	IF (latgeo.GE.50.0 .and. latgeo.LT.60.0) srcbe(i,k)=45.
485	IF (latgeo.GE.40.0 .and. latgeo.LT.50.0) srcbe(i,k)=27.
486	IF (latgeo.GE.30.0 .and. latgeo.LT.40.0) srcbe(i,k)=15.
487	IF (latgeo.GE.20.0 .and. latgeo.LT.30.0) srcbe(i,k)=10.
488	IF (latgeo.GE.10.0 .and. latgeo.LT.20.0) srcbe(i,k)=7.6
489	IF (latgeo.GE.0.0 .and. latgeo.LT.10.0) srcbe(i,k)=7.
490	END IF
491	END DO
492	END DO
493	END IF ! fin de 39 niveaux verticaux
494
495
496	!====================================
497	! Conversion de la source en U/s/kgA
498	!====================================
499	DO k = 1,klev
500	DO i = 1,klon
501	! La source est at/min/m3 -> at/s/kgA
502	! avec une masse volumique de l'air = 1.295 kg/m3
503	! 1/(60*1.295) = 0.01287
504	srcbe(i,k)=srcbe(i,k)*0.01287
505	! La source est at/min/m3 -> at/s/m3
506	! srcbe(i,k)=srcbe(i,k)*0.0166667
507	END DO
508	END DO
509
510	END SUBROUTINE init_be

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

Download in other formats:

Original Format