Context Navigation

phytrac_mod.F90 @ 2236

Last change on this file since 2236 was 2210, checked in by fhourdin, 10 years ago
Control of convective scavenging by .def files Conrôle du lessivage convectif dans les fichiers .def ccntrAA= 1. # fraction d'aerosol en phase condensee dans l'ascendance adiab ccntrENV=1. # fraction d'aerosol en phase condensee dans les melange (Emanuel) coefcoli=1.e-3 # efficacité de collision
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
File size: 36.5 KB

Rev	Line
[1877]	1	!$Id$
[1813]	2	MODULE phytrac_mod
	3	!=================================================================================
	4	! Interface between the LMDZ physical package and tracer computation.
	5	! Chemistry modules (INCA, Reprobus or the more specific traclmdz routine)
	6	! are called from phytrac.
	7	!
	8	!======================================================================
	9	! Auteur(s) FH
	10	! Objet: Moniteur general des tendances traceurs
	11	!
	12	! iflag_vdf_trac : Options for activating transport by vertical diffusion :
	13	! 1. notmal
	14	! 0. emission is injected in the first layer only, without diffusion
	15	! -1 no emission & no diffusion
	16	! Modification 2013/07/22 : transformed into a module to pass tendencies to
	17	! physics outputs. Additional keys for controling activation of sub processes.
	18	! Modification R. Pilon 10 octobre 2012 large scale scavenging incloud_scav + bc_scav
	19	! Modification R. Pilon 01 janvier 2012 transport+scavenging KE scheme : cvltr
	20	!=================================================================================
	21
	22	!
	23	! Tracer tendencies, for outputs
	24	!-------------------------------
	25	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_cl ! Td couche limite/traceur
	26	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_dec !RomP
	27	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_cv ! Td convection/traceur
	28	! RomP >>>
	29	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_insc
	30	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_bcscav
	31	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_evapls
	32	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_ls
	33	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_trsp
	34	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_sscav
	35	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_sat
	36	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_uscav
	37	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: qPr,qDi ! concentration tra dans pluie,air descente insaturee
	38	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: qPa,qMel
	39	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: qTrdi,dtrcvMA ! conc traceur descente air insaturee et td convective MA
	40	! RomP <<<
	41	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_th ! Td thermique
	42	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_lessi_impa ! Td du lessivage par impaction
	43	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_lessi_nucl ! Td du lessivage par nucleation
	44	REAL,DIMENSION(:,:),ALLOCATABLE,SAVE :: qPrls !jyg: concentration tra dans pluie LS a la surf.
	45	REAL,DIMENSION(:,:),ALLOCATABLE,SAVE :: d_tr_dry ! Td depot sec/traceur (1st layer),ALLOCATABLE,SAVE jyg
	46	REAL,DIMENSION(:,:),ALLOCATABLE,SAVE :: flux_tr_dry ! depot sec/traceur (surface),ALLOCATABLE,SAVE jyg
	47
	48	!$OMP THREADPRIVATE(qPa,qMel,qTrdi,dtrcvMA,d_tr_th,d_tr_lessi_impa,d_tr_lessi_nucl)
	49	!$OMP THREADPRIVATE(d_tr_trsp,d_tr_sscav,d_tr_sat,d_tr_uscav,qPr,qDi)
	50	!$OMP THREADPRIVATE(d_tr_insc,d_tr_bcscav,d_tr_evapls,d_tr_ls,qPrls)
	51	!$OMP THREADPRIVATE(d_tr,d_tr_cl,d_tr_dry,flux_tr_dry,d_tr_dec,d_tr_cv)
	52
	53
	54	CONTAINS
	55
[2146]	56	SUBROUTINE phytrac( &
	57	nstep, julien, gmtime, debutphy, &
	58	lafin, pdtphys, u, v, t_seri, &
	59	paprs, pplay, pmfu, pmfd, &
	60	pen_u, pde_u, pen_d, pde_d, &
	61	cdragh, coefh, fm_therm, entr_therm, &
	62	yu1, yv1, ftsol, pctsrf, &
	63	ustar, u10m, v10m, &
	64	wstar, ale_bl, ale_wake, &
	65	xlat, xlon, &
	66	frac_impa,frac_nucl,beta_fisrt,beta_v1, &
	67	presnivs, pphis, pphi, albsol, &
	68	sh, rh, cldfra, rneb, &
	69	diafra, cldliq, itop_con, ibas_con, &
	70	pmflxr, pmflxs, prfl, psfl, &
	71	da, phi, mp, upwd, &
	72	phi2, d1a, dam, sij, wght_cvfd, & ! RomP +RL
	73	wdtrainA, wdtrainM, sigd, clw, elij, & ! RomP
	74	evap, ep, epmlmMm, eplaMm, & ! RomP
	75	dnwd, aerosol_couple, flxmass_w, &
	76	tau_aero, piz_aero, cg_aero, ccm, &
	77	rfname, &
	78	d_tr_dyn, & ! RomP
	79	tr_seri)
	80	!
	81	!======================================================================
	82	! Auteur(s) FH
	83	! Objet: Moniteur general des tendances traceurs
	84	! Modification R. Pilon 01 janvier 2012 transport+scavenging KE scheme : cvltr
	85	! Modification R. Pilon 10 octobre 2012 large scale scavenging incloud_scav + bc_scav
	86	!======================================================================
[1813]	87
[2146]	88	USE ioipsl
	89	USE phys_cal_mod, only : hour
	90	USE dimphy
	91	USE infotrac
	92	USE mod_grid_phy_lmdz
	93	USE mod_phys_lmdz_para
	94	USE comgeomphy
	95	USE iophy
	96	USE traclmdz_mod
	97	USE tracinca_mod
	98	USE tracreprobus_mod
	99	USE control_mod
	100	USE indice_sol_mod
[1813]	101
[2146]	102	USE mod_phys_lmdz_mpi_data, ONLY : is_mpi_root
[1813]	103
[2146]	104	IMPLICIT NONE
[1813]	105
[2146]	106	INCLUDE "YOMCST.h"
	107	INCLUDE "dimensions.h"
	108	INCLUDE "clesphys.h"
	109	INCLUDE "temps.h"
	110	INCLUDE "paramet.h"
	111	INCLUDE "thermcell.h"
	112	INCLUDE "iniprint.h"
	113	!==========================================================================
	114	! -- ARGUMENT DESCRIPTION --
	115	!==========================================================================
