Context Navigation

infotrac.F90 @ 4539

Last change on this file since 4539 was 4523, checked in by evignon, 16 months ago
merge de la branche blowing snow vers la trunk premiere tentative Etienne
Property copyright set to Name of program: LMDZ Creation date: 1984 Version: LMDZ5 License: CeCILL version 2 Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539 See the license file in the root directory Property svn:keywords set to `Id`
File size: 27.5 KB

Rev	Line
[4064]	1	!$Id: infotrac.F90 4523 2023-05-03 16:21:08Z fhourdin $
[1279]	2	!
[1114]	3	MODULE infotrac
	4
[4325]	5	USE strings_mod, ONLY: msg, fmsg, maxlen, cat, dispTable, int2str, bool2str, strStack, strParse
[4389]	6	USE readTracFiles_mod, ONLY: trac_type, readTracersFiles, tracers, setGeneration, itZonIso, nzone, tran0, isoZone, &
	7	delPhase, niso, getKey, isot_type, readIsotopesFile, isotope, maxTableWidth, iqIsoPha, nphas, ixIso, isoPhas, &
[4487]	8	addPhase, iH2O, addKey, isoSelect, testTracersFiles, isoKeys, indexUpdate, isoCheck, nbIso, ntiso, isoName
[4063]	9	IMPLICIT NONE
[1114]	10
[4063]	11	PRIVATE
[1114]	12
[4063]	13	!=== FOR TRACERS:
[4325]	14	PUBLIC :: init_infotrac !--- Initialization of the tracers
[4389]	15	PUBLIC :: tracers, type_trac !--- Full tracers database, tracers type keyword
[4063]	16	PUBLIC :: nqtot, nbtr, nqo, nqCO2, nqtottr !--- Main dimensions
[4124]	17	PUBLIC :: conv_flg, pbl_flg !--- Convection & boundary layer activation keys
[3923]	18
[4063]	19	!=== FOR ISOTOPES: General
[4325]	20	PUBLIC :: isot_type, nbIso !--- Derived type, full isotopes families database + nb of families
[4063]	21	PUBLIC :: isoSelect, ixIso !--- Isotopes family selection tool + selected family index
	22	!=== FOR ISOTOPES: Specific to water
[4325]	23	PUBLIC :: iH2O !--- H2O isotopes class index
[4120]	24	PUBLIC :: min_qParent, min_qMass, min_ratio !--- Min. values for various isotopic quantities
[4063]	25	!=== FOR ISOTOPES: Depending on the selected isotopes family
	26	PUBLIC :: isotope, isoKeys !--- Selected isotopes database + associated keys (cf. getKey)
	27	PUBLIC :: isoName, isoZone, isoPhas !--- Isotopes and tagging zones names, phases
	28	PUBLIC :: niso, nzone, nphas, ntiso !--- " " numbers + isotopes & tagging tracers number
[4143]	29	PUBLIC :: itZonIso !--- idx "it" (in "isoName(1:niso)") = function(tagging idx, isotope idx)
	30	PUBLIC :: iqIsoPha !--- idx "iq" (in "qx") = function(isotope idx, phase idx) + aliases
[4063]	31	PUBLIC :: isoCheck !--- Run isotopes checking routines
	32	!=== FOR BOTH TRACERS AND ISOTOPES
	33	PUBLIC :: getKey !--- Get a key from "tracers" or "isotope"
[3870]	34
[4063]	35	!=== CONVENTIONS FOR TRACERS NUMBERS:
	36	! \|--------------------+-----------------------+-----------------+---------------+----------------------------\|
	37	! \| water in different \| water tagging \| water isotopes \| other tracers \| additional tracers moments \|
