Context Navigation

isotrac_mod.F90 @ 4933

Last change on this file since 4933 was 4368, checked in by lguez, 2 years ago
Sync latest trunk changes to Ocean_skin
Property svn:executable set to ``*
File size: 32.6 KB

Rev	Line
[3927]	1	#ifdef ISO
	2	#ifdef ISOTRAC
	3
	4	MODULE isotrac_mod
[4368]	5	USE infotrac_phy, ONLY: niso, ntiso, nzone
	6	USE readTracFiles_mod, ONLY: delPhase
	7	USE isotopes_mod, ONLY: ridicule, get_in
[3927]	8
[4368]	9	IMPLICIT NONE
	10	SAVE
[3927]	11
[4368]	12	!=== CONTENT: ALL THE ISOTOPIC TRACERS RELATED VARIABLES ===
	13	!
	14	! option 1: on trace evap ocean et continent separement
[3927]	15	! option 2: on trace evap ocean, continent et evap precip
[4368]	16	! option 3: on trace evap differents bassins oceaniques
	17	! + continents + residu
	18	! attention, choisir dans ce cas les bassins oceaniques
[3927]	19	! dans iso_traceurs_opt3F90.h
[4368]	20	! option 4: tracage par temperature minimale
	21	! dans ce cas, on definit des bins dans iso_traceurs_opt4.h
	22	! option 5: pour AMMA: on taggue residu/AEJ/flux mousson/Harmattan
[3927]	23	! option 6: taggage des ddfts
[4368]	24	! option 7: pour Sandrine: taggage de la vapeur a 700hPa pour omega500<-20 TODO
	25	! option 8: pour Sandrine: taggage de la vapeur entre 950 et 800hPa, omega de 0 a 25 hPa et de l'evaoration en omega<-20. TODO
[3927]	26	! option 9: taggage du condensat et de la revap precip
	27	! option 10: taggage evap oce, transpiration et evaporation
	28	! fractionante. A utiliser quand on couple avec ORCHIDEE
	29	! option 11: comme 2, mais on trace juste revap precip et reste
	30	! option 12: taggage evap oce, sol nu, canop et reste evap cont.
	31	! A utiliser quand on couple avec ORCHIDEE
[4368]	32	! option 13: taggage temperature minimale + revap precip
	33	! option 14: taggage lat et altitude de derniere saturation (niveaux de pression) + evap surf
[3927]	34	! otion 15: taggage irrigation
	35	! option 16: taggage precip selon saisons et fonte neige: seulement pour ORCHIDEE
[4368]	36	! option 17: taggage temperature minimum de condensation directement dans la convection et la cond LS, + evap sfc, condensat et precipitation
[3927]	37	! option 18: idem 17 mais on tague qsmin au lieu de Tmin
	38	! option 19: on tag vap residuelle, vap residuelle dans ddfts, sfc, cond, rev
	39	! option 20: on taggue vapeur tropicale vs vapeur extratropicale
	40	! option 21: taggage de 2 boites 3D: extratropiques (>35°) et UT tropicale (15-15°, > 500hPa)
[4368]	41	! option 22: tagage de la vapeur proccessee dans les zones tres convectives
[3927]	42
[4368]	43	!--- nzone_opt (value of nzone for the selected option) must be equal to nzone as defined in onfotrac
	44	REAL, PARAMETER :: ridicule_trac = ridicule * 1e4
	45	INTEGER, SAVE :: option_traceurs, nzone_opt, nzoneOR
	46	!$OMP THREADPRIVATE(option_traceurs,nzone_opt,nzoneOR)
	47	INTEGER, SAVE :: initialisation_isotrac
	48	!$OMP THREADPRIVATE(initialisation_isotrac)
	49	! 1 pour idealise
	50	! 0 pour lecture dans fichier
[3927]	51
[4368]	52	!=== VARIABLES SPECIFIC TO THE SELECTED OPTION, BUT NEEDED FOR THE COMPUTATION OF THE NUMBER OF ZONES ; TO BE INITIALIZED IN traceurs_init
[3927]	53
[4368]	54	!--- option 3
	55	LOGICAL, SAVE :: use_bassin_Austral, use_bassin_Atlantic, use_bassin_MidLats, use_bassin_SouthIndian, use_bassin_MerArabie
	56	!$OMP THREADPRIVATE(use_bassin_Austral, use_bassin_Atlantic, use_bassin_MidLats, use_bassin_SouthIndian, use_bassin_MerArabie)
	57	INTEGER, SAVE :: bassin_Austral, bassin_Atlantic, bassin_MidLats, bassin_SouthIndian, bassin_MerArabie
	58	!$OMP THREADPRIVATE( bassin_Austral, bassin_Atlantic, bassin_MidLats, bassin_SouthIndian, bassin_MerArabie)
	59	LOGICAL, SAVE :: use_bassin_Pacific, use_bassin_Indian, use_bassin_Tropics, use_bassin_BengalGolf, use_bassin_HighLats, use_bassin_Medit
	60	!$OMP THREADPRIVATE(use_bassin_Pacific, use_bassin_Indian, use_bassin_Tropics, use_bassin_BengalGolf, use_bassin_HighLats, use_bassin_Medit)
	61	INTEGER, SAVE :: bassin_Pacific, bassin_Indian, bassin_Tropics, bassin_BengalGolf, bassin_HighLats, bassin_Medit
	62	!$OMP THREADPRIVATE( bassin_Pacific, bassin_Indian, bassin_Tropics, bassin_BengalGolf, bassin_HighLats, bassin_Medit)