[1813]	116
[2146]	117	! Input arguments
	118	!----------------
	119	!Configuration grille,temps:
	120	INTEGER,INTENT(IN) :: nstep ! Appel physique
	121	INTEGER,INTENT(IN) :: julien ! Jour julien
	122	REAL,INTENT(IN) :: gmtime ! Heure courante
	123	REAL,INTENT(IN) :: pdtphys ! Pas d'integration pour la physique (seconde)
	124	LOGICAL,INTENT(IN) :: debutphy ! le flag de l'initialisation de la physique
	125	LOGICAL,INTENT(IN) :: lafin ! le flag de la fin de la physique
[1813]	126
[2146]	127	REAL,DIMENSION(klon),INTENT(IN) :: xlat ! latitudes pour chaque point
	128	REAL,DIMENSION(klon),INTENT(IN) :: xlon ! longitudes pour chaque point
	129	!
	130	!Physique:
	131	!--------
	132	REAL,DIMENSION(klon,klev),INTENT(IN) :: t_seri ! Temperature
	133	REAL,DIMENSION(klon,klev),INTENT(IN) :: u ! variable not used
	134	REAL,DIMENSION(klon,klev),INTENT(IN) :: v ! variable not used
	135	REAL,DIMENSION(klon,klev),INTENT(IN) :: sh ! humidite specifique
	136	REAL,DIMENSION(klon,klev),INTENT(IN) :: rh ! humidite relative
	137	REAL,DIMENSION(klon,klev+1),INTENT(IN) :: paprs ! pression pour chaque inter-couche (en Pa)
	138	REAL,DIMENSION(klon,klev),INTENT(IN) :: pplay ! pression pour le mileu de chaque couche (en Pa)
	139	REAL,DIMENSION(klon,klev),INTENT(IN) :: pphi ! geopotentiel
	140	REAL,DIMENSION(klon),INTENT(IN) :: pphis
	141	REAL,DIMENSION(klev),INTENT(IN) :: presnivs
	142	REAL,DIMENSION(klon,klev),INTENT(IN) :: cldliq ! eau liquide nuageuse
	143	REAL,DIMENSION(klon,klev),INTENT(IN) :: cldfra ! fraction nuageuse (tous les nuages)
	144	REAL,DIMENSION(klon,klev),INTENT(IN) :: diafra ! fraction nuageuse (convection ou stratus artificiels)
	145	REAL,DIMENSION(klon,klev),INTENT(IN) :: rneb ! fraction nuageuse (grande echelle)
	146	!
	147	REAL :: ql_incl ! contenu en eau liquide nuageuse dans le nuage ! ql_incl=oliq/rneb
	148	REAL,DIMENSION(klon,klev),INTENT(IN) :: beta_fisrt ! taux de conversion de l'eau cond (de fisrtilp)
	149	REAL,DIMENSION(klon,klev),INTENT(out) :: beta_v1 ! -- (originale version)
[1813]	150
[2146]	151	!
	152	INTEGER,DIMENSION(klon),INTENT(IN) :: itop_con
	153	INTEGER,DIMENSION(klon),INTENT(IN) :: ibas_con
	154	REAL,DIMENSION(klon),INTENT(IN) :: albsol ! albedo surface
	155	!
	156	!Dynamique
	157	!--------
	158	REAL,DIMENSION(klon,klev,nbtr),INTENT(IN) :: d_tr_dyn
	159	!
	160	!Convection:
	161	!----------
	162	REAL,DIMENSION(klon,klev),INTENT(IN) :: pmfu ! flux de masse dans le panache montant
	163	REAL,DIMENSION(klon,klev),INTENT(IN) :: pmfd ! flux de masse dans le panache descendant
	164	REAL,DIMENSION(klon,klev),INTENT(IN) :: pen_u ! flux entraine dans le panache montant
	165	REAL,DIMENSION(klon,klev),INTENT(IN) :: pde_u ! flux detraine dans le panache montant
	166	REAL,DIMENSION(klon,klev),INTENT(IN) :: pen_d ! flux entraine dans le panache descendant
	167	REAL,DIMENSION(klon,klev),INTENT(IN) :: pde_d ! flux detraine dans le panache descendant
[1813]	168
[2146]	169	!...Tiedke
	170	REAL,DIMENSION(klon,klev+1),INTENT(IN) :: pmflxr, pmflxs ! Flux precipitant de pluie, neige aux interfaces [convection]
	171	REAL,DIMENSION(klon,klev+1),INTENT(IN) :: prfl, psfl ! Flux precipitant de pluie, neige aux interfaces [large-scale]
[1813]	172
[2146]	173	LOGICAL,INTENT(IN) :: aerosol_couple
	174	REAL,DIMENSION(klon,klev),INTENT(IN) :: flxmass_w
	175	REAL,DIMENSION(klon,klev,9,2),INTENT(IN) :: tau_aero
	176	REAL,DIMENSION(klon,klev,9,2),INTENT(IN) :: piz_aero
	177	REAL,DIMENSION(klon,klev,9,2),INTENT(IN) :: cg_aero
	178	CHARACTER(len=4),DIMENSION(9),INTENT(IN) :: rfname
	179	REAL,DIMENSION(klon,klev,2),INTENT(IN) :: ccm
	180	!... K.Emanuel
	181	REAL,DIMENSION(klon,klev),INTENT(IN) :: da
	182	REAL,DIMENSION(klon,klev,klev),INTENT(IN):: phi
	183	! RomP >>>
	184	REAL,DIMENSION(klon,klev),INTENT(IN) :: d1a,dam
	185	REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: phi2
	186	!
	187	REAL,DIMENSION(klon,klev),INTENT(IN) :: wdtrainA
	188	REAL,DIMENSION(klon,klev),INTENT(IN) :: wdtrainM
	189	REAL,DIMENSION(klon),INTENT(IN) :: sigd
	190	! ---- RomP flux entraine, detraine et precipitant kerry Emanuel
	191	REAL,DIMENSION(klon,klev),INTENT(IN) :: evap
	192	REAL,DIMENSION(klon,klev),INTENT(IN) :: ep
	193	REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: sij
	194	REAL,DIMENSION(klon,klev),INTENT(IN) :: wght_cvfd !RL
	195	REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: elij
	196	REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: epmlmMm
	197	REAL,DIMENSION(klon,klev),INTENT(IN) :: eplaMm
	198	REAL,DIMENSION(klon,klev),INTENT(IN) :: clw
	199	! RomP <<<