[4523]	38	! \| phases: H2O_[glsb] \| isotopes \| \| \| for higher order schemes \|
[4063]	39	! \|--------------------+-----------------------+-----------------+---------------+----------------------------\|
	40	! \| \| \| \| \| \|
	41	! \|<-- nqo -->\|<-- nqoniso nzone -->\|<-- nqo*niso -->\|<-- nbtr -->\|<-- (nmom) -->\|
	42	! \| \| \| \|
	43	! \| \|<-- nqoniso(nzone+1) = nqo*ntiso -->\|<-- nqtottr = nbtr + nmom -->\|
	44	! \| = nqtot - nqo*(ntiso+1) \|
	45	! \| \|
	46	! \|<-- nqtrue = nbtr + nqo*(ntiso+1) -->\| \|
	47	! \| \|
	48	! \|<-- nqtot = nqtrue + nmom -->\|
	49	! \| \|
	50	! \|-----------------------------------------------------------------------------------------------------------\|
	51	! NOTES FOR THIS TABLE:
	52	! * Used "niso", "nzone" and "ntiso" are components of "isotopes(ip)" for water (isotopes(ip)%parent == 'H2O'),
	53	! since water is so far the sole tracers family, except passive CO2, removed from the main tracers table.
	54	! * For water, "nqo" is equal to the more general field "isotopes(ip)%nphas".
	55	! * "niso", "nzone", "ntiso", "nphas" are defined for other isotopic tracers families, if any.
	56	!
	57	!=== DERIVED TYPE EMBEDDING MOST OF THE TRACERS-RELATED QUANTITIES (LENGTH: nqtot)
	58	! Each entry is accessible using "%" sign.
	59	! \|-------------+------------------------------------------------------+-------------+------------------------+
	60	! \| entry \| Meaning \| Former name \| Possible values \|
	61	! \|-------------+------------------------------------------------------+-------------+------------------------+
	62	! \| name \| Name (short) \| tname \| \|
	63	! \| gen0Name \| Name of the 1st generation ancestor \| / \| \|
	64	! \| parent \| Name of the parent \| / \| \|
	65	! \| longName \| Long name (with adv. scheme suffix) for outputs \| ttext \| \|
	66	! \| type \| Type (so far: tracer or tag) \| / \| tracer,tag \|
[4523]	67	! \| phase \| Phases list ("g"as / "l"iquid / "s"olid / "b"lowing) \| / \| [g][l][s][b] \|
[4063]	68	! \| component \| Name(s) of the merged/cumulated section(s) \| / \| coma-separated names \|
	69	! \| iGeneration \| Generation (>=1) \| / \| \|
	70	! \| iqParent \| Index of the parent tracer \| iqpere \| 1:nqtot \|
	71	! \| iqDescen \| Indexes of the childs (all generations) \| iqfils \| 1:nqtot \|
	72	! \| nqDescen \| Number of the descendants (all generations) \| nqdesc \| 1:nqtot \|
[4325]	73	! \| nqChildren \| Number of childs (1st generation only) \| nqfils \| 1:nqtot \|
[4358]	74	! \| keys \| key/val pairs accessible with "getKey" routine \| / \| \|
	75	! \| iadv \| Advection scheme number \| iadv \| 1,2,10-20(exc.15,19),30\|
	76	! \| isAdvected \| advected tracers flag (.TRUE. if iadv >= 0) \| / \| nqtrue .TRUE. values \|
	77	! \| isInPhysics \| tracers not extracted from the main table in physics \| / \| nqtottr .TRUE. values \|
[4063]	78	! \| iso_iGroup \| Isotopes group index in isotopes(:) \| / \| 1:nbIso \|
	79	! \| iso_iName \| Isotope name index in isotopes(iso_iGroup)%trac(:) \| iso_indnum \| 1:niso \|
	80	! \| iso_iZone \| Isotope zone index in isotopes(iso_iGroup)%zone(:) \| zone_num \| 1:nzone \|
	81	! \| iso_iPhas \| Isotope phase index in isotopes(iso_iGroup)%phas(:) \| phase_num \| 1:nphas \|
	82	! +-------------+------------------------------------------------------+-------------+------------------------+
	83	!
	84	!=== DERIVED TYPE EMBEDDING MOST OF THE ISOTOPES-RELATED QUANTITIES (LENGTH: nbIso, NUMBER OF ISOTOPES FAMILIES)
	85	! Each entry is accessible using "%" sign.
[4120]	86	! \|-----------------+--------------------------------------------------+--------------------+-----------------+
	87	! \| entry \| length \| Meaning \| Former name \| Possible values \|
	88	! \|-----------------+--------------------------------------------------+--------------------+-----------------+
	89	! \| parent \| Parent tracer (isotopes family name) \| \| \|
	90	! \| keys \| niso \| Isotopes keys/values pairs list + number \| \| \|
	91	! \| trac \| ntiso \| Isotopes + tagging tracers list + number \| / \| ntraciso \| \|
	92	! \| zone \| nzone \| Geographic tagging zones list + number \| / \| ntraceurs_zone \| \|