[3927]	63
[4368]	64	!--- option 4
	65	INTEGER, PARAMETER :: nzone_temp = 1
	66	REAL, SAVE :: zone_temp1, zone_tempf, zone_tempa
	67	!$OMP THREADPRIVATE(zone_temp1, zone_tempf, zone_tempa)
	68	REAL, SAVE :: zone_temp(nzone_temp-1)
[3927]	69	!$OMP THREADPRIVATE(zone_temp)
	70
[4368]	71	!--- option 5
	72	INTEGER, SAVE :: izone_aej, izone_harmattan, izone_mousson
	73	!$OMP THREADPRIVATE(izone_aej, izone_harmattan, izone_mousson)
[3927]	74
[4368]	75	!--- option 6
	76	INTEGER, SAVE :: izone_ddft
	77	!$OMP THREADPRIVATE(izone_ddft)
[3927]	78
[4368]	79	!--- option 10
	80	INTEGER, SAVE :: izone_contfrac
	81	!$OMP THREADPRIVATE(izone_contfrac)
[3927]	82
[4368]	83	!--- option 12
	84	INTEGER, SAVE :: izone_contcanop
	85	!$OMP THREADPRIVATE(izone_contcanop)
[3927]	86
[4368]	87	!--- option 13
	88	INTEGER, PARAMETER :: nzone_pres = 3
	89	REAL, SAVE :: zone_pres(nzone_pres-1)
	90	!$OMP THREADPRIVATE(zone_pres)
[3927]	91
[4368]	92	!--- option 14
	93	INTEGER, PARAMETER :: nzone_lat = 4
	94	REAL, SAVE :: zone_pres1, zone_presf, zone_presa
	95	!$OMP THREADPRIVATE(zone_pres1, zone_presf, zone_presa)
	96	REAL, SAVE :: dlattag, lattag_min, zone_lat(nzone_lat-1)
	97	!$OMP THREADPRIVATE(dlattag, lattag_min, zone_lat)
[3927]	98
[4368]	99	!--- option 15
	100	INTEGER, SAVE :: izone_irrig
	101	!$OMP THREADPRIVATE(izone_irrig)
[3927]	102
[4368]	103	!--- option 17
	104	REAL, SAVE :: seuil_tag_tmin, seuil_tag_tmin_ls
	105	!$OMP THREADPRIVATE(seuil_tag_tmin, seuil_tag_tmin_ls)
	106	INTEGER, SAVE :: option_seuil_tag_tmin
	107	!$OMP THREADPRIVATE(option_seuil_tag_tmin)
[3927]	108
[4368]	109	!--- option 20
	110	INTEGER, SAVE :: izone_trop, izone_extra
	111	!$OMP THREADPRIVATE(izone_trop, izone_extra)
	112	REAL, SAVE :: lim_tag20
	113	!$OMP THREADPRIVATE(lim_tag20)
[3927]	114
[4368]	115	!--- option 21: on garde izone_trop, izone_extra
[3927]	116
[4368]	117	!--- option 22
	118	INTEGER, SAVE :: izone_conv_BT, izone_conv_UT
	119	!$OMP THREADPRIVATE(izone_conv_BT, izone_conv_UT)
	120	REAL, SAVE :: lim_precip_tag22
	121	!$OMP THREADPRIVATE(lim_precip_tag22)
[3927]	122
[4368]	123
	124	INTEGER, ALLOCATABLE, SAVE :: index_iso(:), index_zone(:), itZonIso_loc(:,:)
	125	!$OMP THREADPRIVATE( index_iso, index_zone, itZonIso_loc)
	126	CHARACTER(LEN=3), ALLOCATABLE :: strtrac(:)
	127	!$OMP THREADPRIVATE( strtrac)
	128	INTEGER, ALLOCATABLE, SAVE :: bassin_map(:), boite_map(:,:)
	129	!$OMP THREADPRIVATE( bassin_map, boite_map)
[3927]	130
[4368]	131	!=== RECYCLING AND EVAPORATION TREATMENT
	132	INTEGER, SAVE :: izone_cont, izone_oce !--- For land and ocean recycling
	133	!$OMP THREADPRIVATE(izone_cont, izone_oce)
	134	INTEGER, SAVE :: izone_poubelle !--- For small numerical residues
	135	!$OMP THREADPRIVATE(izone_poubelle)
	136	INTEGER, SAVE :: izone_init !--- For default initialization
	137	!$OMP THREADPRIVATE(izone_init)
	138	INTEGER, SAVE :: izone_revap !--- For droplets evaporation
	139	!$OMP THREADPRIVATE(izone_revap)
	140	INTEGER, SAVE :: option_revap, option_tmin, option_cond, izone_cond
	141	!$OMP THREADPRIVATE(option_revap, option_tmin, option_cond, izone_cond)
	142	REAL, PARAMETER :: evap_franche = 1e-6 !--- In kg/m2/s