[1813]	200
[2146]	201	!
	202	REAL,DIMENSION(klon,klev),INTENT(IN) :: mp
	203	REAL,DIMENSION(klon,klev),INTENT(IN) :: upwd ! saturated updraft mass flux
	204	REAL,DIMENSION(klon,klev),INTENT(IN) :: dnwd ! saturated downdraft mass flux
	205	!
	206	!Thermiques:
	207	!----------
	208	REAL,DIMENSION(klon,klev+1),INTENT(IN) :: fm_therm
	209	REAL,DIMENSION(klon,klev),INTENT(IN) :: entr_therm
	210	!
	211	!Couche limite:
	212	!--------------
	213	!
	214	!
	215	REAL,DIMENSION(:),INTENT(IN) :: cdragh ! (klon) coeff drag pour T et Q
	216	REAL,DIMENSION(:,:),INTENT(IN) :: coefh ! (klon,klev) coeff melange CL (m**2/s)
	217	REAL,DIMENSION(:),INTENT(IN) :: ustar,u10m,v10m ! (klon) u* & vent a 10m (m/s)
	218	REAL,DIMENSION(:),INTENT(IN) :: wstar,ale_bl,ale_wake ! (klon) w* and Avail. Lifting Ener.
	219	REAL,DIMENSION(:),INTENT(IN) :: yu1 ! (klon) vents au premier niveau
	220	REAL,DIMENSION(:),INTENT(IN) :: yv1 ! (klon) vents au premier niveau
[1813]	221
[2146]	222	!
	223	!Lessivage:
	224	!----------
	225	!
	226	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: ccntrAA
	227	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: ccntrENV
	228	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: coefcoli
	229	LOGICAL, DIMENSION(:), ALLOCATABLE, SAVE :: flag_cvltr
	230	!$OMP THREADPRIVATE(ccntrAA,ccntrENV,coefcoli,flag_cvltr)
	231	REAL, DIMENSION(klon,klev) :: ccntrAA_3d
	232	REAL, DIMENSION(klon,klev) :: ccntrENV_3d
	233	REAL, DIMENSION(klon,klev) :: coefcoli_3d
	234	!
	235	! pour le ON-LINE
	236	!
	237	REAL,DIMENSION(klon,klev),INTENT(IN) :: frac_impa ! fraction d'aerosols non impactes
	238	REAL,DIMENSION(klon,klev),INTENT(IN) :: frac_nucl ! fraction d'aerosols non nuclees
[1813]	239
[2146]	240	! Arguments necessaires pour les sources et puits de traceur:
	241	REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: ftsol ! Temperature du sol (surf)(Kelvin)
	242	REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: pctsrf ! Pourcentage de sol (nature du sol)
[1813]	243
	244
[2146]	245	! Output argument
	246	!----------------
	247	REAL,DIMENSION(klon,klev,nbtr),INTENT(INOUT) :: tr_seri ! Concentration Traceur [U/KgA]
	248	REAL,DIMENSION(klon,klev) :: sourceBE
	249	!=======================================================================================
	250	! -- LOCAL VARIABLES --
	251	!=======================================================================================
	252
	253	INTEGER :: i, k, it
	254	INTEGER :: nsplit
	255
	256	!Sources et Reservoirs de traceurs (ex:Radon):
	257	!--------------------------------------------
	258	!
	259	REAL,DIMENSION(:,:),ALLOCATABLE,SAVE :: source ! a voir lorsque le flux de surface est prescrit
[1813]	260	!$OMP THREADPRIVATE(source)
	261
[2146]	262	!
	263	!Entrees/Sorties: (cf ini_histrac.h et write_histrac.h)
	264	!---------------
	265	INTEGER :: iiq, ierr
	266	INTEGER :: nhori, nvert
	267	REAL :: zsto, zout, zjulian
	268	INTEGER,SAVE :: nid_tra ! pointe vers le fichier histrac.nc
[1813]	269	!$OMP THREADPRIVATE(nid_tra)
[2146]	270	REAL,DIMENSION(klon) :: zx_tmp_fi2d ! variable temporaire grille physique
	271	INTEGER :: itau_w ! pas de temps ecriture = nstep + itau_phy
	272	LOGICAL,PARAMETER :: ok_sync=.TRUE.
[1820]	273
[2146]	274	!
	275	! Nature du traceur
	276	!------------------
	277	LOGICAL,DIMENSION(:),ALLOCATABLE,SAVE :: aerosol ! aerosol(it) = true => aerosol => lessivage
[1813]	278	!$OMP THREADPRIVATE(aerosol) ! aerosol(it) = false => gaz
[2146]	279	REAL,DIMENSION(klon,klev) :: delp ! epaisseur de couche (Pa)
	280	!
	281	! Tendances de traceurs (Td) et flux de traceurs:
	282	!------------------------
	283	REAL,DIMENSION(klon,klev) :: d_tr ! Td dans l'atmosphere
	284	REAL,DIMENSION(klon,klev) :: Mint
	285	REAL,DIMENSION(klon,klev,nbtr) :: zmfd1a
	286	REAL,DIMENSION(klon,klev,nbtr) :: zmfdam
	287	REAL,DIMENSION(klon,klev,nbtr) :: zmfphi2
	288	! Physique
	289	!----------
	290	REAL,DIMENSION(klon,klev,nbtr) :: flestottr ! flux de lessivage dans chaque couche
	291	REAL,DIMENSION(klon,klev) :: zmasse ! densité atmosphérique Kg/m2
	292	REAL,DIMENSION(klon,klev) :: ztra_th
	293	!PhH
	294	REAL,DIMENSION(klon,klev) :: zrho
	295	REAL,DIMENSION(klon,klev) :: zdz
	296	REAL :: evaplsc,dx,beta ! variable pour lessivage Genthon
	297	REAL,DIMENSION(klon) :: his_dh ! ---
	298	! in-cloud scav variables
	299	REAL :: ql_incloud_ref ! ref value of in-cloud condensed water content
	300
	301	!Controles:
	302	!---------
	303	INTEGER,SAVE :: iflag_vdf_trac,iflag_con_trac,iflag_the_trac
	304	INTEGER,SAVE :: iflag_con_trac_omp, iflag_vdf_trac_omp,iflag_the_trac_omp