[4523]	93	! \| phase \| nphas \| Phases list + number \| \|[g][l][s][b],1:4 \|
[4143]	94	! \| iqIsoPha \| Index in "qx" = f(name(1:ntiso)),phas) \| iqiso \| 1:nqtot \|
[4120]	95	! \| itZonIso \| Index in "trac(1:ntiso)"= f(zone, name(1:niso)) \| index_trac \| 1:ntiso \|
	96	! +-----------------+--------------------------------------------------+--------------------+-----------------+
[1114]	97
[4063]	98	REAL, PARAMETER :: min_qParent = 1.e-30, min_qMass = 1.e-18, min_ratio = 1.e-16 ! MVals et CRisi
[3870]	99
[4063]	100	!=== DIMENSIONS OF THE TRACERS TABLES AND OTHER SCALAR VARIABLES
[4325]	101	INTEGER, SAVE :: nqtot, & !--- Tracers nb in dynamics (incl. higher moments + H2O)
	102	nbtr, & !--- Tracers nb in physics (excl. higher moments + H2O)
	103	nqo, & !--- Number of water phases
	104	nqtottr, & !--- Number of tracers passed to phytrac (TO BE DELETED ?)
	105	nqCO2 !--- Number of tracers of CO2 (ThL)
[4389]	106	CHARACTER(LEN=maxlen), SAVE :: type_trac !--- Keyword for tracers type
[4063]	107
[4325]	108	!=== VARIABLES FOR INCA
	109	INTEGER, SAVE, ALLOCATABLE :: conv_flg(:), & !--- Convection activation ; needed for INCA (nbtr)
	110	pbl_flg(:) !--- Boundary layer activation ; needed for INCA (nbtr)
[4063]	111
[1114]	112	CONTAINS
	113
[4325]	114	SUBROUTINE init_infotrac
[4358]	115	USE control_mod, ONLY: planet_type
[1565]	116	#ifdef REPROBUS
[4063]	117	USE CHEM_REP, ONLY: Init_chem_rep_trac
[1565]	118	#endif
[4063]	119	IMPLICIT NONE
	120	!==============================================================================================================================
[1114]	121	!
	122	! Auteur: P. Le Van /L. Fairhead/F.Hourdin
	123	! -------
	124	!
[4063]	125	! Modifications:
	126	! --------------
	127	! 05/94: F.Forget Modif special traceur
	128	! 02/02: M-A Filiberti Lecture de traceur.def
[4120]	129	! 01/22: D. Cugnet Nouveaux tracer.def et tracer_*.def + encapsulation (types trac_type et isot_type)
[4063]	130	!
[1114]	131	! Objet:
	132	! ------
	133	! GCM LMD nouvelle grille
	134	!
[4063]	135	!==============================================================================================================================
[1114]	136	! ... modification de l'integration de q ( 26/04/94 ) ....
[4063]	137	!------------------------------------------------------------------------------------------------------------------------------
	138	! Declarations:
	139	INCLUDE "dimensions.h"
	140	INCLUDE "iniprint.h"
[1114]	141
[4063]	142	!------------------------------------------------------------------------------------------------------------------------------
[1114]	143	! Local variables
[4064]	144	INTEGER, ALLOCATABLE :: hadv(:), vadv(:) !--- Horizontal/vertical transport scheme number
	145	#ifdef INCA
	146	INTEGER, ALLOCATABLE :: had (:), hadv_inca(:), conv_flg_inca(:), &!--- Variables specific to INCA
	147	vad (:), vadv_inca(:), pbl_flg_inca(:)
[4063]	148	CHARACTER(LEN=8), ALLOCATABLE :: solsym_inca(:) !--- Tracers names for INCA
[4064]	149	INTEGER :: nqINCA
	150	#endif
[4063]	151	CHARACTER(LEN=2) :: suff(9) !--- Suffixes for schemes of order 3 or 4 (Prather)
	152	CHARACTER(LEN=3) :: descrq(30) !--- Advection scheme description tags