[3927]	143
[4368]	144	CONTAINS
[3927]	145
[4368]	146	SUBROUTINE iso_traceurs_init()
[3927]	147
[4368]	148	USE infotrac_phy, ONLY: itZonIso, isoName, isoZone
	149	USE isotopes_mod, ONLY: iso_eau, ntracisoOR, initialisation_iso
	150	USE dimphy, ONLY: klon, klev
	151	USE strings_mod, ONLY: int2str, strStack, strTail, strHead, fmsg
[3927]	152
[4368]	153	IMPLICIT NONE
	154	! Define which zones and isotopes correspond to isotopic tagging tracers
	155	! Modify traceurs.h variables
	156	INTEGER :: izone, ixt, k
	157	INTEGER :: izone_pres, izone_lat
	158	INTEGER :: nzone_opt
[3927]	159
[4368]	160	IF(fmsg("traceurs_init 18: isotrac ne marche que si on met l'eau comme isotope", 'iso_traceurs_init', iso_eau==0)) STOP
[3927]	161
[4368]	162	!--- Initialize
	163	option_traceurs = 0
	164	initialisation_isotrac = 0
[3927]	165
[4368]	166	!--- Allocate
	167	ALLOCATE(index_iso (ntiso))
	168	ALLOCATE(index_zone(ntiso))
	169	ALLOCATE(itZonIso_loc(nzone,niso))
	170	ALLOCATE(strtrac(nzone))
	171	ALLOCATE(bassin_map(klon))
	172	ALLOCATE( boite_map(klon,klev))
[3927]	173
[4368]	174	IF(initialisation_iso == 0) CALL get_in('initialisation_isotrac', initialisation_isotrac)
[3927]	175
[4368]	176	!--- Read tracing option
	177	CALL get_in('option_traceurs', option_traceurs)
[3927]	178
[4368]	179	!--- Genral case: no traceurs in ORCHIDEE
	180	ntracisoOR=niso
[3927]	181
[4368]	182	! partie a editer ! pour definir les differentes zones
	183	SELECT CASE(option_traceurs)
	184	!========================================================================================================================
	185	CASE(1) !=== TRACING LAND/OCEAN
	186	!========================================================================================================================
	187	nzone_opt=2
	188	izone_cont=1
	189	izone_oce=2
	190	izone_poubelle=2 ! zone ou on met les flux non physiques, de reajustement
	191	izone_init=2 ! zone d'initialisation par defaut
	192	option_revap=0
	193	option_tmin=0
	194	izone_revap=0
	195	option_cond=0
	196	strtrac(izone_cont) = 'con'
	197	strtrac(izone_oce) = 'oce'
	198	!========================================================================================================================
	199	CASE(2) !=== TRACING LAND/OCEAN/DROPLETS REEVAPORATION
	200	!========================================================================================================================
	201	nzone_opt=3
	202	izone_cont=1
	203	izone_oce=2
	204	izone_poubelle=2 ! zone ou on met les flux non physiques, de reajustement
	205	izone_init=2 ! zone d'initialisation par defaut
	206	option_revap=1
	207	option_tmin=0
	208	izone_revap=3
	209	option_cond=0
	210	strtrac(izone_cont) = 'con'
	211	strtrac(izone_oce) = 'oce'
	212	strtrac(izone_revap)= 'rev'
	213	!========================================================================================================================
	214	CASE(3) !=== TRACING OCEANS BASINS + RESIDUE (LAST DIMENSION). NO DROPLETS EVAPORATION TRACING.
	215	!========================================================================================================================
	216	! lire les use_bassin
	217	CALL get_in('use_bassin_Atlantic', use_bassin_Atlantic)
	218	CALL get_in('use_bassin_Medit', use_bassin_Medit)
	219	CALL get_in('use_bassin_Indian', use_bassin_Indian)
	220	CALL get_in('use_bassin_Austral', use_bassin_Austral)
	221	CALL get_in('use_bassin_Pacific', use_bassin_Pacific)
	222	CALL get_in('use_bassin_MerArabie', use_bassin_MerArabie)
	223	CALL get_in('use_bassin_BengalGolf', use_bassin_BengalGolf)
	224	CALL get_in('use_bassin_SouthIndian',use_bassin_SouthIndian)
	225	CALL get_in('use_bassin_Tropics', use_bassin_Tropics)
	226	CALL get_in('use_bassin_Midlats', use_bassin_Midlats)
	227	CALL get_in('use_bassin_HighLats', use_bassin_HighLats)
	228	nzone_opt = 2 + COUNT([use_bassin_Atlantic, use_bassin_Medit, use_bassin_Indian, &
	229	use_bassin_Austral, use_bassin_Pacific, use_bassin_MerArabie, use_bassin_BengalGolf, &
	230	use_bassin_SouthIndian, use_bassin_Tropics, use_bassin_Midlats, use_bassin_HighLats])
	231	izone_cont=nzone
	232	izone_oce=0 ! pas de sens car separee en bassins
	233	izone_poubelle=nzone-1 ! zone ou on met les flux non physiques, de reajustement
	234	izone_init=nzone-1 ! zone d'initialisation par defaut
	235	option_revap=0 ! on ne trace pas les gouttes
	236	option_tmin=0
	237	izone_revap=0 ! pas de sens car on taggue pas les gouttes separemment
	238	option_cond=0