[1820]	305	!$OMP THREADPRIVATE(iflag_vdf_trac,iflag_con_trac,iflag_the_trac)
[1813]	306
[2146]	307	LOGICAL,SAVE :: lessivage
[1813]	308	!$OMP THREADPRIVATE(lessivage)
	309
[2146]	310	!RomP >>>
	311	INTEGER,SAVE :: iflag_lscav_omp,iflag_lscav
[2210]	312	REAL, SAVE :: ccntrAA_in,ccntrAA_omp
	313	REAL, SAVE :: ccntrENV_in,ccntrENV_omp
	314	REAL, SAVE :: coefcoli_in,coefcoli_omp
	315
[2146]	316	LOGICAL,SAVE :: convscav_omp,convscav
[1822]	317	!$OMP THREADPRIVATE(iflag_lscav)
[2210]	318	!$OMP THREADPRIVATE(ccntrAA_in,ccntrENV_in,coefcoli_in)
[1822]	319	!$OMP THREADPRIVATE(convscav)
[2146]	320	!RomP <<<
	321	!######################################################################
	322	! -- INITIALIZATION --
	323	!######################################################################
	324	IF (debutphy) THEN
	325	ALLOCATE(d_tr_cl(klon,klev,nbtr),d_tr_dry(klon,nbtr))
	326	ALLOCATE(flux_tr_dry(klon,nbtr),d_tr_dec(klon,klev,nbtr),d_tr_cv(klon,klev,nbtr))
	327	ALLOCATE(d_tr_insc(klon,klev,nbtr),d_tr_bcscav(klon,klev,nbtr))
	328	ALLOCATE(d_tr_evapls(klon,klev,nbtr),d_tr_ls(klon,klev,nbtr))
	329	ALLOCATE(qPrls(klon,nbtr),d_tr_trsp(klon,klev,nbtr))
	330	ALLOCATE(d_tr_sscav(klon,klev,nbtr),d_tr_sat(klon,klev,nbtr))
	331	ALLOCATE(d_tr_uscav(klon,klev,nbtr),qPr(klon,klev,nbtr),qDi(klon,klev,nbtr))
	332	ALLOCATE(qPa(klon,klev,nbtr),qMel(klon,klev,nbtr))
	333	ALLOCATE(qTrdi(klon,klev,nbtr),dtrcvMA(klon,klev,nbtr))
	334	ALLOCATE(d_tr_th(klon,klev,nbtr))
	335	ALLOCATE(d_tr_lessi_impa(klon,klev,nbtr),d_tr_lessi_nucl(klon,klev,nbtr))
	336	ENDIF
[1813]	337
[2146]	338	DO k=1,klev
	339	DO i=1,klon
	340	sourceBE(i,k)=0.
	341	Mint(i,k)=0.
	342	zrho(i,k)=0.
	343	zdz(i,k)=0.
	344	END DO
	345	END DO
[1813]	346
[2146]	347	DO it=1, nbtr
	348	DO k=1,klev
	349	DO i=1,klon
	350	d_tr_insc(i,k,it)=0.
	351	d_tr_bcscav(i,k,it)=0.
	352	d_tr_evapls(i,k,it)=0.
	353	d_tr_ls(i,k,it)=0.
	354	d_tr_cv(i,k,it)=0.
	355	d_tr_cl(i,k,it)=0.
	356	d_tr_trsp(i,k,it)=0.
	357	d_tr_sscav(i,k,it)=0.
	358	d_tr_sat(i,k,it)=0.
	359	d_tr_uscav(i,k,it)=0.
	360	d_tr_lessi_impa(i,k,it)=0.
	361	d_tr_lessi_nucl(i,k,it)=0.
	362	qDi(i,k,it)=0.
	363	qPr(i,k,it)=0.
	364	qPa(i,k,it)=0.
	365	qMel(i,k,it)=0.
	366	qTrdi(i,k,it)=0.
	367	dtrcvMA(i,k,it)=0.
	368	zmfd1a(i,k,it)=0.
	369	zmfdam(i,k,it)=0.
	370	zmfphi2(i,k,it)=0.
	371	END DO
	372	END DO
[1813]	373	END DO
[2146]	374
	375	DO k = 1, klev
	376	DO i = 1, klon
	377	delp(i,k) = paprs(i,k)-paprs(i,k+1)
	378	END DO
	379	END DO
	380
	381	IF (debutphy) THEN
	382	!!jyg
[1813]	383	!$OMP BARRIER
[2146]	384	ecrit_tra=86400. ! frequence de stokage en dur
	385	! obsolete car remplace par des ecritures dans phys_output_write
	386	!RomP >>>
	387	!
	388	!Config Key = convscav
	389	!Config Desc = Convective scavenging switch: 0=off, 1=on.
	390	!Config Def = .false.
	391	!Config Help =
	392	!
[1813]	393	!$OMP MASTER
[2146]	394	convscav_omp=.false.
	395	call getin('convscav', convscav_omp)
	396	iflag_vdf_trac_omp=1
	397	call getin('iflag_vdf_trac', iflag_vdf_trac_omp)
	398	iflag_con_trac_omp=1
	399	call getin('iflag_con_trac', iflag_con_trac_omp)
	400	iflag_the_trac_omp=1
	401	call getin('iflag_the_trac', iflag_the_trac_omp)
[1813]	402	!$OMP END MASTER
	403	!$OMP BARRIER