[4389]	153	CHARACTER(LEN=maxlen) :: msg1, texp, ttp !--- Strings for messages and expanded tracers type
[4063]	154	INTEGER :: fType !--- Tracers description file type ; 0: none
	155	!--- 1/2/3: "traceur.def"/"tracer.def"/"tracer_*.def"
	156	INTEGER :: nqtrue !--- Tracers nb from tracer.def (no higher order moments)
	157	INTEGER :: iad !--- Advection scheme number
[4325]	158	INTEGER :: ic, iq, jq, it, nt, im, nm, iz, k !--- Indexes and temporary variables
[4389]	159	LOGICAL :: lerr, ll
[4063]	160	CHARACTER(LEN=1) :: p
	161	TYPE(trac_type), ALLOCATABLE, TARGET :: ttr(:)
[4064]	162	TYPE(trac_type), POINTER :: t1, t(:)
	163	INTEGER :: ierr
[2566]	164
[4325]	165	CHARACTER(LEN=*), PARAMETER :: modname="init_infotrac"
[4063]	166	!------------------------------------------------------------------------------------------------------------------------------
	167	! Initialization :
	168	!------------------------------------------------------------------------------------------------------------------------------
	169	suff = ['x ','y ','z ','xx','xy','xz','yy','yz','zz']
[4487]	170	descrq( 1:30) = ' '
[4063]	171	descrq( 1: 2) = ['LMV','BAK']
	172	descrq(10:20) = ['VL1','VLP','FH1','FH2','VLH',' ','PPM','PPS','PPP',' ','SLP']
	173	descrq(30) = 'PRA'
	174
