source: LMDZ4/branches/LMDZ4-dev-20091210/libf/phylmd/init_be.F90 @ 5440

Last change on this file since 5440 was 1191, checked in by jghattas, 16 years ago

Reecriture de phytrac et les routines concernes (Anthony Jamelot)

  • les suffix change de F -> F90 (nflxtr.F90,cltracrn.F90,initrrnpb.F90,cvltr.F90,minmaxqfi.F90,cltrac.F90,phytrac.F90)

Traitement d'un nouveau traceur berelium (optionel, toujours pour des
tests)(Anthony Jamelot)

  • radiornpb.F change du nom pour radio_decay.F90 car il traite maintenant tout les traceurs radioactives
  • ajoute init_be.F90

Nouveau interface dans phytrac pour serparer les calculs et appels
specifique a INCA avec les traitements des traceurs specifiques au LMDZ
(JG)

  • ajoute tracinca_mod.F90 pour les appeles a INCA
  • ajoute traclmdz_mod.F90 pour les calculs des traceurs specifiques a LMDZ
  • enleve fichier restartrac et ajoute la variable trs dans restartphy.nc

La convergence numerique a etait rompue uniquement pour les traceurs
LMDZ RN et PB.
File size: 24.8 KB
Line 
1!$Id $
2
3SUBROUTINE 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
510END SUBROUTINE init_be
Note: See TracBrowser for help on using the repository browser.