[2146]	404	convscav=convscav_omp
	405	iflag_vdf_trac=iflag_vdf_trac_omp
	406	iflag_con_trac=iflag_con_trac_omp
	407	iflag_the_trac=iflag_the_trac_omp
	408	write(lunout,*) 'phytrac passage dans routine conv avec lessivage', convscav
	409	!
	410	!Config Key = iflag_lscav
	411	!Config Desc = Large scale scavenging parametrization: 0=none, 1=old(Genthon92),
	412	! 2=1+PHeinrich, 3=Reddy_Boucher2004, 4=3+RPilon.
	413	!Config Def = 1
	414	!Config Help =
	415	!
[1813]	416	!$OMP MASTER
[2146]	417	iflag_lscav_omp=1
	418	call getin('iflag_lscav', iflag_lscav_omp)
[2210]	419	ccntrAA_omp=1
	420	ccntrENV_omp=1.
	421	coefcoli_omp=0.001
	422	call getin('ccntrAA', ccntrAA_omp)
	423	call getin('ccntrENV', ccntrENV_omp)
	424	call getin('coefcoli', coefcoli_omp)
[1813]	425	!$OMP END MASTER
	426	!$OMP BARRIER
[2146]	427	iflag_lscav=iflag_lscav_omp
[2210]	428	ccntrAA_in=ccntrAA_omp
	429	ccntrENV_in=ccntrENV_omp
	430	coefcoli_in=coefcoli_omp
[2146]	431	!
	432	SELECT CASE(iflag_lscav)
	433	CASE(0)
	434	WRITE(lunout,*) 'Large scale scavenging: none'
	435	CASE(1)
	436	WRITE(lunout,*) 'Large scale scavenging: C. Genthon, Tellus(1992), 44B, 371-389'
	437	CASE(2)
	438	WRITE(lunout,*) 'Large scale scavenging: C. Genthon, modified P. Heinrich'
	439	CASE(3)
	440	WRITE(lunout,*) 'Large scale scavenging: M. Shekkar Reddy and O. Boucher, JGR(2004), 109, D14202'
	441	CASE(4)
	442	WRITE(lunout,*) 'Large scale scavenging: Reddy and Boucher, modified R. Pilon'
	443	END SELECT
	444	!RomP <<<
	445	WRITE(,) 'FIRST TIME IN PHYTRAC : pdtphys(sec) = ',pdtphys,'ecrit_tra (sec) = ',ecrit_tra
	446	ALLOCATE( source(klon,nbtr), stat=ierr)
	447	IF (ierr /= 0) CALL abort_gcm('phytrac', 'pb in allocation 1',1)
[1813]	448
[2146]	449	ALLOCATE( aerosol(nbtr), stat=ierr)
	450	IF (ierr /= 0) CALL abort_gcm('phytrac', 'pb in allocation 2',1)
	451
	452
	453	! Initialize module for specific tracers
	454	SELECT CASE(type_trac)
	455	CASE('lmdz')
	456	CALL traclmdz_init(pctsrf,xlat,xlon,ftsol,tr_seri,t_seri,pplay,sh,pdtphys,aerosol,lessivage)
	457	CASE('inca')
	458	source(:,:)=0.
	459	CALL tracinca_init(aerosol,lessivage)
	460	CASE('repr')
	461	source(:,:)=0.
	462	END SELECT
	463
	464	!
	465	!--initialising coefficients for scavenging in the case of NP
	466	!
	467	ALLOCATE(flag_cvltr(nbtr))
	468	IF (iflag_con.EQ.3) THEN
	469	!
	470	ALLOCATE(ccntrAA(nbtr))
	471	ALLOCATE(ccntrENV(nbtr))
	472	ALLOCATE(coefcoli(nbtr))
	473	!
	474	DO it=1, nbtr
	475	SELECT CASE(type_trac)
	476	CASE('lmdz')
	477	IF (convscav.and.aerosol(it)) THEN
	478	flag_cvltr(it)=.true.
[2210]	479	ccntrAA(it) =ccntrAA_in !--a modifier par JYG a lire depuis fichier
	480	ccntrENV(it)=ccntrENV_in
	481	coefcoli(it)=coefcoli_in
[2146]	482	ELSE
	483	flag_cvltr(it)=.false.
	484	ENDIF
	485
	486	CASE('inca')
	487	! IF ((it.EQ.id_Rn222) .OR. ((it.GE.id_SO2) .AND. (it.LE.id_NH3)) ) THEN
	488	! !--gas-phase species
	489	! flag_cvltr(it)=.false.
[1813]	490	!
[2146]	491	! ELSEIF ( (it.GE.id_CIDUSTM) .AND. (it.LE.id_AIN) ) THEN
	492	! !--insoluble aerosol species
	493	! flag_cvltr(it)=.true.
	494	! ccntrAA(it)=0.7
	495	! ccntrENV(it)=0.7
	496	! coefcoli(it)=0.001
	497	! ELSEIF ( (it.EQ.id_Pb210) .OR. ((it.GE.id_CSSSM) .AND. (it.LE.id_SSN))) THEN
	498	! !--soluble aerosol species
	499	! flag_cvltr(it)=.true.
	500	! ccntrAA(it)=0.9
	501	! ccntrENV(it)=0.9
	502	! coefcoli(it)=0.001
	503	! ELSE
	504	! WRITE(lunout,*) 'pb it=', it
	505	! CALL abort_gcm('phytrac','pb it scavenging',1)
	506	! ENDIF
	507	!--test OB
	508	!--for now we do not scavenge in cvltr
	509	flag_cvltr(it)=.false.
	510	END SELECT
	511	ENDDO
	512	!
	513	ELSE ! iflag_con .ne. 3
	514	flag_cvltr(:) = .false.
	515	ENDIF
	516	!
	517	! Initialize diagnostic output
	518	! ----------------------------