[3927]	239	#ifdef ISOVERIF
[4368]	240	IF(use_bassin_Indian) THEN !=== NON COMPATIBLE WITH A DETAILED INDIAN CUTTING
	241	IF(use_bassin_MerArabie .OR. use_bassin_SouthIndian .OR. use_bassin_BengalGolf) THEN
	242	WRITE(,) 'traceurs_init 73'; STOP
	243	END IF
	244	! CALL iso_verif_egalite(float(use_bassin_MerArabie), 0.0, 'iso_traceurs_init 73: revoir def des bassins')
	245	! CALL iso_verif_egalite(float(use_bassin_BengalGolf), 0.0, 'iso_traceurs_init 73: revoir def des bassins')
	246	! CALL iso_verif_egalite(float(use_bassin_SouthIndian), 0.0, 'iso_traceurs_init 73: revoir def des bassins')
	247	END IF
[3927]	248	#endif
[4368]	249	bassin_Atlantic = 1
	250	bassin_Medit = bassin_Atlantic + COUNT([use_bassin_Medit]); WRITE(,) 'bassin_Atlantic =' ,bassin_Atlantic
	251	bassin_Indian = bassin_Medit + COUNT([use_bassin_Indian]); WRITE(,) 'bassin_Medit =' ,bassin_Medit
	252	bassin_Austral = bassin_Indian + COUNT([use_bassin_Austral]); WRITE(,) 'bassin_Indian =' ,bassin_Indian
	253	bassin_Pacific = bassin_Austral + COUNT([use_bassin_Pacific]); WRITE(,) 'bassin_Austral =' ,bassin_Austral
	254	bassin_MerArabie = bassin_Pacific + COUNT([use_bassin_MerArabie]); WRITE(,) 'bassin_MerArabie =' ,bassin_MerArabie
	255	bassin_BengalGolf = bassin_MerArabie + COUNT([use_bassin_BengalGolf]); WRITE(,) 'bassin_BengalGolf =' ,bassin_BengalGolf
	256	bassin_SouthIndian= bassin_BengalGolf + COUNT([use_bassin_SouthIndian]); WRITE(,) 'bassin_SouthIndian =' ,bassin_SouthIndian
	257	bassin_Tropics = bassin_SouthIndian + COUNT([use_bassin_Tropics]); WRITE(,) 'bassin_Tropics =' ,bassin_Tropics
	258	bassin_MidLats = bassin_Tropics + COUNT([use_bassin_MidLats]); WRITE(,) 'bassin_MidLats =' ,bassin_MidLats
	259	bassin_HighLats = bassin_MidLats + COUNT([use_bassin_HighLats]); WRITE(,) 'bassin_HighLats =' ,bassin_HighLats
	260	IF(use_bassin_atlantic ) strtrac(bassin_atlantic) = 'atl'
	261	IF(use_bassin_medit ) strtrac(bassin_medit) = 'med'
	262	IF(use_bassin_indian ) strtrac(bassin_indian) = 'ind'
	263	IF(use_bassin_austral ) strtrac(bassin_austral) = 'aus'
	264	IF(use_bassin_pacific ) strtrac(bassin_pacific) = 'pac'
	265	IF(use_bassin_merarabie ) strtrac(bassin_merarabie) = 'ara'
	266	IF(use_bassin_BengalGolf ) strtrac(bassin_BengalGolf) = 'ben'
	267	IF(use_bassin_SouthIndian) strtrac(bassin_SouthIndian)= 'ins'
	268	IF(use_bassin_tropics ) strtrac(bassin_tropics) = 'tro'
	269	IF(use_bassin_midlats ) strtrac(bassin_midlats) = 'mid'
	270	IF(use_bassin_HighLats ) strtrac(bassin_HighLats) = 'hau'
	271	strtrac(nzone-1)='res'
	272	strtrac(nzone)='con'
	273	!========================================================================================================================
	274	CASE(4) !=== TRACING MINIMAL EXPERIENCED TEMPERATURE AS IN THE ARTICLE ON LfG, EXCEPT NO REVAPORATION
	275	!========================================================================================================================
	276	zone_temp1 = 293.0 ! en K
	277	! zone_tempf = 223.0 ! en K
	278	zone_tempf = 243.0 ! en K
	279	! courbure de la relation entre l'indice et la temperature: 0 pour lineaire, <0 pour plus de detal en bas
