Context Navigation

isotrac_mod.F90 @ 4113

Last change on this file since 4113 was 3927, checked in by Laurent Fairhead, 3 years ago
Initial import of the physics wih isotopes from Camille Risi CR
Property svn:executable set to ``*
File size: 36.8 KB

Rev	Line
[3927]	1	#ifdef ISO
	2	#ifdef ISOTRAC
	3	! $Id: $
	4
	5	MODULE isotrac_mod
	6	use infotrac_phy, ONLY: niso,ntraciso,ntraceurs_zone
	7	use isotopes_mod, only: ridicule
	8
	9	IMPLICIT NONE
	10	SAVE
	11
	12	! contient toutes les variables traceurs isotopiques + les routines specifiquement
	13	! traceurs isotopiques
	14
	15	real ridicule_trac
	16	parameter (ridicule_trac=ridicule*1e4)
	17
	18	integer, save :: option_traceurs
	19	integer, save :: ntraceurs_zone_opt ! ntraceurs_zone propre à l'option
	20	! on vérifie que ça correspond bien à ntraceurs_zone d'infotrac
	21	integer, save :: ntraceurs_zoneOR
	22	!$OMP THREADPRIVATE(option_traceurs,ntraceurs_zone_opt,ntraceurs_zoneOR)
	23	integer, save :: initialisation_isotrac
	24	! 1 pour idéalisé
	25	! 0 pour lecture dans fichier
	26	!$OMP THREADPRIVATE(initialisation_isotrac)
	27
	28	! variables spécifiques aux différentes options, mais necessaires au
	29	! calcul du nombre de zones de traceurs
	30	! si option=3
	31	integer, save :: use_bassin_atlantic
	32	!$OMP THREADPRIVATE(use_bassin_atlantic)
	33	integer, save :: use_bassin_medit
	34	!$OMP THREADPRIVATE(use_bassin_medit)
	35	integer, save :: use_bassin_indian
	36	!$OMP THREADPRIVATE(use_bassin_indian)
	37	integer, save :: use_bassin_austral
	38	!$OMP THREADPRIVATE(use_bassin_austral)
	39	integer, save :: use_bassin_pacific
	40	!$OMP THREADPRIVATE(use_bassin_pacific)
	41	integer, save :: use_bassin_merarabie
	42	!$OMP THREADPRIVATE(use_bassin_merarabie)
	43	integer, save :: use_bassin_golfebengale
	44	!$OMP THREADPRIVATE(use_bassin_golfebengale)
	45	integer, save :: use_bassin_indiansud
	46	!$OMP THREADPRIVATE(use_bassin_indiansud)
	47	integer, save :: use_bassin_tropics
	48	!$OMP THREADPRIVATE(use_bassin_tropics)
	49	integer, save :: use_bassin_midlats
	50	!$OMP THREADPRIVATE(use_bassin_midlats)
	51	integer, save :: use_bassin_hauteslats
	52	!$OMP THREADPRIVATE(use_bassin_hauteslats)
	53	integer, save :: bassin_atlantic
	54	!$OMP THREADPRIVATE(bassin_atlantic)
	55	integer, save :: bassin_medit
	56	!$OMP THREADPRIVATE(bassin_medit)
	57	integer, save :: bassin_indian
	58	!$OMP THREADPRIVATE(bassin_indian)
	59	integer, save :: bassin_austral
	60	!$OMP THREADPRIVATE(bassin_austral)
	61	integer, save :: bassin_pacific
	62	!$OMP THREADPRIVATE(bassin_pacific)
	63	integer, save :: bassin_merarabie
	64	!$OMP THREADPRIVATE(bassin_merarabie)
	65	integer, save :: bassin_golfebengale
	66	!$OMP THREADPRIVATE(bassin_golfebengale)
	67	integer, save :: bassin_indiansud
	68	!$OMP THREADPRIVATE(bassin_indiansud)
	69	integer, save :: bassin_tropics
	70	!$OMP THREADPRIVATE(bassin_tropics)
	71	integer, save :: bassin_midlats
	72	!$OMP THREADPRIVATE(bassin_midlats)
	73	integer, save :: bassin_hauteslats
	74	!$OMP THREADPRIVATE(bassin_hauteslats)
	75	! si option=4
	76	integer nzone_temp
	77	parameter (nzone_temp=1)
	78	real, save :: zone_temp1,zone_tempf,zone_tempa
	79	!$OMP THREADPRIVATE(zone_temp1,zone_tempf,zone_tempa)
	80	! si option 14
	81	integer nzone_lat
	82	parameter (nzone_lat=4)
	83	integer nzone_pres
	84	parameter (nzone_pres=3)
	85	real, save :: zone_pres1,zone_presf,zone_presa
	86	!$OMP THREADPRIVATE(zone_pres1,zone_presf,zone_presa)
	87	real, save :: dlattag,lattag_min
	88	!$OMP THREADPRIVATE(dlattag,lattag_min)
	89
	90
	91	! option 1: on trace evap ocean et continent séparement
	92	! option 2: on trace evap ocean, continent et evap precip
	93	! option 3: on trace evap différents bassins océaniques
	94	! + continents + résidu
	95	! attention, choisir dans ce cas les bassins océaniques
	96	! dans iso_traceurs_opt3F90.h
	97	! option 4: tracage par température minimale
	98	! dans ce cas, on définit des bins dans iso_traceurs_opt4.h
	99	! option 5: pour AMMA: on taggue résidu/AEJ/flux mousson/Harmattan
	100	! option 6: taggage des ddfts
	101	! option 7: pour Sandrine: taggage de la vapeur à 700hPa pour omega500<-20 TODO
	102	! option 8: pour Sandrine: taggage de la vapeur entre 950 et 800hPa, omega de 0 à 25 hPa et de l'évaoration en omega<-20. TODO
	103	! option 9: taggage du condensat et de la revap precip
	104	! option 10: taggage evap oce, transpiration et evaporation
	105	! fractionante. A utiliser quand on couple avec ORCHIDEE
	106	! option 11: comme 2, mais on trace juste revap precip et reste
	107	! option 12: taggage evap oce, sol nu, canop et reste evap cont.
	108	! A utiliser quand on couple avec ORCHIDEE
	109	! option 13: taggage température minimale + revap precip
	110	! option 14: taggage lat et altitude de dernière saturation (niveaux de pression) + evap surf
	111	! otion 15: taggage irrigation
	112	! option 16: taggage precip selon saisons et fonte neige: seulement pour ORCHIDEE
	113	! option 17: taggage température minimum de condensation directement dans la convection et la cond LS, + evap sfc, condensat et precipitation
	114	! option 18: idem 17 mais on tague qsmin au lieu de Tmin
	115	! option 19: on tag vap residuelle, vap residuelle dans ddfts, sfc, cond, rev
	116	! option 20: on taggue vapeur tropicale vs vapeur extratropicale
	117	! option 21: taggage de 2 boites 3D: extratropiques (>35°) et UT tropicale (15-15°, > 500hPa)
	118	! option 22: tagage de la vapeur proccessée dans les zones très convectives
	119
	120	! ces variables sont initialisées dans traceurs_init
	121
	122	!integer ntraciso
	123	!parameter (ntraciso=(ntraceurs_zone+1)*niso)
	124	!integer ntracisoOR ! défini dans traceurs_init
	125	integer, ALLOCATABLE, DIMENSION(:), save :: index_iso
	126	!$OMP THREADPRIVATE(index_iso)
	127	integer, ALLOCATABLE, DIMENSION(:), save :: index_zone
	128	!$OMP THREADPRIVATE(index_zone)
	129	integer, ALLOCATABLE, DIMENSION(:,:), save :: index_trac_loc ! il y a déjà un index_trac dans infotrac: vérifier que c'est le même
	130	!$OMP THREADPRIVATE(index_trac_loc)
	131	character*3, ALLOCATABLE, DIMENSION(:), save :: strtrac
	132	!$OMP THREADPRIVATE(strtrac)
	133	! -> tout ça passe maintenant par infotrac
	134
	135	integer, ALLOCATABLE, DIMENSION(:), save :: bassin_map
	136	integer, ALLOCATABLE, DIMENSION(:,:), save :: boite_map
	137	!$OMP THREADPRIVATE(bassin_map,boite_map)
	138
	139
	140	! traitement recyclage et evap
	141	integer, save :: izone_cont ! pour le recyclage continental
	142	!$OMP THREADPRIVATE(izone_cont)
	143	integer, save :: izone_oce ! pour l'océan
	144	!$OMP THREADPRIVATE(izone_oce)
	145	integer, save :: izone_poubelle ! pour les petits résidus numériques
	146	!$OMP THREADPRIVATE(izone_poubelle)
	147	integer, save :: izone_init ! pour l'initialisation par défaut
	148	!$OMP THREADPRIVATE(izone_init)
	149	integer, save :: izone_revap ! pour l'évap des gouttes
	150	!$OMP THREADPRIVATE(izone_revap)
	151	integer, save :: option_revap
	152	!$OMP THREADPRIVATE(option_revap)
	153	integer, save :: option_tmin
	154	!$OMP THREADPRIVATE(option_tmin)
	155	integer, save :: option_cond
	156	!$OMP THREADPRIVATE(option_cond)
	157	integer, save :: izone_cond
	158	!$OMP THREADPRIVATE(izone_cond)
	159	real evap_franche
	160	parameter (evap_franche=1e-6) ! en kg/m2/s
	161
	162	! specifique à option 4:
	163	real, save :: zone_temp(nzone_temp-1)
	164	!$OMP THREADPRIVATE(zone_temp)
	165	! si option 5
	166	integer, save :: izone_aej
	167	!$OMP THREADPRIVATE(izone_aej)
	168	integer, save :: izone_harmattan
	169	!$OMP THREADPRIVATE(izone_harmattan)
	170	integer, save :: izone_mousson
	171	!$OMP THREADPRIVATE(izone_mousson)
	172	! si option 6
	173	integer, save :: izone_ddft
	174	!$OMP THREADPRIVATE(izone_ddft)
	175	! si option 10
	176	integer, save :: izone_contfrac
	177	!$OMP THREADPRIVATE(izone_contfrac)
	178	! si option 12
	179	integer, save :: izone_contcanop
	180	!$OMP THREADPRIVATE(izone_contcanop)
	181	! specifique à option 13:
	182	real, save :: zone_pres(nzone_pres-1)
	183	!$OMP THREADPRIVATE(zone_pres)
	184	! si option 14
	185	real, save :: zone_lat(nzone_lat-1)
	186	!$OMP THREADPRIVATE(zone_lat)
	187	! si option 15
	188	integer, save :: izone_irrig
	189	!$OMP THREADPRIVATE(izone_irrig)
	190	! si option 17
	191	real, save :: seuil_tag_tmin
	192	!$OMP THREADPRIVATE(seuil_tag_tmin)
	193	real, save :: seuil_tag_tmin_ls
	194	!$OMP THREADPRIVATE(seuil_tag_tmin_ls)
	195	integer, save :: option_seuil_tag_tmin
	196	!$OMP THREADPRIVATE(option_seuil_tag_tmin)
	197	! si option 20
	198	integer, save :: izone_trop,izone_extra
	199	real, save :: lim_tag20
	200	!$OMP THREADPRIVATE(izone_trop,izone_extra,lim_tag20)
	201	! si option 21: on garde izone_trop,izone_extra
	202	! si opt 22
	203	integer, save :: izone_conv_BT,izone_conv_UT
	204	real, save :: lim_precip_tag22
	205	!$OMP THREADPRIVATE(izone_conv_BT,izone_conv_UT,lim_precip_tag22)
	206
	207
	208	CONTAINS
	209
	210	subroutine iso_traceurs_init()
	211
	212	use IOIPSL ! getin
	213	USE infotrac_phy, ONLY: ntraciso,niso,ntraceurs_zone,index_trac
	214	USE isotopes_mod, ONLY: iso_eau,ntracisoOR,initialisation_iso, &
	215	& iso_eau_possible
	216	USE dimphy, only: klon,klev
	217
	218	implicit none
	219
	220
	221	! définition de quelles zones et quelles isotopes représentent
	222	! les traceurs
	223
	224	! inputs, outputs
	225	! ! c'est les variables dans traceurs.h qui sont modifiées
	226
	227	! locals
	228	integer itrac,izone,ixt,k
	229	integer izone_pres,izone_lat
	230	character*2 strz,strz_preslat
	231	character*1 strz_pres,strz_lat
	232	integer ntraceurs_zone_opt
	233
	234	! vérifier que on a bien l'eau comme traceurs
	235	if (iso_eau.eq.0) then
	236	write(,) 'traceurs_init 18: isotrac ne marche que si ', &
	237	& 'on met l''eau comme isotope'
	238	stop
	239	endif
	240
	241	! initialiser
	242	option_traceurs=0
	243	initialisation_isotrac=0
	244
	245	! allouer
	246	allocate (index_iso(ntraciso))
	247	allocate (index_zone(ntraciso))
	248	allocate (index_trac_loc(ntraceurs_zone,niso))
	249	allocate (strtrac(ntraceurs_zone))
	250	allocate (bassin_map(klon))
	251	allocate (boite_map(klon,klev))
	252
	253	if (initialisation_iso.eq.0) then
	254	call getin('initialisation_isotrac',initialisation_isotrac)
	255	write(,) 'initialisation_isotrac=',initialisation_isotrac
	256	endif !if (initialisation_iso.eq.0) then
	257
	258	! lire l'option de traçage
	259	call getin('option_traceurs',option_traceurs)
	260	write(,) 'option_traceurs=',option_traceurs
	261
	262	! cas général: pas de traceurs dans ORCHIDEE
	263	ntracisoOR=niso
	264
	265	! partie à éditer ! pour définir les différentes zones
	266	if (option_traceurs.eq.1) then
	267	! on trace continents/ocean
	268
	269	ntraceurs_zone_opt=2
	270	izone_cont=1
	271	izone_oce=2
	272	izone_poubelle=2 ! zone où on met les flux non physiques, de
	273	! réajustement
	274	izone_init=2 ! zone d'initialisation par défaut
	275	option_revap=0
	276	option_tmin=0
	277	izone_revap=0
	278	option_cond=0
	279
	280	strtrac(izone_cont)='con'
	281	strtrac(izone_oce)='oce'
	282
	283	elseif (option_traceurs.eq.2) then
	284	! on trace continent/ ocean/reevap des gouttes
	285
	286	ntraceurs_zone_opt=3
	287	izone_cont=1
	288	izone_oce=2
	289	izone_poubelle=2 ! zone où on met les flux non physiques, de
	290	! réajustement
	291	izone_init=2 ! zone d'initialisation par défaut
	292	option_revap=1
	293	option_tmin=0
	294	izone_revap=3
	295	option_cond=0
	296
	297	strtrac(izone_cont)='con'
	298	strtrac(izone_oce)='oce'
	299	strtrac(izone_revap)='rev'
	300
	301
	302	else if (option_traceurs.eq.3) then
	303	! on trace des bassins océaniques + un résidu. On ne trace
	304	! pas l'évap des gouttes à part
	305	! le résidu est la dernère dimension
	306
	307	! lire les use_bassin
	308	call getin('use_bassin_atlantic',use_bassin_atlantic)
	309	call getin('use_bassin_medit',use_bassin_medit)
	310	call getin('use_bassin_indian',use_bassin_indian)
	311	call getin('use_bassin_austral',use_bassin_austral)
	312	call getin('use_bassin_pacific',use_bassin_pacific)
	313	call getin('use_bassin_merarabie',use_bassin_merarabie)
	314	call getin('use_bassin_golfebengale',use_bassin_golfebengale)
	315	call getin('use_bassin_indiansud',use_bassin_indiansud)
	316	call getin('use_bassin_tropics',use_bassin_tropics)
	317	call getin('use_bassin_midlats',use_bassin_midlats)
	318	call getin('use_bassin_hauteslats',use_bassin_hauteslats)
	319
	320	write(,) 'use_bassin_atlantic=' ,use_bassin_atlantic
	321	write(,) 'use_bassin_medit=' ,use_bassin_medit
	322	write(,) 'use_bassin_indian=' ,use_bassin_indian
	323	write(,) 'use_bassin_austral=' ,use_bassin_austral
	324	write(,) 'use_bassin_merarabie=' ,use_bassin_merarabie
	325	write(,) 'use_bassin_golfebengale=' ,use_bassin_golfebengale
	326	write(,) 'use_bassin_indiansud=' ,use_bassin_indiansud
	327	write(,) 'use_bassin_tropics=' ,use_bassin_tropics
	328	write(,) 'use_bassin_midlats=' ,use_bassin_midlats
	329	write(,) 'use_bassin_hauteslats=' ,use_bassin_hauteslats
	330
	331
	332	ntraceurs_zone_opt=2 &
	333	& +use_bassin_atlantic &
	334	& +use_bassin_medit &
	335	& +use_bassin_indian &
	336	& +use_bassin_austral &
	337	& +use_bassin_pacific &
	338	& +use_bassin_merarabie &
	339	& +use_bassin_golfebengale &
	340	& +use_bassin_indiansud &
	341	& +use_bassin_tropics &
	342	& +use_bassin_midlats &
	343	& +use_bassin_hauteslats
	344
	345	izone_cont=ntraceurs_zone
	346	izone_oce=0 ! pas de sens car séparée en bassins
	347	izone_poubelle=ntraceurs_zone-1 ! zone où on met les flux non physiques, de
	348	! réajustement
	349	izone_init=ntraceurs_zone-1 ! zone d'initialisation par défaut
	350	option_revap=0 ! on ne trace pas les gouttes
	351	option_tmin=0
	352	izone_revap=0 ! pas de sens car on taggue pas les gouttes séparemment
	353	option_cond=0
	354
	355	! si on a use_bassin_indian, on n'a pas le découpage détaillé
	356	! de l'indian:
	357	#ifdef ISOVERIF
	358	if (use_bassin_indian.eq.1) then
	359	! call iso_verif_egalite(float(use_bassin_merarabie), &
	360	! & 0.0,'iso_traceurs_init 73: revoir def des bassins')
	361	if ((use_bassin_merarabie.ne.0).or. &
	362	& (use_bassin_indiansud.ne.0).or. &
	363	& (use_bassin_golfebengale.ne.0)) then
	364	write(,) 'traceurs_init 73'
	365	stop
	366	endif
	367	! call iso_verif_egalite(float(use_bassin_golfebengale), &
	368	! & 0.0,'iso_traceurs_init 73: revoir def des bassins')
	369	! call iso_verif_egalite(float(use_bassin_indiansud), &
	370	! & 0.0,'iso_traceurs_init 73: revoir def des bassins')
	371	endif
	372	#endif
	373
	374	bassin_atlantic= max(use_bassin_atlantic,1)
	375	bassin_medit=max(use_bassin_atlantic &
	376	& +use_bassin_medit,1)
	377	bassin_indian=max(use_bassin_atlantic &
	378	& +use_bassin_medit &
	379	& +use_bassin_indian,1)
	380	bassin_austral=max(use_bassin_atlantic &
	381	& +use_bassin_medit &
	382	& +use_bassin_indian &
	383	& +use_bassin_austral,1)
	384	bassin_pacific=max(use_bassin_atlantic &
	385	& +use_bassin_medit &
	386	& +use_bassin_indian &
	387	& +use_bassin_austral &
	388	& +use_bassin_pacific,1)
	389	bassin_merarabie=max(use_bassin_atlantic &
	390	& +use_bassin_medit &
	391	& +use_bassin_indian &
	392	& +use_bassin_austral &
	393	& +use_bassin_pacific &
	394	& +use_bassin_merarabie,1)
	395	bassin_golfebengale=max(use_bassin_atlantic&
	396	& +use_bassin_medit &
	397	& +use_bassin_indian &
	398	& +use_bassin_austral &
	399	& +use_bassin_pacific &
	400	& +use_bassin_merarabie &
	401	& +use_bassin_golfebengale,1)
	402	bassin_indiansud=max(use_bassin_atlantic &
	403	& +use_bassin_medit &
	404	& +use_bassin_indian &
	405	& +use_bassin_austral &
	406	& +use_bassin_pacific &
	407	& +use_bassin_merarabie &
	408	& +use_bassin_golfebengale &
	409	& +use_bassin_indiansud,1)
	410	bassin_tropics=max(use_bassin_atlantic &
	411	& +use_bassin_medit &
	412	& +use_bassin_indian &
	413	& +use_bassin_austral &
	414	& +use_bassin_pacific &
	415	& +use_bassin_merarabie &
	416	& +use_bassin_golfebengale &
	417	& +use_bassin_indiansud, &
	418	& +use_bassin_tropics,1)
	419	bassin_midlats=max(use_bassin_atlantic &
	420	& +use_bassin_medit &
	421	& +use_bassin_indian &
	422	& +use_bassin_austral &
	423	& +use_bassin_pacific &
	424	& +use_bassin_merarabie &
	425	& +use_bassin_golfebengale &
	426	& +use_bassin_indiansud &
	427	& +use_bassin_tropics &
	428	& +use_bassin_midlats,1)
	429	bassin_hauteslats=max(use_bassin_atlantic &
	430	& +use_bassin_medit &
	431	& +use_bassin_indian &
	432	& +use_bassin_austral &
	433	& +use_bassin_pacific &
	434	& +use_bassin_merarabie &
	435	& +use_bassin_golfebengale &
	436	& +use_bassin_indiansud &
	437	& +use_bassin_tropics &
	438	& +use_bassin_midlats &
	439	& +use_bassin_hauteslats,1)
	440
	441	write(,) 'bassin_atlantic=' ,bassin_atlantic
	442	write(,) 'bassin_medit=' ,bassin_medit
	443	write(,) 'bassin_indian=' ,bassin_indian
	444	write(,) 'bassin_austral=' ,bassin_austral
	445	write(,) 'bassin_merarabie=' ,bassin_merarabie
	446	write(,) 'bassin_golfebengale=' ,bassin_golfebengale
	447	write(,) 'bassin_indiansud=' ,bassin_indiansud
	448	write(,) 'bassin_tropics=' ,bassin_tropics
	449	write(,) 'bassin_midlats=' ,bassin_midlats
	450	write(,) 'bassin_hauteslats=' ,bassin_hauteslats
	451
	452	if (use_bassin_atlantic.eq.1) then
	453	strtrac(bassin_atlantic)='atl'
	454	endif
	455	if (use_bassin_medit.eq.1) then
	456	strtrac(bassin_medit)='med'
	457	endif
	458	if (use_bassin_indian.eq.1) then
	459	strtrac(bassin_indian)='ind'
	460	endif
	461	if (use_bassin_austral.eq.1) then
	462	strtrac(bassin_austral)='aus'
	463	endif
	464	if (use_bassin_pacific.eq.1) then
	465	strtrac(bassin_pacific)='pac'
	466	endif
	467	if (use_bassin_merarabie.eq.1) then
	468	strtrac(bassin_merarabie)='ara'
	469	endif
	470	if (use_bassin_golfebengale.eq.1) then
	471	strtrac(bassin_golfebengale)='ben'
	472	endif
	473	if (use_bassin_indiansud.eq.1) then
	474	strtrac(bassin_indiansud)='ins'
	475	endif
	476	if (use_bassin_tropics.eq.1) then
	477	strtrac(bassin_tropics)='tro'
	478	endif
	479	if (use_bassin_midlats.eq.1) then
	480	strtrac(bassin_midlats)='mid'
	481	endif
	482	if (use_bassin_hauteslats.eq.1) then
	483	strtrac(bassin_hauteslats)='hau'
	484	endif
	485	strtrac(ntraceurs_zone-1)='res'
	486	strtrac(ntraceurs_zone)='con'
	487
	488	else if (option_traceurs.eq.4) then
	489	! on trace les température minimales vécues
	490	! comme dans article sur LdG sauf pas de revap
	491
	492	zone_temp1=293.0 ! en K
	493	! zone_tempf=223.0 ! en K
	494	zone_tempf=243.0 ! en K
	495	! courbure de la relation entre l'indice et la température: 0 pour linéaire, <0 pour plus de détal en bas
	496
	497	! zone 1: >= zone_temp1
	498	! zone 2 à 4: intermédiaire,
	499	! zone 5: <zone_tempf
	500
	501	ntraceurs_zone_opt=nzone_temp+1
	502
	503	zone_tempa=-4.0 ! en K
	504	izone_cont=ntraceurs_zone
	505	izone_oce=ntraceurs_zone
	506	izone_poubelle=ntraceurs_zone
	507	izone_init=ntraceurs_zone ! zone d'initialisation par défaut
	508	option_revap=0
	509	option_tmin=0
	510	izone_revap=0
	511	option_cond=0
	512	do izone=1,nzone_temp
	513	write(strz,'(i2.2)') izone
	514	strtrac(izone)='t'//strz
	515	write(,) 'izone,strz,strtrac=',izone,strz,strtrac(izone)
	516	enddo
	517	strtrac(izone_poubelle)='pou'
	518
	519	! initialisation des zones de tempéarture
	520	do izone=1,nzone_temp-1
	521	zone_temp(izone)=zone_temp1+float(izone-1) &
	522	& (zone_tempafloat(izone-nzone_temp+1) &
	523	& +(zone_tempf-zone_temp1)/float(nzone_temp-2))
	524	enddo
	525	write(,) 'iso_trac_init 183: zone_temp=',zone_temp
	526
	527	elseif (option_traceurs.eq.5) then
	528	! on trace AEJ/flux de mousson/Harmattan
	529	! write(,) 'iso_traceurs_init 129'
	530
	531	ntraceurs_zone_opt=4
	532	izone_cont=1
	533	izone_oce=1
	534	izone_poubelle=1 ! zone où on met les flux non physiques, de
	535	! réajustement
	536	izone_init=1 ! zone d'initialisation par défaut
	537	option_revap=0
	538	option_tmin=0
	539	izone_revap=0
	540	izone_aej=2
	541	izone_mousson=3
	542	izone_harmattan=4
	543	option_cond=0
	544
	545	strtrac(izone_poubelle)='res'
	546	strtrac(izone_aej)='aej'
	547	strtrac(izone_mousson)='mou'
	548	strtrac(izone_harmattan)='sah'
	549
	550	elseif (option_traceurs.eq.6) then
	551	! on trace les ddfts
	552
	553	ntraceurs_zone_opt=2
	554	izone_cont=1
	555	izone_oce=1
	556	izone_poubelle=1 ! zone où on met les flux non physiques, de
	557	! réajustement
	558	izone_init=1 ! zone d'initialisation par défaut
	559	option_revap=0
	560	option_tmin=0
	561	izone_revap=0
	562	izone_ddft=2
	563	option_cond=0
	564
	565	strtrac(izone_poubelle)='res'
	566	strtrac(izone_ddft)='dft'
	567
	568	elseif (option_traceurs.eq.9) then
	569	! on trace le condensat
	570
	571	ntraceurs_zone_opt=3
	572	izone_cont=1
	573	izone_oce=1
	574	izone_poubelle=1 ! zone où on met les flux non physiques, de
	575	! réajustement
	576	izone_init=1 ! zone d'initialisation par défaut
	577	option_revap=1
	578	option_tmin=0
	579	izone_revap=2
	580	izone_cond=3
	581	option_cond=1
	582
	583	! 1 par défaut pour colorier à la fois condensat LS et
	584	! condensat convectif. Mais on peut mettre 2 si on ne veut que
	585	! collorier que le condensat convectif.
	586	call getin('option_cond',option_cond)
	587	write(,) 'option_cond=',option_cond
	588
	589	strtrac(izone_poubelle)='res'
	590	strtrac(izone_cond)='con'
	591	strtrac(izone_revap)='rev'
	592
	593	elseif (option_traceurs.eq.10) then
	594	! on trace l'évap venant de ocean/continent no frac/continent frac
	595	! utilse seulement si couplé avec ORCHIDEE
	596	#ifdef CPP_VEGET
	597	#else
	598	write(,) 'iso_traceurs_init 219: option_traceurs=10 ', &
	599	& 'inutile si on ne couple pas avec ORCHIDEE'
	600	stop
	601	#endif
	602
	603	ntraceurs_zone_opt=3
	604	izone_cont=1 ! sous-entendu non fractionnant
	605	izone_oce=2
	606	izone_poubelle=2 ! zone où on met les flux non physiques, de
	607	! réajustement
	608	izone_init=2 ! zone d'initialisation par défaut
	609	option_revap=0
	610	option_tmin=0
	611	izone_revap=0
	612	izone_contfrac=3
	613	izone_contcanop=3
	614	izone_irrig=0
	615	option_cond=0
	616
	617	strtrac(izone_oce)='oce'
	618	strtrac(izone_cont)='con'
	619	strtrac(izone_contfrac)='enu' ! evap sol nu
	620
	621	elseif (option_traceurs.eq.11) then
	622	! on trace reevap des gouttes et le reste
	623
	624	ntraceurs_zone_opt=2
	625	izone_cont=1
	626	izone_oce=1
	627	izone_poubelle=1 ! zone où on met les flux non physiques, de
	628	! réajustement
	629	izone_init=1 ! zone d'initialisation par défaut
	630	option_revap=1
	631	option_tmin=0
	632	izone_revap=2
	633	izone_irrig=0
	634	option_cond=0
	635
	636	strtrac(izone_poubelle)='res'
	637	strtrac(izone_revap)='rev'
	638
	639	elseif (option_traceurs.eq.12) then
	640	! on trace evap du sol nu, evap de la canopée, reste de l'evap cont et
	641	! evap oce
	642	#ifdef CPP_VEGET
	643	#else
	644	write(,) 'iso_traceurs_init 257: option_traceurs=10 ', &
	645	& 'inutile si on ne couple pas avec ORCHIDEE'
	646	stop
	647	#endif
	648
	649	ntraceurs_zone_opt=2
	650	izone_cont=1
	651	izone_oce=2
	652	izone_poubelle=2 ! zone où on met les flux non physiques, de
	653	! réajustement
	654	izone_init=2 ! zone d'initialisation par défaut
	655	option_revap=0
	656	option_tmin=0
	657	izone_revap=0
	658	izone_contfrac=3
	659	izone_contcanop=4
	660	izone_irrig=0
	661	option_cond=0
	662
	663	strtrac(izone_oce)='oce'
	664	strtrac(izone_cont)='con'
	665	strtrac(izone_contfrac)='enu' ! evap sol nu
	666	strtrac(izone_contcanop)='eca' ! evap canop
	667
	668	else if (option_traceurs.eq.13) then
	669	! on trace les température minimales vécues + la revap
	670	! comme dans article sur LdG
	671
	672	zone_temp1=293.0 ! en K
	673	! parameter (zone_tempf=223.0) ! en K
	674	zone_tempf=243.0 ! en K
	675	zone_tempa=-4.0 ! courbure de la relation entre l'indice et la température: 0 pour linéaire, <0 pour plus de détal en bas
	676
	677	! zone 1: >= zone_temp1
	678	! zone 2 à 4: intermédiaire,
	679	! zone 5: <zone_tempf
	680
	681	ntraceurs_zone_opt=nzone_temp+1
	682
	683	izone_cont=1
	684	izone_oce=1
	685	izone_poubelle=1
	686	izone_init=1 ! zone d'initialisation par défaut
	687	option_revap=1
	688	option_tmin=0
	689	izone_revap=ntraceurs_zone
	690	izone_irrig=0
	691	option_cond=0
	692	do izone=1,nzone_temp
	693	write(strz,'(i2.2)') izone
	694	strtrac(izone)='t'//strz
	695	write(,) 'izone,strz,strtrac=',izone,strz,strtrac(izone)
	696	enddo
	697	strtrac(izone_revap)='rev'
	698
	699	! initialisation des zones de tempéarture
	700	do izone=1,nzone_temp-1
	701	zone_temp(izone)=zone_temp1+float(izone-1) &
	702	& (zone_tempafloat(izone-nzone_temp+1) &
	703	& +(zone_tempf-zone_temp1)/float(nzone_temp-2))
	704	enddo
	705	write(,) 'zone_temp=',zone_temp
	706
	707	else if (option_traceurs.eq.14) then
	708	! on trace les pres et lat de dernière saturation définies
	709	! comme rh>90%
	710
	711	zone_pres1=600.0*100.0 ! en Pa
	712	zone_presf=300.0*100.0 ! en Pa
	713	zone_presa=0.0 ! courbure de la relation entre l'indice et la température: 0 pour linéaire, <0 pour plus de détal en bas
	714
	715	lattag_min=10.0 ! en degrès
	716	dlattag=15.0
	717
	718	! zone 1: >= zone_pres1
	719	! zone 2 à 4: intermédiaire,
	720	! zone 5: <zone_presf
	721
	722	ntraceurs_zone_opt=nzone_pres*nzone_lat+1
	723	izone_cont=ntraceurs_zone
	724	izone_oce=ntraceurs_zone
	725	izone_poubelle=ntraceurs_zone
	726	izone_init=ntraceurs_zone ! zone d'initialisation par défaut
	727	option_revap=0
	728	option_tmin=0
	729	izone_revap=0
	730	izone_irrig=0
	731	option_cond=0
	732	do izone_pres=1,nzone_pres
	733	do izone_lat=1,nzone_lat
	734	write(strz_pres,'(i1.1)') izone_pres
	735	write(strz_lat,'(i1.1)') izone_lat
	736	strz_preslat=strz_pres//strz_lat
	737	izone=izone_lat+(izone_pres-1)*nzone_lat
	738	strtrac(izone)='t'//strz_preslat
	739	write(,) 'izone_pres,izone_lat,strtrac=', &
	740	& izone_pres,izone_lat,izone,strtrac(izone)
	741	enddo !do izone_lat=1,nzone_lat
	742	enddo !do izone_pres=1,nzone_pres
	743	strtrac(ntraceurs_zone)='sfc'
	744
	745	! initialisation des zones de tempéarture
	746	do izone=1,nzone_pres-1
	747	zone_pres(izone)=zone_pres1+float(izone-1) &
	748	& (zone_presafloat(izone-nzone_pres+1) &
	749	& +(zone_presf-zone_pres1)/float(nzone_pres-2))
	750	enddo !do izone=1,nzone_pres-1
	751	write(,) 'traceurs_init 332: zone_pres=',zone_pres
	752	! stop
	753	!
	754	elseif (option_traceurs.eq.15) then
	755	! on trace l'irrigation dans ORCHIDEE
	756	#ifdef CPP_VEGET
	757	#else
	758	write(,) 'iso_traceurs_init 257: option_traceurs=15 ', &
	759	& 'inutile si on ne couple pas avec ORCHIDEE'
	760	stop
	761	#endif
	762
	763	ntraceurs_zone_opt=1
	764	izone_cont=1
	765	izone_oce=1
	766	izone_poubelle=1 ! zone où on met les flux non physiques, de
	767	! réajustement
	768	izone_init=1 ! zone d'initialisation par défaut
	769	option_revap=0
	770	option_tmin=0
	771	izone_revap=0
	772	izone_contfrac=0
	773	izone_contcanop=0
	774	izone_irrig=2
	775	option_cond=0
	776
	777	strtrac(izone_poubelle)='res'
	778	strtrac(izone_irrig)='irrig'
	779
	780	! dans ce cas particulier, il y a des traceurs dans ORCHIDEE
	781	ntracisoOR=ntraciso
	782
	783	else if ((option_traceurs.eq.17).or. &
	784	& (option_traceurs.eq.18)) then
	785	! on trace les température minimales vécues
	786	! comme dans article sur LdG sauf pas de revap
	787
	788	zone_temp1=12.0e-3 ! en kg/kg
	789	zone_tempf=0.2e-3 ! en kg/kg
	790	zone_tempa=1.2e-3 ! courbure de la relation entre l'indice et la température: 0 pour linéaire, <0 pour plus de détail en bas
	791
	792	! parameter (zone_temp1=14.0e-3) ! en kg/kg
	793	! parameter (zone_tempf=0.2e-3) ! en kg/kg
	794	! parameter (zone_tempa=0.5e-3)
	795
	796	! parameter (zone_temp1=10.0e-3) ! en kg/kg
	797	! parameter (zone_tempf=0.5e-3) ! en kg/kg
	798	! parameter (zone_tempa=0.5e-3)
	799
	800	! zone 1: >= zone_temp1
	801	! zone 2 à 4: intermédiaire,
	802	! zone 5: <zone_tempf
	803
	804	ntraceurs_zone_opt=nzone_temp+3
	805
	806	izone_cont=nzone_temp+1
	807	izone_oce=nzone_temp+1
	808	izone_poubelle=nzone_temp+1
	809	izone_init=nzone_temp+1 ! zone d'initialisation par défaut
	810	option_revap=1
	811	option_tmin=1
	812	option_cond=1
	813
	814	izone_revap=nzone_temp+3
	815	izone_cond=nzone_temp+2
	816	do izone=1,nzone_temp
	817	write(strz,'(i2.2)') izone
	818	strtrac(izone)='t'//strz
	819	write(,) 'izone,strz,strtrac=',izone,strz,strtrac(izone)
	820	enddo !do izone=1,nzone_temp
	821	strtrac(izone_poubelle)='sfc'
	822	strtrac(izone_cond)='con'
	823	strtrac(izone_revap)='rev'
	824
	825	! initialisation des zones de tempéarture
	826	do izone=1,nzone_temp-1
	827	zone_temp(izone)=zone_temp1+float(izone-1) &
	828	& (zone_tempafloat(izone-nzone_temp+1) &
	829	& +(zone_tempf-zone_temp1)/float(nzone_temp-2))
	830	enddo
	831	write(,) 'zone_temp1,zone_tempf,zone_tempa=', &
	832	& zone_temp1,zone_tempf,zone_tempa
	833	write(,) 'zone_temp=',zone_temp
	834	! stop
	835
	836	else if (option_traceurs.eq.19) then
	837
	838	zone_temp1=12.0e-3 ! en kg/kg
	839	zone_tempf=0.2e-3 ! en kg/kg
	840	zone_tempa=1.2e-3 ! courbure de la relation entre l'indice et la température: 0 pour linéaire, <0 pour plus de détail en bas
	841
	842	! parameter (zone_temp1=14.0e-3) ! en kg/kg
	843	! parameter (zone_tempf=0.2e-3) ! en kg/kg
	844	! parameter (zone_tempa=0.5e-3)
	845
	846	! parameter (zone_temp1=10.0e-3) ! en kg/kg
	847	! parameter (zone_tempf=0.5e-3) ! en kg/kg
	848	! parameter (zone_tempa=0.5e-3)
	849
	850	! zone 1: >= zone_temp1
	851	! zone 2 à 4: intermédiaire,
	852	! zone 5: <zone_tempf
	853
	854	ntraceurs_zone_opt=nzone_temp+4
	855
	856	izone_cont=nzone_temp+1
	857	izone_oce=nzone_temp+1
	858	izone_poubelle=nzone_temp+1
	859	if (option_seuil_tag_tmin.eq.1) then
	860	izone_init=nzone_temp+1 ! zone d'initialisation par défaut
	861	else
	862	izone_init=nzone_temp
	863	endif
	864	option_revap=1
	865	izone_revap=nzone_temp+3
	866	izone_cond=nzone_temp+2
	867	izone_ddft=nzone_temp+4
	868	option_tmin=1
	869	option_cond=1
	870	do izone=1,nzone_temp
	871	write(strz,'(i2.2)') izone
	872	strtrac(izone)='t'//strz
	873	write(,) 'izone,strz,strtrac=',izone,strz,strtrac(izone)
	874	enddo !do izone=1,nzone_temp
	875	strtrac(izone_poubelle)='sfc'
	876	strtrac(izone_cond)='con'
	877	strtrac(izone_revap)='rev'
	878	strtrac(izone_ddft)='dft'
	879
	880	elseif (option_traceurs.eq.20) then
	881	! on vapeur tropical/extractropicale/recyclage extractropical
	882	! pour comprendre controles humidité et isotopes subtropicaux.
	883
	884	lim_tag20=35.0
	885	call getin('lim_tag20',lim_tag20)
	886	write(,) 'lim_tag20=',lim_tag20
	887
	888	ntraceurs_zone_opt=3
	889	izone_cont=1
	890	izone_oce=1
	891	izone_poubelle=2 ! zone où on met les flux non physiques, de
	892	! réajustement
	893	izone_init=2 ! zone d'initialisation par défaut
	894	option_revap=0
	895	option_tmin=0
	896	izone_revap=0
	897	izone_trop=2
	898	izone_extra=3
	899
	900	strtrac(izone_trop)='tro' ! vapeur tropicale
	901	strtrac(izone_extra)='ext' ! vapeur extractropicale evaporée
	902	! dans les tropiques
	903	strtrac(izone_cont)='rec' ! recyclage
	904
	905	elseif (option_traceurs.eq.21) then
	906	! on trace 2 boites 3D: UT tropicale et extratropiques
	907	! fonctionnement similaire à option 5 pour taggage des zones
	908	! AMMA
	909	! write(,) 'iso_traceurs_init 129'
	910
	911	ntraceurs_zone_opt=3
	912	izone_cont=1
	913	izone_oce=1
	914	izone_poubelle=1 ! zone où on met les flux non physiques, de
	915	! réajustement
	916	izone_init=1 ! zone d'initialisation par défaut
	917	option_revap=0
	918	option_tmin=0
	919	izone_revap=0
	920	izone_trop=2
	921	izone_extra=3
	922	option_cond=0
	923
	924	strtrac(izone_poubelle)='res'
	925	strtrac(izone_trop)='tro'
	926	strtrac(izone_extra)='ext'
	927
	928	elseif (option_traceurs.eq.22) then
	929	! on trace la vapeur qui a été processée dans zones de
	930	! convections à 3 niveaux: BT, MT et UT
	931
	932	lim_precip_tag22=20.0
	933	call getin('lim_precip_tag22',lim_precip_tag22)
	934	write(,) 'lim_precip_tag22=',lim_precip_tag22
	935
	936	ntraceurs_zone_opt=3
	937	izone_cont=1
	938	izone_oce=1
	939	izone_poubelle=1 ! zone où on met les flux non physiques, de
	940	! réajustement
	941	izone_init=1 ! zone d'initialisation par défaut
	942	option_revap=0
	943	option_tmin=0
	944	izone_revap=0
	945	izone_conv_BT=2
	946	izone_conv_UT=3
	947	option_cond=0
	948
	949	strtrac(izone_poubelle)='res'
	950	strtrac(izone_conv_BT)='cbt'
	951	strtrac(izone_conv_UT)='cut'
	952
	953	else
	954	write(,) 'traceurs_init 36: option pas encore prévue'
	955	stop
	956	endif
	957
	958
	959	if (ntraceurs_zone_opt.ne.ntraceurs_zone) then
	960	write(,) 'ntraceurs_zone_opt,ntraceurs_zone=', &
	961	& ntraceurs_zone_opt,ntraceurs_zone
	962	call abort_physic ('isotrac_mod','ntraceurs_zone incoherent',1)
	963	endif
	964
	965
	966	! seuil sur le taux de condensation
	967	if (option_tmin.eq.1) then
	968	seuil_tag_tmin=0.01
	969	call getin('seuil_tag_tmin',seuil_tag_tmin)
	970	write(,) 'seuil_tag_tmin=',seuil_tag_tmin
	971
	972	seuil_tag_tmin_ls=seuil_tag_tmin
	973	call getin('seuil_tag_tmin_ls',seuil_tag_tmin_ls)
	974	write(,) 'seuil_tag_tmin_ls=',seuil_tag_tmin_ls
	975
	976	option_seuil_tag_tmin=1
	977	call getin('option_seuil_tag_tmin',option_seuil_tag_tmin)
	978	write(,) 'option_seuil_tag_tmin=',option_seuil_tag_tmin
	979	endif
	980
	981
	982	do ixt=1,niso
	983	index_zone(ixt)=0
	984	index_iso(ixt)=ixt
	985	enddo
	986	itrac=niso
	987	do izone=1,ntraceurs_zone
	988	do ixt=1,niso
	989	itrac=itrac+1
	990	index_zone(itrac)=izone
	991	index_iso(itrac)=ixt
	992	index_trac_loc(izone,ixt)=itrac
	993	if (index_trac(izone,ixt).ne.index_trac_loc(izone,ixt)) then
	994	write(,) 'isotrac 989: izone,ixt,itrac=',izone,ixt,itrac
	995	CALL abort_physic ('isotrac','isotrac 989',1)
	996	endif
	997	enddo
	998	enddo
	999	#ifdef ISOVERIF
	1000	! call iso_verif_egalite(float(itrac),float(ntraciso), &
	1001	! & 'traceurs_init 50')
	1002	if (itrac.ne.ntraciso) then
	1003	write(,) 'traceurs_init 50'
	1004	stop
	1005	endif
	1006
	1007	write(,) 'traceurs_init 65: bilan de l''init:'
	1008	write(,) 'index_zone=',index_zone(1:ntraciso)
	1009	write(,) 'index_iso=',index_iso(1:ntraciso)
	1010	write(,) 'index_trac=',index_trac(1:ntraceurs_zone,1:niso)
	1011	do izone=1,ntraceurs_zone
	1012	write(,) 'strtrac(',izone,')=',strtrac(izone)
	1013	enddo !do izone=1,ntraceurs_zone
	1014	write(,) 'ntracisoOR=',ntracisoOR
	1015	#endif
	1016
	1017	end subroutine iso_traceurs_init
	1018
	1019
	1020	END MODULE isotrac_mod
	1021	#endif
	1022	#endif

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

Download in other formats:

Original Format