Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

init_be.F90 @ 4530

Last change on this file since 4530 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

Rev	Line
[1191]	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