[3927]	280	! zone 1: >= zone_temp1
[4368]	281	! zone 2 a 4: intermediaire,
[3927]	282	! zone 5: <zone_tempf
[4368]	283	nzone_opt=nzone_temp+1
	284	zone_tempa=-4.0 ! en K
	285	izone_cont=nzone
	286	izone_oce=nzone
	287	izone_poubelle=nzone
	288	izone_init=nzone ! zone d'initialisation par defaut
	289	option_revap=0
	290	option_tmin=0
	291	izone_revap=0
	292	option_cond=0
	293	DO izone=1,nzone_temp
	294	strtrac(izone) = 't'//TRIM(int2str(izone))
	295	WRITE(,) 'izone, strtrac=', izone, strtrac(izone)
	296	END DO
	297	strtrac(izone_poubelle)='pou'
	298	! Initialization of temperatures zones
	299	DO izone=1,nzone_temp-1
	300	zone_temp(izone) = zone_temp1+float(izone-1) &
	301	* (zone_tempa*float(izone-nzone_temp+1) &
	302	+ (zone_tempf-zone_temp1)/float(nzone_temp-2))
	303	END DO
	304	WRITE(,) 'iso_trac_init 183: zone_temp=', zone_temp
	305	!========================================================================================================================
	306	CASE(5) !=== TRACING AEJ/MOONSOON FLUX/Harmattan
	307	!========================================================================================================================
	308	! WRITE,) 'iso_traceurs_init 129'
	309	nzone_opt=4
	310	izone_cont=1
	311	izone_oce=1
	312	izone_poubelle=1 ! zone ou on met les flux non physiques, de reajustement
	313	izone_init=1 ! zone d'initialisation par defaut
	314	option_revap=0
	315	option_tmin=0
	316	izone_revap=0
	317	izone_aej=2
	318	izone_mousson=3
	319	izone_harmattan=4
	320	option_cond=0
	321	strtrac(izone_poubelle) = 'res'
	322	strtrac(izone_aej) = 'aej'
	323	strtrac(izone_mousson) = 'mou'
	324	strtrac(izone_harmattan)= 'sah'
	325	!========================================================================================================================
	326	CASE(6) !=== TRACING DDFTS
	327	!========================================================================================================================
	328	nzone_opt=2
	329	izone_cont=1
	330	izone_oce=1
	331	izone_poubelle=1 ! zone ou on met les flux non physiques, de reajustement
	332	izone_init=1 ! zone d'initialisation par defaut
	333	option_revap=0
	334	option_tmin=0
	335	izone_revap=0
	336	izone_ddft=2
	337	option_cond=0
	338	strtrac(izone_poubelle)='res'
	339	strtrac(izone_ddft)='dft'
	340	!========================================================================================================================
	341	CASE(9) !=== TRACING CONDENSATION
	342	!========================================================================================================================
	343	nzone_opt=3
	344	izone_cont=1
	345	izone_oce=1
	346	izone_poubelle=1 ! zone ou on met les flux non physiques, de reajustement
	347	izone_init=1 ! zone d'initialisation par defaut
	348	option_revap=1
	349	option_tmin=0
	350	izone_revap=2
	351	izone_cond=3
	352	option_cond=1
	353	! 1 par defaut pour colorier a la fois condensat LS et condensat convectif.
	354	! Mais on peut mettre 2 si on ne veut que colorier que le condensat convectif.
	355	CALL get_in('option_cond',option_cond)
	356	strtrac(izone_poubelle)='res'
	357	strtrac(izone_cond)='con'
	358	strtrac(izone_revap)='rev'
	359	!========================================================================================================================
	360	CASE(10) !=== TRACING EVAPORATION FROM OCEAN/LAND, NON FRAC/LAND FRAC ; ONLY WHEN COUPLED WITH ORCHIDEE
	361	!========================================================================================================================
	362	#ifndef CPP_VEGET
	363	WRITE(,) 'iso_traceurs_init 219: option_traceurs=10 inutile si on ne couple pas avec ORCHIDEE'; STOP
[3927]	364	#endif
[4368]	365	nzone_opt=3
	366	izone_cont=1 ! sous-entendu non fractionnant
	367	izone_oce=2
	368	izone_poubelle=2 ! zone ou on met les flux non physiques, de reajustement
	369	izone_init=2 ! zone d'initialisation par defaut
	370	option_revap=0
	371	option_tmin=0
	372	izone_revap=0
	373	izone_contfrac=3
	374	izone_contcanop=3
	375	izone_irrig=0
	376	option_cond=0
	377	strtrac(izone_oce)='oce'
	378	strtrac(izone_cont)='con'
	379	strtrac(izone_contfrac)='enu' ! evap sol nu
	380	!========================================================================================================================
	381	CASE(11) !=== TRACING DROPLETS REEVAPORATION + REST
	382	!========================================================================================================================
	383	nzone_opt=2
	384	izone_cont=1
	385	izone_oce=1
	386	izone_poubelle=1 ! zone ou on met les flux non physiques, de reajustement
	387	izone_init=1 ! zone d'initialisation par defaut
	388	option_revap=1
	389	option_tmin=0
	390	izone_revap=2
	391	izone_irrig=0
	392	option_cond=0
	393	strtrac(izone_poubelle)='res'
	394	strtrac(izone_revap)='rev'
	395	!========================================================================================================================
	396	CASE(12) !=== TRACING NAKED GROUND EVAPORATION, CANOPY EVAPORATION, REST OF LAND EVAPORATION AND OCEAN EVAPORATION
	397	!========================================================================================================================
	398	#ifndef CPP_VEGET
	399	WRITE(,) 'iso_traceurs_init 257: option_traceurs=10 inutile si on ne couple pas avec ORCHIDEE'; STOP