[4067]	175	CALL msg('type_trac = "'//TRIM(type_trac)//'"', modname)
[3870]	176
[4389]	177	!--- MESSAGE ABOUT THE CHOSEN CONFIGURATION
	178	msg1 = 'For type_trac = "'//TRIM(type_trac)//'":'
	179	SELECT CASE(type_trac)
	180	CASE('inca'); CALL msg(TRIM(msg1)//' coupling with INCA chemistry model', modname)
	181	CASE('inco'); CALL msg(TRIM(msg1)//' coupling jointly with INCA and CO2 cycle', modname)
	182	CASE('repr'); CALL msg(TRIM(msg1)//' coupling with REPROBUS chemistry model', modname)
	183	CASE('co2i'); CALL msg(TRIM(msg1)//' you have chosen to run with CO2 cycle', modname)
	184	CASE('coag'); CALL msg(TRIM(msg1)//' tracers are treated for COAGULATION tests', modname)
	185	CASE('lmdz'); CALL msg(TRIM(msg1)//' tracers are treated in LMDZ only', modname)
	186	CASE DEFAULT; CALL abort_gcm(modname,'type_trac='//TRIM(type_trac)//' not possible yet.',1)
	187	END SELECT
[1454]	188
[4389]	189	!--- COHERENCE TEST BETWEEN "type_trac" AND PREPROCESSING KEYS
	190	SELECT CASE(type_trac)
	191	CASE('inca', 'inco')
[1569]	192	#ifndef INCA
[4389]	193	CALL abort_gcm(modname, 'You must add cpp key INCA and compile with INCA code', 1)
[1569]	194	#endif
[4389]	195	CASE('repr')
[1569]	196	#ifndef REPROBUS
[4389]	197	CALL abort_gcm(modname, 'You must add cpp key REPROBUS and compile with REPROBUS code', 1)
[1569]	198	#endif
[4389]	199	CASE('coag')
[2690]	200	#ifndef CPP_StratAer
[4389]	201	CALL abort_gcm(modname, 'You must add cpp key StratAer and compile with StratAer code', 1)
[2690]	202	#endif
[4389]	203	END SELECT
[1279]	204
[4389]	205	nqCO2 = COUNT( [type_trac == 'inco', type_trac == 'co2i'] )
[4063]	206
	207	!==============================================================================================================================
	208	! 1) Get the numbers of: true (first order only) tracers "nqtrue", water tracers "nqo" (vapor/liquid/solid)
	209	!==============================================================================================================================
[4389]	210	texp = type_trac !=== EXPANDED VERSION OF "type_trac", WITH "\|" SEPARATOR
	211	IF(texp == 'inco') texp = 'co2i\|inca'
	212	IF(texp /= 'lmdz') texp = 'lmdz\|'//TRIM(texp)
	213
	214	!=== DETERMINE THE TYPE OF THE INPUT TRACERS DESCRIPTION FILE
	215	IF(testTracersFiles(modname, texp, fType, .TRUE.)) CALL abort_gcm(modname, 'problem with tracers file(s)',1)
	216	ttp = type_trac; IF(fType /= 1) ttp = texp
	217
	218	IF(readTracersFiles(ttp, type_trac == 'repr')) CALL abort_gcm(modname, 'problem with tracers file(s)',1)
[4301]	219	!---------------------------------------------------------------------------------------------------------------------------
[4063]	220	IF(fType == 0) CALL abort_gcm(modname, 'Missing tracers file: "traceur.def", "tracer.def" or "tracer_<keyword>.def file.',1)
	221	!---------------------------------------------------------------------------------------------------------------------------
[4389]	222	IF(fType == 1 .AND. ANY(['inca','inco']==type_trac)) THEN !=== FOUND OLD STYLE INCA "traceur.def"
[4063]	223	!---------------------------------------------------------------------------------------------------------------------------
	224	#ifdef INCA
[4301]	225	nqo = SIZE(tracers) - nqCO2
[4063]	226	CALL Init_chem_inca_trac(nqINCA) !--- Get nqINCA from INCA
	227	nbtr = nqINCA + nqCO2 !--- Number of tracers passed to phytrac
	228	nqtrue = nbtr + nqo !--- Total number of "true" tracers
[4523]	229	IF(ALL([2,3,4,5] /= nqo)) CALL abort_gcm(modname, 'Only 2, 3, 4 , 5 water phases allowed ; found nqo='//TRIM(int2str(nqo)), 1)
[4143]	230	ALLOCATE(hadv(nqtrue), hadv_inca(nqINCA), conv_flg_inca(nqINCA), solsym_inca(nqINCA))
	231	ALLOCATE(vadv(nqtrue), vadv_inca(nqINCA), pbl_flg_inca(nqINCA))
[4127]	232	CALL init_transport(solsym_inca, conv_flg_inca, pbl_flg_inca, hadv_inca, vadv_inca)
[4063]	233	ALLOCATE(ttr(nqtrue))
[4389]	234	ttr(1:nqo+nqCO2) = tracers
	235	ttr(1 : nqo )%component = 'lmdz'
	236	ttr(1+nqo:nqCO2+nqo )%component = 'co2i'
	237	ttr(1+nqo+nqCO2:nqtrue)%component = 'inca'
	238	ttr(1+nqo :nqtrue)%name = [('CO2 ', k=1, nqCO2), solsym_inca]
	239	ttr(1+nqo+nqCO2:nqtrue)%parent = tran0
	240	ttr(1+nqo+nqCO2:nqtrue)%phase = 'g'
[4189]	241	lerr = getKey('hadv', had, ky=tracers(:)%keys)
	242	lerr = getKey('vadv', vad, ky=tracers(:)%keys)
[4389]	243	hadv(1:nqo+nqCO2) = had(:); hadv(1+nqo+nqCO2:nqtrue) = hadv_inca
	244	vadv(1:nqo+nqCO2) = vad(:); vadv(1+nqo+nqCO2:nqtrue) = vadv_inca
[4063]	245	CALL MOVE_ALLOC(FROM=ttr, TO=tracers)
[4429]	246	DO iq = 1, nqtrue
	247	t1 => tracers(iq)
[4431]	248	CALL addKey('name', t1%name, t1%keys)
	249	CALL addKey('component', t1%component, t1%keys)
	250	CALL addKey('parent', t1%parent, t1%keys)
	251	CALL addKey('phase', t1%phase, t1%keys)
[4429]	252	END DO
[4389]	253	IF(setGeneration(tracers)) CALL abort_gcm(modname,'See above',1) !- SET FIELDS %iGeneration, %gen0Name
[4120]	254	DEALLOCATE(had, hadv_inca, vad, vadv_inca, conv_flg_inca, pbl_flg_inca, solsym_inca)
[4063]	255	#endif
	256	!---------------------------------------------------------------------------------------------------------------------------
[4301]	257	ELSE !=== OTHER CASES (OLD OR NEW FORMAT, NO INCA MODULE)
[4063]	258	!---------------------------------------------------------------------------------------------------------------------------
[4120]	259	nqo = COUNT(delPhase(tracers(:)%name) == 'H2O' &
	260	.AND. tracers(:)%component == 'lmdz') !--- Number of water phases
[4063]	261	nqtrue = SIZE(tracers) !--- Total number of "true" tracers
[4120]	262	nbtr = nqtrue-COUNT(delPhase(tracers(:)%gen0Name) == 'H2O' &
	263	.AND. tracers(:)%component == 'lmdz') !--- Number of tracers passed to phytrac
[4172]	264	#ifdef INCA
	265	nqINCA = COUNT(tracers(:)%component == 'inca')
	266	#endif
[4063]	267	lerr = getKey('hadv', hadv, ky=tracers(:)%keys)
	268	lerr = getKey('vadv', vadv, ky=tracers(:)%keys)
	269	!---------------------------------------------------------------------------------------------------------------------------
	270	END IF
	271	!---------------------------------------------------------------------------------------------------------------------------
[2262]	272
[4429]	273	#ifdef REPROBUS
[4063]	274	!--- Transfert the number of tracers to Reprobus
	275	CALL Init_chem_rep_trac(nbtr, nqo, tracers(:)%name)
[4429]	276
[3945]	277	#endif
[4063]	278	!==============================================================================================================================
	279	! 2) Calculate nqtot, number of tracers needed (greater if advection schemes 20 or 30 have been chosen).
	280	!==============================================================================================================================
	281	DO iq = 1, nqtrue
	282	IF( hadv(iq)<20 .OR. (ANY(hadv(iq)==[20,30]) .AND. hadv(iq)==vadv(iq)) ) CYCLE
	283	WRITE(msg1,'("The choice hadv=",i0,", vadv=",i0,a)')hadv(iq),vadv(iq),' for "'//TRIM(tracers(iq)%name)//'" is not available'
	284	CALL abort_gcm(modname, TRIM(msg1), 1)
	285	END DO
	286	nqtot = COUNT( hadv< 20 .AND. vadv< 20 ) & !--- No additional tracer
	287	+ 4*COUNT( hadv==20 .AND. vadv==20 ) & !--- 3 additional tracers
	288	+ 10*COUNT( hadv==30 .AND. vadv==30 ) !--- 9 additional tracers
	289
	290	!--- More tracers due to the choice of advection scheme => assign total number of tracers
	291	IF( nqtot /= nqtrue ) THEN
	292	CALL msg('The choice of advection scheme for one or more tracers makes it necessary to add tracers')
	293	CALL msg('The number of true tracers is '//TRIM(int2str(nqtrue)))
	294	CALL msg('The total number of tracers needed is '//TRIM(int2str(nqtot)))
	295	END IF
	296
	297	!==============================================================================================================================
	298	! 3) Determine the advection scheme choice for water and tracers "iadv" and the fields long name, isAdvected.
	299	! iadv = 1 "LMDZ-specific humidity transport" (for H2O vapour) LMV
	300	! iadv = 2 backward (for H2O liquid) BAK
	301	! iadv = 14 Van-Leer + specific humidity, modified by Francis Codron VLH
	302	! iadv = 10 Van-Leer (chosen for vapour and liquid water) VL1
	303	! iadv = 11 Van-Leer for hadv and PPM version (Monotonic) for vadv VLP
	304	! iadv = 12 Frederic Hourdin I FH1
	305	! iadv = 13 Frederic Hourdin II FH2
	306	! iadv = 16 Monotonic PPM (Collela & Woodward 1984) PPM
	307	! iadv = 17 Semi-monotonic PPM (overshoots allowed) PPS
	308	! iadv = 18 Definite positive PPM (overshoots and undershoots allowed) PPP
	309	! iadv = 20 Slopes SLP
	310	! iadv = 30 Prather PRA
[3945]	311	!
[4063]	312	! In array q(ij,l,iq) : iq = 1/2[/3] for vapour/liquid[/ice] water
	313	! And optionaly: iq = 3[4],nqtot for other tracers
	314	!==============================================================================================================================
	315	ALLOCATE(ttr(nqtot))
	316	jq = nqtrue+1; tracers(:)%iadv = -1
	317	DO iq = 1, nqtrue
	318	t1 => tracers(iq)
[3945]	319
[4063]	320	!--- VERIFY THE CHOICE OF ADVECTION SCHEME
	321	iad = -1
	322	IF(hadv(iq) == vadv(iq) ) iad = hadv(iq)
	323	IF(hadv(iq)==10 .AND. vadv(iq)==16) iad = 11
	324	WRITE(msg1,'("Bad choice of advection scheme for ",a,": hadv = ",i0,", vadv = ",i0)')TRIM(t1%name), hadv(iq), vadv(iq)
	325	IF(iad == -1) CALL abort_gcm(modname, msg1, 1)
[3945]	326
[4071]	327	!--- SET FIELDS %longName, %iadv, %isAdvected, %isInPhysics
[4063]	328	t1%longName = t1%name; IF(iad > 0) t1%longName=TRIM(t1%name)//descrq(iad)
	329	t1%iadv = iad
[4071]	330	t1%isAdvected = iad >= 0
[4120]	331	t1%isInPhysics= delPhase(t1%gen0Name) /= 'H2O' &
	332	.OR. t1%component /= 'lmdz' !=== OTHER EXCEPTIONS TO BE ADDED: CO2i, SURSATURATED WATER CLOUD...
[4063]	333	ttr(iq) = t1
[3945]	334
[4063]	335	!--- DEFINE THE HIGHER ORDER TRACERS, IF ANY
	336	nm = 0
[4120]	337	IF(iad == 20) nm = 3 !--- 2nd order scheme
	338	IF(iad == 30) nm = 9 !--- 3rd order scheme
	339	IF(nm == 0) CYCLE !--- No higher moments
[4063]	340	ttr(jq+1:jq+nm) = t1
	341	ttr(jq+1:jq+nm)%name = [(TRIM(t1%name) //'-'//TRIM(suff(im)), im=1, nm) ]
	342	ttr(jq+1:jq+nm)%parent = [(TRIM(t1%parent) //'-'//TRIM(suff(im)), im=1, nm) ]
	343	ttr(jq+1:jq+nm)%longName = [(TRIM(t1%longName)//'-'//TRIM(suff(im)), im=1, nm) ]
	344	ttr(jq+1:jq+nm)%iadv = [(-iad, im=1, nm) ]
	345	ttr(jq+1:jq+nm)%isAdvected = [(.FALSE., im=1, nm) ]
	346	jq = jq + nm
	347	END DO
	348	DEALLOCATE(hadv, vadv)
	349	CALL MOVE_ALLOC(FROM=ttr, TO=tracers)
[3945]	350
[4325]	351	!--- SET FIELDS %iqParent, %nqChildren, %iGeneration, %iqDescen, %nqDescen
[4063]	352	CALL indexUpdate(tracers)
[1114]	353
[4063]	354	!=== TEST ADVECTION SCHEME
	355	DO iq=1,nqtot ; t1 => tracers(iq); iad = t1%iadv
[1114]	356
[4063]	357	!--- ONLY TESTED VALUES FOR TRACERS FOR NOW: iadv = 14, 10 (and 0 for non-transported tracers)
	358	IF(ALL([10,14,0] /= iad)) &
	359	CALL abort_gcm(modname, 'Not tested for iadv='//TRIM(int2str(iad))//' ; 10 or 14 only are allowed !', 1)
[1114]	360
[4063]	361	!--- ONLY TESTED VALUES FOR PARENTS HAVING CHILDS FOR NOW: iadv = 14, 10 (PARENTS: TRACERS OF GENERATION 1)
	362	IF(ALL([10,14] /= iad) .AND. t1%iGeneration == 1 .AND. ANY(tracers(:)%iGeneration > 1)) &
	363	CALL abort_gcm(modname, 'iadv='//TRIM(int2str(iad))//' not implemented for parents ; 10 or 14 only are allowed !', 1)
[1114]	364
[4063]	365	!--- ONLY TESTED VALUES FOR CHILDS FOR NOW: iadv = 10 (CHILDS: TRACERS OF GENERATION GREATER THAN 1)
	366	IF(fmsg('WARNING ! iadv='//TRIM(int2str(iad))//' not implemented for childs. Setting iadv=10 for "'//TRIM(t1%name)//'"',&
	367	modname, iad /= 10 .AND. t1%iGeneration > 1)) t1%iadv = 10
[1454]	368
[4063]	369	!--- ONLY VALID SCHEME NUMBER FOR WATER VAPOUR: iadv = 14
[4120]	370	ll = t1%name /= addPhase('H2O','g')
[4063]	371	IF(fmsg('WARNING ! iadv=14 is valid for water vapour only. Setting iadv=10 for "'//TRIM(t1%name)//'".', &
	372	modname, iad == 14 .AND. ll)) t1%iadv = 10
	373	END DO
[1114]	374
[4325]	375	!=== READ PHYSICAL PARAMETERS FOR ISOTOPES ; DONE HERE BECAUSE dynetat0 AND iniacademic NEED "tnat" AND "alpha_ideal"
	376	niso = 0; nzone = 0; nphas = nqo; ntiso = 0; isoCheck = .FALSE.
	377	IF(readIsotopesFile()) CALL abort_gcm(modname, 'Problem when reading isotopes parameters', 1)
[2270]	378
[4143]	379	!--- Convection / boundary layer activation for all tracers
	380	ALLOCATE(conv_flg(nbtr)); conv_flg(1:nbtr) = 1
	381	ALLOCATE( pbl_flg(nbtr)); pbl_flg(1:nbtr) = 1
	382
[4120]	383	!--- Note: nqtottr can differ from nbtr when nmom/=0
[4325]	384	nqtottr = nqtot - COUNT(delPhase(tracers%gen0Name) == 'H2O' .AND. tracers%component == 'lmdz')
[4301]	385	IF(COUNT(tracers%iso_iName == 0) - COUNT(delPhase(tracers%name) == 'H2O' .AND. tracers%component == 'lmdz') /= nqtottr) &
[4325]	386	CALL abort_gcm(modname, 'pb dans le calcul de nqtottr', 1)
[3923]	387
[4120]	388	!=== DISPLAY THE RESULTS
[4403]	389	CALL msg('nqo = '//TRIM(int2str(nqo)), modname)
	390	CALL msg('nbtr = '//TRIM(int2str(nbtr)), modname)
	391	CALL msg('nqtrue = '//TRIM(int2str(nqtrue)), modname)
	392	CALL msg('nqtot = '//TRIM(int2str(nqtot)), modname)
	393	CALL msg('niso = '//TRIM(int2str(niso)), modname)
	394	CALL msg('ntiso = '//TRIM(int2str(ntiso)), modname)
[4120]	395	#ifdef INCA
[4403]	396	CALL msg('nqCO2 = '//TRIM(int2str(nqCO2)), modname)
	397	CALL msg('nqINCA = '//TRIM(int2str(nqINCA)), modname)
[4120]	398	#endif
	399	t => tracers
	400	CALL msg('Information stored in infotrac :', modname)
	401	IF(dispTable('isssssssssiiiiiiiii', &
[4325]	402	['iq ', 'name ', 'lName ', 'gen0N ', 'parent', 'type ', 'phase ', 'compon', 'isPhy ', 'isAdv ', &
[4120]	403	'iadv ', 'iGen ', 'iqPar ', 'nqDes ', 'nqChld', 'iGroup', 'iName ', 'iZone ', 'iPhase'], &
[4325]	404	cat(t%name, t%longName, t%gen0Name, t%parent, t%type, t%phase, t%component, bool2str(t%isInPhysics), &
	405	bool2str(t%isAdvected)), &
	406	cat([(iq, iq=1, nqtot)], t%iadv, t%iGeneration, t%iqParent, t%nqDescen, t%nqChildren, t%iso_iGroup, &
[4193]	407	t%iso_iName, t%iso_iZone, t%iso_iPhase), nColMax=maxTableWidth, nHead=2, sub=modname)) &
[4120]	408	CALL abort_gcm(modname, "problem with the tracers table content", 1)
	409	IF(niso > 0) THEN
	410	CALL msg('Where, for isotopes family "'//TRIM(isotope%parent)//'":', modname)
[4325]	411	CALL msg(' isoKeys%name = '//strStack(isoKeys%name), modname)
[4120]	412	CALL msg(' isoName = '//strStack(isoName), modname)
	413	CALL msg(' isoZone = '//strStack(isoZone), modname)
	414	CALL msg(' isoPhas = '//TRIM(isoPhas), modname)
	415	ELSE
	416	CALL msg('No isotopes identified.', modname)
	417	END IF
[4063]	418	CALL msg('end', modname)
[3923]	419
[4325]	420	END SUBROUTINE init_infotrac
[3923]	421
[1114]	422	END MODULE infotrac

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

Download in other formats:

Original Format