[1813]	519	#ifdef CPP_IOIPSL
[2146]	520	! INCLUDE "ini_histrac.h"
[1813]	521	#endif
[2146]	522	!
	523	! print out all tracer flags
	524	!
	525	WRITE(lunout,*) 'print out all tracer flags'
	526	WRITE(lunout,*) 'type_trac =', type_trac
	527	WRITE(lunout,*) 'config_inca =', config_inca
	528	WRITE(lunout,*) 'iflag_con_trac =', iflag_con_trac
	529	WRITE(lunout,*) 'iflag_con =', iflag_con
	530	WRITE(lunout,*) 'convscav =', convscav
	531	WRITE(lunout,*) 'iflag_lscav =', iflag_lscav
	532	WRITE(lunout,*) 'aerosol =', aerosol
	533	WRITE(lunout,*) 'iflag_the_trac =', iflag_the_trac
	534	WRITE(lunout,*) 'iflag_thermals =', iflag_thermals
	535	WRITE(lunout,*) 'iflag_vdf_trac =', iflag_vdf_trac
	536	WRITE(lunout,*) 'pbl_flg =', pbl_flg
	537	WRITE(lunout,*) 'lessivage =', lessivage
	538	write(lunout,*) 'flag_cvltr = ', flag_cvltr
[1813]	539
[2146]	540	IF (lessivage.AND.config_inca.EQ.'inca') THEN
	541	CALL abort_gcm('phytrac', 'lessivage=T config_inca=inca impossible',1)
	542	STOP
	543	ENDIF
	544	!
	545	END IF ! of IF (debutphy)
	546	!############################################ END INITIALIZATION #######
[1813]	547
[2146]	548	DO k=1,klev
	549	DO i=1,klon
	550	zmasse(i,k)=(paprs(i,k)-paprs(i,k+1))/rg
	551	END DO
	552	END DO
	553	!
	554	IF (id_be .GT. 0) THEN
	555	DO k=1,klev
	556	DO i=1,klon
	557	sourceBE(i,k)=srcbe(i,k) !RomP -> pour sortie histrac
	558	END DO
	559	END DO
	560	ENDIF
[1813]	561
[2146]	562	!===============================================================================
	563	! -- Do specific treatment according to chemestry model or local LMDZ tracers
	564	!
	565	!===============================================================================
	566	SELECT CASE(type_trac)
	567	CASE('lmdz')
	568	! -- Traitement des traceurs avec traclmdz
	569	CALL traclmdz(nstep, julien, gmtime, pdtphys, t_seri, paprs, pplay, &
	570	cdragh, coefh, yu1, yv1, ftsol, pctsrf, xlat, xlon,iflag_vdf_trac>=0,sh, &
	571	rh, pphi, ustar, wstar, ale_bl, ale_wake, u10m, v10m, &
[2180]	572	tr_seri, source, d_tr_cl,d_tr_dec, zmasse) !RomP
[1813]	573
[2146]	574	CASE('inca')
	575	! -- CHIMIE INCA config_inca = aero or chem --
[1813]	576
[2146]	577	CALL tracinca(&
	578	nstep, julien, gmtime, lafin, &
	579	pdtphys, t_seri, paprs, pplay, &
	580	pmfu, upwd, ftsol, pctsrf, pphis, &
	581	pphi, albsol, sh, rh, &
	582	cldfra, rneb, diafra, cldliq, &
	583	itop_con, ibas_con, pmflxr, pmflxs, &
	584	prfl, psfl, aerosol_couple, flxmass_w, &
	585	tau_aero, piz_aero, cg_aero, ccm, &
	586	rfname, &
[2180]	587	tr_seri, source)
[1813]	588
[2146]	589	CASE('repr')
	590	! -- CHIMIE REPROBUS --
[1813]	591
[2146]	592	CALL tracreprobus(pdtphys, gmtime, debutphy, julien, &
	593	presnivs, xlat, xlon, pphis, pphi, &
	594	t_seri, pplay, paprs, sh , &
[2180]	595	tr_seri)
[1813]	596
[2146]	597	END SELECT
	598	!======================================================================
	599	! -- Calcul de l'effet de la convection --
	600	!======================================================================
[1813]	601
[2146]	602	IF (iflag_con_trac==1) THEN
[1813]	603
[2146]	604	DO it=1, nbtr
	605	IF ( conv_flg(it) == 0 ) CYCLE
	606	IF (iflag_con.LT.2) THEN
	607	!--pas de transport convectif
[1813]	608
[2146]	609	d_tr_cv(:,:,it)=0.
	610	ELSE IF (iflag_con.EQ.2) THEN
	611	!--ancien transport convectif de Tiedtke
[1813]	612
[2146]	613	CALL nflxtr(pdtphys, pmfu, pmfd, pen_u, pde_u, pen_d, pde_d, &
	614	pplay, paprs, tr_seri(:,:,it), d_tr_cv(:,:,it))
	615	ELSE
	616	!--nouveau transport convectif de Emanuel
[1813]	617
[2146]	618	IF (flag_cvltr(it)) THEN
	619	!--nouveau transport convectif de Emanuel avec lessivage convectif
	620	!
	621	!
	622	ccntrAA_3d(:,:) =ccntrAA(it)
	623	ccntrENV_3d(:,:)=ccntrENV(it)
	624	coefcoli_3d(:,:)=coefcoli(it)
[1813]	625
[2146]	626	!--beware this interface is a bit weird because it is called for each tracer
	627	!--with the full array tr_seri even if only item it is processed
[1813]	628
[2210]	629	print*,'CV SCAV ',it,ccntrAA(it),ccntrENV(it)
	630
[2146]	631	CALL cvltr_scav(pdtphys, da, phi,phi2,d1a,dam, mp,ep, &
	632	sigd,sij,wght_cvfd,clw,elij,epmlmMm,eplaMm, &
	633	pmflxr,pmflxs,evap,t_seri,wdtrainA,wdtrainM, &
	634	paprs,it,tr_seri,upwd,dnwd,itop_con,ibas_con, &
	635	ccntrAA_3d,ccntrENV_3d,coefcoli_3d, &
	636	d_tr_cv,d_tr_trsp,d_tr_sscav,d_tr_sat,d_tr_uscav,qDi,qPr,&
	637	qPa,qMel,qTrdi,dtrcvMA,Mint, &
	638	zmfd1a,zmfphi2,zmfdam)