[3927]	400	#endif
[4368]	401	nzone_opt=2
	402	izone_cont=1
	403	izone_oce=2
	404	izone_poubelle=2 ! zone ou on met les flux non physiques, de reajustement
	405	izone_init=2 ! zone d'initialisation par defaut
	406	option_revap=0
	407	option_tmin=0
	408	izone_revap=0
	409	izone_contfrac=3
	410	izone_contcanop=4
	411	izone_irrig=0
	412	option_cond=0
	413	strtrac(izone_oce)='oce'
	414	strtrac(izone_cont)='con'
	415	strtrac(izone_contfrac)='enu' ! evap sol nu
	416	strtrac(izone_contcanop)='eca'! evap canop
	417	!========================================================================================================================
	418	CASE(13) !=== TRACING MINIMUM EXPERIENCED TEMPERATIRES + REEVAPORATION AS IN THE ARTICLE ON LdG
	419	!========================================================================================================================
	420	zone_temp1=293.0 ! en K
	421	! zone_tempf=223.0 ! en K
	422	zone_tempf=243.0 ! en K
	423	zone_tempa=-4.0 ! courbure de la relation entre l'indice et la temperature: 0 pour lineaire, <0 pour plus de detal en bas
	424	! zone 1: >= zone_temp1
	425	! zone 2 a 4: intermediaire,
	426	! zone 5: <zone_tempf
	427	nzone_opt=nzone_temp+1
	428	izone_cont=1
	429	izone_oce=1
	430	izone_poubelle=1
	431	izone_init=1 ! zone d'initialisation par defaut
	432	option_revap=1
	433	option_tmin=0
	434	izone_revap=nzone
	435	izone_irrig=0
	436	option_cond=0
	437	DO izone=1,nzone_temp
	438	strtrac(izone) = 't'//TRIM(int2str(izone))
	439	WRITE(,) 'izone, strtrac = ', izone, strtrac(izone)
	440	END DO
	441	strtrac(izone_revap)='rev'
	442	! initialisation des zones de tempearture
	443	DO izone=1,nzone_temp-1
	444	zone_temp(izone) = zone_temp1+float(izone-1) &
	445	(zone_tempafloat(izone-nzone_temp+1) &
	446	+(zone_tempf-zone_temp1)/float(nzone_temp-2))
	447	END DO
	448	WRITE(,) 'zone_temp=',zone_temp
	449	!========================================================================================================================
	450	CASE(14) !=== TRACING PRES AND LAT OF LAST SATURATION DEFINED AS rh>90%
	451	!========================================================================================================================
	452	zone_pres1=600.0*100.0 ! en Pa
	453	zone_presf=300.0*100.0 ! en Pa
	454	zone_presa=0.0 ! courbure de la relation entre l'indice et la temperature: 0 pour lineaire
	455	lattag_min=10.0 ! en degres
	456	dlattag=15.0
	457	! zone 1: >= zone_pres1
	458	! zone 2 a 4: intermediaire,
	459	! zone 5: <zone_presf
	460	nzone_opt=nzone_pres*nzone_lat+1
	461	izone_cont=nzone
	462	izone_oce=nzone
	463	izone_poubelle=nzone
	464	izone_init=nzone ! zone d'initialisation par defaut
	465	option_revap=0
	466	option_tmin=0
	467	izone_revap=0
	468	izone_irrig=0
	469	option_cond=0
	470	DO izone_pres=1,nzone_pres
	471	DO izone_lat=1,nzone_lat
	472	izone=izone_lat+(izone_pres-1)*nzone_lat
	473	strtrac(izone) = 't'//TRIM(int2str(izone_pres))//TRIM(int2str(izone_lat))
	474	write(,) 'izone_pres, izone_lat, izone, strtrac = ',izone_pres, izone_lat, izone, strtrac(izone)
	475	END DO
	476	END DO
	477	strtrac(nzone)='sfc'
	478	! initialisation des zones de temperature
	479	DO izone=1,nzone_pres-1
	480	zone_pres(izone) = zone_pres1+float(izone-1) &
	481	(zone_presafloat(izone-nzone_pres+1) &
	482	+(zone_presf-zone_pres1)/float(nzone_pres-2))
	483	END DO
	484	WRITE(,) 'traceurs_init 332: zone_pres=',zone_pres
	485	!========================================================================================================================
	486	CASE(15) !=== TRACING IRRIGATION IN ORCHIDEE
	487	!========================================================================================================================
	488	#ifndef CPP_VEGET
	489	WRITE(,) 'iso_traceurs_init 257: option_traceurs=15 inutile si on ne couple pas avec ORCHIDEE'; STOP
[3927]	490	#endif
[4368]	491	nzone_opt=1
	492	izone_cont=1
	493	izone_oce=1
	494	izone_poubelle=1 ! zone ou on met les flux non physiques, de reajustement
	495	izone_init=1 ! zone d'initialisation par defaut
	496	option_revap=0
	497	option_tmin=0
	498	izone_revap=0
	499	izone_contfrac=0
	500	izone_contcanop=0
	501	izone_irrig=2
	502	option_cond=0
	503	strtrac(izone_poubelle)='res'
	504	strtrac(izone_irrig)='irrig'
	505	! dans ce cas particulier, il y a des traceurs dans ORCHIDEE
	506	ntracisoOR=ntiso
	507	!========================================================================================================================
	508	CASE(17,18) !=== TRACING MINIMAL EXPERIENCES TEMPERATURES AS IN THE ARTICLE ABOUT LdG, BUT NO EVAPORATION
	509	!========================================================================================================================
	510	zone_temp1=12.0e-3 ! en kg/kg
	511	zone_tempf=0.2e-3 ! en kg/kg
	512	zone_tempa=1.2e-3 ! courbure de la relation entre l'indice et la temperature: 0 pour lineaire
	513	! zone_temp1=14.0e-3 ! en kg/kg
	514	! zone_tempf=0.2e-3 ! en kg/kg
	515	! zone_tempa=0.5e-3
	516	! zone_temp1=10.0e-3 ! en kg/kg
	517	! zone_tempf=0.5e-3 ! en kg/kg
	518	! zone_tempa=0.5e-3
	519	! zone 1: >= zone_temp1
	520	! zone 2 a 4: intermediaire,
	521	! zone 5: <zone_tempf
	522	nzone_opt=nzone_temp+3
	523	izone_cont=nzone_temp+1
	524	izone_oce=nzone_temp+1
	525	izone_poubelle=nzone_temp+1
	526	izone_init=nzone_temp+1 ! zone d'initialisation par defaut
	527	option_revap=1
	528	option_tmin=1
	529	option_cond=1
	530	izone_revap=nzone_temp+3
	531	izone_cond=nzone_temp+2
	532	DO izone=1,nzone_temp
	533	strtrac(izone) = 't'//TRIM(int2str(izone))
	534	WRITE(,) 'izone, strtrac = ', izone, strtrac(izone)
	535	END DO !do izone=1,nzone_temp
	536	strtrac(izone_poubelle)='sfc'
	537	strtrac(izone_cond)='con'
	538	strtrac(izone_revap)='rev'
	539	! initialisation des zones de tempearture
	540	DO izone=1,nzone_temp-1
	541	zone_temp(izone) = zone_temp1+float(izone-1) &
	542	(zone_tempafloat(izone-nzone_temp+1) &
	543	+(zone_tempf-zone_temp1)/float(nzone_temp-2))
	544	END DO
	545	WRITE(,) 'zone_temp1,zone_tempf,zone_tempa=',zone_temp1,zone_tempf,zone_tempa
	546	WRITE(,) 'zone_temp=',zone_temp
	547	! STOP
	548	!========================================================================================================================
	549	CASE(19) !=== TRACING TROPICAL AND EXTRATROPICAL VAPOUR
	550	!========================================================================================================================
	551	zone_temp1=12.0e-3 ! en kg/kg
	552	zone_tempf=0.2e-3 ! en kg/kg
	553	zone_tempa=1.2e-3 ! courbure de la relation entre l'indice et la temperature: 0 pour lineaire, <0 pour plus de detail en bas
	554	! zone_temp1=14.0e-3 ! en kg/kg
	555	! zone_tempf=0.2e-3 ! en kg/kg
	556	! zone_tempa=0.5e-3
	557	! zone_temp1=10.0e-3 ! en kg/kg
	558	! zone_tempf=0.5e-3 ! en kg/kg
	559	! zone_tempa=0.5e-3
	560	! zone 1: >= zone_temp1
	561	! zone 2 a 4: intermediaire,
	562	! zone 5: <zone_tempf
	563	nzone_opt=nzone_temp+4
	564	izone_cont=nzone_temp+1
	565	izone_oce=nzone_temp+1
	566	izone_poubelle=nzone_temp+1
	567	IF(option_seuil_tag_tmin == 1) THEN
	568	izone_init=nzone_temp+1 ! zone d'initialisation par defaut
	569	ELSE