[1813]	639
[2146]	640
	641	ELSE !---flag_cvltr(it).EQ.FALSE
	642	!--nouveau transport convectif de Emanuel mais pas de lessivage convectif
	643
	644	!--beware this interface is a bit weird because it is called for each tracer
	645	!--with the full array tr_seri even if only item it is processed
	646	!
	647	CALL cvltr_noscav(it,pdtphys, da, phi,mp,wght_cvfd,paprs,pplay, & !jyg
	648	tr_seri,upwd,dnwd,d_tr_cv) !jyg
	649
	650	ENDIF
	651
	652	ENDIF !--iflag
	653
	654	!--on ajoute les tendances
	655
	656	DO k = 1, klev
	657	DO i = 1, klon
	658	tr_seri(i,k,it) = tr_seri(i,k,it) + d_tr_cv(i,k,it)
	659	END DO
	660	END DO
	661
	662	CALL minmaxqfi(tr_seri(:,:,it),0.,1.e33,'convection it = '//solsym(it))
	663
	664	END DO ! nbtr
	665
	666	END IF ! convection
	667
	668	!======================================================================
	669	! -- Calcul de l'effet des thermiques --
	670	!======================================================================
	671
[1813]	672	DO it=1,nbtr
[2146]	673	DO k=1,klev
	674	DO i=1,klon
	675	d_tr_th(i,k,it)=0.
	676	tr_seri(i,k,it)=MAX(tr_seri(i,k,it),0.)
	677	tr_seri(i,k,it)=MIN(tr_seri(i,k,it),1.e10)
	678	END DO
	679	END DO
	680	END DO
[1813]	681
[2146]	682	IF (iflag_thermals.GT.0.AND.iflag_the_trac>0) THEN
[1813]	683
[2146]	684	DO it=1, nbtr
[1813]	685
[2146]	686	CALL thermcell_dq(klon,klev,1,pdtphys,fm_therm,entr_therm, &
	687	zmasse,tr_seri(1:klon,1:klev,it), &
	688	d_tr_th(1:klon,1:klev,it),ztra_th,0 )
[1813]	689
[2146]	690	DO k=1,klev
	691	DO i=1,klon
	692	d_tr_th(i,k,it)=pdtphys*d_tr_th(i,k,it)
	693	tr_seri(i,k,it)=MAX(tr_seri(i,k,it)+d_tr_th(i,k,it),0.)
	694	END DO
	695	END DO
[1813]	696
[2146]	697	END DO ! it
[1813]	698
[2146]	699	END IF ! Thermiques
[1813]	700
[2146]	701	!======================================================================
	702	! -- Calcul de l'effet de la couche limite --
	703	!======================================================================
	704
	705	IF (iflag_vdf_trac==1) THEN
	706
	707	! Injection during BL mixing
	708	!
	709	DO it=1, nbtr
	710	!
	711	IF( pbl_flg(it) /= 0 ) THEN
	712	!
	713	CALL cltrac(pdtphys, coefh,t_seri, &
	714	tr_seri(:,:,it), source(:,it), &
	715	paprs, pplay, delp, &
	716	d_tr_cl(:,:,it),d_tr_dry(:,it),flux_tr_dry(:,it))
	717	!
	718	tr_seri(:,:,it)=tr_seri(:,:,it)+d_tr_cl(:,:,it)
	719	!
	720	END IF
	721	!
	722	END DO
	723	!
	724	ELSE IF (iflag_vdf_trac==0) THEN
	725	!
	726	! Injection of source in the first model layer
	727	!
	728	DO it=1,nbtr
	729	d_tr_cl(:,1,it)=source(:,it)RG/delp(:,1)pdtphys
	730	tr_seri(:,1,it)=tr_seri(:,1,it)+d_tr_cl(:,1,it)
	731	ENDDO
	732	d_tr_cl(:,2:klev,1:nbtr)=0.
	733	!
	734	ELSE IF (iflag_vdf_trac==-1) THEN
	735	!
	736	! Nothing happens
	737	!
	738	d_tr_cl=0.
	739	!
	740	ELSE
	741	!
	742	CALL abort_gcm('iflag_vdf_trac', 'cas non prevu',1)
	743	!
	744	END IF ! couche limite
	745
	746	!======================================================================
	747	! Calcul de l'effet de la precipitation grande echelle
	748	! POUR INCA le lessivage est fait directement dans INCA
	749	!======================================================================
	750
	751	IF (lessivage) THEN
	752
	753	ql_incloud_ref = 10.e-4
	754	ql_incloud_ref = 5.e-4
	755
	756
	757	! calcul du contenu en eau liquide au sein du nuage
	758	ql_incl = ql_incloud_ref
	759	! choix du lessivage
	760	!
	761	IF (iflag_lscav .EQ. 3 .OR. iflag_lscav .EQ. 4) THEN
	762	! ******** Olivier Boucher version (3) possibly with modified ql_incl (4)
	763	!
	764	DO it = 1, nbtr
[2210]	765
	766	IF (aerosol(it)) THEN
[2146]	767	! incloud scavenging and removal by large scale rain ! orig : ql_incl was replaced by 0.5e-3 kg/kg
	768	! the value 0.5e-3 kg/kg is from Giorgi and Chameides (1986), JGR
	769	! Liu (2001) proposed to use 1.5e-3 kg/kg
	770
	771	CALL lsc_scav(pdtphys,it,iflag_lscav,ql_incl,prfl,psfl,rneb,beta_fisrt, &
[1813]	772	beta_v1,pplay,paprs,t_seri,tr_seri,d_tr_insc, &
	773	d_tr_bcscav,d_tr_evapls,qPrls)
	774
[2146]	775	!large scale scavenging tendency
	776	DO k = 1, klev
	777	DO i = 1, klon
	778	d_tr_ls(i,k,it)=d_tr_insc(i,k,it)+d_tr_bcscav(i,k,it)+d_tr_evapls(i,k,it)
	779	tr_seri(i,k,it)=tr_seri(i,k,it)+d_tr_ls(i,k,it)
	780	ENDDO
	781	ENDDO
	782	CALL minmaxqfi(tr_seri(:,:,it),0.,1.e33,'lsc scav it = '//solsym(it))
[2210]	783	ENDIF
	784
[2146]	785	END DO !tr
[1813]	786
[2146]	787	ELSE IF (iflag_lscav .EQ. 2) THEN ! frac_impa, frac_nucl
	788	! ********* modified old version
[1813]	789
[2146]	790	d_tr_lessi_nucl(:,:,:) = 0.
	791	d_tr_lessi_impa(:,:,:) = 0.
	792	flestottr(:,:,:) = 0.
	793	! Tendance des aerosols nuclees et impactes
	794	DO it = 1, nbtr
	795	IF (aerosol(it)) THEN
	796	his_dh(:)=0.
	797	DO k = 1, klev
	798	DO i = 1, klon
	799	!PhH
	800	zrho(i,k)=pplay(i,k)/t_seri(i,k)/RD
	801	zdz(i,k)=(paprs(i,k)-paprs(i,k+1))/zrho(i,k)/RG
	802	!
	803	END DO
	804	END DO
[1813]	805
[2146]	806	DO k=klev-1, 1, -1
	807	DO i=1, klon
	808	! d_tr_ls(i,k,it)=tr_seri(i,k,it)(frac_impa(i,k)frac_nucl(i,k)-1.)
	809	dx=d_tr_ls(i,k,it)
	810	his_dh(i)=his_dh(i)-dxzrho(i,k)zdz(i,k)/pdtphys ! kg/m2/s
	811	evaplsc = prfl(i,k) - prfl(i,k+1) + psfl(i,k) - psfl(i,k+1)
	812	! Evaporation Partielle -> Liberation Partielle 0.5*evap
	813	IF ( evaplsc .LT.0..and.abs(prfl(i,k+1)+psfl(i,k+1)).gt.1.e-10) THEN
	814	evaplsc = (-evaplsc)/(prfl(i,k+1)+psfl(i,k+1))
	815	! evaplsc est donc positif, his_dh(i) est positif
	816	!--------------
	817	d_tr_evapls(i,k,it)=0.5evaplsc(d_tr_lessi_nucl(i,k+1,it) &
	818	+d_tr_lessi_impa(i,k+1,it))
	819	!------------- d_tr_evapls(i,k,it)=-0.5evaplsc(d_tr_lsc(i,k+1,it))
	820	beta=0.5*evaplsc
	821	if ((prfl(i,k)+psfl(i,k)).lt.1.e-10) THEN
	822	beta=1.0*evaplsc
	823	endif
	824	dx=betahis_dh(i)/zrho(i,k)/zdz(i,k)pdtphys
	825	his_dh(i)=(1.-beta)*his_dh(i) ! tracer from
	826	d_tr_evapls(i,k,it)=dx
	827	ENDIF
	828	d_tr_ls(i,k,it)=tr_seri(i,k,it)(frac_impa(i,k)frac_nucl(i,k)-1.) &
	829	+d_tr_evapls(i,k,it)
[1813]	830
[2146]	831	!--------------
	832	d_tr_lessi_nucl(i,k,it) = d_tr_lessi_nucl(i,k,it) + &
	833	( 1 - frac_nucl(i,k) )*tr_seri(i,k,it)
	834	d_tr_lessi_impa(i,k,it) = d_tr_lessi_impa(i,k,it) + &
	835	( 1 - frac_impa(i,k) )*tr_seri(i,k,it)
	836	!
	837	! Flux lessivage total
	838	flestottr(i,k,it) = flestottr(i,k,it) - &
	839	( d_tr_lessi_nucl(i,k,it) + &
	840	d_tr_lessi_impa(i,k,it) ) * &
	841	( paprs(i,k)-paprs(i,k+1) ) / &
	842	(RG * pdtphys)
	843	!! Mise a jour des traceurs due a l'impaction,nucleation
	844	! tr_seri(i,k,it)=tr_seri(i,k,it)frac_impa(i,k)frac_nucl(i,k)
	845	!! calcul de la tendance liee au lessivage stratiforme
	846	! d_tr_ls(i,k,it)=tr_seri(i,k,it)*&
	847	! (1.-1./(frac_impa(i,k)*frac_nucl(i,k)))
	848	!--------------
	849	END DO
	850	END DO
	851	END IF
	852	END DO
	853	! ********* end modified old version
[1813]	854
[2146]	855	ELSE IF (iflag_lscav .EQ. 1) THEN ! frac_impa, frac_nucl
	856	! ********* old version
[1813]	857
[2146]	858	d_tr_lessi_nucl(:,:,:) = 0.
	859	d_tr_lessi_impa(:,:,:) = 0.
	860	flestottr(:,:,:) = 0.
	861	!=========================
	862	! LESSIVAGE LARGE SCALE :
	863	!=========================
[1813]	864
[2146]	865	! Tendance des aerosols nuclees et impactes
	866	! -----------------------------------------
	867	DO it = 1, nbtr
	868	IF (aerosol(it)) THEN
	869	DO k = 1, klev
	870	DO i = 1, klon
	871	d_tr_lessi_nucl(i,k,it) = d_tr_lessi_nucl(i,k,it) + &
	872	( 1 - frac_nucl(i,k) )*tr_seri(i,k,it)
	873	d_tr_lessi_impa(i,k,it) = d_tr_lessi_impa(i,k,it) + &
	874	( 1 - frac_impa(i,k) )*tr_seri(i,k,it)
[1813]	875
[2146]	876	!
	877	! Flux lessivage total
	878	! ------------------------------------------------------------
	879	flestottr(i,k,it) = flestottr(i,k,it) - &
	880	( d_tr_lessi_nucl(i,k,it) + &
	881	d_tr_lessi_impa(i,k,it) ) * &
	882	( paprs(i,k)-paprs(i,k+1) ) / &
	883	(RG * pdtphys)
	884	!
	885	! Mise a jour des traceurs due a l'impaction,nucleation
	886	! ----------------------------------------------------------------------
	887	tr_seri(i,k,it)=tr_seri(i,k,it)frac_impa(i,k)frac_nucl(i,k)
	888	END DO
	889	END DO
	890	END IF
	891	END DO
[1813]	892
[2146]	893	! ********* end old version
	894	ENDIF ! iflag_lscav . EQ. 1, 2, 3 or 4
	895	!
	896	END IF ! lessivage
	897
	898	!=============================================================
	899	! Ecriture des sorties
	900	!=============================================================
[1813]	901	#ifdef CPP_IOIPSL
[2146]	902	! INCLUDE "write_histrac.h"
[1813]	903	#endif
	904
[2146]	905	END SUBROUTINE phytrac
[1813]	906
	907	END MODULE

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

Download in other formats:

Original Format