[3927]	570	izone_init=nzone_temp
[4368]	571	END IF
	572	option_revap=1
	573	izone_revap=nzone_temp+3
	574	izone_cond=nzone_temp+2
	575	izone_ddft=nzone_temp+4
	576	option_tmin=1
	577	option_cond=1
	578	DO izone=1,nzone_temp
	579	strtrac(izone) = 't'//TRIM(int2str(izone))
	580	WRITE(,) 'izone, strtrac = ', izone, strtrac(izone)
	581	END DO
	582	strtrac(izone_poubelle)='sfc'
	583	strtrac(izone_cond)='con'
	584	strtrac(izone_revap)='rev'
	585	strtrac(izone_ddft)='dft'
	586	!========================================================================================================================
	587	CASE(20) !=== TRACING TROPICAL/EXTRATROPICAL/EXTRATROPICAL RECYCLING TO STUDY HUMIDITY AND SUBTROPICAL ISOTOPES CONTROL
	588	!========================================================================================================================
	589	CALL get_in('lim_tag20', lim_tag20, 35.0)
	590	nzone_opt=3
	591	izone_cont=1
	592	izone_oce=1
	593	izone_poubelle=2 ! zone ou on met les flux non physiques, de reajustement
	594	izone_init=2 ! zone d'initialisation par defaut
	595	option_revap=0
	596	option_tmin=0
	597	izone_revap=0
	598	izone_trop=2
	599	izone_extra=3
	600	strtrac(izone_trop)='tro' ! tropical vapour
	601	strtrac(izone_extra)='ext' ! extratropical vapour evaporated in the tropics
	602	strtrac(izone_cont)='rec' ! recycling
	603	!========================================================================================================================
	604	CASE(21) !=== TRACING TWO 3D BOXES: TROPICAL UT AND EXTRATROPICS ; SIMILAR TO 5 FOR AMMA ZONES TAGGING
	605	!========================================================================================================================
	606	! WRITE(,) 'iso_traceurs_init 129'
	607	nzone_opt=3
	608	izone_cont=1
	609	izone_oce=1
	610	izone_poubelle=1 ! zone ou on met les flux non physiques, de reajustement
	611	izone_init=1 ! zone d'initialisation par defaut
	612	option_revap=0
	613	option_tmin=0
	614	izone_revap=0
	615	izone_trop=2
	616	izone_extra=3
	617	option_cond=0
	618	strtrac(izone_poubelle)='res'
	619	strtrac(izone_trop)='tro'
	620	strtrac(izone_extra)='ext'
	621	!========================================================================================================================
	622	CASE(22) !=== TRACING WATER VAPOUR PROCESSED IN THE 3-LEVELS SCONVECTION ZONES BT, MT AND UT
	623	!========================================================================================================================
	624	CALL get_in('lim_precip_tag22', lim_precip_tag22, 20.0)
	625	nzone_opt=3
	626	izone_cont=1
	627	izone_oce=1
	628	izone_poubelle=1 ! zone ou on met les flux non physiques, de reajustement
	629	izone_init=1 ! zone d'initialisation par defaut
	630	option_revap=0
	631	option_tmin=0
	632	izone_revap=0
	633	izone_conv_BT=2
	634	izone_conv_UT=3
	635	option_cond=0
	636	strtrac(izone_poubelle)='res'
	637	strtrac(izone_conv_BT)='cbt'
	638	strtrac(izone_conv_UT)='cut'
	639	CASE DEFAULT
	640	WRITE(,) 'traceurs_init 36: option pas encore prevue' ; STOP
	641	END SELECT
[3927]	642
[4368]	643	IF(nzone_opt /= nzone) THEN
	644	WRITE(,) 'nzone_opt, nzone=', nzone_opt, nzone
	645	CALL abort_physic ('isotrac_mod','nzone incoherent',1)
	646	END IF
[3927]	647
[4368]	648	!--- Condensation rate threshold
	649	IF(option_tmin == 1) THEN
	650	seuil_tag_tmin = 0.01
	651	CALL get_in('seuil_tag_tmin', seuil_tag_tmin, 0.01)
	652	CALL get_in('seuil_tag_tmin_ls', seuil_tag_tmin_ls, seuil_tag_tmin)
	653	CALL get_in('option_seuil_tag_tmin', option_seuil_tag_tmin, 1)
	654	END IF
	655	DO ixt=1,niso
	656	index_zone(ixt)=0
	657	index_iso(ixt)=ixt
	658	END DO
[3927]	659
[4368]	660	index_zone = [(INDEX(isoZone, strTail( isoName(ixt) ,'_')), ixt=1, ntiso)]
	661	index_iso = [(INDEX(isoName, strHead(delPhase(isoName(ixt)),'_')), ixt=1, ntiso)]
	662	itZonIso_loc = itZonIso(:,:)
[3927]	663	#ifdef ISOVERIF
[4368]	664	WRITE(,) 'traceurs_init 65: bilan de l''init:'
	665	WRITE(,) 'index_zone = '//TRIM(strStack(int2str(index_zone(1:ntiso))))
	666	WRITE(,) 'index_iso = '//TRIM(strStack(int2str(index_iso (1:ntiso))))
	667	DO izone=1,nzone
	668	WRITE(,)'itZonIso('//TRIM(int2str(izone))//',:) = '//strStack(int2str(itZonIso(izone,:)))
	669	END DO
	670	DO izone=1,nzone
	671	WRITE(,)'strtrac('//TRIM(int2str(izone))//',:) = '//TRIM(strtrac(izone))
	672	END DO
	673	WRITE(,) 'ntracisoOR=',ntracisoOR
[3927]	674	#endif
	675
[4368]	676	END SUBROUTINE iso_traceurs_init
[3927]	677
	678	END MODULE isotrac_mod
	679	#endif
	680	#endif

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

Download in other formats:

Original Format