Context Navigation

phytracr_spl_mod.F90 @ 3817

Last change on this file since 3817 was 3814, checked in by asima, 4 years ago
SPLA : unplugging Jeronimo Escribano’s assimilation emission coefficient files, in a minimalist, non-distructive way. Each file reading is put under an IF("ASSIM"=="YES"), always false. Thorough removal of assimilation coefficients will follow.
File size: 174.0 KB

Rev	Line
[2630]	1	! $Id: physiq.F90 2298 2015-06-14 19:13:32Z fairhead $
	2	!#define IO_DEBUG
	3
	4	MODULE phytracr_spl_mod
	5
	6	! Recuperation des morceaux de la physique de Jeronimo specifiques
	7	! du modele d'aerosols d'Olivier n'co.
	8	!
	9	INCLUDE "chem.h"
	10	INCLUDE "chem_spla.h"
	11
	12	REAL,SAVE :: scale_param_ssacc !Scaling parameter for Fine Sea Salt
	13	REAL,SAVE :: scale_param_sscoa !Scaling parameter for Coarse Sea Salt
	14
	15
	16
	17	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: scale_param_ind !Scaling parameter for industrial emissions of SO2
	18	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: scale_param_bb !Scaling parameter for biomas burning (SO2,BC & OM)
	19	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: scale_param_ff !Scaling parameter for industrial emissions (fossil fuel)
	20	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: scale_param_dustacc !Scaling parameter for Fine Dust
	21	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: scale_param_dustcoa !Scaling parameter for Coarse Dust
	22	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: scale_param_dustsco !Scaling parameter for SCoarse Dust
	23	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: param_wstarBLperregion !parameter for ..
	24	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: param_wstarWAKEperregion !parameter for ..
	25	!$OMP THREADPRIVATE(scale_param_ind,scale_param_bb,scale_param_ff)
	26	!$OMP THREADPRIVATE(scale_param_dustacc,scale_param_dustcoa,scale_param_dustsco)
	27	!$OMP THREADPRIVATE(scale_param_ssacc,scale_param_sscoa)
	28	!$OMP THREADPRIVATE(param_wstarBLperregion,param_wstarWAKEperregion)
	29	REAL, DIMENSION(:),ALLOCATABLE,SAVE ::dust_ec, u10m_ec, v10m_ec
	30	!$OMP THREADPRIVATE(dust_ec, u10m_ec, v10m_ec)
	31
	32	CHARACTER*800 fileregionsdimsind
	33	CHARACTER*800 fileregionsdimsdust
	34	CHARACTER*800 fileregionsdimsbb
	35	CHARACTER*800 fileregionsdimswstar
	36	CHARACTER*100 paramname_ind
	37	CHARACTER*100 paramname_bb
	38	CHARACTER*100 paramname_ff
	39	CHARACTER*100 paramname_dustacc
	40	CHARACTER*100 paramname_dustcoa
	41	CHARACTER*100 paramname_dustsco
	42	CHARACTER*100 paramname_ssacc
	43	CHARACTER*100 paramname_sscoa
	44	CHARACTER*100 paramname_wstarBL
	45	CHARACTER*100 paramname_wstarWAKE
	46
	47
	48	CHARACTER*800 filescaleparams
	49	CHARACTER*800 paramsname
	50
	51
	52	!!------------------------ SULFUR emissions ----------------------------
	53	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_so2volc_cont ! emissions so2 volcan continuous
	54	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_altvolc_cont ! altitude so2 volcan continuous
	55	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_so2volc_expl ! emissions so2 volcan explosive
	56	!$OMP THREADPRIVATE( lmt_so2volc_cont,lmt_altvolc_cont,lmt_so2volc_expl )
	57	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_altvolc_expl ! altitude so2 volcan explosive
	58	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_so2ff_l ! emissions so2 fossil fuel (low)
	59	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_so2ff_h ! emissions so2 fossil fuel (high)
	60	!$OMP THREADPRIVATE( lmt_altvolc_expl,lmt_so2ff_l,lmt_so2ff_h )
	61	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_so2nff ! emissions so2 non-fossil fuel
	62	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_so2ba ! emissions de so2 bateau
	63	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_so2bb_l ! emissions de so2 biomass burning (low)
	64	!$OMP THREADPRIVATE( lmt_so2nff,lmt_so2ba,lmt_so2bb_l )
	65	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_so2bb_h ! emissions de so2 biomass burning (high)
	66	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_dmsconc ! concentration de dms oceanique
	67	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_dmsbio ! emissions de dms bio
	68	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_h2sbio ! emissions de h2s bio
	69	!$OMP THREADPRIVATE(lmt_so2bb_h,lmt_dmsconc,lmt_dmsbio,lmt_h2sbio )
	70	!------------------------- BLACK CARBON emissions ----------------------
	71	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_bcff ! emissions de BC fossil fuels
	72	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_bcnff ! emissions de BC non-fossil fuels
	73	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_bcbb_l ! emissions de BC biomass basses
	74	!$OMP THREADPRIVATE( lmt_bcff,lmt_bcnff,lmt_bcbb_l)
	75	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_bcbb_h ! emissions de BC biomass hautes
	76	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_bcba ! emissions de BC bateau
	77	!$OMP THREADPRIVATE(lmt_bcbb_h,lmt_bcba)
	78	!------------------------ ORGANIC MATTER emissions ---------------------
	79	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_omff ! emissions de OM fossil fuels
	80	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_omnff ! emissions de OM non-fossil fuels
	81	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_ombb_l ! emissions de OM biomass basses
	82	!$OMP THREADPRIVATE( lmt_omff,lmt_omnff,lmt_ombb_l)
	83	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_ombb_h ! emissions de OM biomass hautes
	84	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_omnat ! emissions de OM Natural
	85	REAL , DIMENSION(:),ALLOCATABLE,SAVE :: lmt_omba ! emissions de OM bateau
	86	REAL , DIMENSION(:,:),ALLOCATABLE,SAVE :: lmt_sea_salt ! emissions de OM Natural
	87	!$OMP THREADPRIVATE(lmt_ombb_h,lmt_omnat,lmt_omba,lmt_sea_salt)
	88
	89	!JE20141224 >>
	90	! others
	91	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: tsol
	92	!$OMP THREADPRIVATE(tsol)
	93	INTEGER :: ijulday
	94	LOGICAL , parameter :: edgar = .true.
	95	INTEGER , parameter :: flag_dms=4
	96	INTEGER*4 nbjour
	97
	98	!
	99	! Tracer tendencies, for outputs
	100	!-------------------------------
	101	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_cl ! Td couche
	102	!. limite/traceur
	103	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_dec
	104	!RomP
	105	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_cv ! Td
	106	!onvection/traceur
	107	! RomP >>>
	108	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_insc
	109	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_bcscav
	110	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_evapls
	111	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_ls
	112	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_trsp
	113	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_sscav
	114	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_sat
	115	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_uscav
	116	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: qPr,qDi ! concentration tra
	117	!dans pluie,air descente insaturee
	118	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: qPa,qMel
	119	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: qTrdi,dtrcvMA ! conc traceur
	120	!descente air insaturee et td convective MA
	121	!! RomP <<<
	122	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_th ! Td thermique
	123	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_lessi_impa ! Td du
	124	!lessivage par impaction
	125	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_lessi_nucl ! Td du
	126	!lessivage par nucleation
	127	REAL,DIMENSION(:,:),ALLOCATABLE,SAVE :: qPrls !jyg:
	128	!oncentration tra dans pluie LS a la surf.
	129	REAL,DIMENSION(:,:),ALLOCATABLE,SAVE :: d_tr_dry ! Td depot
	130	!sec/traceur (1st layer),ALLOCATABLE,SAVE jyg
	131	REAL,DIMENSION(:,:),ALLOCATABLE,SAVE :: flux_tr_dry ! depot
	132	!sec/traceur (surface),ALLOCATABLE,SAVE jyg
	133
	134	! Index of each traceur
	135	INTEGER,SAVE :: id_prec, id_fine, id_coss, id_codu, id_scdu
	136
	137	!$OMP THREADPRIVATE(d_tr_cl,d_tr_dec,d_tr_cv,d_tr_insc,d_tr_bcscav,d_tr_evapls)
	138	!$OMP THREADPRIVATE(d_tr_ls,d_tr_trsp,d_tr_sscav,d_tr_sat,d_tr_uscav)
	139	!$OMP THREADPRIVATE(qPr,qDi,qPa,qMel,qTrdi,dtrcvMA,d_tr_th,d_tr_lessi_impa)
	140	!$OMP THREADPRIVATE(d_tr_lessi_nucl,qPrls,d_tr_dry,flux_tr_dry)
	141	!$OMP THREADPRIVATE(id_prec,id_fine,id_coss,id_codu,id_scdu)
	142
	143	! JE20141224 <<
	144
	145	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diff_aod550_tot ! epaisseur optique total aerosol 550 nm
	146	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diag_aod670_tot ! epaisseur optique total aerosol 670 nm
	147	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diag_aod865_tot ! epaisseur optique total aerosol 865 nm
	148	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diff_aod550_tr2 ! epaisseur optique Traceur 2 aerosol 550 nm, diagnostic
	149	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diag_aod670_tr2 ! epaisseur optique Traceur 2 aerosol 670 nm, diagnostic
	150	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diag_aod865_tr2 ! epaisseur optique Traceur 2 aerosol 865 nm, diagnostic
	151	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diag_aod550_ss ! epaisseur optique Sels marins aerosol 550 nm, diagnostic
	152	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diag_aod670_ss ! epaisseur optique Sels marins aerosol 670 nm, diagnostic
	153	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diag_aod865_ss ! epaisseur optique Sels marins aerosol 865 nm, diagnostic
	154	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diag_aod550_dust ! epaisseur optique Dust aerosol 550 nm, diagnostic
	155	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diag_aod670_dust ! epaisseur optique Dust aerosol 670 nm, diagnostic
	156	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diag_aod865_dust ! epaisseur optique Dust aerosol 865 nm, diagnostic
	157	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diag_aod550_dustsco ! epaisseur optique Dust SCOarse aerosol 550 nm, diagnostic
	158	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diag_aod670_dustsco ! epaisseur optique Dust SCOarse aerosol 670 nm, diagnostic
	159	REAL,DIMENSION(:),ALLOCATABLE,SAVE :: diag_aod865_dustsco ! epaisseur optique Dust SCOarse aerosol 865 nm, diagnostic
	160
	161	!$OMP THREADPRIVATE(diff_aod550_tot,diag_aod670_tot,diag_aod865_tot)
	162	!$OMP THREADPRIVATE(diff_aod550_tr2,diag_aod670_tr2,diag_aod865_tr2)
	163	!$OMP THREADPRIVATE(diag_aod550_ss,diag_aod670_ss,diag_aod865_ss,diag_aod550_dust)
	164	!$OMP THREADPRIVATE(diag_aod670_dust,diag_aod865_dust,diag_aod550_dustsco)
	165	!$OMP THREADPRIVATE(diag_aod670_dustsco,diag_aod865_dustsco)
	166
	167
	168	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_terra ! AOD at terra overpass time ( 10.30 local hour)
	169	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_tr2_terra ! AOD at terra overpass time ( 10.30 local hour)
	170	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_ss_terra ! AOD at terra overpass time ( 10.30 local hour)
	171	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_dust_terra ! AOD at terra overpass time ( 10.30 local hour)
	172	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_dustsco_terra ! AOD at terra overpass time ( 10.30 local hour)
	173	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_terra ! AOD at terra overpass time ( 10.30 local hour)
	174	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_tr2_terra ! AOD at terra overpass time ( 10.30 local hour)
	175	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_ss_terra ! AOD at terra overpass time ( 10.30 local hour)
	176	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_dust_terra ! AOD at terra overpass time ( 10.30 local hour)
	177	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_dustsco_terra ! AOD at terra overpass time ( 10.30 local hour)
	178	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_terra ! AOD at terra overpass time ( 10.30 local hour)
	179	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_tr2_terra ! AOD at terra overpass time ( 10.30 local hour)
	180	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_ss_terra ! AOD at terra overpass time ( 10.30 local hour)
	181	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_dust_terra ! AOD at terra overpass time ( 10.30 local hour)
	182	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_dustsco_terra ! AOD at terra overpass time ( 10.30 local hour)
	183
	184
	185	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	186	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_tr2_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	187	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_ss_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	188	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_dust_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	189	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_dustsco_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	190	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	191	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_tr2_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	192	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_ss_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	193	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_dust_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	194	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_dustsco_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	195	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	196	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_tr2_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	197	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_ss_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	198	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_dust_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	199	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_dustsco_aqua ! AOD at aqua overpass time ( 13.30 local hour)
	200
	201	!$OMP THREADPRIVATE(aod550_aqua,aod550_tr2_aqua,aod550_ss_aqua,aod550_dust_aqua,aod550_dustsco_aqua)
	202	!$OMP THREADPRIVATE(aod670_aqua,aod670_tr2_aqua,aod670_ss_aqua,aod670_dust_aqua,aod670_dustsco_aqua)
	203	!$OMP THREADPRIVATE(aod865_aqua,aod865_tr2_aqua,aod865_ss_aqua,aod865_dust_aqua,aod865_dustsco_aqua)
	204	!$OMP THREADPRIVATE(aod550_terra,aod550_tr2_terra,aod550_ss_terra,aod550_dust_terra,aod550_dustsco_terra)
	205	!$OMP THREADPRIVATE(aod670_terra,aod670_tr2_terra,aod670_ss_terra,aod670_dust_terra,aod670_dustsco_terra)
	206	!$OMP THREADPRIVATE(aod865_terra,aod865_tr2_terra,aod865_ss_terra,aod865_dust_terra,aod865_dustsco_terra)
	207
	208
	209	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: sconc01 ! surface concentration
	210	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: trm01 ! burden
	211	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: sconc02 ! surface concentration
	212	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: trm02 ! burden
	213	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: sconc03 ! surface concentration
	214	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: trm03 ! burden
	215	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: sconc04 ! surface concentration
	216	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: trm04 ! burden
	217	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: sconc05 ! surface concentration
	218	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: trm05 ! burden
	219	!$OMP THREADPRIVATE(sconc01,sconc02,sconc03,sconc04,sconc05)
	220	!$OMP THREADPRIVATE(trm01,trm02,trm03,trm04,trm05)
	221	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: flux01
	222	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: flux02
	223	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: flux03
	224	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: flux04
	225	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: flux05
	226	!$OMP THREADPRIVATE(flux01,flux02,flux03,flux04,flux05)
	227	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: ds01
	228	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: ds02
	229	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: ds03
	230	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: ds04
	231	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: ds05
	232	!$OMP THREADPRIVATE(ds01,ds02,ds03,ds04,ds05)
	233	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dh01
	234	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dh02
	235	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dh03
	236	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dh04
	237	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dh05
	238	!$OMP THREADPRIVATE(dh01,dh02,dh03,dh04,dh05)
	239	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dtrconv01
	240	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dtrconv02
	241	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dtrconv03
	242	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dtrconv04
	243	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dtrconv05
	244	!$OMP THREADPRIVATE(dtrconv01,dtrconv02,dtrconv03,dtrconv04,dtrconv05)
	245	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dtherm01
	246	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dtherm02
	247	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dtherm03
	248	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dtherm04
	249	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dtherm05
	250	!$OMP THREADPRIVATE(dtherm01,dtherm02,dtherm03,dtherm04,dtherm05)
	251	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dhkecv01
	252	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dhkecv02
	253	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dhkecv03
	254	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dhkecv04
	255	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dhkecv05
	256	!$OMP THREADPRIVATE(dhkecv01,dhkecv02,dhkecv03,dhkecv04,dhkecv05)
	257	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: d_tr_ds01
	258	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: d_tr_ds02
	259	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: d_tr_ds03
	260	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: d_tr_ds04
	261	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: d_tr_ds05
	262	!$OMP THREADPRIVATE(d_tr_ds01,d_tr_ds02,d_tr_ds03,d_tr_ds04,d_tr_ds05)
	263	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dhkelsc01
	264	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dhkelsc02
	265	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dhkelsc03
	266	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dhkelsc04
	267	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: dhkelsc05
	268	!$OMP THREADPRIVATE(dhkelsc01,dhkelsc02,dhkelsc03,dhkelsc04,dhkelsc05)
	269	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_cv01
	270	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_cv02
	271	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_cv03
	272	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_cv04
	273	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_cv05
	274	!$OMP THREADPRIVATE(d_tr_cv01,d_tr_cv02,d_tr_cv03,d_tr_cv04,d_tr_cv05)
	275	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_trsp01
	276	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_trsp02
	277	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_trsp03
	278	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_trsp04
	279	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_trsp05
	280	!$OMP THREADPRIVATE(d_tr_trsp01,d_tr_trsp02,d_tr_trsp03,d_tr_trsp04,d_tr_trsp05)
	281	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_sscav01
	282	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_sscav02
	283	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_sscav03
	284	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_sscav04
	285	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_sscav05
	286	!$OMP THREADPRIVATE(d_tr_sscav01,d_tr_sscav02,d_tr_sscav03,d_tr_sscav04,d_tr_sscav05)
	287	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_sat01
	288	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_sat02
	289	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_sat03
	290	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_sat04
	291	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_sat05
	292	!$OMP THREADPRIVATE(d_tr_sat01,d_tr_sat02,d_tr_sat03,d_tr_sat04,d_tr_sat05)
	293	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_uscav01
	294	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_uscav02
	295	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_uscav03
	296	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_uscav04
	297	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_uscav05
	298	!$OMP THREADPRIVATE(d_tr_uscav01,d_tr_uscav02,d_tr_uscav03,d_tr_uscav04,d_tr_uscav05)
	299
	300	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	301	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	302	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	303
	304	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_insc01
	305	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_insc02
	306	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_insc03
	307	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_insc04
	308	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_insc05
	309	!$OMP THREADPRIVATE(d_tr_insc01,d_tr_insc02,d_tr_insc03,d_tr_insc04,d_tr_insc05)
	310	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_bcscav01
	311	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_bcscav02
	312	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_bcscav03
	313	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_bcscav04
	314	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_bcscav05
	315	!$OMP THREADPRIVATE(d_tr_bcscav01,d_tr_bcscav02,d_tr_bcscav03,d_tr_bcscav04,d_tr_bcscav05)
	316	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_evapls01
	317	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_evapls02
	318	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_evapls03
	319	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_evapls04
	320	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_evapls05
	321	!$OMP THREADPRIVATE(d_tr_evapls01,d_tr_evapls02,d_tr_evapls03,d_tr_evapls04,d_tr_evapls05)
	322	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_ls01
	323	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_ls02
	324	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_ls03
	325	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_ls04
	326	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_ls05
	327	!$OMP THREADPRIVATE(d_tr_ls01,d_tr_ls02,d_tr_ls03,d_tr_ls04,d_tr_ls05)
	328
	329	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_dyn01
	330	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_dyn02
	331	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_dyn03
	332	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_dyn04
	333	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_dyn05
	334	!$OMP THREADPRIVATE(d_tr_dyn01,d_tr_dyn02,d_tr_dyn03,d_tr_dyn04,d_tr_dyn05)
	335
	336	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_cl01
	337	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_cl02
	338	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_cl03
	339	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_cl04
	340	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_cl05
	341	!$OMP THREADPRIVATE(d_tr_cl01,d_tr_cl02,d_tr_cl03,d_tr_cl04,d_tr_cl05)
	342
	343	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_th01
	344	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_th02
	345	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_th03
	346	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_th04
	347	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: d_tr_th05
	348	!$OMP THREADPRIVATE(d_tr_th01,d_tr_th02,d_tr_th03,d_tr_th04,d_tr_th05)
	349
	350	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: sed_ss3D ! corresponds to tracer 3
	351	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: sed_dust3D ! corresponds to tracer 4
	352	REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: sed_dustsco3D ! corresponds to tracer 4
	353	!$OMP THREADPRIVATE(sed_ss3D,sed_dust3D,sed_dustsco3D)
	354
	355	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	356	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	357	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: sed_ss ! corresponds to tracer 3
	358	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: sed_dust ! corresponds to tracer 4
	359	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: sed_dustsco ! corresponds to tracer 4
	360	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: his_g2pgas ! corresponds to tracer 4
	361	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: his_g2paer ! corresponds to tracer 4
	362	!$OMP THREADPRIVATE(sed_ss,sed_dust,sed_dustsco,his_g2pgas,his_g2paer)
	363
	364	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxbb
	365	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxff
	366	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxbcbb
	367	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxbcff
	368	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxbcnff
	369	!$OMP THREADPRIVATE(fluxbb,fluxff,fluxbcbb,fluxbcff,fluxbcnff)
	370	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxbcba
	371	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxbc
	372	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxombb
	373	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxomff
	374	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxomnff
	375	!$OMP THREADPRIVATE(fluxbcba,fluxbc,fluxombb,fluxomff,fluxomnff)
	376	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxomba
	377	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxomnat
	378	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxom
	379	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxh2sff
	380	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxh2snff
	381	!$OMP THREADPRIVATE(fluxomba,fluxomnat,fluxom,fluxh2sff,fluxh2snff)
	382	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxso2ff
	383	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxso2nff
	384	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxso2bb
	385	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxso2vol
	386	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxso2ba
	387	!$OMP THREADPRIVATE(fluxso2ff,fluxso2nff,fluxso2bb,fluxso2vol,fluxso2ba)
	388	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxso2
	389	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxso4ff
	390	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxso4nff
	391	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxso4bb
	392	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxso4ba
	393	!$OMP THREADPRIVATE(fluxso2,fluxso4ff,fluxso4nff,fluxso4ba,fluxso4bb)
	394	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxso4
	395	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxdms
	396	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxh2sbio
	397	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxdustec
	398	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxddfine
	399	!$OMP THREADPRIVATE(fluxso4,fluxdms,fluxh2sbio,fluxdustec,fluxddfine)
	400	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxddcoa
	401	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxddsco
	402	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxdd
	403	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxssfine
	404	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxsscoa
	405	!$OMP THREADPRIVATE(fluxddcoa,fluxddsco,fluxdd,fluxssfine,fluxsscoa)
	406	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: fluxss
	407	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: flux_sparam_ind
	408	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: flux_sparam_bb
	409	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: flux_sparam_ff
	410	!$OMP THREADPRIVATE(fluxss,flux_sparam_ind,flux_sparam_bb,flux_sparam_ff)
	411	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: flux_sparam_ddfine
	412	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: flux_sparam_ddcoa
	413	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: flux_sparam_ddsco
	414	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: flux_sparam_ssfine
	415	!$OMP THREADPRIVATE(flux_sparam_ddfine,flux_sparam_ddcoa)
	416	!$OMP THREADPRIVATE(flux_sparam_ddsco,flux_sparam_ssfine)
	417	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: flux_sparam_sscoa
	418	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: u10m_ss
	419	REAL, DIMENSION(:), ALLOCATABLE, SAVE :: v10m_ss
	420	!$OMP THREADPRIVATE(flux_sparam_sscoa,u10m_ss,v10m_ss)
	421
[3806]	422	! Select dust emission scheme for the Sahara:
[2630]	423	! LOGICAL,PARAMETER,SAVE :: ok_chimeredust=.FALSE.
	424	LOGICAL,PARAMETER :: ok_chimeredust=.TRUE.
	425	!!!!!! !$OMP THREADPRIVATE(ok_chimeredust)
	426
	427
	428	CONTAINS
	429	!
	430	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[3786]	431	SUBROUTINE phytracr_spl_out_init()
	432	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[3806]	433	!AS : This subroutine centralises the ALLOCATE needed for the 1st call of
	434	! phys_output_write_spl in physiq
[3786]	435
	436	USE dimphy
	437	USE infotrac_phy, ONLY: nbtr
	438	USE dustemission_mod, ONLY : dustemis_out_init
	439
	440	! pour les variables m[1-3]dflux
	441	CALL dustemis_out_init()
	442
	443	!traceur
	444	ALLOCATE( diff_aod550_tot(klon) )
	445	ALLOCATE( diag_aod670_tot(klon) )
	446	ALLOCATE( diag_aod865_tot(klon) )
	447	ALLOCATE( diff_aod550_tr2(klon) )
	448	ALLOCATE( diag_aod670_tr2(klon) )
	449	ALLOCATE( diag_aod865_tr2(klon) )
	450	ALLOCATE( diag_aod550_ss(klon) )
	451	ALLOCATE( diag_aod670_ss(klon) )
	452	ALLOCATE( diag_aod865_ss(klon) )
	453	ALLOCATE( diag_aod550_dust(klon) )
	454	ALLOCATE( diag_aod670_dust(klon) )
	455	ALLOCATE( diag_aod865_dust(klon) )
	456	ALLOCATE( diag_aod550_dustsco(klon) )
	457	ALLOCATE( diag_aod670_dustsco(klon) )
	458	ALLOCATE( diag_aod865_dustsco(klon) )
	459	!AS: les 15 vars _terra et 15 _aqua suivantes sont groupees differemment dans spla_output_write.h
	460	ALLOCATE( aod550_terra(klon))
	461	ALLOCATE( aod550_tr2_terra(klon))
	462	ALLOCATE( aod550_ss_terra(klon))
	463	ALLOCATE( aod550_dust_terra(klon))
	464	ALLOCATE( aod550_dustsco_terra(klon))
	465	ALLOCATE( aod670_terra(klon))
	466	ALLOCATE( aod670_tr2_terra(klon))
	467	ALLOCATE( aod670_ss_terra(klon))
	468	ALLOCATE( aod670_dust_terra(klon))
	469	ALLOCATE( aod670_dustsco_terra(klon))
	470	ALLOCATE( aod865_terra(klon))
	471	ALLOCATE( aod865_tr2_terra(klon))
	472	ALLOCATE( aod865_ss_terra(klon))
	473	ALLOCATE( aod865_dust_terra(klon))
	474	ALLOCATE( aod865_dustsco_terra(klon))
	475
	476	ALLOCATE( aod550_aqua(klon))
	477	ALLOCATE( aod550_tr2_aqua(klon))
	478	ALLOCATE( aod550_ss_aqua(klon))
	479	ALLOCATE( aod550_dust_aqua(klon))
	480	ALLOCATE( aod550_dustsco_aqua(klon))
	481	ALLOCATE( aod670_aqua(klon))
	482	ALLOCATE( aod670_tr2_aqua(klon))
	483	ALLOCATE( aod670_ss_aqua(klon))
	484	ALLOCATE( aod670_dust_aqua(klon))
	485	ALLOCATE( aod670_dustsco_aqua(klon))
	486	ALLOCATE( aod865_aqua(klon))
	487	ALLOCATE( aod865_tr2_aqua(klon))
	488	ALLOCATE( aod865_ss_aqua(klon))
	489	ALLOCATE( aod865_dust_aqua(klon))
	490	ALLOCATE( aod865_dustsco_aqua(klon))
	491
	492	ALLOCATE( sconc01(klon) )
	493	ALLOCATE( trm01(klon) )
	494	ALLOCATE( sconc02(klon) )
	495	ALLOCATE( trm02(klon) )
	496	ALLOCATE( sconc03(klon) )
	497	ALLOCATE( trm03(klon) )
	498	ALLOCATE( sconc04(klon) )
	499	ALLOCATE( trm04(klon) )
	500	ALLOCATE( sconc05(klon) )
	501	ALLOCATE( trm05(klon) )
	502
	503	! Lessivage
	504	ALLOCATE( flux01(klon) )
	505	ALLOCATE( flux02(klon) )
	506	ALLOCATE( flux03(klon) )
	507	ALLOCATE( flux04(klon) )
	508	ALLOCATE( flux05(klon) )
	509	ALLOCATE( ds01(klon) )
	510	ALLOCATE( ds02(klon) )
	511	ALLOCATE( ds03(klon) )
	512	ALLOCATE( ds04(klon) )
	513	ALLOCATE( ds05(klon) )
	514	ALLOCATE( dh01(klon) )
	515	ALLOCATE( dh02(klon) )
	516	ALLOCATE( dh03(klon) )
	517	ALLOCATE( dh04(klon) )
	518	ALLOCATE( dh05(klon) )
	519	ALLOCATE( dtrconv01(klon) )
	520	ALLOCATE( dtrconv02(klon) )
	521	ALLOCATE( dtrconv03(klon) )
	522	ALLOCATE( dtrconv04(klon) )
	523	ALLOCATE( dtrconv05(klon) )
	524	ALLOCATE( dtherm01(klon) )
	525	ALLOCATE( dtherm02(klon) )
	526	ALLOCATE( dtherm03(klon) )
	527	ALLOCATE( dtherm04(klon) )
	528	ALLOCATE( dtherm05(klon) )
	529	ALLOCATE( dhkecv01(klon) )
	530	ALLOCATE( dhkecv02(klon) )
	531	ALLOCATE( dhkecv03(klon) )
	532	ALLOCATE( dhkecv04(klon) )
	533	ALLOCATE( dhkecv05(klon) )
	534	ALLOCATE( d_tr_ds01(klon) )
	535	ALLOCATE( d_tr_ds02(klon) )
	536	ALLOCATE( d_tr_ds03(klon) )
	537	ALLOCATE( d_tr_ds04(klon) )
	538	ALLOCATE( d_tr_ds05(klon) )
	539	ALLOCATE( dhkelsc01(klon) )
	540	ALLOCATE( dhkelsc02(klon) )
	541	ALLOCATE( dhkelsc03(klon) )
	542	ALLOCATE( dhkelsc04(klon) )
	543	ALLOCATE( dhkelsc05(klon) )
	544	ALLOCATE( d_tr_cv01(klon,klev))
	545	ALLOCATE( d_tr_cv02(klon,klev))
	546	ALLOCATE( d_tr_cv03(klon,klev))
	547	ALLOCATE( d_tr_cv04(klon,klev))
	548	ALLOCATE( d_tr_cv05(klon,klev))
	549	ALLOCATE( d_tr_trsp01(klon,klev))
	550	ALLOCATE( d_tr_trsp02(klon,klev))
	551	ALLOCATE( d_tr_trsp03(klon,klev))
	552	ALLOCATE( d_tr_trsp04(klon,klev))
	553	ALLOCATE( d_tr_trsp05(klon,klev))
	554	ALLOCATE( d_tr_sscav01(klon,klev))
	555	ALLOCATE( d_tr_sscav02(klon,klev))
	556	ALLOCATE( d_tr_sscav03(klon,klev))
	557	ALLOCATE( d_tr_sscav04(klon,klev))
	558	ALLOCATE( d_tr_sscav05(klon,klev))
	559	ALLOCATE( d_tr_sat01(klon,klev))
	560	ALLOCATE( d_tr_sat02(klon,klev))
	561	ALLOCATE( d_tr_sat03(klon,klev))
	562	ALLOCATE( d_tr_sat04(klon,klev))
	563	ALLOCATE( d_tr_sat05(klon,klev))
	564	ALLOCATE( d_tr_uscav01(klon,klev))
	565	ALLOCATE( d_tr_uscav02(klon,klev))
	566	ALLOCATE( d_tr_uscav03(klon,klev))
	567	ALLOCATE( d_tr_uscav04(klon,klev))
	568	ALLOCATE( d_tr_uscav05(klon,klev))
	569	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	570	ALLOCATE( d_tr_insc01(klon,klev))
	571	ALLOCATE( d_tr_insc02(klon,klev))
	572	ALLOCATE( d_tr_insc03(klon,klev))
	573	ALLOCATE( d_tr_insc04(klon,klev))
	574	ALLOCATE( d_tr_insc05(klon,klev))
	575	ALLOCATE( d_tr_bcscav01(klon,klev))
	576	ALLOCATE( d_tr_bcscav02(klon,klev))
	577	ALLOCATE( d_tr_bcscav03(klon,klev))
	578	ALLOCATE( d_tr_bcscav04(klon,klev))
	579	ALLOCATE( d_tr_bcscav05(klon,klev))
	580	ALLOCATE( d_tr_evapls01(klon,klev))
	581	ALLOCATE( d_tr_evapls02(klon,klev))
	582	ALLOCATE( d_tr_evapls03(klon,klev))
	583	ALLOCATE( d_tr_evapls04(klon,klev))
	584	ALLOCATE( d_tr_evapls05(klon,klev))
	585	ALLOCATE( d_tr_ls01(klon,klev))
	586	ALLOCATE( d_tr_ls02(klon,klev))
	587	ALLOCATE( d_tr_ls03(klon,klev))
	588	ALLOCATE( d_tr_ls04(klon,klev))
	589	ALLOCATE( d_tr_ls05(klon,klev))
	590
	591	ALLOCATE( d_tr_dyn01(klon,klev))
	592	ALLOCATE( d_tr_dyn02(klon,klev))
	593	ALLOCATE( d_tr_dyn03(klon,klev))
	594	ALLOCATE( d_tr_dyn04(klon,klev))
	595	ALLOCATE( d_tr_dyn05(klon,klev))
	596
	597	ALLOCATE( d_tr_cl01(klon,klev))
	598	ALLOCATE( d_tr_cl02(klon,klev))
	599	ALLOCATE( d_tr_cl03(klon,klev))
	600	ALLOCATE( d_tr_cl04(klon,klev))
	601	ALLOCATE( d_tr_cl05(klon,klev))
	602	ALLOCATE( d_tr_th01(klon,klev))
	603	ALLOCATE( d_tr_th02(klon,klev))
	604	ALLOCATE( d_tr_th03(klon,klev))
	605	ALLOCATE( d_tr_th04(klon,klev))
	606	ALLOCATE( d_tr_th05(klon,klev))
	607
	608	ALLOCATE( sed_ss(klon))
	609	ALLOCATE( sed_dust(klon))
	610	ALLOCATE( sed_dustsco(klon))
	611	ALLOCATE( his_g2pgas(klon))
	612	ALLOCATE( his_g2paer(klon))
	613
	614	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	615	ALLOCATE( sed_ss3D(klon,klev))
	616	ALLOCATE( sed_dust3D(klon,klev))
	617	ALLOCATE( sed_dustsco3D(klon,klev))
	618	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	619
	620	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	621	! histrac_spl
	622	!
	623	ALLOCATE( fluxbb(klon))
	624	ALLOCATE( fluxff(klon))
	625	ALLOCATE( fluxbcbb(klon))
	626	ALLOCATE( fluxbcff(klon))
	627	ALLOCATE( fluxbcnff(klon))
	628	ALLOCATE( fluxbcba(klon))
	629	ALLOCATE( fluxbc(klon))
	630	ALLOCATE( fluxombb(klon))
	631	ALLOCATE( fluxomff(klon))
	632	ALLOCATE( fluxomnff(klon))
	633	ALLOCATE( fluxomba(klon))
	634	ALLOCATE( fluxomnat(klon))
	635	ALLOCATE( fluxom(klon))
	636	ALLOCATE( fluxh2sff(klon))
	637	ALLOCATE( fluxh2snff(klon))
	638	ALLOCATE( fluxso2ff(klon))
	639	ALLOCATE( fluxso2nff(klon))
	640	ALLOCATE( fluxso2bb(klon))
	641	ALLOCATE( fluxso2vol(klon))
	642	ALLOCATE( fluxso2ba(klon))
	643	ALLOCATE( fluxso2(klon))
	644	ALLOCATE( fluxso4ff(klon))
	645	ALLOCATE( fluxso4nff(klon))
	646	ALLOCATE( fluxso4bb(klon))
	647	ALLOCATE( fluxso4ba(klon))
	648	ALLOCATE( fluxso4(klon))
	649	ALLOCATE( fluxdms(klon))
	650	ALLOCATE( fluxh2sbio(klon))
	651	ALLOCATE( fluxdustec(klon))
	652	ALLOCATE( fluxddfine(klon))
	653	ALLOCATE( fluxddcoa(klon))
	654	ALLOCATE( fluxddsco(klon))
	655	ALLOCATE( fluxdd(klon))
	656	ALLOCATE( fluxssfine(klon))
	657	ALLOCATE( fluxsscoa(klon))
	658	ALLOCATE( fluxss(klon))
	659	ALLOCATE( flux_sparam_ind(klon))
	660	ALLOCATE( flux_sparam_bb(klon))
	661	ALLOCATE( flux_sparam_ff(klon))
	662	ALLOCATE( flux_sparam_ddfine(klon))
	663	ALLOCATE( flux_sparam_ddcoa(klon))
	664	ALLOCATE( flux_sparam_ddsco(klon))
	665	ALLOCATE( flux_sparam_ssfine(klon))
	666	ALLOCATE( flux_sparam_sscoa(klon))
	667	ALLOCATE( u10m_ss(klon))
	668	ALLOCATE( v10m_ss(klon))
	669
[3806]	670	!AS: in phys_output_write_spl, but not in spla_output_write.h
[3786]	671	!------------------------------------------------------
	672	ALLOCATE(d_tr_cl(klon,klev,nbtr))
	673	ALLOCATE(d_tr_th(klon,klev,nbtr))
	674	ALLOCATE(d_tr_cv(klon,klev,nbtr))
	675	ALLOCATE(d_tr_lessi_impa(klon,klev,nbtr))
	676	ALLOCATE(d_tr_lessi_nucl(klon,klev,nbtr))
	677	ALLOCATE(d_tr_insc(klon,klev,nbtr))
	678	ALLOCATE(d_tr_bcscav(klon,klev,nbtr))
	679	ALLOCATE(d_tr_evapls(klon,klev,nbtr))
	680	ALLOCATE(d_tr_ls(klon,klev,nbtr))
	681	ALLOCATE(d_tr_trsp(klon,klev,nbtr))
	682	ALLOCATE(d_tr_sscav(klon,klev,nbtr))
	683	ALLOCATE(d_tr_sat(klon,klev,nbtr))
	684	ALLOCATE(d_tr_uscav(klon,klev,nbtr))
	685
	686	END SUBROUTINE phytracr_spl_out_init
	687
	688	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[2630]	689	SUBROUTINE phytracr_spl_ini(klon,nbreg_ind,nbreg_bb,nbreg_dust,nbreg_wstardust)
	690	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	691
	692
	693	IMPLICIT NONE
	694	INTEGER klon,nbreg_ind,nbreg_bb,nbreg_dust,nbreg_wstardust
	695
	696	ALLOCATE( tsol(klon) )
[3814]	697
	698	!AS: IF permettant le debranchage des coefs de Jeronimo Escribano: fichiers *_meta
	699	! nbreg_* sont initialisés à 1 dans phytracr_spl, if debutphy,
	700	! avant d'appeler la subroutine presente, phytracr_spl_ini
	701	! (phytracr_spl_ini appele readregionsdims2_spl,
	702	! qui lit et fait "bcast" de nbreg_ind,_bb,_dust,_wstardust dans fichiers regions_*_meta)
	703	IF("ASSIM"=="YES") THEN
[2630]	704	fileregionsdimsind='regions_ind_meta'
	705	fileregionsdimsdust='regions_dustacc_meta'
	706	! fileregionsdimsdust='regions_dust_meta'
	707	fileregionsdimsbb='regions_bb_meta'
	708	fileregionsdimswstar='regions_pwstarwake_meta'
	709	call readregionsdims2_spl(nbreg_ind,fileregionsdimsind)
	710	call readregionsdims2_spl(nbreg_dust,fileregionsdimsdust)
	711	call readregionsdims2_spl(nbreg_bb,fileregionsdimsbb)
	712	call readregionsdims2_spl(nbreg_wstardust,fileregionsdimswstar)
[3814]	713	ENDIF ! ASSIM
	714	! fin debranchage
[2630]	715
	716	!readregions_spl()
	717
	718	ALLOCATE(scale_param_ind(nbreg_ind))
	719	ALLOCATE(scale_param_bb(nbreg_bb))
	720	ALLOCATE(scale_param_ff(nbreg_ind))
	721	ALLOCATE(scale_param_dustacc(nbreg_dust))
	722	ALLOCATE(scale_param_dustcoa(nbreg_dust))
	723	ALLOCATE(scale_param_dustsco(nbreg_dust))
	724	ALLOCATE(param_wstarBLperregion(nbreg_wstardust))
	725	ALLOCATE(param_wstarWAKEperregion(nbreg_wstardust))
	726	ALLOCATE( dust_ec(klon) )
	727	ALLOCATE( u10m_ec(klon) )
	728	ALLOCATE( v10m_ec(klon) )
	729	ALLOCATE( lmt_so2volc_cont(klon) )
	730	ALLOCATE( lmt_altvolc_cont(klon) )
	731	ALLOCATE( lmt_so2volc_expl(klon) )
	732	ALLOCATE( lmt_altvolc_expl(klon) )
	733	ALLOCATE( lmt_so2ff_l(klon) )
	734	ALLOCATE( lmt_so2ff_h(klon) )
	735	ALLOCATE( lmt_so2nff(klon) )
	736	ALLOCATE( lmt_so2ba(klon) )
	737	ALLOCATE( lmt_so2bb_l(klon) )
	738	ALLOCATE( lmt_so2bb_h(klon) )
	739	ALLOCATE( lmt_dmsconc(klon) )
	740	ALLOCATE( lmt_dmsbio(klon) )
	741	ALLOCATE( lmt_h2sbio(klon) )
	742	ALLOCATE( lmt_bcff(klon) )
	743	ALLOCATE( lmt_bcnff(klon) )
	744	ALLOCATE( lmt_bcbb_l(klon) )
	745	ALLOCATE( lmt_bcbb_h(klon) )
	746	ALLOCATE( lmt_bcba(klon) )
	747	ALLOCATE( lmt_omff(klon) )
	748	ALLOCATE( lmt_omnff(klon) )
	749	ALLOCATE( lmt_ombb_l(klon) )
	750	ALLOCATE( lmt_ombb_h(klon) )
	751	ALLOCATE( lmt_omnat(klon) )
	752	ALLOCATE( lmt_omba(klon) )
	753	ALLOCATE(lmt_sea_salt(klon,ss_bins))
	754
	755
	756
	757
	758	!temporal hardcoded null inicialization of assimilation emmision factors
[3814]	759	!AS: scale_param sont ensuite lus dans modvalues.nc
	760	! par la subroutine read_scalenc, appelee par readscaleparamsnc_spl
[2630]	761	scale_param_ssacc=1.
	762	scale_param_sscoa=1.
	763	scale_param_ind(:)=1.
	764	scale_param_bb(:)=1.
	765	scale_param_ff(:)=1.
	766	scale_param_dustacc(:)=1.
	767	scale_param_dustcoa(:)=1.
	768	scale_param_dustsco(:)=1.
	769	param_wstarBLperregion(:)=0.
	770	param_wstarWAKEperregion(:)=0.
	771
	772
	773	RETURN
	774	END SUBROUTINE phytracr_spl_ini
	775
	776
	777
	778
	779	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	780	SUBROUTINE phytracr_spl ( debutphy,lafin,jD_cur,jH_cur,iflag_conv, & ! I
	781	pdtphys,ftsol, & ! I
	782	t_seri,q_seri,paprs,pplay,RHcl, & ! I
	783	pmfu, pmfd, pen_u, pde_u, pen_d, pde_d, & ! I
	784	coefh, cdragh, cdragm, yu1, yv1, & ! I
	785	u_seri, v_seri, rlat,rlon, & ! I
	786	pphis,pctsrf,pmflxr,pmflxs,prfl,psfl, & ! I
	787	da,phi,phi2,d1a,dam,mp,ep,sigd,sij,clw,elij, & ! I
	788	epmlmMm,eplaMm,upwd,dnwd,itop_con,ibas_con, & ! I
	789	evapls,wdtrainA, wdtrainM,wght_cvfd, & ! I
	790	fm_therm, entr_therm, rneb, & ! I
	791	beta_fisrt,beta_v1, & ! I
	792	zu10m,zv10m,wstar,ale_bl,ale_wake, & ! I
	793	d_tr_dyn,tr_seri) ! O
	794	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	795
	796	USE mod_grid_phy_lmdz
	797	USE mod_phys_lmdz_para
	798	USE IOIPSL
	799	USE dimphy
	800	USE infotrac
	801	USE indice_sol_mod
	802	USE write_field_phy
	803
	804
	805	USE mod_phys_lmdz_transfert_para
	806
	807	USE phys_cal_mod, only: jD_1jan,year_len, mth_len, days_elapsed, jh_1jan, year_cur, &
	808	mth_cur, phys_cal_update
	809
	810	!
	811	IMPLICIT none
	812	!
	813
	814	!======================================================================
	815	! Auteur(s) FH
	816	! Objet: Moniteur general des tendances traceurs
	817	!
	818	! Remarques en vrac:
	819	! ------------------
	820	! 1/ le call phytrac se fait avec nqmax-2 donc nous avons bien
	821	! les vrais traceurs (nbtr) dans phytrac (pas la vapeur ni eau liquide)
	822	!======================================================================
	823	#include "dimensions.h"
	824	#include "chem.h"
	825	#include "chem_spla.h"
	826	#include "YOMCST.h"
	827	#include "YOETHF.h"
	828	#include "paramet.h"
	829	#include "thermcell.h"
	830
	831	!======================================================================
	832
	833	! Arguments:
	834	!
	835	! EN ENTREE:
	836	! ==========
	837	!
	838	! divers:
	839	! -------
	840	!
	841	real,intent(in) :: pdtphys ! pas d'integration pour la physique (seconde)
	842	REAL, intent(in):: jD_cur, jH_cur
	843	real, intent(in) :: ftsol(klon,nbsrf) ! temperature du sol par type
	844	real, intent(in) :: t_seri(klon,klev) ! temperature
	845	real, intent(in) :: u_seri(klon,klev) ! vent
	846	real , intent(in) :: v_seri(klon,klev) ! vent
	847	real , intent(in) :: q_seri(klon,klev) ! vapeur d eau kg/kg
	848
	849	LOGICAL, INTENT(IN) :: lafin
	850
	851	real tr_seri(klon,klev,nbtr) ! traceur
	852	real tmp_var(klon,klev) ! auxiliary variable to replace traceur
	853	real tmp_var2(klon,nbtr) ! auxiliary variable to replace source
	854	real tmp_var3(klon,klev,nbtr) ! auxiliary variable 3D
	855	real dummy1d ! JE auxiliary variable
	856	real aux_var2(klon) ! auxiliary variable to replace traceur
	857	real aux_var3(klon,klev) ! auxiliary variable to replace traceur
	858	real d_tr(klon,klev,nbtr) ! traceur tendance
	859	real sconc_seri(klon,nbtr) ! surface concentration of traceur
	860	!
	861	integer nbjour
	862	save nbjour
	863	!$OMP THREADPRIVATE(nbjour)
	864	!
	865	INTEGER masque_aqua_cur(klon)
	866	INTEGER masque_terra_cur(klon)
	867	INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: masque_aqua !mask for 1 day
	868	INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: masque_terra !
	869	!$OMP THREADPRIVATE(masque_aqua,masque_terra)
	870
	871	INTEGER, SAVE :: nbreg_dust, nbreg_ind, nbreg_bb, nbreg_ss,nbreg_wstardust
	872	!$OMP THREADPRIVATE(nbreg_dust, nbreg_ind, nbreg_bb,nbreg_ss,nbreg_wstardust)
	873
	874
	875
	876	REAL lmt_dms(klon) ! emissions de dms
	877
	878	!JE20150518<<
	879	REAL, DIMENSION(klon_glo) :: aod550_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	880	REAL, DIMENSION(klon_glo) :: aod550_tr2_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	881	REAL, DIMENSION(klon_glo) :: aod550_ss_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	882	REAL, DIMENSION(klon_glo) :: aod550_dust_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	883	REAL, DIMENSION(klon_glo) :: aod550_dustsco_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	884	REAL, DIMENSION(klon_glo) :: aod670_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	885	REAL, DIMENSION(klon_glo) :: aod670_tr2_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	886	REAL, DIMENSION(klon_glo) :: aod670_ss_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	887	REAL, DIMENSION(klon_glo) :: aod670_dust_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	888	REAL, DIMENSION(klon_glo) :: aod670_dustsco_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	889	REAL, DIMENSION(klon_glo) :: aod865_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	890	REAL, DIMENSION(klon_glo) :: aod865_tr2_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	891	REAL, DIMENSION(klon_glo) :: aod865_ss_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	892	REAL, DIMENSION(klon_glo) :: aod865_dust_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	893	REAL, DIMENSION(klon_glo) :: aod865_dustsco_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
	894
	895	REAL, DIMENSION(klon_glo) :: aod550_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	896	REAL, DIMENSION(klon_glo) :: aod550_tr2_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	897	REAL, DIMENSION(klon_glo) :: aod550_ss_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	898	REAL, DIMENSION(klon_glo) :: aod550_dust_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	899	REAL, DIMENSION(klon_glo) :: aod550_dustsco_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	900	REAL, DIMENSION(klon_glo) :: aod670_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	901	REAL, DIMENSION(klon_glo) :: aod670_tr2_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	902	REAL, DIMENSION(klon_glo) :: aod670_ss_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	903	REAL, DIMENSION(klon_glo) :: aod670_dust_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	904	REAL, DIMENSION(klon_glo) :: aod670_dustsco_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	905	REAL, DIMENSION(klon_glo) :: aod865_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	906	REAL, DIMENSION(klon_glo) :: aod865_tr2_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	907	REAL, DIMENSION(klon_glo) :: aod865_ss_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	908	REAL, DIMENSION(klon_glo) :: aod865_dust_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	909	REAL, DIMENSION(klon_glo) :: aod865_dustsco_aqua_glo ! AOD at aqua overpass time ( 13.30 local hour)
	910	!!!!!!!!!!!!!
	911	!JE20150518>>
	912
	913
	914
	915
	916	real , intent(in) :: paprs(klon,klev+1) ! pression pour chaque inter-couche (en Pa)
	917	real , intent(in) :: pplay(klon,klev) ! pression pour le mileu de chaque couche (en Pa)
	918	real , intent(in) :: RHcl(klon,klev) ! humidite relativen ciel clair
	919	real znivsig(klev) ! indice des couches
	920	real paire(klon)
	921	real, intent(in) :: pphis(klon)
	922	real, intent(in) :: pctsrf(klon,nbsrf)
	923	logical , intent(in) :: debutphy ! le flag de l'initialisation de la physique
	924	!
	925	! Scaling Parameters:
	926	! ----------------------
	927	!
	928	CHARACTER*50 c_Directory
	929	CHARACTER*80 c_FileName1
	930	CHARACTER*80 c_FileName2
	931	CHARACTER*130 c_FullName1
	932	CHARACTER*130 c_FullName2
	933	INTEGER :: xidx, yidx
	934	INTEGER,DIMENSION(klon) :: mask_bbreg
	935	INTEGER,DIMENSION(klon) :: mask_ffso2reg
	936	INTEGER :: aux_mask1
	937	INTEGER :: aux_mask2
[3814]	938	INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: iregion_so4 !Defines regions for SO4 ; AS: PAS UTILISE!
[2630]	939	INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: iregion_ind !Defines regions for SO2, BC & OM
	940	INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: iregion_bb !Defines regions for SO2, BC & OM
	941	INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: iregion_dust !Defines dust regions
	942	INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: iregion_wstardust !Defines dust regions
	943	!$OMP THREADPRIVATE(iregion_so4,iregion_ind,iregion_bb,iregion_dust,iregion_wstardust)
	944
	945	! Emissions:
	946
	947	!
	948	!---------------------------- SEA SALT & DUST emissions ------------------------
	949	REAL lmt_sea_salt(klon,ss_bins) !Sea salt 0.03-8.0 um
	950	REAL u10m_ec1(klon),v10m_ec1(klon)
	951	REAL u10m_ec2(klon),v10m_ec2(klon),dust_ec2(klon)
	952	REAL dust_ec(klon)
	953	! new dust emission chimere je20140522
	954	REAL,DIMENSION(klon),INTENT(IN) :: zu10m
	955	REAL,DIMENSION(klon),INTENT(IN) :: zv10m
	956	REAL,DIMENSION(klon),INTENT(IN) :: wstar,ale_bl,ale_wake
	957
	958
	959	!
	960	! Rem : nbtr : nombre de vrais traceurs est defini dans dimphy.h
	961
	962	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	963	!Dynamique
	964	!--------
	965	REAL,DIMENSION(klon,klev,nbtr),INTENT(IN) :: d_tr_dyn
	966
	967	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	968	! convection:
	969	! -----------
	970	!
	971	REAL , intent(in) :: pmfu(klon,klev) ! flux de masse dans le panache montant
	972	REAL , intent(in) :: pmfd(klon,klev) ! flux de masse dans le panache descendant
	973	REAL, intent(in) :: pen_u(klon,klev) ! flux entraine dans le panache montant
	974	REAL, intent(in) :: pde_u(klon,klev) ! flux detraine dans le panache montant
	975	REAL, intent(in) :: pen_d(klon,klev) ! flux entraine dans le panache descendant
	976	REAL, intent(in) :: pde_d(klon,klev) ! flux detraine dans le panache descendant
	977	!
	978	! Convection KE scheme:
	979	! ---------------------
	980	!
	981	!! Variables pour le lessivage convectif
	982	REAL,DIMENSION(klon,klev),INTENT(IN) :: da
	983	REAL,DIMENSION(klon,klev,klev),INTENT(IN):: phi
	984	REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: phi2
	985	REAL,DIMENSION(klon,klev),INTENT(IN) :: d1a,dam
	986	REAL,DIMENSION(klon,klev),INTENT(IN) :: mp
	987	REAL,DIMENSION(klon,klev),INTENT(IN) :: upwd ! saturated
	988	! updraft mass flux
	989	REAL,DIMENSION(klon,klev),INTENT(IN) :: dnwd ! saturated
	990	! downdraft mass flux
	991	INTEGER,DIMENSION(klon),INTENT(IN) :: itop_con
	992	INTEGER,DIMENSION(klon),INTENT(IN) :: ibas_con
	993	REAL,DIMENSION(klon,klev) :: evapls
	994	REAL,DIMENSION(klon,klev),INTENT(IN) :: wdtrainA
	995	REAL,DIMENSION(klon,klev),INTENT(IN) :: wdtrainM
	996
	997
	998	REAL,DIMENSION(klon,klev),INTENT(IN) :: ep
	999	REAL,DIMENSION(klon),INTENT(IN) :: sigd
	1000	REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: sij
	1001	REAL,DIMENSION(klon,klev),INTENT(IN) :: clw
	1002	REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: elij
	1003	REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: epmlmMm
	1004	REAL,DIMENSION(klon,klev),INTENT(IN) :: eplaMm
	1005	REAL,DIMENSION(klon,klev),INTENT(IN) :: wght_cvfd !RL
	1006
	1007
	1008	! KE: Tendances de traceurs (Td) et flux de traceurs:
	1009	! ------------------------
	1010	REAL,DIMENSION(klon,klev) :: Mint
	1011	REAL,DIMENSION(klon,klev,nbtr) :: zmfd1a
	1012	REAL,DIMENSION(klon,klev,nbtr) :: zmfdam
	1013	REAL,DIMENSION(klon,klev,nbtr) :: zmfphi2
	1014
	1015	! !tra dans pluie LS a la surf.
	1016	! outputs for cvltr_spl
	1017	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_cv_o
	1018	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_trsp_o
	1019	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_sscav_o
	1020	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_sat_o
	1021	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_uscav_o
	1022	!!!!!!!!!!!!!!!!!
	1023	!!!!!!!!!!!!!!!!!
	1024	!!!!!!!!!!!!!!!!!
	1025	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_insc_o
	1026	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_bcscav_o
	1027	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_evapls_o
	1028	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_ls_o
	1029	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_dyn_o
	1030	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_cl_o
	1031	REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_th_o
	1032	!!!!!!!!!!!!!!!!!
	1033	!!!!!!!!!!!!!!!!!
	1034	!!!!!!!!!!!!!!!!!
	1035
	1036	!$OMP THREADPRIVATE(d_tr_cv_o,d_tr_trsp_o,d_tr_sscav_o,d_tr_sat_o,d_tr_uscav_o)
	1037	!$OMP THREADPRIVATE(d_tr_insc_o,d_tr_bcscav_o,d_tr_evapls_o,d_tr_ls_o)
	1038	!$OMP THREADPRIVATE(d_tr_dyn_o,d_tr_cl_o,d_tr_th_o)
	1039
	1040
	1041	INTEGER :: nsplit
	1042	!
	1043
	1044
	1045
	1046	!
	1047	! Lessivage
	1048	! ---------
	1049	!
	1050	REAL, intent(in) :: pmflxr(klon,klev+1), pmflxs(klon,klev+1) !--convection
	1051	REAL, intent(in) :: prfl(klon,klev+1), psfl(klon,klev+1) !--large-scale
	1052	REAL :: ql_incl ! contenu en eau liquide nuageuse dans le nuage ! ql_incl=oliq/rneb
	1053	REAL :: ql_incloud_ref ! ref value of in-cloud condensed water content
	1054
	1055	REAL,DIMENSION(klon,klev),INTENT(IN) :: rneb ! fraction nuageuse (grande echelle)
	1056	!
	1057
	1058	REAL,DIMENSION(klon,klev) :: beta_fisrt ! taux de conversion
	1059	! ! de l'eau cond (de fisrtilp)
	1060	REAL,DIMENSION(klon,klev) :: beta_v1 ! -- (originale version)
	1061	INTEGER,SAVE :: iflag_lscav_omp,iflag_lscav
	1062	!$OMP THREADPRIVATE(iflag_lscav_omp,iflag_lscav)
	1063
	1064
	1065
	1066
	1067	!Thermiques:
	1068	!----------
	1069	REAL,DIMENSION(klon,klev+1),INTENT(IN) :: fm_therm
	1070	REAL,DIMENSION(klon,klev),INTENT(IN) :: entr_therm
	1071
	1072
	1073	!
	1074	! Couche limite:
	1075	! --------------
	1076	!
	1077	REAL , intent(in) :: coefh(klon,klev) ! coeff melange CL
	1078	REAL , intent(in) :: cdragh(klon), cdragm(klon)
	1079	REAL, intent(in) :: yu1(klon) ! vent dans la 1iere couche
	1080	REAL, intent(in) :: yv1(klon) ! vent dans la 1iere couche
	1081	!
	1082	!
	1083	!----------------------------------------------------------------------
	1084	REAL his_ds(klon,nbtr)
	1085	REAL his_dh(klon,nbtr)
	1086	REAL his_dhlsc(klon,nbtr) ! in-cloud scavenging lsc
	1087	REAL his_dhcon(klon,nbtr) ! in-cloud scavenging con
	1088	REAL his_dhbclsc(klon,nbtr) ! below-cloud scavenging lsc
	1089	REAL his_dhbccon(klon,nbtr) ! below-cloud scavenging con
	1090	REAL trm(klon,nbtr)
	1091	!
	1092	REAL u10m_ec(klon), v10m_ec(klon)
	1093	!
	1094	REAL his_th(klon,nbtr)
	1095	REAL his_dhkecv(klon,nbtr)
	1096	REAL his_dhkelsc(klon,nbtr)
	1097
	1098
	1099	!
	1100	! Coordonnees
	1101	! -----------
	1102	!
	1103	REAL, intent(in) :: rlat(klon) ! latitudes pour chaque point
	1104	REAL, intent(in) :: rlon(klon) ! longitudes pour chaque point
	1105	!
	1106	INTEGER i, k, it, j, ig
	1107	!
	1108	! DEFINITION OF DIAGNOSTIC VARIABLES
	1109	!
	1110	REAL diag_trm(nbtr), diag_drydep(nbtr)
	1111	REAL diag_wetdep(nbtr), diag_cvtdep(nbtr)
	1112	REAL diag_emissn(nbtr), diag_g2part
	1113	REAL diag_sedimt
	1114	REAL trm_aux(nbtr), src_aux(nbtr)
	1115	!
	1116	! Variables locales pour effectuer les appels en serie
	1117	!----------------------------------------------------
	1118	REAL source_tr(klon,nbtr)
	1119	REAL flux_tr(klon,nbtr)
	1120	REAL m_conc(klon,klev)
	1121	REAL henry(nbtr) !--cste de Henry mol/l/atm
	1122	REAL kk(nbtr) !--coefficient de var avec T (K)
	1123	REAL alpha_r(nbtr)!--coefficient d'impaction pour la pluie
	1124	REAL alpha_s(nbtr)!--coefficient d'impaction pour la neige
	1125	REAL vdep_oce(nbtr), vdep_sic(nbtr)
	1126	REAL vdep_ter(nbtr), vdep_lic(nbtr)
	1127	REAL ccntrAA_spla(nbtr)
	1128	REAL ccntrENV_spla(nbtr)
	1129	REAL coefcoli_spla(nbtr)
	1130	REAL dtrconv(klon,nbtr)
	1131	REAL zrho(klon,klev), zdz(klon,klev)
	1132	REAL zalt(klon,klev)
	1133	REAL,DIMENSION(klon,klev) :: zmasse ! densité atmosphérique
	1134	! . Kg/m2
	1135	REAL,DIMENSION(klon,klev) :: ztra_th
	1136	REAL qmin, qmax, aux
	1137	! PARAMETER (qmin=0.0, qmax=1.e33)
	1138	PARAMETER (qmin=1.e33, qmax=-1.e33)
	1139
	1140	! Variables to save data into file
	1141	!----------------------------------
	1142
	1143	CHARACTER*2 str2
[3806]	1144	!!AS: LOGICAL ok_histrac
	1145	!!!JE2014124 PARAMETER (ok_histrac=.true.)
	1146	!! PARAMETER (ok_histrac=.false.)
[2630]	1147	INTEGER ndex2d(iim(jjm+1)), ndex3d(iim(jjm+1)*klev)
	1148	INTEGER nhori1, nhori2, nhori3, nhori4, nhori5, nvert
	1149	INTEGER nid_tra1, nid_tra2, nid_tra3, nid_tra4, nid_tra5
	1150	SAVE nid_tra1, nid_tra2, nid_tra3, nid_tra4, nid_tra5
	1151	!$OMP THREADPRIVATE(nid_tra1, nid_tra2, nid_tra3, nid_tra4, nid_tra5)
	1152	INTEGER itra
	1153	SAVE itra ! compteur pour la physique
	1154	!$OMP THREADPRIVATE(itra)
	1155	INTEGER ecrit_tra, ecrit_tra_h, ecrit_tra_m
	1156	SAVE ecrit_tra, ecrit_tra_h, ecrit_tra_m
	1157	!$OMP THREADPRIVATE(ecrit_tra, ecrit_tra_h, ecrit_tra_m)
	1158	REAL presnivs(klev) ! pressions approximat. des milieux couches ( en PA)
	1159	REAL zx_tmp_2d(iim,jjm+1), zx_tmp_3d(iim,jjm+1,klev)
	1160	REAL zx_tmp_fi2d(klon), zx_tmp_fi3d(klon, klev)
	1161	REAL zx_lon_glo(nbp_lon,nbp_lat), zx_lat_glo(nbp_lon,nbp_lat)
	1162	REAL zsto, zout, zout_h, zout_m, zjulian
	1163
	1164	!------Molar Masses
	1165	REAL masse(nbtr)
	1166	!
	1167	REAL fracso2emis !--fraction so2 emis en so2
	1168	PARAMETER (fracso2emis=0.95)
	1169	REAL frach2sofso2 !--fraction h2s from so2
	1170	PARAMETER (frach2sofso2=0.0426)
	1171	!
	1172	! Controles
	1173	!-------------
	1174	LOGICAL convection,lessivage,lminmax,lcheckmass
	1175	DATA convection,lessivage,lminmax,lcheckmass &
	1176	/.true.,.true.,.true.,.false./
	1177	!
	1178	REAL xconv(nbtr)
	1179	!
	1180	LOGICAL anthropo, bateau, edgar
	1181	DATA anthropo,bateau,edgar/.true.,.true.,.true./
	1182	!
	1183	!c bc_source
	1184	INTEGER kminbc, kmaxbc
	1185	!JE20150715 PARAMETER (kminbc=3, kmaxbc=5)
	1186	PARAMETER (kminbc=4, kmaxbc=7)
	1187	!
	1188	REAL tr1_cont, tr2_cont, tr3_cont, tr4_cont
	1189	!
	1190	! JE for updating in cltrac
	1191	REAL,DIMENSION(klon,klev) :: delp ! epaisseur de couche (Pa)
	1192	!! JE for include gas to particle conversion in output
	1193	! REAL his_g2pgas(klon) ! gastoparticle in gas units (check!)
	1194	! REAL his_g2paer(klon) ! gastoparticle in aerosol units (check!)
	1195	!
	1196	INTEGER ,intent(in) :: iflag_conv
	1197	LOGICAL iscm3 ! debug variable. for checkmass ! JE
	1198
	1199	!------------------------------------------------------------------------
	1200	! only to compute time consumption of each process
	1201	!----
	1202	INTEGER clock_start,clock_end,clock_rate,clock_start_spla
	1203	INTEGER clock_end_outphytracr,clock_start_outphytracr
	1204	INTEGER ti_init,dife,ti_inittype,ti_inittwrite
	1205	INTEGER ti_spla,ti_emis,ti_depo,ti_cltr,ti_ther
	1206	INTEGER ti_sedi,ti_gasp,ti_wetap,ti_cvltr,ti_lscs,ti_brop,ti_outs
	1207	INTEGER ti_nophytracr,clock_per_max
	1208	REAL tia_init,tia_inittype,tia_inittwrite
	1209	REAL tia_spla,tia_emis,tia_depo,tia_cltr,tia_ther
	1210	REAL tia_sedi,tia_gasp,tia_wetap,tia_cvltr,tia_lscs
	1211	REAL tia_brop,tia_outs
	1212	REAL tia_nophytracr
	1213
	1214	SAVE tia_init,tia_inittype,tia_inittwrite
	1215	SAVE tia_spla,tia_emis,tia_depo,tia_cltr,tia_ther
	1216	SAVE tia_sedi,tia_gasp,tia_wetap,tia_cvltr,tia_lscs
	1217	SAVE tia_brop,tia_outs
	1218	SAVE ti_nophytracr
	1219	SAVE tia_nophytracr
	1220	SAVE clock_end_outphytracr,clock_start_outphytracr
	1221	SAVE clock_per_max
	1222	LOGICAL logitime
	1223	!$OMP THREADPRIVATE(tia_init,tia_inittype,tia_inittwrite)
	1224	!$OMP THREADPRIVATE(tia_spla,tia_emis,tia_depo,tia_cltr,tia_ther)
	1225	!$OMP THREADPRIVATE(tia_sedi,tia_gasp,tia_wetap,tia_cvltr,tia_lscs)
	1226	!$OMP THREADPRIVATE(tia_brop,tia_outs)
	1227	!$OMP THREADPRIVATE(ti_nophytracr)
	1228	!$OMP THREADPRIVATE(tia_nophytracr)
	1229	!$OMP THREADPRIVATE(clock_end_outphytracr,clock_start_outphytracr)
	1230	!$OMP THREADPRIVATE(clock_per_max)
	1231
	1232	! utils parallelization
	1233	REAL :: auxklon_glo(klon_glo)
	1234	INTEGER :: iauxklon_glo(klon_glo)
	1235	REAL, DIMENSION(klon_glo,nbp_lev) :: auxklonnbp_lev
	1236	REAL, DIMENSION(klon_glo,nbp_lev,nbtr) :: auxklonklevnbtr_glo
	1237	REAL,DIMENSION(nbp_lon,nbp_lat) :: zx_tmp_2d_glo
	1238	REAL,DIMENSION(nbp_lon,nbp_lat,nbp_lev) :: zx_tmp_3d_glo
	1239	REAL,DIMENSION(klon_glo) :: zx_tmp_fi2d_glo
	1240	REAL,DIMENSION(klon_glo , nbp_lev) :: zx_tmp_fi3d_glo
	1241	REAL,DIMENSION(klon_glo,nbtr) :: auxklonnbtr_glo
	1242
	1243
	1244
	1245	source_tr=0.
	1246
	1247
	1248
	1249	if (debutphy) then
[2632]	1250	#ifdef IOPHYS_DUST
[2630]	1251	CALL iophys_ini
[2632]	1252	#endif
[2630]	1253	nbreg_ind=1
	1254	nbreg_bb=1
	1255	nbreg_dust=1
	1256	nbreg_wstardust=1
	1257	CALL phytracr_spl_ini(klon,nbreg_ind,nbreg_bb,nbreg_dust,nbreg_wstardust)
	1258	endif
	1259
	1260
[2632]	1261	#ifdef IOPHYS_DUST
	1262	do it=1,nbtr
	1263	write(str2,'(i2.2)') it
	1264	call iophys_ecrit('TRA'//str2,klev,'SOURCE','',tr_seri(:,:,it))
	1265	enddo
	1266	#endif
	1267
[2630]	1268
	1269
	1270
	1271	ijulday=jD_cur-jD_1jan+1
	1272	nbjour = 1
	1273
	1274	paramname_ind='ind'
	1275	paramname_bb='bb'
	1276	paramname_ff='ind'
	1277	paramname_dustacc='dustacc'
	1278	paramname_dustcoa='dustcoasco'
	1279	paramname_dustsco='dustcoasco'
	1280	! paramname_dustacc='dust'
	1281	! paramname_dustcoa='dust'
	1282	! paramname_dustsco='dust'
	1283	paramname_wstarBL='pwstarbl'
	1284	paramname_wstarWAKE='pwstarwake'
	1285	paramname_ssacc='ssacc'
	1286	paramname_sscoa='sscoa'
	1287
	1288	filescaleparams='modvalues.nc'
[3814]	1289	!AS: debranchage de lecture des coefs d'assmilation de Jeronimo Escribano
	1290	IF("ASSIM"=="YES") THEN
	1291	CALL readscaleparamsnc_spl(scale_param_ind, &
[2630]	1292	nbreg_ind, paramname_ind, &
	1293	scale_param_ff, nbreg_ind,paramname_ff, &
	1294	scale_param_bb, nbreg_bb,paramname_bb, &
	1295	scale_param_dustacc, nbreg_dust,paramname_dustacc, &
	1296	scale_param_dustcoa, nbreg_dust,paramname_dustcoa, &
	1297	scale_param_dustsco, nbreg_dust,paramname_dustsco, &
	1298	param_wstarBLperregion, nbreg_wstardust, paramname_wstarBL, &
	1299	param_wstarWAKEperregion, nbreg_wstardust, paramname_wstarWAKE, &
	1300	scale_param_ssacc , paramname_ssacc, &
	1301	scale_param_sscoa , paramname_sscoa, &
	1302	filescaleparams,ijulday,jH_cur, pdtphys,debutphy)
[3814]	1303	ENDIF ! ASSIM
	1304	!AS: le commentaire suivant "add seasalt" ne semble pas avoir ete mis en pratique.
	1305	! Des fichiers regions_ssacc et _sscoa existent mais ne semblent pas lus.
	1306	! Ca reste donc aux valeurs initialisées: nbreg_ss=1, scale_param_ss*=1, cf fichiers ss et modvalues
	1307	!! add seasalt
[2630]	1308
	1309	print *,'JE : check scale_params'
	1310
	1311	print *, 'nbreg_ind', nbreg_ind
	1312	print *, 'nbreg_dust', nbreg_dust
	1313	print *, 'nbreg_bb', nbreg_bb
	1314	print *, 'ind', scale_param_ind
	1315	print *, 'dustacc', scale_param_dustacc
	1316	print *, 'dustcoa', scale_param_dustcoa
	1317	print *, 'dustsco', scale_param_dustsco
	1318	print *, 'wstardustBL', param_wstarBLperregion
	1319	print *, 'wstardustWAKE', param_wstarWAKEperregion
	1320	print *, 'ff', scale_param_ff
	1321	print *, 'bb', scale_param_bb
	1322	print *, 'ssacc', scale_param_ssacc
	1323	print *, 'sscoa', scale_param_sscoa
	1324
	1325	print *,'JE: before read_newemissions '
	1326	print *,'JE: jD_cur:',jD_cur,' ijulday:',ijulday,' jH_cur:',jH_cur,' pdtphys:',pdtphys
	1327	print *,'JE: now read_newemissions:'
[3786]	1328	!AS: La ligne suivante fait planter a l'execution : lmt_so2ff_l pas initialise
	1329	! print *,'lmt_so2ff_l AVANT' , MINVAL(lmt_so2ff_l), MAXVAL(lmt_so2ff_l)
[2630]	1330	call read_newemissions(ijulday,jH_cur ,edgar, flag_dms,debutphy, & !I
	1331	pdtphys, lafin, nbjour, pctsrf, & !I
	1332	t_seri, rlat, rlon, & !I
	1333	pmflxr, pmflxs, prfl, psfl, & !I
	1334	u10m_ec, v10m_ec, dust_ec, & !O
	1335	lmt_sea_salt, lmt_so2ff_l, & !O
	1336	lmt_so2ff_h, lmt_so2nff, & !O
	1337	lmt_so2ba, lmt_so2bb_l, lmt_so2bb_h, & !O
	1338	lmt_so2volc_cont, lmt_altvolc_cont, & !O
	1339	lmt_so2volc_expl, lmt_altvolc_expl, & !O
	1340	lmt_dmsbio, lmt_h2sbio, lmt_dmsconc, & !O
	1341	lmt_bcff, lmt_bcnff, lmt_bcbb_l, & !O
	1342	lmt_bcbb_h, lmt_bcba, lmt_omff, & !O
	1343	lmt_omnff, lmt_ombb_l, lmt_ombb_h, & !O
	1344	lmt_omnat, lmt_omba) !O
	1345
	1346
	1347	print *,'Check emissions'
	1348	print *,'lmt_so2ff_l' , MINVAL(lmt_so2ff_l), MAXVAL(lmt_so2ff_l)
	1349	print *,'lmt_so2ff_h' , MINVAL(lmt_so2ff_h), MAXVAL(lmt_so2ff_h)
	1350	print *,'lmt_so2nff' , MINVAL(lmt_so2nff), MAXVAL(lmt_so2nff)
	1351	print *,'lmt_so2ba' , MINVAL(lmt_so2ba), MAXVAL(lmt_so2ba)
	1352	print *,'lmt_so2bb_l' , MINVAL(lmt_so2bb_l), MAXVAL(lmt_so2bb_l)
	1353	print *,'lmt_so2bb_h' , MINVAL(lmt_so2bb_h), MAXVAL(lmt_so2bb_h)
	1354	print *,'lmt_so2volc_cont' , MINVAL(lmt_so2volc_cont), MAXVAL(lmt_so2volc_cont)
	1355	print *,'lmt_altvolc_cont' , MINVAL(lmt_altvolc_cont), MAXVAL(lmt_altvolc_cont)
	1356	print *,'lmt_so2volc_expl' , MINVAL(lmt_so2volc_expl), MAXVAL(lmt_so2volc_expl)
	1357	print *,'lmt_altvolc_expl' , MINVAL(lmt_altvolc_expl), MAXVAL(lmt_altvolc_expl)
	1358	print *,'lmt_dmsbio' , MINVAL(lmt_dmsbio), MAXVAL(lmt_dmsbio)
	1359	print *,'lmt_h2sbio' , MINVAL(lmt_h2sbio), MAXVAL(lmt_h2sbio)
	1360	print *,'lmt_dmsconc' , MINVAL(lmt_dmsconc), MAXVAL(lmt_dmsconc)
	1361	print *,'lmt_bcff' , MINVAL(lmt_bcff), MAXVAL(lmt_bcff)
	1362	print *,'lmt_bcnff' , MINVAL(lmt_bcnff), MAXVAL(lmt_bcnff)
	1363	print *,'lmt_bcbb_l' , MINVAL(lmt_bcbb_l), MAXVAL(lmt_bcbb_l)
	1364	print *,'lmt_bcbb_h' , MINVAL(lmt_bcbb_h), MAXVAL(lmt_bcbb_h)
	1365	print *,'lmt_bcba' , MINVAL(lmt_bcba), MAXVAL(lmt_bcba)
	1366	print *,'lmt_omff' , MINVAL(lmt_omff), MAXVAL(lmt_omff)
	1367	print *,'lmt_omnff' , MINVAL(lmt_omnff), MAXVAL(lmt_omnff)
	1368	print *,'lmt_ombb_l' , MINVAL(lmt_ombb_l), MAXVAL(lmt_ombb_l)
	1369	print *,'lmt_ombb_h' , MINVAL(lmt_ombb_h), MAXVAL(lmt_ombb_h)
	1370	print *,'lmt_omnat' , MINVAL(lmt_omnat), MAXVAL(lmt_omnat)
	1371	print *,'lmt_omba' , MINVAL(lmt_omba), MAXVAL(lmt_omba)
	1372	print *,'JE iflag_con',iflag_conv
	1373
	1374
	1375	!JE_dbg
	1376	do i=1,klon
	1377	tsol(i)=0.0
	1378	do j=1,nbsrf
	1379	tsol(i)=tsol(i)+ftsol(i,j)*pctsrf(i,j)
	1380	enddo
	1381	enddo
	1382
	1383
	1384	!======================================================================
	1385	! INITIALISATIONS
	1386	!======================================================================
	1387	! CALL checknanqfi(da(:,:),1.,-1.,' da_ before
	1388	! . phytracr_inphytracr')
	1389
	1390	!
	1391	! computing time
	1392	! logitime=.true.
	1393	logitime=.false.
	1394	IF (logitime) THEN
	1395	clock_start=0
	1396	clock_end=0
	1397	clock_rate=0
	1398	CALL SYSTEM_CLOCK(COUNT_RATE=clock_rate,COUNT_MAX=clock_per_max)
	1399	CALL SYSTEM_CLOCK(COUNT=clock_start_spla)
	1400	clock_start=clock_start_spla
	1401	clock_end_outphytracr=clock_start_spla
	1402	ENDIF
	1403
	1404
	1405	! Definition of tracers index.
	1406	print*,'OK ON PASSSE BIEN LA'
	1407	CALL minmaxsource(source_tr,qmin,qmax,'A1 maxsource init phytracr')
	1408
	1409
	1410	IF (debutphy) THEN
	1411	id_prec=-1
	1412	id_fine=-1
	1413	id_coss=-1
	1414	id_codu=-1
	1415	id_scdu=-1
	1416	!print *,nbtr
	1417	do it=1,nbtr
[2647]	1418	print *, it, tname(it+nqo)
	1419	if (tname(it+nqo) == 'PREC' ) then
[2630]	1420	id_prec=it
	1421	endif
[2647]	1422	if (tname(it+nqo) == 'FINE' ) then
[2630]	1423	id_fine=it
	1424	endif
[2647]	1425	if (tname(it+nqo) == 'COSS' ) then
[2630]	1426	id_coss=it
	1427	endif
[2647]	1428	if (tname(it+nqo) == 'CODU' ) then
[2630]	1429	id_codu=it
	1430	endif
[2647]	1431	if (tname(it+nqo) == 'SCDU' ) then
[2630]	1432	id_scdu=it
	1433	endif
	1434	enddo
	1435	! check consistency with dust emission scheme:
	1436	if (ok_chimeredust) then
	1437	if (.not.( id_scdu>0 .and. id_codu>0 .and. id_fine>0)) then
	1438	call abort_gcm('phytracr_mod', 'pb in ok_chimdust 0',1)
	1439	endif
	1440	else
	1441	if (id_scdu>0) then
	1442	call abort_gcm('phytracr_mod', 'pb in ok_chimdust 1 SCDU',1)
	1443	endif
	1444	if ( (id_codu .le. 0) .or. ( id_fine.le.0) ) then
	1445	call abort_gcm('phytracr_mod', 'pb in ok_chimdust 1',1)
	1446	endif
	1447	endif
	1448
	1449
	1450	!print *,id_prec,id_fine,id_coss,id_codu,id_scdu
	1451	ENDIF
	1452
	1453
	1454
	1455
	1456
	1457
	1458	!---fraction of tracer that is convected (Tiedke)
	1459	xconv(:)=0.
	1460	if(id_prec>0) xconv(id_prec)=0.8
	1461	if(id_fine>0) xconv(id_fine)=0.5
	1462	if(id_coss>0) xconv(id_coss)=0.5
	1463	if(id_codu>0) xconv(id_codu)=0.6
	1464	if(id_scdu>0) xconv(id_scdu)=0.6 !!JE fix
	1465
	1466	masse(:)=1.
	1467	if(id_prec>0) masse(id_prec)=32.
	1468	if(id_fine>0) masse(id_fine)=6.02e23
	1469	if(id_coss>0) masse(id_coss)=6.02e23
	1470	if(id_codu>0) masse(id_codu)=6.02e23
	1471	if(id_scdu>0) masse(id_scdu)=6.02e23
	1472
	1473	henry(:)=0.
	1474	if(id_prec>0) henry(id_prec)=1.4
	1475	if(id_fine>0) henry(id_fine)=0.0
	1476	if(id_coss>0) henry(id_coss)=0.0
	1477	if(id_codu>0) henry(id_codu)=0.0
	1478	if(id_scdu>0) henry(id_scdu)=0.0
	1479	!henry= (/1.4, 0.0, 0.0, 0.0/)
	1480	kk(:)=0.
	1481	if(id_prec>0) kk(id_prec)=2900.
	1482	if(id_fine>0) kk(id_fine)=0.0
	1483	if(id_coss>0) kk(id_coss)=0.0
	1484	if(id_codu>0) kk(id_codu)=0.0
	1485	if(id_scdu>0) kk(id_scdu)=0.0
	1486	!kk = (/2900., 0., 0., 0./)
	1487	alpha_r(:)=0.
	1488	if(id_prec>0) alpha_r(id_prec)=0.0
	1489	if(id_fine>0) alpha_r(id_fine)=0.001
	1490	if(id_coss>0) alpha_r(id_coss)=0.001
	1491	if(id_codu>0) alpha_r(id_codu)=0.001
	1492	if(id_scdu>0) alpha_r(id_scdu)=0.001 !JE fix
	1493	alpha_s(:)=0.
	1494	if(id_prec>0) alpha_s(id_prec)=0.0
	1495	if(id_fine>0) alpha_s(id_fine)=0.01
	1496	if(id_coss>0) alpha_s(id_coss)=0.01
	1497	if(id_codu>0) alpha_s(id_codu)=0.01
	1498	if(id_scdu>0) alpha_s(id_scdu)=0.01 !JE fix
	1499
	1500	! alpha_r = (/0., 0.001, 0.001, 0.001/)
	1501	! alpha_s = (/0., 0.01, 0.01, 0.01/)
	1502
	1503	! nhl DATA vdep_oce /0.7, 0.05, 1.2, 1.2/
	1504	! nhl vdep_oce for tr1 is a weighted average of dms and so2 dep velocities
	1505	!vdep_oce = (/0.28, 0.28, 1.2, 1.2/)
	1506	vdep_oce(:)=0.
	1507	if(id_prec>0) vdep_oce(id_prec) = 0.28
	1508	if(id_fine>0) vdep_oce(id_fine) = 0.28
	1509	if(id_coss>0) vdep_oce(id_coss) = 1.2
	1510	if(id_codu>0) vdep_oce(id_codu) = 1.2
	1511	if(id_scdu>0) vdep_oce(id_scdu) = 1.2
	1512	vdep_sic(:)=0.
	1513	if(id_prec>0) vdep_sic(id_prec) = 0.2
	1514	if(id_fine>0) vdep_sic(id_fine) = 0.17
	1515	if(id_coss>0) vdep_sic(id_coss) = 1.2
	1516	if(id_codu>0) vdep_sic(id_codu) = 1.2
	1517	if(id_scdu>0) vdep_sic(id_scdu) = 1.2
	1518
	1519	!vdep_sic = (/0.2, 0.17, 1.2, 1.2/)
	1520	!vdep_ter = (/0.3, 0.14, 1.2, 1.2/)
	1521	vdep_ter(:)=0.
	1522	if(id_prec>0) vdep_ter(id_prec) = 0.3
	1523	if(id_fine>0) vdep_ter(id_fine) = 0.14
	1524	if(id_coss>0) vdep_ter(id_coss) = 1.2
	1525	if(id_codu>0) vdep_ter(id_codu) = 1.2
	1526	if(id_scdu>0) vdep_ter(id_scdu) = 1.2
	1527
	1528	vdep_lic(:)=0.
	1529	if(id_prec>0) vdep_lic(id_prec) = 0.2
	1530	if(id_fine>0) vdep_lic(id_fine) = 0.17
	1531	if(id_coss>0) vdep_lic(id_coss) = 1.2
	1532	if(id_codu>0) vdep_lic(id_codu) = 1.2
	1533	if(id_scdu>0) vdep_lic(id_scdu) = 1.2
	1534
	1535
	1536	! convective KE lessivage aer params:
[3786]	1537	! AS: #DFB (Binta) a aussi teste ccntrAA_spla=ccntrENV_spla=0.9/1.0/0.9/0.9
	1538	! mais effet negligeable sur l'AOD
[2630]	1539	ccntrAA_spla(:)=0.
	1540	if(id_prec>0) ccntrAA_spla(id_prec)=-9999.
	1541	if(id_fine>0) ccntrAA_spla(id_fine)=0.7
	1542	if(id_coss>0) ccntrAA_spla(id_coss)=1.0
	1543	if(id_codu>0) ccntrAA_spla(id_codu)=0.7
	1544	if(id_scdu>0) ccntrAA_spla(id_scdu)=0.7
	1545
	1546	ccntrENV_spla(:)=0.
	1547	if(id_prec>0) ccntrENV_spla(id_prec)=-9999.
	1548	if(id_fine>0) ccntrENV_spla(id_fine)=0.7
	1549	if(id_coss>0) ccntrENV_spla(id_coss)=1.0
	1550	if(id_codu>0) ccntrENV_spla(id_codu)=0.7
	1551	if(id_scdu>0) ccntrENV_spla(id_scdu)=0.7
[3786]	1552	! #DFB
[2630]	1553	coefcoli_spla(:)=0.
	1554	if(id_prec>0) coefcoli_spla(id_prec)=-9999.
	1555	if(id_fine>0) coefcoli_spla(id_fine)=0.001
	1556	if(id_coss>0) coefcoli_spla(id_coss)=0.001
	1557	if(id_codu>0) coefcoli_spla(id_codu)=0.001
	1558	if(id_scdu>0) coefcoli_spla(id_scdu)=0.001
	1559
	1560	!vdep_lic = (/0.2, 0.17, 1.2, 1.2/)
	1561	!
	1562
	1563	iscm3=.false.
	1564	if (debutphy) then
	1565	!$OMP MASTER
	1566	CALL suphel
	1567	print *, 'let s check nbtr=', nbtr
	1568	! JE before put in zero
	1569	IF (lminmax) THEN
	1570	DO it=1,nbtr
	1571	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan init phytracr')
	1572	ENDDO
	1573	DO it=1,nbtr
	1574	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'minmax init phytracr')
	1575	ENDDO
	1576	CALL minmaxsource(source_tr,qmin,qmax,'maxsource init phytracr')
	1577	ENDIF
	1578	! JE initializon to cero the tracers
	1579	! DO it=1, nbtr
	1580	! tr_seri(:,:,it)=0.0
	1581	! ENDDO
	1582	! JE end
	1583	! Initializing to zero tr_seri for comparison purposes
	1584	! tr_seri(:,:,:)=0.0
	1585	!
	1586	! DO it=1,nbtr
	1587	! trm_aux(it)=0.0
	1588	! src_aux(it)=0.0
	1589	! diag_trm(it)=0.0
	1590	! diag_drydep(it)=0.0
	1591	! diag_wetdep(it)=0.0
	1592	! diag_cvtdep(it)=0.0
	1593	! diag_emissn(it)=0.0
	1594	! ENDDO
	1595	! diag_g2part=0.0
	1596	print *,'PREPARE FILES TO SAVE VARIABLES'
	1597	!
	1598	nbjour=30
	1599	ecrit_tra = NINT(86400./pdtphys) !--1-day average
	1600	ecrit_tra_h = NINT(86400./pdtphys*0.25) !--6-hour average
	1601	ecrit_tra_m = NINT(86400./pdtphys*FLOAT(nbjour)) !--1-mth average
	1602	print *,'ecrit_tra=', pdtphys, ecrit_tra
	1603
[3806]	1604	!!AS deleting lines
	1605	!! IF (ok_histrac) THEN
	1606	!! IF (is_mpi_root .AND. is_omp_root) THEN
	1607	!!-----many deleted lines----
	1608	!!! nbjour=1
	1609	!! ENDIF ! mpi root
	1610	!! ENDIF !--ok_histrac
[2630]	1611
	1612	!$OMP END MASTER
	1613	!$OMP BARRIER
	1614	endif ! debutphy
	1615	!
	1616	!======================================================================
	1617	! Initialisations
	1618	!======================================================================
	1619	!
	1620	!
	1621	! je KE init
	1622	IF (debutphy) THEN
	1623	!$OMP MASTER
	1624
[3786]	1625	ALLOCATE(d_tr_dry(klon,nbtr))
[2630]	1626	ALLOCATE(flux_tr_dry(klon,nbtr),d_tr_dec(klon,klev,nbtr))
[3786]	1627	ALLOCATE(qPrls(klon,nbtr),qPr(klon,klev,nbtr))
[2630]	1628	ALLOCATE(qDi(klon,klev,nbtr))
	1629	ALLOCATE(qPa(klon,klev,nbtr),qMel(klon,klev,nbtr))
	1630	ALLOCATE(qTrdi(klon,klev,nbtr),dtrcvMA(klon,klev,nbtr))
	1631
[3806]	1632	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	1633	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[2630]	1634	ALLOCATE(d_tr_cv_o(klon,klev,nbtr))
	1635	ALLOCATE(d_tr_trsp_o(klon,klev,nbtr))
	1636	ALLOCATE(d_tr_sscav_o(klon,klev,nbtr), &
	1637	d_tr_sat_o(klon,klev,nbtr))
	1638	ALLOCATE(d_tr_uscav_o(klon,klev,nbtr))
	1639
	1640	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	1641	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	1642	ALLOCATE(d_tr_insc_o(klon,klev,nbtr))
	1643	ALLOCATE(d_tr_bcscav_o(klon,klev,nbtr))
	1644	ALLOCATE(d_tr_evapls_o(klon,klev,nbtr))
	1645	ALLOCATE(d_tr_ls_o(klon,klev,nbtr))
	1646	ALLOCATE(d_tr_dyn_o(klon,klev,nbtr))
	1647	ALLOCATE(d_tr_cl_o(klon,klev,nbtr))
	1648	ALLOCATE(d_tr_th_o(klon,klev,nbtr))
	1649	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	1650	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	1651	ALLOCATE(iregion_so4(klon))
	1652	ALLOCATE(iregion_bb(klon))
	1653	ALLOCATE(iregion_ind(klon))
	1654	ALLOCATE(iregion_dust(klon))
	1655	ALLOCATE(iregion_wstardust(klon))
	1656
	1657	!JE20150518<<
	1658	ALLOCATE(masque_aqua(klon))
	1659	ALLOCATE(masque_terra(klon))
	1660
	1661
	1662	masque_aqua(:)=0
	1663	masque_terra(:)=0
	1664
	1665	aod550_terra(:)=0.
	1666	aod550_tr2_terra(:)=0.
	1667	aod550_ss_terra(:)=0.
	1668	aod550_dust_terra(:)=0.
	1669	aod550_dustsco_terra(:)=0.
	1670	aod670_terra(:)=0.
	1671	aod670_tr2_terra(:)=0.
	1672	aod670_ss_terra(:)=0.
	1673	aod670_dust_terra(:)=0.
	1674	aod670_dustsco_terra(:)=0.
	1675	aod865_terra(:)=0.
	1676	aod865_tr2_terra(:)=0.
	1677	aod865_ss_terra(:)=0.
	1678	aod865_dust_terra(:)=0.
	1679	aod865_dustsco_terra(:)=0.
	1680	aod550_aqua(:)=0.
	1681	aod550_tr2_aqua(:)=0.
	1682	aod550_ss_aqua(:)=0.
	1683	aod550_dust_aqua(:)=0.
	1684	aod550_dustsco_aqua(:)=0.
	1685	aod670_aqua(:)=0.
	1686	aod670_tr2_aqua(:)=0.
	1687	aod670_ss_aqua(:)=0.
	1688	aod670_dust_aqua(:)=0.
	1689	aod670_dustsco_aqua(:)=0.
	1690	aod865_aqua(:)=0.
	1691	aod865_tr2_aqua(:)=0.
	1692	aod865_ss_aqua(:)=0.
	1693	aod865_dust_aqua(:)=0.
	1694	aod865_dustsco_aqua(:)=0.
	1695	!JE20150518>>
	1696
	1697
	1698
	1699
	1700
	1701	!
	1702	!Config Key = iflag_lscav
	1703	!Config Desc = Large scale scavenging parametrization: 0=none,
	1704	!1=old(Genthon92),
	1705	! 2=1+PHeinrich, 3=Reddy_Boucher2004, 4=3+RPilon.
	1706	!Config Def = 4
	1707	!Config
	1708	iflag_lscav_omp=4
	1709	call getin('iflag_lscav', iflag_lscav_omp)
	1710	iflag_lscav=iflag_lscav_omp
	1711	! initialiation for time computation
	1712
	1713	tia_spla=0.
	1714	tia_emis=0.
	1715	tia_depo=0.
	1716	tia_cltr=0.
	1717	tia_ther=0.
	1718	tia_sedi=0.
	1719	tia_gasp=0.
	1720	tia_wetap=0.
	1721	tia_cvltr=0.
	1722	tia_lscs=0.
	1723	tia_brop=0.
	1724	tia_outs=0.
	1725	tia_nophytracr=0.
	1726	clock_start_outphytracr=clock_end_outphytracr+1
	1727	!$OMP END MASTER
	1728	!$OMP BARRIER
	1729	ENDIF ! debutphy
	1730
	1731	lmt_dms(:)=0.0
	1732	aux_var2(:)=0.0
	1733	aux_var3(:,:)=0.0
	1734	source_tr(:,:)=0.0
	1735	flux_tr(:,:)=0.0
	1736	flux_sparam_bb(:)=0.0
	1737	flux_sparam_ff(:)=0.0
	1738	flux_sparam_ind(:)=0.0
	1739	flux_sparam_ddfine(:)=0.0
	1740	flux_sparam_ddcoa(:)=0.0
	1741	flux_sparam_ddsco(:)=0.0
	1742	flux_sparam_ssfine(:)=0.0
	1743	flux_sparam_sscoa(:)=0.0
	1744
	1745	! initialiation for time computation
	1746
	1747	ti_spla=0
	1748	ti_emis=0
	1749	ti_depo=0
	1750	ti_cltr=0
	1751	ti_ther=0
	1752	ti_sedi=0
	1753	ti_gasp=0
	1754	ti_wetap=0
	1755	ti_cvltr=0
	1756	ti_lscs=0
	1757	ti_brop=0
	1758	ti_outs=0
	1759
	1760
	1761	DO k=1,klev
	1762	DO i=1,klon
	1763	Mint(i,k)=0.
	1764	END DO
	1765	END DO
	1766
	1767
	1768	!
	1769	DO it=1, nbtr
	1770	DO k=1,klev
	1771	DO i=1,klon
	1772	d_tr_cv(i,k,it)=0.
	1773	d_tr_trsp(i,k,it)=0.
	1774	d_tr_sscav(i,k,it)=0.
	1775	d_tr_sat(i,k,it)=0.
	1776	d_tr_uscav(i,k,it)=0.
	1777	d_tr(i,k,it)=0.
	1778	d_tr_insc(i,k,it)=0.
	1779	d_tr_bcscav(i,k,it)=0.
	1780	d_tr_evapls(i,k,it)=0.
	1781	d_tr_ls(i,k,it)=0.
	1782	d_tr_cl(i,k,it)=0.
	1783	d_tr_th(i,k,it)=0.
	1784
	1785	d_tr_cv_o(i,k,it)=0.
	1786	d_tr_trsp_o(i,k,it)=0.
	1787	d_tr_sscav_o(i,k,it)=0.
	1788	d_tr_sat_o(i,k,it)=0.
	1789	d_tr_uscav_o(i,k,it)=0.
	1790
	1791
	1792	qDi(i,k,it)=0.
	1793	qPr(i,k,it)=0.
	1794	qPa(i,k,it)=0.
	1795	qMel(i,k,it)=0.
	1796	qTrdi(i,k,it)=0.
	1797	dtrcvMA(i,k,it)=0.
	1798	zmfd1a(i,k,it)=0.
	1799	zmfdam(i,k,it)=0.
	1800	zmfphi2(i,k,it)=0.
	1801	END DO
	1802	END DO
	1803	END DO
	1804
	1805
	1806	DO it=1, nbtr
	1807	DO i=1,klon
	1808	qPrls(i,it)=0.0
	1809	dtrconv(i,it)=0.0
	1810	!JE20140507<<
	1811	d_tr_dry(i,it)=0.0
	1812	flux_tr_dry(i,it)=0.0
	1813	!JE20140507>>
	1814	ENDDO
	1815	ENDDO
	1816
	1817	DO it=1, nbtr
	1818	DO i=1, klon
	1819	his_dh(i,it)=0.0
	1820	his_dhlsc(i,it)=0.0
	1821	his_dhcon(i,it)=0.0
	1822	his_dhbclsc(i,it)=0.0
	1823	his_dhbccon(i,it)=0.0
	1824	trm(i,it)=0.0
	1825	his_th(i,it)=0.0
	1826	his_dhkecv(i,it)=0.0
	1827	his_ds(i,it)=0.0
	1828	his_dhkelsc(i,it)=0.0
	1829
	1830	ENDDO
	1831	ENDDO
	1832	!JE:
	1833	DO i=1, klon
	1834	his_g2pgas(i) = 0.0
	1835	his_g2paer(i) = 0.0
	1836	ENDDO
	1837	! endJE
	1838	!
	1839
	1840	DO k=1, klev
	1841	DO i = 1, klon
	1842	zrho(i,k)=pplay(i,k)/t_seri(i,k)/RD
	1843	zdz(i,k)=(paprs(i,k)-paprs(i,k+1))/zrho(i,k)/RG
	1844	zmasse(i,k)=(paprs(i,k)-paprs(i,k+1))/RG
	1845	ENDDO
	1846	ENDDO
	1847	!
	1848	DO i = 1, klon
	1849	zalt(i,1)=pphis(i)/RG
	1850	ENDDO
	1851	DO k=1, klev-1
	1852	DO i = 1, klon
	1853	zalt(i,k+1)=zalt(i,k)+zdz(i,k)
	1854	ENDDO
	1855	ENDDO
	1856
	1857
	1858
	1859	IF (logitime) THEN
	1860	CALL SYSTEM_CLOCK(COUNT=clock_end)
	1861	dife=clock_end-clock_start
	1862	ti_init=dife*MAX(0,SIGN(1,dife)) &
	1863	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	1864	tia_init=tia_init+REAL(ti_init)/REAL(clock_rate)
	1865	ENDIF
	1866	IF (logitime) THEN
	1867	CALL SYSTEM_CLOCK(COUNT=clock_start)
	1868	ENDIF
	1869
	1870
[3814]	1871	IF (debutphy) then
[2630]	1872
[3814]	1873	! AS: initialisation des indices par point de grille physique iregion_*
	1874	! (variables tenant de l'assimilation, a eliminer dans un 2eme temps)
	1875	iregion_dust(:)=1
	1876	iregion_ind(:)=1
	1877	iregion_bb(:)=1
	1878	iregion_wstardust(:)=1
[2630]	1879
[3814]	1880	!AS: lecture des indices dans fichiers "regions_*" eliminee par IF("ASSIM"="YES") (faux donc)
	1881	IF("ASSIM"=="YES") THEN
[2630]	1882	c_FullName1='regions_dustacc'
	1883	!c_FullName1='regions_dust'
	1884	call readregions_spl(iregion_dust,c_FullName1)
	1885	c_FullName1='regions_ind'
	1886	call readregions_spl(iregion_ind,c_FullName1)
	1887	c_FullName1='regions_bb'
	1888	call readregions_spl(iregion_bb,c_FullName1)
	1889	c_FullName1='regions_pwstarwake'
	1890	call readregions_spl(iregion_wstardust,c_FullName1)
	1891
	1892	!$OMP MASTER
	1893	IF (is_mpi_root .AND. is_omp_root) THEN
	1894
	1895	OPEN(25,FILE='dustregions_pyvar_je.data')
	1896	OPEN(55,FILE='indregions_pyvar_je.data')
	1897	OPEN(75,FILE='bbregions_pyvar_je.data')
	1898	OPEN(95,FILE='wstardustregions_pyvar_je.data')
	1899	OPEN(76,FILE='xlat.data')
	1900	OPEN(77,FILE='xlon.data')
	1901	ENDIF ! mpi root
	1902	!$OMP END MASTER
	1903	!$OMP BARRIER
	1904
	1905	CALL gather(iregion_dust,iauxklon_glo)
	1906	!$OMP MASTER
	1907	IF (is_mpi_root .AND. is_omp_root) THEN
	1908	DO k=1,klon_glo
	1909	WRITE(25,'(i10)') iauxklon_glo(k)
	1910	ENDDO
	1911	ENDIF ! mpi root
	1912	!$OMP END MASTER
	1913	!$OMP BARRIER
	1914	CALL gather(iregion_ind,iauxklon_glo)
	1915	!$OMP MASTER
	1916	IF (is_mpi_root .AND. is_omp_root) THEN
	1917	DO k=1,klon_glo
	1918	WRITE(55,'(i10)') iauxklon_glo(k)
	1919	ENDDO
	1920	ENDIF ! mpi root
	1921	!$OMP END MASTER
	1922	!$OMP BARRIER
	1923	CALL gather(iregion_bb,iauxklon_glo)
	1924	!$OMP MASTER
	1925	IF (is_mpi_root .AND. is_omp_root) THEN
	1926	DO k=1,klon_glo
	1927	WRITE(75,'(i10)') iauxklon_glo(k)
	1928	ENDDO
	1929	ENDIF ! mpi root
	1930	!$OMP END MASTER
	1931	!$OMP BARRIER
	1932	CALL gather(iregion_wstardust,iauxklon_glo)
	1933	!$OMP MASTER
	1934	IF (is_mpi_root .AND. is_omp_root) THEN
	1935	DO k=1,klon_glo
	1936	WRITE(95,'(i10)') iauxklon_glo(k)
	1937	ENDDO
	1938	ENDIF ! mpi root
	1939	!$OMP END MASTER
	1940	!$OMP BARRIER
	1941
	1942
	1943	CALL gather(rlat,auxklon_glo)
	1944	!$OMP MASTER
	1945	IF (is_mpi_root .AND. is_omp_root) THEN
	1946	DO k=1,klon_glo
	1947	WRITE(76,*) auxklon_glo(k)
	1948	ENDDO
	1949	ENDIF ! mpi root
	1950	!$OMP END MASTER
	1951	!$OMP BARRIER
	1952	CALL gather(rlon,auxklon_glo)
	1953	!$OMP MASTER
	1954	IF (is_mpi_root .AND. is_omp_root) THEN
	1955	DO k=1,klon_glo
	1956	WRITE(77,*) auxklon_glo(k)
	1957	ENDDO
	1958
	1959	CLOSE(25)
	1960	CLOSE(55)
	1961	CLOSE(75)
	1962	CLOSE(76)
	1963	CLOSE(77)
[3806]	1964	CLOSE(95)
[2630]	1965
	1966	ENDIF ! mpi root
	1967	!$OMP END MASTER
	1968	!$OMP BARRIER
	1969
[3814]	1970	ENDIF ! ASSIM
	1971
[2630]	1972	ENDIF ! debutphy
	1973
	1974	IF (logitime) THEN
	1975	CALL SYSTEM_CLOCK(COUNT=clock_end)
	1976	dife=clock_end-clock_start
	1977	ti_inittype=dife*MAX(0,SIGN(1,dife)) &
	1978	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	1979	tia_inittype=tia_inittype+REAL(ti_inittype)/REAL(clock_rate)
	1980	ENDIF
	1981
	1982	IF (logitime) THEN
	1983	CALL SYSTEM_CLOCK(COUNT=clock_start)
	1984	ENDIF
	1985
	1986	!
	1987	!=======================================================================
	1988	! SAVING SURFACE TYPE
	1989	!=======================================================================
	1990	IF (debutphy) THEN
	1991	!$OMP MASTER
	1992	IF (is_mpi_root .AND. is_omp_root) THEN
	1993
	1994	OPEN(35,FILE='surface_ocean.data')
	1995	OPEN(45,FILE='surface_seaice.data')
	1996	OPEN(65,FILE='surface_land.data')
	1997	OPEN(85,FILE='surface_landice.data')
	1998	ENDIF ! mpi root
	1999	!$OMP END MASTER
	2000	!$OMP BARRIER
	2001	do i = 1, klon
	2002	aux_var2(i) = pctsrf(i,is_oce)
	2003	enddo
	2004	call gather(aux_var2,auxklon_glo)
	2005	!$OMP MASTER
	2006	IF (is_mpi_root .AND. is_omp_root) THEN
	2007	DO i = 1, klon_glo
	2008	WRITE (35,103) auxklon_glo(i)
	2009	ENDDO
	2010	ENDIF ! mpi root
	2011	!$OMP END MASTER
	2012	!$OMP BARRIER
	2013
	2014	do i = 1, klon
	2015	aux_var2(i) = pctsrf(i,is_sic)
	2016	enddo
	2017	call gather(aux_var2,auxklon_glo)
	2018	!$OMP MASTER
	2019	IF (is_mpi_root .AND. is_omp_root) THEN
	2020	DO i = 1, klon_glo
	2021	WRITE (45,103) auxklon_glo(i)
	2022	ENDDO
	2023	ENDIF ! mpi root
	2024	!$OMP END MASTER
	2025	!$OMP BARRIER
	2026
	2027	do i = 1, klon
	2028	aux_var2(i) = pctsrf(i,is_ter)
	2029	enddo
	2030	call gather(aux_var2,auxklon_glo)
	2031	!$OMP MASTER
	2032	IF (is_mpi_root .AND. is_omp_root) THEN
	2033	DO i = 1, klon_glo
	2034	WRITE (65,103) auxklon_glo(i)
	2035	ENDDO
	2036	ENDIF ! mpi root
	2037	!$OMP END MASTER
	2038	!$OMP BARRIER
	2039
	2040	do i = 1, klon
	2041	aux_var2(i) = pctsrf(i,is_lic)
	2042	enddo
	2043	call gather(aux_var2,auxklon_glo)
	2044	!$OMP MASTER
	2045	IF (is_mpi_root .AND. is_omp_root) THEN
	2046	DO i = 1, klon_glo
	2047	WRITE (85,103) auxklon_glo(i)
	2048	ENDDO
	2049	!
	2050	! DO i = 1, klon
	2051	! WRITE (35,103) pctsrf(i,is_oce)
	2052	! WRITE (45,103) pctsrf(i,is_sic)
	2053	! WRITE (65,103) pctsrf(i,is_ter)
	2054	! WRITE (85,103) pctsrf(i,is_lic)
	2055	! ENDDO
	2056	CLOSE(35)
	2057	CLOSE(45)
	2058	CLOSE(65)
	2059	CLOSE(85)
	2060	103 FORMAT (f6.2)
	2061	ENDIF ! mpi root
	2062	!$OMP END MASTER
	2063	!$OMP BARRIER
	2064	ENDIF ! debutphy
	2065
	2066	! stop
	2067	!
	2068	!=======================================================================
	2069	!
	2070	DO it=1, nbtr
	2071	DO j=1,klev
	2072	DO i=1,klon
	2073	tmp_var(i,j)=tr_seri(i,j,it)
	2074	ENDDO
	2075	ENDDO
	2076	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	2077	DO j=1,klev
	2078	DO i=1,klon
	2079	tr_seri(i,j,it)=tmp_var(i,j)
	2080	ENDDO
	2081	ENDDO
	2082	ENDDO
	2083	iscm3=.true.
	2084
	2085	!=======================================================================
	2086	!
	2087	DO k=1, klev
	2088	DO i=1, klon
	2089	m_conc(i,k)=pplay(i,k)/t_seri(i,k)/RKBOL*1.e-6
	2090	ENDDO
	2091	ENDDO
	2092
	2093	!
	2094	!
	2095	IF (lminmax) THEN
	2096	DO it=1,nbtr
	2097	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_avt_coarem')
	2098	ENDDO
	2099	DO it=1,nbtr
	2100	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'avt coarem')
	2101	ENDDO
	2102	CALL minmaxsource(source_tr,qmin,qmax,'src: avt coarem')
	2103	ENDIF
	2104
	2105	IF (logitime) THEN
	2106	CALL SYSTEM_CLOCK(COUNT=clock_end)
	2107	dife=clock_end-clock_start
	2108	ti_inittwrite=dife*MAX(0,SIGN(1,dife)) &
	2109	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	2110	tia_inittwrite=tia_inittwrite+REAL(ti_inittwrite)/REAL(clock_rate)
	2111	ENDIF
	2112
	2113	!
	2114	!
	2115	!=======================================================================
	2116	! EMISSIONS OF COARSE AEROSOLS
	2117	!=======================================================================
	2118
	2119
	2120	IF (logitime) THEN
	2121	CALL SYSTEM_CLOCK(COUNT=clock_start)
	2122	ENDIF
	2123
	2124
	2125
	2126	!
	2127	print *,'Number of tracers = ',nbtr
	2128
	2129	print *,'AT BEGINNING OF PHYTRACR_SPL'
	2130	! print *,'tr_seri = ',SUM(tr_seri(:,:,3)),MINVAL(tr_seri(:,:,3)),
	2131	! . MAXVAL(tr_seri(:,:,3))
[2632]	2132	#ifdef IOPHYS_DUST
[2630]	2133	do it=1,nbtr
	2134	write(str2,'(i2.2)') it
	2135	call iophys_ecrit('sav'//str2,1,'SOURCE','',source_tr(:,it))
	2136	call iophys_ecrit('fav'//str2,1,'SOURCE','',source_tr(:,it))
	2137	enddo
	2138	do it=1,nbtr
	2139	write(str2,'(i2.2)') it
	2140	call iophys_ecrit('TRB'//str2,klev,'SOURCE','',tr_seri(:,:,it))
	2141	enddo
[2648]	2142	#endif
[2630]	2143
	2144
	2145	CALL coarsemission(pctsrf,pdtphys,t_seri, &
	2146	pmflxr,pmflxs,prfl,psfl, &
	2147	rlat,rlon,debutphy, &
	2148	zu10m,zv10m,wstar,ale_bl,ale_wake, &
	2149	scale_param_ssacc,scale_param_sscoa, &
	2150	scale_param_dustacc,scale_param_dustcoa, &
	2151	scale_param_dustsco, &
	2152	nbreg_dust, &
	2153	iregion_dust,dust_ec, &
	2154	param_wstarBLperregion,param_wstarWAKEperregion, &
	2155	nbreg_wstardust, &
	2156	iregion_wstardust, &
	2157	lmt_sea_salt,qmin,qmax, &
	2158	flux_sparam_ddfine,flux_sparam_ddcoa, &
	2159	flux_sparam_ddsco, &
	2160	flux_sparam_ssfine,flux_sparam_sscoa, &
	2161	id_prec,id_fine,id_coss,id_codu,id_scdu, &
	2162	ok_chimeredust, &
	2163	source_tr,flux_tr)
	2164
[2632]	2165	#ifdef IOPHYS_DUST
[2630]	2166	do it=1,nbtr
	2167	write(str2,'(i2.2)') it
	2168	call iophys_ecrit('sap'//str2,1,'SOURCE','',source_tr(:,it))
	2169	call iophys_ecrit('fap'//str2,1,'SOURCE','',source_tr(:,it))
	2170	enddo
[2632]	2171	#endif
[2630]	2172
	2173	IF (lminmax) THEN
	2174	DO it=1,nbtr
	2175	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_after_coarem')
	2176	ENDDO
	2177	DO it=1,nbtr
	2178	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'after coarem')
	2179	ENDDO
	2180	CALL minmaxsource(source_tr,qmin,qmax,'src: after coarem')
	2181	ENDIF
	2182
	2183	!
	2184	!
	2185	!
	2186	!======================================================================
	2187	! EMISSIONS OF AEROSOL PRECURSORS
	2188	!======================================================================
	2189	!
[2632]	2190	#ifdef IOPHYS_DUST
[2630]	2191	print *,'INPUT TO PRECUREMISSION'
	2192	call iophys_ecrit('ftsol',4,'ftsol','',ftsol)
	2193	call iophys_ecrit('u10m_ec',1,'u10m_ec','',u10m_ec)
	2194	call iophys_ecrit('v10m_ec',1,'v10m_ec','',v10m_ec)
	2195	call iophys_ecrit('pctsrf',4,'pctsrf','',pctsrf)
	2196	call iophys_ecrit('u_seri',klev,'u_seri','',u_seri)
	2197	call iophys_ecrit('v_seri',klev,'v_seri','',v_seri)
	2198	call iophys_ecrit('paprs',klev,'paprs','',paprs)
	2199	call iophys_ecrit('pplay',klev,'pplay','',pplay)
	2200	call iophys_ecrit('cdragh',1,'cdragh','',cdragh)
	2201	call iophys_ecrit('cdragm',1,'cdragm','',cdragm)
	2202	call iophys_ecrit('t_seri',klev,'t_seri','',t_seri)
	2203	call iophys_ecrit('q_seri',klev,'q_seri','',q_seri)
	2204	call iophys_ecrit('tsol',1,'tsol','',tsol)
	2205	print*,'fracso2emis,frach2sofso2,bateau',fracso2emis,frach2sofso2,bateau
	2206	print*,'kminbc,kmaxbc,pdtphys',kminbc,kmaxbc,pdtphys
	2207	print*,'scale_param_bb,scale_param_ind',scale_param_bb,scale_param_ind
	2208	print*,'iregion_ind,iregion_bb,nbreg_ind, nbreg_bb',iregion_ind,iregion_bb,nbreg_ind, nbreg_bb
	2209	print*,'id_prec,id_fine',id_prec,id_fine
	2210	call iophys_ecrit('zdz',klev,'zdz','',zdz)
	2211	call iophys_ecrit('zalt',klev,'zalt','',zalt)
	2212	call iophys_ecrit('lmt_so2ff_l',1,'lmt_so2ff_l','',lmt_so2ff_l)
	2213	call iophys_ecrit('lmt_so2ff_h',1,'lmt_so2ff_h','',lmt_so2ff_h)
	2214	call iophys_ecrit('lmt_so2nff',1,'lmt_so2nff','',lmt_so2nff)
	2215	call iophys_ecrit('lmt_so2ba',1,'lmt_so2ba','',lmt_so2ba)
	2216	call iophys_ecrit('lmt_so2bb_l',1,'lmt_so2bb_l','',lmt_so2bb_l)
	2217	call iophys_ecrit('lmt_so2bb_h',1,'lmt_so2bb_h','',lmt_so2bb_h)
	2218	call iophys_ecrit('lmt_so2volc_cont',1,'lmt_so2volc_cont','',lmt_so2volc_cont)
	2219	call iophys_ecrit('lmt_altvolc_cont',1,'lmt_altvolc_cont','',lmt_altvolc_cont)
	2220	call iophys_ecrit('lmt_so2volc_expl',1,'lmt_so2volc_expl','',lmt_so2volc_expl)
	2221	call iophys_ecrit('lmt_altvolc_expl',1,'lmt_altvolc_expl','',lmt_altvolc_expl)
	2222	call iophys_ecrit('lmt_dmsbio',1,'lmt_dmsbio','',lmt_dmsbio)
	2223	call iophys_ecrit('lmt_h2sbio',1,'lmt_h2sbio','',lmt_h2sbio)
	2224	call iophys_ecrit('lmt_dmsconc',1,'lmt_dmsconc','',lmt_dmsconc)
	2225	call iophys_ecrit('lmt_dms',1,'lmt_dms','',lmt_dms)
	2226	call iophys_ecrit('flux_sparam_ind',1,'flux_sparam_ind','',flux_sparam_ind)
	2227	call iophys_ecrit('flux_sparam_bb',1,'flux_sparam_bb','',flux_sparam_bb)
[2632]	2228	#endif
[2630]	2229
	2230
	2231
	2232	print*,'ON PASSE DANS precuremission'
	2233	CALL precuremission(ftsol,u10m_ec,v10m_ec,pctsrf, &
	2234	u_seri,v_seri,paprs,pplay,cdragh,cdragm, &
	2235	t_seri,q_seri,tsol,fracso2emis,frach2sofso2, &
	2236	bateau,zdz,zalt,kminbc,kmaxbc,pdtphys, &
	2237	scale_param_bb,scale_param_ind, &
	2238	iregion_ind, iregion_bb, &
	2239	nbreg_ind, nbreg_bb, &
	2240	lmt_so2ff_l,lmt_so2ff_h, lmt_so2nff,lmt_so2ba, &
	2241	lmt_so2bb_l,lmt_so2bb_h, &
	2242	lmt_so2volc_cont,lmt_altvolc_cont, &
	2243	lmt_so2volc_expl,lmt_altvolc_expl, &
	2244	lmt_dmsbio,lmt_h2sbio, lmt_dmsconc, lmt_dms, &
	2245	id_prec,id_fine, &
	2246	flux_sparam_ind, flux_sparam_bb, &
	2247	source_tr,flux_tr,tr_seri)
	2248	!
	2249	IF (lminmax) THEN
	2250	DO it=1,nbtr
	2251	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_after precur')
	2252	ENDDO
	2253	DO it=1,nbtr
	2254	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'after precur')
	2255	ENDDO
	2256	CALL minmaxsource(source_tr,qmin,qmax,'src: after precur')
	2257	ENDIF
	2258
	2259	!=======================================================================
	2260	! EMISSIONS OF FINE AEROSOLS
	2261	!=======================================================================
[2632]	2262	#ifdef IOPHYS_DUST
[2630]	2263	!
	2264	do it=1,nbtr
	2265	write(str2,'(i2.2)') it
	2266	call iophys_ecrit('tpr'//str2,1,'SOURCE','',source_tr(:,it))
	2267	call iophys_ecrit('fpr'//str2,1,'SOURCE','',flux_tr(:,it))
	2268	enddo
[2632]	2269	#endif
[2630]	2270
	2271	CALL finemission(zdz,pdtphys,zalt,kminbc,kmaxbc, &
	2272	scale_param_bb,scale_param_ff, &
	2273	iregion_ind,iregion_bb, &
	2274	nbreg_ind,nbreg_bb, &
	2275	lmt_bcff, lmt_bcnff, lmt_bcbb_l,lmt_bcbb_h, &
	2276	lmt_bcba, lmt_omff, lmt_omnff, &
	2277	lmt_ombb_l, lmt_ombb_h, lmt_omnat, lmt_omba, &
	2278	id_fine, &
	2279	flux_sparam_bb, flux_sparam_ff, &
	2280	source_tr,flux_tr,tr_seri)
	2281	!
	2282	!
	2283	IF (lminmax) THEN
	2284	DO it=1,nbtr
	2285	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_after_fineem')
	2286	ENDDO
	2287	DO it=1,nbtr
	2288	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'after fineem')
	2289	ENDDO
	2290	IF (lcheckmass) THEN
	2291	DO it=1,nbtr
	2292	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	2293	pplay,t_seri,iscm3,'after fineem')
	2294	ENDDO
	2295	ENDIF
	2296	CALL minmaxsource(source_tr,qmin,qmax,'src: after fineem')
	2297	ENDIF
	2298
	2299	!
	2300
	2301	IF (logitime) THEN
	2302	CALL SYSTEM_CLOCK(COUNT=clock_end)
	2303	dife=clock_end-clock_start
	2304	ti_emis=dife*MAX(0,SIGN(1,dife)) &
	2305	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	2306	tia_emis=tia_emis+REAL(ti_emis)/REAL(clock_rate)
	2307	ENDIF
	2308
	2309
[2632]	2310	#ifdef IOPHYS_DUST
[2630]	2311	do it=1,nbtr
	2312	write(str2,'(i2.2)') it
	2313	call iophys_ecrit('t'//str2,1,'SOURCE','',source_tr(:,it))
	2314	call iophys_ecrit('f'//str2,1,'SOURCE','',flux_tr(:,it))
	2315	enddo
[2632]	2316	#endif
[2630]	2317	!
	2318	!
	2319
	2320
	2321
	2322	!
	2323	!=======================================================================
	2324	! DRY DEPOSITION AND BOUNDARY LAYER MIXING
	2325	!=======================================================================
	2326	!
	2327	! DO it=1,nbtr
	2328	! CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz,
	2329	! . pplay,t_seri,iscm3,'')
	2330	! ENDDO
	2331
	2332	!======================================================================
	2333	! -- Dry deposition --
	2334	!======================================================================
	2335	IF (logitime) THEN
	2336	CALL SYSTEM_CLOCK(COUNT=clock_start)
	2337	ENDIF
	2338
	2339	DO it=1, nbtr
	2340	DO j=1,klev
	2341	DO i=1,klon
	2342	tmp_var(i,j)=tr_seri(i,j,it)
	2343	ENDDO
	2344	ENDDO
	2345	CALL cm3_to_kg(pplay,t_seri,tmp_var)
	2346	DO j=1,klev
	2347	DO i=1,klon
	2348	tr_seri(i,j,it)=tmp_var(i,j)
	2349	ENDDO
	2350	ENDDO
	2351	ENDDO
	2352	iscm3=.false.
	2353	!----------------------------
	2354	IF (lminmax) THEN
	2355	DO it=1,nbtr
	2356	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_before_depo')
	2357	ENDDO
	2358	DO it=1,nbtr
	2359	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'before depo')
	2360	ENDDO
	2361	IF (lcheckmass) THEN
	2362	DO it=1,nbtr
	2363	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	2364	pplay,t_seri,iscm3,'before depo')
	2365	ENDDO
	2366	ENDIF
	2367	CALL minmaxsource(source_tr,qmin,qmax,'src: before depo')
	2368	ENDIF
	2369
[2632]	2370	#ifdef IOPHYS_DUST
[2630]	2371	do it=1,nbtr
	2372	write(str2,'(i2.2)') it
	2373	call iophys_ecrit('TRC'//str2,klev,'SOURCE','',tr_seri(:,:,it))
	2374	enddo
[2632]	2375	#endif
[2630]	2376
	2377	CALL deposition(vdep_oce,vdep_sic,vdep_ter,vdep_lic,pctsrf, &
	2378	zrho,zdz,pdtphys,RHcl,masse,t_seri,pplay,paprs, &
	2379	lminmax,qmin,qmax, &
	2380	his_ds,source_tr,tr_seri)
	2381	!
	2382	IF (lminmax) THEN
	2383	DO it=1,nbtr
	2384	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_after_depo')
	2385	ENDDO
	2386	DO it=1,nbtr
	2387	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'after depo')
	2388	ENDDO
	2389	IF (lcheckmass) THEN
	2390	DO it=1,nbtr
	2391	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	2392	pplay,t_seri,iscm3,'after depo')
	2393	ENDDO
	2394	ENDIF
	2395	CALL minmaxsource(source_tr,qmin,qmax,'src: after depo')
	2396	ENDIF
	2397
	2398	IF (logitime) THEN
	2399	CALL SYSTEM_CLOCK(COUNT=clock_end)
	2400	dife=clock_end-clock_start
	2401	ti_depo=dife*MAX(0,SIGN(1,dife)) &
	2402	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	2403	tia_depo=tia_depo+REAL(ti_depo)/REAL(clock_rate)
	2404	ENDIF
	2405
	2406
	2407	!
	2408	!======================================================================
	2409	! -- Boundary layer mixing --
	2410	!======================================================================
	2411
[2632]	2412	#ifdef IOPHYS_DUST
[2630]	2413	do it=1,nbtr
	2414	write(str2,'(i2.2)') it
	2415	call iophys_ecrit('TRD'//str2,klev,'SOURCE','',tr_seri(:,:,it))
	2416	enddo
[2632]	2417	#endif
[2630]	2418
	2419
	2420
	2421	IF (logitime) THEN
	2422	CALL SYSTEM_CLOCK(COUNT=clock_start)
	2423	ENDIF
	2424
	2425	!
	2426
	2427	DO k = 1, klev
	2428	DO i = 1, klon
	2429	delp(i,k) = paprs(i,k)-paprs(i,k+1)
	2430	END DO
	2431	END DO
	2432	!
	2433	DO it=1, nbtr
	2434	DO j=1, klev
	2435	DO i=1, klon
	2436	tmp_var(i,j)=tr_seri(i,j,it)
	2437	aux_var2(i)=source_tr(i,it)
	2438	ENDDO
	2439	ENDDO
	2440	IF (iflag_conv.EQ.2) THEN
	2441	! Tiedke
	2442	CALL cltrac_spl(pdtphys,coefh,yu1,yv1,t_seri,tmp_var, &
	2443	aux_var2,paprs,pplay,aux_var3)
	2444
	2445	ELSE IF (iflag_conv.GE.3) THEN
	2446	!KE
	2447	CALL cltrac(pdtphys, coefh,t_seri,tmp_var,aux_var2,paprs,pplay, &
	2448	delp,aux_var3,d_tr_dry,flux_tr_dry(:,it))
	2449	ENDIF
	2450
	2451	DO i=1, klon
	2452	DO j=1, klev
	2453	tr_seri(i,j,it)=tmp_var(i,j)
	2454	d_tr(i,j,it)=aux_var3(i,j)
	2455	d_tr_cl(i,j,it)=d_tr(i,j,it)
	2456	ENDDO
	2457	ENDDO
	2458	DO k = 1, klev
	2459	DO i = 1, klon
	2460	tr_seri(i,k,it) = tr_seri(i,k,it) + d_tr(i,k,it)
	2461	ENDDO
	2462	ENDDO
	2463	print *,' AFTER Cltrac'
	2464	IF (lminmax) THEN
	2465	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'after cltrac')
	2466	ENDIF
	2467	ENDDO !--end itr loop
	2468
	2469	IF (logitime) THEN
	2470	CALL SYSTEM_CLOCK(COUNT=clock_end)
	2471	dife=clock_end-clock_start
	2472	ti_cltr=dife*MAX(0,SIGN(1,dife)) &
	2473	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	2474	tia_cltr=tia_cltr+REAL(ti_cltr)/REAL(clock_rate)
	2475	ENDIF
	2476
	2477
	2478
	2479	!======================================================================
	2480	! -- Calcul de l'effet des thermiques for KE--
	2481	!======================================================================
	2482
[2632]	2483	#ifdef IOPHYS_DUST
[2630]	2484	print*,'iflag_conv=',iflag_conv
	2485	call iophys_ecrit('coefh',klev,'coefh','',coefh)
	2486	call iophys_ecrit('yu1',1,'yu1','',yu1)
	2487	call iophys_ecrit('yv1',1,'yv1','',yv1)
	2488	call iophys_ecrit('delp',klev,'delp','',delp)
	2489	do it=1,nbtr
	2490	write(str2,'(i2.2)') it
	2491	call iophys_ecrit('TRE'//str2,klev,'SOURCE','',tr_seri(:,:,it))
	2492	enddo
[2632]	2493	#endif
[2630]	2494
	2495
	2496
	2497	IF (iflag_conv.GE.3) THEN
	2498
	2499	IF (logitime) THEN
	2500	CALL SYSTEM_CLOCK(COUNT=clock_start)
	2501	ENDIF
	2502
	2503
	2504
	2505
	2506
	2507	IF (lminmax) THEN
	2508	DO it=1,nbtr
	2509	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_before therm')
	2510	ENDDO
	2511	DO it=1,nbtr
	2512	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'before therm')
	2513	ENDDO
	2514	IF (lcheckmass) THEN
	2515	DO it=1,nbtr
	2516	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	2517	pplay,t_seri,iscm3,'before therm')
	2518	ENDDO
	2519	ENDIF
	2520	CALL minmaxsource(source_tr,qmin,qmax,'before therm')
	2521	ENDIF
	2522
	2523	DO it=1,nbtr
	2524	DO k=1,klev
	2525	DO i=1,klon
	2526	tmp_var3(i,k,it)=tr_seri(i,k,it)
	2527	d_tr_th(i,k,it)=0.
	2528	tr_seri(i,k,it)=MAX(tr_seri(i,k,it),0.)
	2529	!JE: precursor >>1e10 tr_seri(i,k,it)=MIN(tr_seri(i,k,it),1.e10)
	2530	END DO
	2531	END DO
	2532	END DO
	2533
	2534	!JE new implicit scheme 20140323
	2535	DO it=1,nbtr
	2536	CALL thermcell_dq(klon,klev,1,pdtphys,fm_therm,entr_therm, &
	2537	zmasse,tr_seri(1:klon,1:klev,it), &
	2538	d_tr(1:klon,1:klev,it),ztra_th,0 )
	2539
	2540	DO k=1,klev
	2541	DO i=1,klon
	2542	d_tr(i,k,it)=pdtphys*d_tr(i,k,it)
	2543	d_tr_th(i,k,it)=d_tr_th(i,k,it)+d_tr(i,k,it)
	2544	tr_seri(i,k,it)=MAX(tr_seri(i,k,it)+d_tr(i,k,it),0.)
	2545	END DO
	2546	END DO
	2547
	2548	ENDDO
	2549
	2550	! old scheme explicit
	2551	! nsplit=10
	2552	! DO it=1,nbtr
	2553	! DO isplit=1,nsplit
	2554	! CALL dqthermcell(klon,klev,pdtphys/nsplit,
	2555	! . fm_therm,entr_therm,zmasse,
	2556	! . tr_seri(1:klon,1:klev,it),
	2557	! . d_tr(1:klon,1:klev,it),ztra_th)
	2558	! DO k=1,klev
	2559	! DO i=1,klon
	2560	! d_tr(i,k,it)=pdtphys*d_tr(i,k,it)/nsplit
	2561	! d_tr_th(i,k,it)=d_tr_th(i,k,it)+d_tr(i,k,it)
	2562	! tr_seri(i,k,it)=MAX(tr_seri(i,k,it)+d_tr(i,k,it),0.)
	2563	! END DO
	2564	! END DO
	2565	! END DO ! nsplit1
	2566	! END DO ! it
	2567	!JE end modif 20140323
	2568
	2569	DO it=1,nbtr
	2570	DO k=1,klev
	2571	DO i=1,klon
	2572	tmp_var(i,k)=tr_seri(i,k,it)-tmp_var3(i,k,it)
	2573	ENDDO
	2574	ENDDO
	2575	IF (lminmax) THEN
	2576	IF (lcheckmass) THEN
	2577	CALL checkmass(tmp_var(:,:),RNAVO,masse(it),zdz, &
	2578	pplay,t_seri,iscm3,'dtr therm ')
	2579	ENDIF
	2580	ENDIF
	2581	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	2582
	2583	DO k=1,klev
	2584	DO i=1,klon
	2585	his_th(i,it)=his_th(i,it)+ &
	2586	(tmp_var(i,k))/RNAVO* &
	2587	masse(it)1.e31.e6*zdz(i,k)/pdtphys
	2588	END DO !klon
	2589	END DO !klev
	2590
	2591	END DO !it
	2592	IF (lminmax) THEN
	2593	DO it=1,nbtr
	2594	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_after therm')
	2595	ENDDO
	2596	DO it=1,nbtr
	2597	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'after therm')
	2598	ENDDO
	2599	IF (lcheckmass) THEN
	2600	DO it=1,nbtr
	2601	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	2602	pplay,t_seri,iscm3,'after therm')
	2603	ENDDO
	2604	ENDIF
	2605	CALL minmaxsource(source_tr,qmin,qmax,'after therm')
	2606	ENDIF
	2607
	2608	IF (logitime) THEN
	2609	CALL SYSTEM_CLOCK(COUNT=clock_end)
	2610	dife=clock_end-clock_start
	2611	ti_ther=dife*MAX(0,SIGN(1,dife)) &
	2612	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	2613	tia_ther=tia_ther+REAL(ti_ther)/REAL(clock_rate)
	2614	ENDIF
	2615
	2616
	2617	ENDIF ! iflag_conv KE
	2618	!------------------------------------
	2619	! Sedimentation
	2620	!-----------------------------------
	2621	IF (logitime) THEN
	2622	CALL SYSTEM_CLOCK(COUNT=clock_start)
	2623	ENDIF
	2624
	2625
	2626	DO it=1,nbtr
	2627	DO j=1,klev
	2628	DO i=1,klon
	2629	tmp_var(i,j)=tr_seri(i,j,it)
	2630	ENDDO
	2631	ENDDO
	2632	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	2633	DO j=1,klev
	2634	DO i=1,klon
	2635	tr_seri(i,j,it)=tmp_var(i,j)
	2636	ENDDO
	2637	ENDDO
	2638	ENDDO !--end itr loop
	2639	iscm3=.true.
	2640	!--------------------------------------
	2641	print *,' BEFORE Sediment'
	2642
	2643	IF (lminmax) THEN
	2644	DO it=1,nbtr
	2645	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_before_sedi')
	2646	ENDDO
	2647	DO it=1,nbtr
	2648	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'before sedi')
	2649	ENDDO
	2650	IF (lcheckmass) THEN
	2651	DO it=1,nbtr
	2652	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	2653	pplay,t_seri,iscm3,'before sedi')
	2654	ENDDO
	2655	ENDIF
	2656	CALL minmaxsource(source_tr,qmin,qmax,'src: before sedi')
	2657	ENDIF
	2658
	2659	print *,'SPLA VERSION OF SEDIMENTATION IS USED'
	2660	CALL sediment_mod(t_seri,pplay,zrho,paprs,pdtphys,RHcl, &
	2661	id_coss,id_codu,id_scdu, &
	2662	ok_chimeredust, &
	2663	sed_ss,sed_dust,sed_dustsco, &
	2664	sed_ss3D,sed_dust3D,sed_dustsco3D,tr_seri)
	2665	CALL cm3_to_kg(pplay,t_seri,sed_ss3D)
	2666	CALL cm3_to_kg(pplay,t_seri,sed_dust3D)
	2667	CALL cm3_to_kg(pplay,t_seri,sed_dustsco3D)
	2668
	2669	IF (lminmax) THEN
	2670	DO it=1,nbtr
	2671	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_after_sedi')
	2672	ENDDO
	2673	DO it=1,nbtr
	2674	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'after sedi')
	2675	ENDDO
	2676	IF (lcheckmass) THEN
	2677	DO it=1,nbtr
	2678	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	2679	pplay,t_seri,iscm3,'after sedi')
	2680	ENDDO
	2681	ENDIF
	2682	CALL minmaxsource(source_tr,qmin,qmax,'src: after sedi')
	2683	ENDIF
	2684
	2685	!
	2686	!=======================================================================
[2632]	2687	#ifdef IOPHYS_DUST
[2630]	2688	do it=1,nbtr
	2689	write(str2,'(i2.2)') it
	2690	call iophys_ecrit('TRF'//str2,klev,'SOURCE','',tr_seri(:,:,it))
	2691	enddo
[2632]	2692	#endif
[2630]	2693
	2694
	2695
	2696	!
	2697	IF (logitime) THEN
	2698	CALL SYSTEM_CLOCK(COUNT=clock_end)
	2699	dife=clock_end-clock_start
	2700	ti_sedi=dife*MAX(0,SIGN(1,dife)) &
	2701	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	2702	tia_sedi=tia_sedi+REAL(ti_sedi)/REAL(clock_rate)
	2703	ENDIF
	2704
	2705	DO it=1, nbtr
	2706	DO j=1,klev
	2707	DO i=1,klon
	2708	tmp_var(i,j)=tr_seri(i,j,it)
	2709	ENDDO
	2710	ENDDO
	2711	CALL cm3_to_kg(pplay,t_seri,tmp_var)
	2712	DO j=1,klev
	2713	DO i=1,klon
	2714	tr_seri(i,j,it)=tmp_var(i,j)
	2715	ENDDO
	2716	ENDDO
	2717	ENDDO
	2718	iscm3=.false.
	2719	!
	2720	!
	2721	!======================================================================
	2722	! GAS TO PARTICLE CONVERSION
	2723	!======================================================================
	2724	!
	2725
	2726	IF (logitime) THEN
	2727	CALL SYSTEM_CLOCK(COUNT=clock_start)
	2728	ENDIF
	2729
	2730	IF (lminmax) THEN
	2731	DO it=1,nbtr
	2732	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_beforegastopar')
	2733	ENDDO
	2734	DO it=1,nbtr
	2735	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'before gastopar')
	2736	ENDDO
	2737	IF (lcheckmass) THEN
	2738	DO it=1,nbtr
	2739	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	2740	pplay,t_seri,iscm3,'before gastopar')
	2741	ENDDO
	2742	ENDIF
	2743	CALL minmaxsource(source_tr,qmin,qmax,'src: before gastopar')
	2744	ENDIF
	2745
	2746	CALL gastoparticle(pdtphys,zdz,zrho,rlat, &
	2747	pplay,t_seri,id_prec,id_fine, &
	2748	tr_seri,his_g2pgas ,his_g2paer)
	2749	!
	2750	IF (lminmax) THEN
	2751	DO it=1,nbtr
	2752	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_after_gastopar')
	2753	ENDDO
	2754	DO it=1,nbtr
	2755	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'after gastopar')
	2756	ENDDO
	2757	IF (lcheckmass) THEN
	2758	DO it=1,nbtr
	2759	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	2760	pplay,t_seri,iscm3,'after gastopar')
	2761	ENDDO
	2762	ENDIF
	2763	CALL minmaxsource(source_tr,qmin,qmax,'src: after gastopar')
	2764	ENDIF
	2765
	2766	IF (logitime) THEN
	2767	CALL SYSTEM_CLOCK(COUNT=clock_end)
	2768	dife=clock_end-clock_start
	2769	ti_gasp=dife*MAX(0,SIGN(1,dife)) &
	2770	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	2771	tia_gasp=tia_gasp+REAL(ti_gasp)/REAL(clock_rate)
	2772	ENDIF
	2773
	2774
	2775	!
	2776	!======================================================================
	2777	! EFFECT OF PRECIPITATION: iflag_conv=2
	2778	!======================================================================
	2779	!
	2780
[2632]	2781	#ifdef IOPHYS_DUST
[2630]	2782	do it=1,nbtr
	2783	write(str2,'(i2.2)') it
	2784	call iophys_ecrit('TRG'//str2,klev,'SOURCE','',tr_seri(:,:,it))
	2785	enddo
[2632]	2786	#endif
[2630]	2787
	2788
	2789	IF (iflag_conv.EQ.2) THEN
	2790
	2791	IF (logitime) THEN
	2792	CALL SYSTEM_CLOCK(COUNT=clock_start)
	2793	ENDIF
	2794
	2795
	2796
	2797
	2798	DO it=1, nbtr
	2799	DO j=1,klev
	2800	DO i=1,klon
	2801	tmp_var(i,j)=tr_seri(i,j,it)
	2802	ENDDO
	2803	ENDDO
	2804	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	2805	DO j=1,klev
	2806	DO i=1,klon
	2807	tr_seri(i,j,it)=tmp_var(i,j)
	2808	ENDDO
	2809	ENDDO
	2810	ENDDO
	2811	iscm3=.true.
	2812	!------------------------------
	2813
	2814	print *,'iflag_conv bef lessiv',iflag_conv
	2815	IF (lessivage) THEN
	2816	!
	2817	print *,' BEFORE Incloud'
	2818
	2819	IF (lminmax) THEN
	2820	DO it=1,nbtr
	2821	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_before_incloud')
	2822	ENDDO
	2823	DO it=1,nbtr
	2824	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'before incloud')
	2825	ENDDO
	2826	IF (lcheckmass) THEN
	2827	DO it=1,nbtr
	2828	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	2829	pplay,t_seri,iscm3,'before incloud')
	2830	ENDDO
	2831	ENDIF
	2832	CALL minmaxsource(source_tr,qmin,qmax,'src: before incloud')
	2833	ENDIF
	2834
	2835
	2836	! CALL incloud_scav(lminmax,qmin,qmax,masse,henry,kk,prfl,
	2837	! . psfl,pmflxr,pmflxs,zrho,zdz,t_seri,pdtphys,
	2838
	2839	! . his_dhlsc,his_dhcon,tr_seri)
	2840	print *,'iflag_conv bef incloud',iflag_conv
	2841
	2842	IF (iflag_conv.EQ.2) THEN
	2843	! Tiedke
	2844	CALL incloud_scav(.false.,qmin,qmax,masse,henry,kk,prfl, &
	2845	psfl,pmflxr,pmflxs,zrho,zdz,t_seri,pdtphys, &
	2846	his_dhlsc,his_dhcon,tr_seri)
	2847
	2848	!---------- to use this option please comment lsc_scav at the end
	2849	! ELSE IF (iflag_conv.GE.3) THEN
	2850	!
	2851	! CALL incloud_scav_lsc(.false.,qmin,qmax,masse,henry,kk,prfl,
	2852	! . psfl,pmflxr,pmflxs,zrho,zdz,t_seri,pdtphys,
	2853	! . his_dhlsc,his_dhcon,tr_seri)
	2854	!--------------------------------------------------------------
	2855
	2856	ENDIF
	2857	!
	2858	!
	2859	print *,' BEFORE blcloud (after incloud)'
	2860	IF (lminmax) THEN
	2861	DO it=1,nbtr
	2862	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_before_blcloud')
	2863	ENDDO
	2864	DO it=1,nbtr
	2865	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'before blcloud')
	2866	ENDDO
	2867	IF (lcheckmass) THEN
	2868	DO it=1,nbtr
	2869	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	2870	pplay,t_seri,iscm3,'before blcloud')
	2871	ENDDO
	2872	ENDIF
	2873	CALL minmaxsource(source_tr,qmin,qmax,'src: before blcloud')
	2874	ENDIF
	2875
	2876	! CALL blcloud_scav(lminmax,qmin,qmax,pdtphys,prfl,psfl,
	2877	! . pmflxr,pmflxs,zdz,alpha_r,alpha_s,masse,
	2878	! . his_dhbclsc,his_dhbccon,tr_seri)
	2879
	2880	IF (iflag_conv.EQ.2) THEN
	2881	! Tiedke
	2882
	2883	CALL blcloud_scav(.false.,qmin,qmax,pdtphys,prfl,psfl, &
	2884	pmflxr,pmflxs,zdz,alpha_r,alpha_s,masse, &
	2885	his_dhbclsc,his_dhbccon,tr_seri)
	2886
	2887	!---------- to use this option please comment lsc_scav at the end
	2888	! and comment IF iflag=2 after "EFFECT OF PRECIPITATION:"
	2889	!
	2890	!
	2891	! ELSE IF (iflag_conv.GE.3) THEN
	2892	!
	2893	! CALL blcloud_scav_lsc(.false.,qmin,qmax,pdtphys,prfl,psfl,
	2894	! . pmflxr,pmflxs,zdz,alpha_r,alpha_s,masse,
	2895	! . his_dhbclsc,his_dhbccon,tr_seri)
	2896	!
	2897	!----------------------------------------------------------------------
	2898	ENDIF
	2899
	2900
	2901	print *,' AFTER blcloud '
	2902
	2903	IF (lminmax) THEN
	2904	DO it=1,nbtr
	2905	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_after_blcloud')
	2906	ENDDO
	2907	DO it=1,nbtr
	2908	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'after blcloud')
	2909	ENDDO
	2910	IF (lcheckmass) THEN
	2911	DO it=1,nbtr
	2912	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	2913	pplay,t_seri,iscm3,'after blcloud')
	2914	ENDDO
	2915	ENDIF
	2916	CALL minmaxsource(source_tr,qmin,qmax,'src: after blcloud')
	2917	ENDIF
	2918
	2919
	2920	ENDIF !--lessivage
	2921
	2922	DO it=1, nbtr
	2923	DO j=1,klev
	2924	DO i=1,klon
	2925	tmp_var(i,j)=tr_seri(i,j,it)
	2926	ENDDO
	2927	ENDDO
	2928	CALL cm3_to_kg(pplay,t_seri,tmp_var)
	2929	DO j=1,klev
	2930	DO i=1,klon
	2931	tr_seri(i,j,it)=tmp_var(i,j)
	2932	ENDDO
	2933	ENDDO
	2934	ENDDO
	2935	iscm3=.false.
	2936	!
	2937	IF (logitime) THEN
	2938	CALL SYSTEM_CLOCK(COUNT=clock_end)
	2939	dife=clock_end-clock_start
	2940	ti_wetap=dife*MAX(0,SIGN(1,dife)) &
	2941	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	2942	tia_wetap=tia_wetap+REAL(ti_wetap)/REAL(clock_rate)
	2943	ENDIF
	2944
	2945
	2946
	2947
	2948	ENDIF ! iflag_conv=2
	2949
	2950	!
	2951	!
	2952	!======================================================================
	2953	! EFFECT OF CONVECTION
	2954	!======================================================================
	2955	!
[2632]	2956	#ifdef IOPHYS_DUST
[2630]	2957	do it=1,nbtr
	2958	write(str2,'(i2.2)') it
	2959	call iophys_ecrit('TRH'//str2,klev,'SOURCE','',tr_seri(:,:,it))
	2960	enddo
[2632]	2961	#endif
[2630]	2962
	2963
	2964	IF (logitime) THEN
	2965	CALL SYSTEM_CLOCK(COUNT=clock_start)
	2966	ENDIF
	2967
	2968
	2969	IF (convection) THEN
	2970	!
	2971	print *,' BEFORE trconvect'
	2972
	2973	IF (lminmax) THEN
	2974	DO it=1,nbtr
	2975	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_before_trconve')
	2976	ENDDO
	2977	DO it=1,nbtr
	2978	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'before trconve')
	2979	ENDDO
	2980	IF (lcheckmass) THEN
	2981	DO it=1,nbtr
	2982	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	2983	pplay,t_seri,iscm3,'before trconve')
	2984	ENDDO
	2985	ENDIF
	2986	CALL minmaxsource(source_tr,qmin,qmax,'src: before trconve')
	2987	ENDIF
	2988
	2989
	2990	! JE CALL trconvect(pplay,t_seri,pdtphys,pmfu,pmfd,pen_u,pde_u,
	2991	! . pen_d,pde_d,paprs,zdz,xconv,qmin,qmax,lminmax,masse,
	2992	! . dtrconv,tr_seri)
	2993	! -------------------------------------------------------------
	2994	IF (iflag_conv.EQ.2) THEN
	2995	! Tiedke
	2996	CALL trconvect(pplay,t_seri,pdtphys,pmfu,pmfd,pen_u,pde_u, &
	2997	pen_d,pde_d,paprs,zdz,xconv,qmin,qmax,.false.,masse, &
	2998	dtrconv,tr_seri)
	2999	DO it=1, nbtr
	3000	d_tr_cv(:,:,it)=0.
	3001	ENDDO
	3002
	3003	ELSE IF (iflag_conv.GE.3) THEN
	3004	! KE
	3005	print *,'JE: KE in phytracr_spl'
	3006	DO it=1, nbtr
	3007	DO k = 1, klev
	3008	DO i = 1, klon
	3009	tmp_var3(i,k,it)=tr_seri(i,k,it)
	3010	END DO
	3011	END DO
	3012	ENDDO
	3013
	3014	DO it=1, nbtr
	3015	! routine for aerosols . otherwise, check cvltrorig
	3016	print *,'Check sum before cvltr it',it,SUM(tr_seri(:,:,it))
	3017	! IF (.FALSE.) THEN
	3018	CALL cvltr_spl(pdtphys, da, phi,phi2,d1a,dam, mp,ep, &
	3019	sigd,sij,wght_cvfd,clw,elij,epmlmMm,eplaMm, &
	3020	pmflxr,pmflxs,evapls,t_seri,wdtrainA,wdtrainM, &
	3021	! paprs,it,tr_seri,upwd,dnwd,itop_con,ibas_con, &
	3022	paprs,it,tmp_var3,upwd,dnwd,itop_con,ibas_con, &
	3023	henry,kk,zrho,ccntrAA_spla,ccntrENV_spla,coefcoli_spla, &
	3024	id_prec,id_fine,id_coss, id_codu, id_scdu, &
	3025	d_tr_cv,d_tr_trsp,d_tr_sscav,d_tr_sat,d_tr_uscav,qDi,qPr, &
	3026	qPa,qMel,qTrdi,dtrcvMA,Mint, &
	3027	zmfd1a,zmfphi2,zmfdam)
	3028	! ENDIF
	3029	!
	3030	! IF (.FALSE.) THEN
	3031	! CALL cvltr(pdtphys, da, phi,phi2,d1a,dam, mp,ep,
	3032	! . sigd,sij,wght_cvfd,clw,elij,epmlmMm,eplaMm,
	3033	! . pmflxr,pmflxs,evapls,t_seri,wdtrainA,wdtrainM,
	3034	! . paprs,it,tmp_var3,upwd,dnwd,itop_con,ibas_con,
	3035	! . d_tr_cv,d_tr_trsp,d_tr_sscav,d_tr_sat,d_tr_uscav,qDi,qPr,
	3036	! . qPa,qMel,qTrdi,dtrcvMA,Mint,
	3037	! . zmfd1a,zmfphi2,zmfdam)
	3038	!! pas lessivage convectif pou n'est pas un aerosol (i/else with cvltr)
	3039	! ENDIF
	3040
	3041
	3042
	3043	!!!!!!! CALL cvltrorig(it,pdtphys, da, phi,mp,paprs,pplay,tr_seri,
	3044	!!! CALL cvltrorig(it,pdtphys, da, phi,mp,paprs,pplay,tmp_var3,
	3045	!!! . upwd,dnwd,d_tr_cv)
	3046	! print *,'justbefore cvltrnoscav it= ',it
	3047	! CALL checknanqfi(da(:,:),1.,-1.,' da')
	3048	! CALL checknanqfi(wght_cvfd(:,:),1.,-1.,'weigth ')
	3049	! CALL checknanqfi(mp(:,:),1.,-1.,'mp ')
	3050	! CALL checknanqfi(paprs(:,:),1.,-1.,'paprs ')
	3051	! CALL checknanqfi(pplay(:,:),1.,-1.,'pplay ')
	3052	! CALL checknanqfi(tmp_var3(:,:,it),1.,-1.,'tmp_var3 ')
	3053	! CALL checknanqfi(upwd(:,:),1.,-1.,'upwd ')
	3054	! CALL checknanqfi(dnwd(:,:),1.,-1.,'dnwd ')
	3055	! CALL checknanqfi(d_tr_cv(:,:,it),1.,-1.,'d_tr_cv ')
	3056	! IF (.TRUE.) THEN
	3057	! CALL cvltr_noscav(it,pdtphys, da, phi,mp,wght_cvfd,paprs,
	3058	! . pplay,tmp_var3,upwd,dnwd,d_tr_cv)
	3059	! ENDIF
	3060	DO k = 1, klev
	3061	DO i = 1, klon
	3062	! tr_seri(i,k,it) = tr_seri(i,k,it) + d_tr_cv(i,k,it)
	3063	tr_seri(i,k,it)=(tmp_var3(i,k,it)+d_tr_cv(i,k,it))
	3064	tmp_var(i,k)=d_tr_cv(i,k,it)
	3065
	3066	END DO
	3067	END DO
	3068
	3069	CALL kg_to_cm3(pplay,t_seri,tmp_var) !just for his_* computation
	3070
	3071	DO k = 1, klev
	3072	DO i = 1, klon
	3073	dtrconv(i,it)=0.0
	3074	his_dhkecv(i,it)=his_dhkecv(i,it)-tmp_var(i,k) &
	3075	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3076	END DO
	3077	END DO
	3078
	3079	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	3080	CALL kg_to_cm3(pplay,t_seri,tmp_var) !just for his_* computation
	3081
	3082	DO k = 1, klev
	3083	DO i = 1, klon
	3084	dtrconv(i,it)=0.0
	3085	his_ds(i,it)=his_ds(i,it)-tmp_var(i,k) &
	3086	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3087	END DO
	3088	END DO
	3089	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	3090	IF (lminmax) THEN
	3091
	3092	print *,'Check sum after cvltr it',it,SUM(tr_seri(:,:,it))
	3093	CALL minmaxqfi2(d_tr_cv(:,:,it),qmin,qmax,'d_tr_cv:')
	3094	CALL minmaxqfi2(d_tr_trsp(:,:,it),qmin,qmax,'d_tr_trsp:')
	3095	CALL minmaxqfi2(d_tr_sscav(:,:,it),qmin,qmax,'d_tr_sscav:')
	3096	CALL minmaxqfi2(d_tr_sat(:,:,it),qmin,qmax,'d_tr_sat:')
	3097	CALL minmaxqfi2(d_tr_uscav(:,:,it),qmin,qmax,'d_tr_uscav:')
	3098	IF (lcheckmass) THEN
	3099	CALL checkmass(d_tr_cv(:,:,it),RNAVO,masse(it),zdz, &
	3100	pplay,t_seri,.false.,'d_tr_cv:')
	3101	ENDIF
	3102	ENDIF
	3103	ENDDO ! it=1,nbtr
	3104
	3105	ENDIF ! iflag_conv
	3106	IF (lminmax) THEN
	3107	DO it=1,nbtr
	3108	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_after_trcon')
	3109	ENDDO
	3110	DO it=1,nbtr
	3111	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'after trconv')
	3112	ENDDO
	3113	IF (lcheckmass) THEN
	3114	DO it=1,nbtr
	3115	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	3116	pplay,t_seri,iscm3,'after trconv')
	3117	ENDDO
	3118	ENDIF
	3119	CALL minmaxsource(source_tr,qmin,qmax,'src: after trconv')
	3120	ENDIF
	3121	ENDIF ! convection
	3122
	3123	IF (logitime) THEN
	3124	CALL SYSTEM_CLOCK(COUNT=clock_end)
	3125	dife=clock_end-clock_start
	3126	ti_cvltr=dife*MAX(0,SIGN(1,dife)) &
	3127	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	3128	tia_cvltr=tia_cvltr+REAL(ti_cvltr)/REAL(clock_rate)
	3129	ENDIF
	3130
	3131
	3132
	3133	!
	3134	!
	3135	!=======================================================================
	3136	! LARGE SCALE SCAVENGING KE
	3137	!=======================================================================
	3138	!
[2632]	3139	#ifdef IOPHYS_DUST
[2630]	3140	call iophys_ecrit('da',klev,'da','',da)
	3141	call iophys_ecrit('phi',klev,'phi','',phi)
	3142	call iophys_ecrit('phi2',klev,'phi2','',phi2)
	3143	call iophys_ecrit('d1a',klev,'d1a','',d1a)
	3144	call iophys_ecrit('dam',klev,'dam','',dam)
	3145	call iophys_ecrit('mp',klev,'mp','',mp)
	3146	call iophys_ecrit('ep',klev,'ep','',ep)
	3147	call iophys_ecrit('sigd',klev,'sigd','',sigd)
	3148	call iophys_ecrit('sij',klev,'sij','',sij)
	3149	call iophys_ecrit('wght_cvfd',klev,'wght_cvfd','',wght_cvfd)
	3150	call iophys_ecrit('clw',klev,'clw','',clw)
	3151	call iophys_ecrit('elij',klev,'elij','',elij)
	3152	call iophys_ecrit('epmlmMm',klev,'epmlmMm','',epmlmMm)
	3153	call iophys_ecrit('eplaMm',klev,'eplaMm','',eplaMm)
	3154	call iophys_ecrit('pmflxr',klev,'pmflxr','',pmflxr)
	3155	call iophys_ecrit('pmflxs',klev,'pmflxs','',pmflxs)
	3156	call iophys_ecrit('evapls',klev,'evapls','',evapls)
	3157	call iophys_ecrit('wdtrainA',klev,'wdtrainA','',wdtrainA)
	3158	call iophys_ecrit('wdtrainM',klev,'wdtrainM','',wdtrainM)
	3159
	3160	do it=1,nbtr
	3161	write(str2,'(i2.2)') it
	3162	call iophys_ecrit('TRI'//str2,klev,'SOURCE','',tr_seri(:,:,it))
	3163	enddo
[2632]	3164	#endif
[2630]	3165
	3166
	3167	IF (iflag_conv.GE.3) THEN
	3168	IF (logitime) THEN
	3169	CALL SYSTEM_CLOCK(COUNT=clock_start)
	3170	ENDIF
	3171
	3172
	3173	IF (lessivage) THEN
	3174	print *,' BEFORE lsc_scav '
	3175	IF (lminmax) THEN
	3176	DO it=1,nbtr
	3177	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_before_lsc_scav')
	3178	ENDDO
	3179	DO it=1,nbtr
	3180	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'before lsc_scav')
	3181	ENDDO
	3182	IF (lcheckmass) THEN
	3183	DO it=1,nbtr
	3184	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	3185	pplay,t_seri,iscm3,'before lsc_scav')
	3186	ENDDO
	3187	ENDIF
	3188	CALL minmaxsource(source_tr,qmin,qmax,'src: before lsc_scav')
	3189	ENDIF
	3190
	3191
	3192
	3193	ql_incloud_ref = 10.e-4
	3194	ql_incloud_ref = 5.e-4
	3195	! calcul du contenu en eau liquide au sein du nuage
	3196	ql_incl = ql_incloud_ref
	3197	! choix du lessivage
	3198	IF (iflag_lscav .EQ. 3 .OR. iflag_lscav .EQ. 4) THEN
[3786]	3199	!IF (.false.) THEN ! test #DFB (Binta) sans lsc_scav_spl
	3200	print *,'JE iflag_lscav',iflag_lscav
	3201	DO it = 1, nbtr
[2630]	3202
	3203	! incloud scavenging and removal by large scale rain ! orig : ql_incl
	3204	! was replaced by 0.5e-3 kg/kg
	3205	! the value 0.5e-3 kg/kg is from Giorgi and Chameides (1986), JGR
	3206	! Liu (2001) proposed to use 1.5e-3 kg/kg
	3207
	3208	! CALL lsc_scav_orig(pdtphys,it,iflag_lscav,ql_incl,prfl,psfl,
	3209	! . rneb,beta_fisrt, beta_v1,pplay,paprs,
	3210	! . t_seri,tr_seri,d_tr_insc,
	3211	! . d_tr_bcscav,d_tr_evapls,qPrls)
[3786]	3212	CALL lsc_scav_spl(pdtphys,it,iflag_lscav,ql_incl,prfl,psfl, &
[2630]	3213	rneb,beta_fisrt, beta_v1,pplay,paprs, &
	3214	t_seri,tr_seri,d_tr_insc, &
	3215	alpha_r,alpha_s,kk, henry, &
	3216	id_prec,id_fine,id_coss, id_codu, id_scdu, &
	3217	d_tr_bcscav,d_tr_evapls,qPrls)
	3218
	3219	!large scale scavenging tendency
[3786]	3220	DO k = 1, klev
	3221	DO i = 1, klon
	3222	d_tr_ls(i,k,it)=d_tr_insc(i,k,it)+d_tr_bcscav(i,k,it) &
[2630]	3223	+d_tr_evapls(i,k,it)
[3786]	3224	tr_seri(i,k,it)=tr_seri(i,k,it)+d_tr_ls(i,k,it)
	3225	tmp_var(i,k)=d_tr_ls(i,k,it)
	3226	ENDDO
	3227	ENDDO
[2630]	3228
[3786]	3229	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	3230
	3231	DO k=1,klev
	3232	DO i=1,klon
[2630]	3233	his_dhkelsc(i,it)=his_dhkelsc(i,it)-tmp_var(i,k) &
	3234	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3235
[3786]	3236	END DO
	3237	END DO
[2630]	3238
[3786]	3239	END DO !it=1,nbtr
	3240
[2630]	3241	ELSE
	3242	print *,'WARNING: NO lsc_scav, Please choose iflag_lscav=3 or 4'
[3786]	3243	DO it = 1, nbtr
	3244	DO i=1,klon
	3245	his_dhkelsc(i,it)=0.0
	3246	END DO ! klon
	3247	END DO !it=1,nbtr
	3248	ENDIF !iflag_lscav
[2630]	3249
	3250	print *,' AFTER lsc_scav '
	3251	IF (lminmax) THEN
	3252	DO it=1,nbtr
	3253	CALL checknanqfi(tr_seri(:,:,it),qmin,qmax,'nan_after_lsc_scav')
	3254	ENDDO
	3255	DO it=1,nbtr
	3256	CALL minmaxqfi2(tr_seri(:,:,it),qmin,qmax,'after lsc_scav')
	3257	ENDDO
	3258	IF (lcheckmass) THEN
	3259	DO it=1,nbtr
	3260	CALL checkmass(tr_seri(:,:,it),RNAVO,masse(it),zdz, &
	3261	pplay,t_seri,iscm3,'after lsc_scav')
	3262	ENDDO
	3263	ENDIF
	3264	CALL minmaxsource(source_tr,qmin,qmax,'src: after lsc_scav')
	3265	ENDIF
	3266
	3267	ENDIF ! lessivage
	3268
	3269	IF (logitime) THEN
	3270	CALL SYSTEM_CLOCK(COUNT=clock_end)
	3271	dife=clock_end-clock_start
	3272	ti_lscs=dife*MAX(0,SIGN(1,dife)) &
	3273	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	3274	tia_lscs=tia_lscs+REAL(ti_lscs)/REAL(clock_rate)
	3275	ENDIF
	3276
	3277
	3278
	3279	ENDIF !iflag_conv
	3280
	3281
	3282	!=======================================================================
	3283	! COMPUTING THE BURDEN
	3284	!=======================================================================
[2632]	3285	#ifdef IOPHYS_DUST
[2630]	3286	do it=1,nbtr
	3287	write(str2,'(i2.2)') it
	3288	call iophys_ecrit('TRJ'//str2,klev,'SOURCE','',tr_seri(:,:,it))
	3289	enddo
[2632]	3290	#endif
[2630]	3291
	3292	!
	3293	IF (logitime) THEN
	3294	CALL SYSTEM_CLOCK(COUNT=clock_start)
	3295	ENDIF
	3296
	3297
	3298	DO it=1, nbtr
	3299	DO j=1,klev
	3300	DO i=1,klon
	3301	tmp_var(i,j)=tr_seri(i,j,it)
	3302	ENDDO
	3303	ENDDO
	3304	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	3305	DO j=1,klev
	3306	DO i=1,klon
	3307	tr_seri(i,j,it)=tmp_var(i,j)
	3308	ENDDO
	3309	ENDDO
	3310	ENDDO
	3311	iscm3=.true.
	3312
	3313	!
	3314	! Computing burden in mg/m2
	3315	DO it=1, nbtr
	3316	DO k=1, klev
	3317	DO i=1, klon
	3318	trm(i,it)=trm(i,it)+tr_seri(i,k,it)1.e6zdz(i,k)* &
	3319	masse(it)*1.e3/RNAVO !--mg S/m2
	3320	ENDDO
	3321	ENDDO
	3322	ENDDO
	3323	!
	3324	! Computing Surface concentration in ug/m3
	3325	!
	3326	DO it=1, nbtr
	3327	DO i=1, klon
	3328	sconc_seri(i,it)=tr_seri(i,1,it)1.e6 &
	3329	masse(it)*1.e3/RNAVO !--mg/m3 (tr_seri ist in g/cm3)
	3330	ENDDO
	3331	ENDDO
	3332	!
	3333	!=======================================================================
	3334	! CALCULATION OF OPTICAL PROPERTIES
	3335	!=======================================================================
	3336	!
	3337	CALL aeropt_spl(zdz, tr_seri, RHcl, &
	3338	id_prec, id_fine, id_coss, id_codu, id_scdu, &
	3339	ok_chimeredust, &
	3340	diff_aod550_tot, diag_aod670_tot, diag_aod865_tot, &
	3341	diff_aod550_tr2, diag_aod670_tr2, diag_aod865_tr2, &
	3342	diag_aod550_ss, diag_aod670_ss, diag_aod865_ss, &
	3343	diag_aod550_dust,diag_aod670_dust,diag_aod865_dust, &
	3344	diag_aod550_dustsco,diag_aod670_dustsco,diag_aod865_dustsco)
	3345
	3346
	3347
	3348	IF (logitime) THEN
	3349	CALL SYSTEM_CLOCK(COUNT=clock_end)
	3350	dife=clock_end-clock_start
	3351	ti_brop=dife*MAX(0,SIGN(1,dife)) &
	3352	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	3353	tia_brop=tia_brop+REAL(ti_brop)/REAL(clock_rate)
	3354	ENDIF
	3355
	3356
	3357	!=======================================================================
	3358	! MODIS terra/aqua simulation output
	3359	!=======================================================================
	3360	masque_aqua_cur(:)=0
	3361	masque_terra_cur(:)=0
	3362
	3363	CALL satellite_out_spla(jD_cur,jH_cur,pdtphys,rlat,rlon, &
	3364	masque_aqua_cur, masque_terra_cur )
	3365	IF (jH_cur-pdtphys/86400. .LT. 0.) THEN
	3366	!new utc day: put in 0 everything
	3367	!JE20150518<<
	3368	masque_aqua(:) =0
	3369	masque_terra(:) =0
	3370	aod550_terra(:)=0.
	3371	aod550_tr2_terra(:)=0.
	3372	aod550_ss_terra(:)=0.
	3373	aod550_dust_terra(:)=0.
	3374	aod550_dustsco_terra(:)=0.
	3375	aod670_terra(:)=0.
	3376	aod670_tr2_terra(:)=0.
	3377	aod670_ss_terra(:)=0.
	3378	aod670_dust_terra(:)=0.
	3379	aod670_dustsco_terra(:)=0.
	3380	aod865_terra(:)=0.
	3381	aod865_tr2_terra(:)=0.
	3382	aod865_ss_terra(:)=0.
	3383	aod865_dust_terra(:)=0.
	3384	aod865_dustsco_terra(:)=0.
	3385	aod550_aqua(:)=0.
	3386	aod550_tr2_aqua(:)=0.
	3387	aod550_ss_aqua(:)=0.
	3388	aod550_dust_aqua(:)=0.
	3389	aod550_dustsco_aqua(:)=0.
	3390	aod670_aqua(:)=0.
	3391	aod670_tr2_aqua(:)=0.
	3392	aod670_ss_aqua(:)=0.
	3393	aod670_dust_aqua(:)=0.
	3394	aod670_dustsco_aqua(:)=0.
	3395	aod865_aqua(:)=0.
	3396	aod865_tr2_aqua(:)=0.
	3397	aod865_ss_aqua(:)=0.
	3398	aod865_dust_aqua(:)=0.
	3399	aod865_dustsco_aqua(:)=0.
	3400	!JE20150518>>
	3401	ENDIF
	3402
	3403	DO i=1,klon
	3404
	3405	aod550_terra(i)=aod550_terra(i)+ &
	3406	masque_terra_cur(i)*diff_aod550_tot(i)
	3407	aod550_tr2_terra(i)= aod550_tr2_terra(i)+ &
	3408	masque_terra_cur(i)*diff_aod550_tr2(i)
	3409	aod550_ss_terra(i)=aod550_ss_terra(i) + &
	3410	masque_terra_cur(i)*diag_aod550_ss(i)
	3411	aod550_dust_terra(i)= aod550_dust_terra(i) + &
	3412	masque_terra_cur(i)*diag_aod550_dust(i)
	3413	aod550_dustsco_terra(i)= aod550_dustsco_terra(i) + &
	3414	masque_terra_cur(i)*diag_aod550_dustsco(i)
	3415	aod670_terra(i)=aod670_terra(i)+ &
	3416	masque_terra_cur(i)*diag_aod670_tot(i)
	3417	aod670_tr2_terra(i)= aod670_tr2_terra(i)+ &
	3418	masque_terra_cur(i)*diag_aod670_tr2(i)
	3419	aod670_ss_terra(i)=aod670_ss_terra(i) + &
	3420	masque_terra_cur(i)*diag_aod670_ss(i)
	3421	aod670_dust_terra(i)= aod670_dust_terra(i) + &
	3422	masque_terra_cur(i)*diag_aod670_dust(i)
	3423	aod670_dustsco_terra(i)= aod670_dustsco_terra(i) + &
	3424	masque_terra_cur(i)*diag_aod670_dustsco(i)
	3425	aod865_terra(i)=aod865_terra(i)+ &
	3426	masque_terra_cur(i)*diag_aod865_tot(i)
	3427	aod865_tr2_terra(i)= aod865_tr2_terra(i)+ &
	3428	masque_terra_cur(i)*diag_aod865_tr2(i)
	3429	aod865_ss_terra(i)=aod865_ss_terra(i) + &
	3430	masque_terra_cur(i)*diag_aod865_ss(i)
	3431	aod865_dust_terra(i)= aod865_dust_terra(i) + &
	3432	masque_terra_cur(i)*diag_aod865_dust(i)
	3433	aod865_dustsco_terra(i)= aod865_dustsco_terra(i) + &
	3434	masque_terra_cur(i)*diag_aod865_dustsco(i)
	3435
	3436
	3437
	3438	aod550_aqua(i)=aod550_aqua(i)+ &
	3439	masque_aqua_cur(i)*diff_aod550_tot(i)
	3440	aod550_tr2_aqua(i)= aod550_tr2_aqua(i)+ &
	3441	masque_aqua_cur(i)*diff_aod550_tr2(i)
	3442	aod550_ss_aqua(i)=aod550_ss_aqua(i) + &
	3443	masque_aqua_cur(i)*diag_aod550_ss(i)
	3444	aod550_dust_aqua(i)= aod550_dust_aqua(i) + &
	3445	masque_aqua_cur(i)*diag_aod550_dust(i)
	3446	aod550_dustsco_aqua(i)= aod550_dustsco_aqua(i) + &
	3447	masque_aqua_cur(i)*diag_aod550_dustsco(i)
	3448	aod670_aqua(i)=aod670_aqua(i)+ &
	3449	masque_aqua_cur(i)*diag_aod670_tot(i)
	3450	aod670_tr2_aqua(i)= aod670_tr2_aqua(i)+ &
	3451	masque_aqua_cur(i)*diag_aod670_tr2(i)
	3452	aod670_ss_aqua(i)=aod670_ss_aqua(i) + &
	3453	masque_aqua_cur(i)*diag_aod670_ss(i)
	3454	aod670_dust_aqua(i)= aod670_dust_aqua(i) + &
	3455	masque_aqua_cur(i)*diag_aod670_dust(i)
	3456	aod670_dustsco_aqua(i)= aod670_dustsco_aqua(i) + &
	3457	masque_aqua_cur(i)*diag_aod670_dustsco(i)
	3458	aod865_aqua(i)=aod865_aqua(i)+ &
	3459	masque_aqua_cur(i)*diag_aod865_tot(i)
	3460	aod865_tr2_aqua(i)= aod865_tr2_aqua(i)+ &
	3461	masque_aqua_cur(i)*diag_aod865_tr2(i)
	3462	aod865_ss_aqua(i)=aod865_ss_aqua(i) + &
	3463	masque_aqua_cur(i)*diag_aod865_ss(i)
	3464	aod865_dust_aqua(i)= aod865_dust_aqua(i) + &
	3465	masque_aqua_cur(i)*diag_aod865_dust(i)
	3466	aod865_dustsco_aqua(i)= aod865_dustsco_aqua(i) + &
	3467	masque_aqua_cur(i)*diag_aod865_dustsco(i)
[3806]	3468
	3469	masque_aqua(i)=masque_aqua(i)+masque_aqua_cur(i)
[2630]	3470	masque_terra(i)=masque_terra(i)+masque_terra_cur(i)
	3471	ENDDO
	3472
	3473	IF (jH_cur+pdtphys/86400. .GE. 1.) THEN
	3474	! print *,'last step of the day'
	3475	DO i=1,klon
	3476	IF (masque_aqua(i).GT. 0) THEN
	3477	aod550_aqua(i)=aod550_aqua(i)/masque_aqua(i)
	3478	aod670_aqua(i)=aod670_aqua(i)/masque_aqua(i)
	3479	aod865_aqua(i)=aod865_aqua(i)/masque_aqua(i)
	3480	aod550_tr2_aqua(i)=aod550_tr2_aqua(i)/masque_aqua(i)
	3481	aod670_tr2_aqua(i)=aod670_tr2_aqua(i)/masque_aqua(i)
	3482	aod865_tr2_aqua(i)=aod865_tr2_aqua(i)/masque_aqua(i)
	3483	aod550_ss_aqua(i)=aod550_ss_aqua(i)/masque_aqua(i)
	3484	aod670_ss_aqua(i)=aod670_ss_aqua(i)/masque_aqua(i)
	3485	aod865_ss_aqua(i)=aod865_ss_aqua(i)/masque_aqua(i)
	3486	aod550_dust_aqua(i)=aod550_dust_aqua(i)/masque_aqua(i)
	3487	aod670_dust_aqua(i)=aod670_dust_aqua(i)/masque_aqua(i)
	3488	aod865_dust_aqua(i)=aod865_dust_aqua(i)/masque_aqua(i)
	3489	aod550_dustsco_aqua(i)=aod550_dustsco_aqua(i)/masque_aqua(i)
	3490	aod670_dustsco_aqua(i)=aod670_dustsco_aqua(i)/masque_aqua(i)
	3491	aod865_dustsco_aqua(i)=aod865_dustsco_aqua(i)/masque_aqua(i)
	3492	ELSE
	3493	aod550_aqua(i) = -999.
	3494	aod670_aqua(i) = -999.
	3495	aod865_aqua(i) = -999.
	3496	aod550_tr2_aqua(i)= -999.
	3497	aod670_tr2_aqua(i)= -999.
	3498	aod865_tr2_aqua(i)= -999.
	3499	aod550_ss_aqua(i)= -999.
	3500	aod670_ss_aqua(i)= -999.
	3501	aod865_ss_aqua(i)= -999.
	3502	aod550_dust_aqua(i)= -999.
	3503	aod670_dust_aqua(i)= -999.
	3504	aod865_dust_aqua(i)= -999.
	3505	aod550_dustsco_aqua(i)= -999.
	3506	aod670_dustsco_aqua(i)= -999.
	3507	aod865_dustsco_aqua(i)= -999.
	3508	ENDIF
	3509	IF (masque_terra(i).GT. 0) THEN
	3510	aod550_terra(i)=aod550_terra(i)/masque_terra(i)
	3511	aod670_terra(i)=aod670_terra(i)/masque_terra(i)
	3512	aod865_terra(i)=aod865_terra(i)/masque_terra(i)
	3513	aod550_tr2_terra(i)=aod550_tr2_terra(i)/masque_terra(i)
	3514	aod670_tr2_terra(i)=aod670_tr2_terra(i)/masque_terra(i)
	3515	aod865_tr2_terra(i)=aod865_tr2_terra(i)/masque_terra(i)
	3516	aod550_ss_terra(i)=aod550_ss_terra(i)/masque_terra(i)
	3517	aod670_ss_terra(i)=aod670_ss_terra(i)/masque_terra(i)
	3518	aod865_ss_terra(i)=aod865_ss_terra(i)/masque_terra(i)
	3519	aod550_dust_terra(i)=aod550_dust_terra(i)/masque_terra(i)
	3520	aod670_dust_terra(i)=aod670_dust_terra(i)/masque_terra(i)
	3521	aod865_dust_terra(i)=aod865_dust_terra(i)/masque_terra(i)
	3522	aod550_dustsco_terra(i)=aod550_dustsco_terra(i)/masque_terra(i)
	3523	aod670_dustsco_terra(i)=aod670_dustsco_terra(i)/masque_terra(i)
	3524	aod865_dustsco_terra(i)=aod865_dustsco_terra(i)/masque_terra(i)
	3525	ELSE
	3526	aod550_terra(i) = -999.
	3527	aod670_terra(i) = -999.
	3528	aod865_terra(i) = -999.
	3529	aod550_tr2_terra(i)= -999.
	3530	aod670_tr2_terra(i)= -999.
	3531	aod865_tr2_terra(i)= -999.
	3532	aod550_ss_terra(i)= -999.
	3533	aod670_ss_terra(i)= -999.
	3534	aod865_ss_terra(i)= -999.
	3535	aod550_dust_terra(i)= -999.
	3536	aod670_dust_terra(i)= -999.
	3537	aod865_dust_terra(i)= -999.
	3538	aod550_dustsco_terra(i)= -999.
	3539	aod670_dustsco_terra(i)= -999.
	3540	aod865_dustsco_terra(i)= -999.
	3541	ENDIF
	3542	ENDDO
	3543
[3806]	3544	!!AS deleting lines
	3545	!! IF (ok_histrac) THEN
	3546	!!!! write in output file
	3547	!!----many deleted lines
	3548	!! ENDIF !mpi_root
	3549	!!!$OMP END MASTER
	3550	!!!$OMP BARRIER
	3551	!! ENDIF !--ok_histrac
[2630]	3552
[3806]	3553	ENDIF ! jH_cur...
[2630]	3554
	3555
	3556	!
	3557	!======================================================================
	3558	! Stockage sur bande histoire
	3559	!======================================================================
[2632]	3560	#ifdef IOPHYS_DUST
[2630]	3561	do it=1,nbtr
	3562	write(str2,'(i2.2)') it
	3563	call iophys_ecrit('TRK'//str2,klev,'SOURCE','',tr_seri(:,:,it))
	3564	enddo
[2632]	3565	#endif
[2630]	3566
	3567
	3568	!
	3569	IF (logitime) THEN
	3570	CALL SYSTEM_CLOCK(COUNT=clock_start)
	3571	ENDIF
	3572
	3573	DO it=1, nbtr
	3574	DO j=1,klev
	3575	DO i=1,klon
	3576	tmp_var(i,j)=tr_seri(i,j,it)
	3577	ENDDO
	3578	ENDDO
	3579	CALL cm3_to_kg(pplay,t_seri,tmp_var)
	3580	DO j=1,klev
	3581	DO i=1,klon
	3582	tr_seri(i,j,it)=tmp_var(i,j)
	3583	ENDDO
	3584	ENDDO
	3585	ENDDO
	3586	iscm3=.false.
	3587
	3588	!
	3589	!
	3590	!======================================================================
	3591	! SAVING AEROSOL RELATED VARIABLES INTO FILE
	3592	!======================================================================
	3593	!
	3594	ndex2d = 0
	3595	ndex3d = 0
	3596	!
	3597	itra=itra+1
	3598
	3599	print *,'SAVING VARIABLES FOR DAY ',itra
	3600	!
	3601	fluxbb(:)=0.0
	3602	fluxff(:)=0.0
	3603	fluxbcbb(:)=0.0
	3604	fluxbcff(:)=0.0
	3605	fluxbcnff(:)=0.0
	3606	fluxbcba(:)=0.0
	3607	fluxbc(:)=0.0
	3608	fluxombb(:)=0.0
	3609	fluxomff(:)=0.0
	3610	fluxomnat(:)=0.0
	3611	fluxomba(:)=0.0
	3612	fluxomnff(:)=0.0
	3613	fluxom(:)=0.0
	3614	fluxh2sff(:)=0.0
	3615	fluxh2snff(:)=0.0
	3616	fluxh2sbio(:)=0.0
	3617	fluxso2ff(:)=0.0
	3618	fluxso2nff(:)=0.0
	3619	fluxso2bb(:)=0.0
	3620	fluxso2vol(:)=0.0
	3621	fluxso2ba(:)=0.0
	3622	fluxso2(:)=0.0
	3623	fluxso4ff(:)=0.0
	3624	fluxso4nff(:)=0.0
	3625	fluxso4bb(:)=0.0
	3626	fluxso4ba(:)=0.0
	3627	fluxso4(:)=0.0
	3628	fluxdms(:)=0.0
	3629	fluxdustec(:)=0.0
	3630	fluxddfine(:)=0.0
	3631	fluxddcoa(:)=0.0
	3632	fluxddsco(:)=0.0
	3633	fluxdd(:)=0.0
	3634	fluxssfine(:)=0.0
	3635	fluxsscoa(:)=0.0
	3636	fluxss(:)=0.0
	3637	DO i=1, klon
	3638	IF (iregion_ind(i).GT.0) THEN ! LAND
	3639	! SULFUR EMISSIONS
	3640	fluxh2sff(i)= (lmt_so2ff_l(i)+lmt_so2ff_h(i))frach2sofso2 &
	3641	scale_param_ind(iregion_ind(i))* &
	3642	1.e4/RNAVOmasse_s1.e3 ! mgS/m2/s
	3643	fluxso2ff(i)=scale_param_ind(iregion_ind(i)) * fracso2emis * &
	3644	(lmt_so2ff_l(i)+lmt_so2ff_h(i)) * 1.e4/RNAVO * &
	3645	masse_s * 1.e3 ! mgS/m2/s
	3646	! SULPHATE EMISSIONS
	3647	fluxso4ff(i)=scale_param_ind(iregion_ind(i))(1-fracso2emis) &
	3648	(lmt_so2ff_l(i)+lmt_so2ff_h(i)) * 1.e4/RNAVO * &
	3649	masse_s * 1.e3 ! mgS/m2/s
	3650	! BLACK CARBON EMISSIONS
	3651	fluxbcff(i)=scale_param_ff(iregion_ind(i))* &
	3652	lmt_bcff(i)1.e41.e3 !/g/m2/s
	3653	! ORGANIC MATTER EMISSIONS
	3654	fluxomff(i)=scale_param_ff(iregion_ind(i))* &
	3655	(lmt_omff(i))1.e41.e3 !/g/m2/s
	3656	! FOSSIL FUEL EMISSIONS
	3657	fluxff(i)=fluxbcff(i)+fluxomff(i)
	3658	ENDIF
	3659	IF (iregion_bb(i).GT.0) THEN ! LAND
	3660	! SULFUR EMISSIONS
	3661	fluxso2bb(i) =scale_param_bb(iregion_bb(i)) * fracso2emis * &
	3662	(lmt_so2bb_l(i)+lmt_so2bb_h(i))* &
	3663	(1.-pctsrf(i,is_oce))1.e4/RNAVOmasse_s*1.e3 ! mgS/m2/s
	3664	! SULPHATE EMISSIONS
	3665	fluxso4bb(i) =scale_param_bb(iregion_bb(i))(1-fracso2emis) &
	3666	(lmt_so2bb_l(i)+lmt_so2bb_h(i))* &
	3667	(1.-pctsrf(i,is_oce))1.e4/RNAVOmasse_s*1.e3 ! mgS/m2/s
	3668	! BLACK CARBON EMISSIONS
	3669	fluxbcbb(i)=scale_param_bb(iregion_bb(i))* &
	3670	(lmt_bcbb_l(i)+lmt_bcbb_h(i))1.e41.e3 !mg/m2/s
	3671	! ORGANIC MATTER EMISSIONS
	3672	fluxombb(i)=scale_param_bb(iregion_bb(i))* &
	3673	(lmt_ombb_l(i)+lmt_ombb_h(i))1.e41.e3 !mg/m2/s
	3674	! BIOMASS BURNING EMISSIONS
	3675	fluxbb(i)=fluxbcbb(i)+fluxombb(i)
	3676	ENDIF
	3677	! H2S EMISSIONS
	3678	fluxh2sbio(i)=lmt_h2sbio(i)1.e4/RNAVOmasse_s*1.e3 ! mgS/m2/s
	3679	fluxh2snff(i)= lmt_so2nff(i)frach2sofso2 &
	3680	1.e4/RNAVOmasse_s1.e3 ! mgS/m2/s
	3681	! SULFUR DIOXIDE EMISSIONS
	3682	fluxso2nff(i)=fracso2emis * lmt_so2nff(i) * 1.e4/RNAVO * &
	3683	masse_s * 1.e3 ! mgS/m2/s
	3684	fluxso2vol(i)=(lmt_so2volc_cont(i)+lmt_so2volc_expl(i)) &
	3685	1.e4/RNAVOmasse_s*1.e3 ! mgS/m2/s
	3686	fluxso2ba(i) =lmt_so2ba(i)1.e4/RNAVOmasse_s1.e3 &
	3687	fracso2emis ! mgS/m2/s
	3688	fluxso2(i)=fluxso2ff(i)+fluxso2bb(i)+fluxso2nff(i)+ &
	3689	fluxso2vol(i)+fluxso2ba(i)
	3690	! DMS EMISSIONS
	3691	fluxdms(i)=( lmt_dms(i)+lmt_dmsbio(i) ) &
	3692	1.e4/RNAVOmasse_s*1.e3 ! mgS/m2/s
	3693	! SULPHATE EMISSIONS
	3694	fluxso4ba(i) =lmt_so2ba(i)1.e4/RNAVOmasse_s*1.e3 &
	3695	*(1-fracso2emis) ! mgS/m2/s
	3696	fluxso4nff(i)=(1-fracso2emis)lmt_so2nff(i) 1.e4/RNAVO * &
	3697	masse_s * 1.e3 ! mgS/m2/s
	3698	fluxso4(i)=fluxso4ff(i)+fluxso4bb(i)+fluxso4ba(i)+fluxso4nff(i)
	3699	! BLACK CARBON EMISSIONS
	3700
	3701	fluxbcnff(i)=lmt_bcnff(i)1.e41.e3 !mg/m2/s
	3702	fluxbcba(i)=lmt_bcba(i)1.e41.e3 !mg/m2/s
	3703	fluxbc(i)=fluxbcbb(i)+fluxbcff(i)+fluxbcnff(i)+fluxbcba(i)
	3704	! ORGANIC MATTER EMISSIONS
	3705	fluxomnat(i)=lmt_omnat(i)1.e41.e3 !mg/m2/s
	3706	fluxomba(i)=lmt_omba(i)1.e41.e3 !mg/m2/s
	3707	fluxomnff(i)=lmt_omnff(i)1.e41.e3 !mg/m2/s
	3708	fluxom(i)=fluxombb(i)+fluxomff(i)+fluxomnat(i)+fluxomba(i)+ &
	3709	fluxomnff(i)
	3710	! DUST EMISSIONS
	3711	fluxdustec(i)=dust_ec(i)*1.e6 ! old dust emission scheme
	3712	!JE20140605<< old dust emission version
	3713	! fluxddfine(i)=scale_param_dustacc(iregion_dust(i))
	3714	! . * dust_ec(i)0.0931.e6
	3715	! fluxddcoa(i)=scale_param_dustcoa(iregion_dust(i))
	3716	! . * dust_ec(i)0.9051.e6
	3717	! fluxdd(i)=fluxddfine(i)+fluxddcoa(i)
	3718	!JE20140605>>
	3719	fluxddfine(i)=flux_sparam_ddfine(i)
	3720	fluxddcoa(i)=flux_sparam_ddcoa(i)
	3721	fluxddsco(i)=flux_sparam_ddsco(i)
	3722	fluxdd(i)=fluxddfine(i)+fluxddcoa(i)+fluxddsco(i)
	3723	! SEA SALT EMISSIONS
	3724	fluxssfine(i)=scale_param_ssacclmt_sea_salt(i,1)1.e4*1.e3
	3725	fluxsscoa(i)=scale_param_sscoalmt_sea_salt(i,2)1.e4*1.e3
	3726	fluxss(i)=fluxssfine(i)+fluxsscoa(i)
	3727	ENDDO
[3814]	3728
[2630]	3729	! prepare outputs cvltr
	3730
	3731	DO it=1, nbtr
	3732	DO k=1,klev
	3733	DO i=1,klon
	3734	tmp_var(i,k)=d_tr_cv(i,k,it)
	3735	ENDDO
	3736	ENDDO
	3737	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	3738	DO k=1,klev
	3739	DO i=1,klon
	3740	d_tr_cv_o(i,k,it)=tmp_var(i,k) &
	3741	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3742	ENDDO
	3743	ENDDO
	3744	ENDDO
	3745	DO it=1, nbtr
	3746	DO k=1,klev
	3747	DO i=1,klon
	3748	tmp_var(i,k)=d_tr_trsp(i,k,it)
	3749	ENDDO
	3750	ENDDO
	3751	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	3752	DO k=1,klev
	3753	DO i=1,klon
	3754	d_tr_trsp_o(i,k,it)=tmp_var(i,k) &
	3755	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3756	ENDDO
	3757	ENDDO
	3758	ENDDO
	3759	DO it=1, nbtr
	3760	DO k=1,klev
	3761	DO i=1,klon
	3762	tmp_var(i,k)=d_tr_sscav(i,k,it)
	3763	ENDDO
	3764	ENDDO
	3765	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	3766	DO k=1,klev
	3767	DO i=1,klon
	3768	d_tr_sscav_o(i,k,it)=tmp_var(i,k) &
	3769	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3770	ENDDO
	3771	ENDDO
	3772	ENDDO
	3773	DO it=1, nbtr
	3774	DO k=1,klev
	3775	DO i=1,klon
	3776	tmp_var(i,k)=d_tr_sat(i,k,it)
	3777	ENDDO
	3778	ENDDO
	3779	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	3780	DO k=1,klev
	3781	DO i=1,klon
	3782	d_tr_sat_o(i,k,it)=tmp_var(i,k) &
	3783	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3784	ENDDO
	3785	ENDDO
	3786	ENDDO
	3787	DO it=1, nbtr
	3788	DO k=1,klev
	3789	DO i=1,klon
	3790	tmp_var(i,k)=d_tr_uscav(i,k,it)
	3791	ENDDO
	3792	ENDDO
	3793	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	3794	DO k=1,klev
	3795	DO i=1,klon
	3796	d_tr_uscav_o(i,k,it)=tmp_var(i,k) &
	3797	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3798	ENDDO
	3799	ENDDO
	3800	ENDDO
	3801	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	3802	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	3803	DO it=1, nbtr
	3804	DO k=1,klev
	3805	DO i=1,klon
	3806	tmp_var(i,k)=d_tr_insc(i,k,it)
	3807	ENDDO
	3808	ENDDO
	3809	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	3810	DO k=1,klev
	3811	DO i=1,klon
	3812	d_tr_insc_o(i,k,it)=tmp_var(i,k) &
	3813	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3814	ENDDO
	3815	ENDDO
	3816	ENDDO
	3817
	3818
	3819	DO it=1, nbtr
	3820	DO k=1,klev
	3821	DO i=1,klon
	3822	tmp_var(i,k)=d_tr_bcscav(i,k,it)
	3823	ENDDO
	3824	ENDDO
	3825	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	3826	DO k=1,klev
	3827	DO i=1,klon
	3828	d_tr_bcscav_o(i,k,it)=tmp_var(i,k) &
	3829	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3830	ENDDO
	3831	ENDDO
	3832	ENDDO
	3833
	3834
	3835	DO it=1, nbtr
	3836	DO k=1,klev
	3837	DO i=1,klon
	3838	tmp_var(i,k)=d_tr_evapls(i,k,it)
	3839	ENDDO
	3840	ENDDO
	3841	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	3842	DO k=1,klev
	3843	DO i=1,klon
	3844	d_tr_evapls_o(i,k,it)=tmp_var(i,k) &
	3845	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3846	ENDDO
	3847	ENDDO
	3848	ENDDO
	3849
	3850
	3851	DO it=1, nbtr
	3852	DO k=1,klev
	3853	DO i=1,klon
	3854	tmp_var(i,k)=d_tr_ls(i,k,it)
	3855	ENDDO
	3856	ENDDO
	3857	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	3858	DO k=1,klev
	3859	DO i=1,klon
	3860	d_tr_ls_o(i,k,it)=tmp_var(i,k) &
	3861	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3862	ENDDO
	3863	ENDDO
	3864	ENDDO
	3865
	3866
	3867	DO it=1, nbtr
	3868	DO k=1,klev
	3869	DO i=1,klon
	3870	tmp_var(i,k)=d_tr_dyn(i,k,it)
	3871	ENDDO
	3872	ENDDO
	3873	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	3874	DO k=1,klev
	3875	DO i=1,klon
	3876	d_tr_dyn_o(i,k,it)=tmp_var(i,k) &
	3877	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3878	ENDDO
	3879	ENDDO
	3880	ENDDO
	3881
	3882
	3883	DO it=1, nbtr
	3884	DO k=1,klev
	3885	DO i=1,klon
	3886	tmp_var(i,k)=d_tr_cl(i,k,it)
	3887	ENDDO
	3888	ENDDO
	3889	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	3890	DO k=1,klev
	3891	DO i=1,klon
	3892	d_tr_cl_o(i,k,it)=tmp_var(i,k) &
	3893	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3894	ENDDO
	3895	ENDDO
	3896	ENDDO
	3897
	3898
	3899	DO it=1, nbtr
	3900	DO k=1,klev
	3901	DO i=1,klon
	3902	tmp_var(i,k)=d_tr_th(i,k,it)
	3903	ENDDO
	3904	ENDDO
	3905	CALL kg_to_cm3(pplay,t_seri,tmp_var)
	3906	DO k=1,klev
	3907	DO i=1,klon
	3908	d_tr_th_o(i,k,it)=tmp_var(i,k) &
	3909	/RNAVOmasse(it)1.e31.e6zdz(i,k)/pdtphys
	3910	ENDDO
	3911	ENDDO
	3912	ENDDO
	3913	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	3914	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	3915
	3916	DO it=1,nbtr
	3917	WRITE(str2,'(i2.2)') it
	3918	DO i=1, klon
	3919	his_dh(i,it)= his_dhlsc(i,it)+his_dhcon(i,it)+ &
	3920	his_dhbclsc(i,it)+his_dhbccon(i,it)
	3921
	3922	ENDDO
	3923	ENDDO
	3924
[3806]	3925	!AS: commenting out and deleting lines
	3926	!! IF (ok_histrac) THEN
	3927	!!
	3928	!! SAVING VARIABLES IN TRACEUR
	3929	!!----- many lines deleted----
	3930	!! ENDIF ! ok_histrac
[2630]	3931
	3932
	3933
	3934
	3935	!JE20141224
	3936	! saving variables for output
	3937	! 2D outputs
	3938	DO i=1, klon
	3939	trm01(i)=0.
	3940	trm02(i)=0.
	3941	trm03(i)=0.
	3942	trm04(i)=0.
	3943	trm05(i)=0.
	3944	sconc01(i)=0.
	3945	sconc02(i)=0.
	3946	sconc03(i)=0.
	3947	sconc04(i)=0.
	3948	sconc05(i)=0.
	3949	flux01(i)=0.
	3950	flux02(i)=0.
	3951	flux03(i)=0.
	3952	flux04(i)=0.
	3953	flux05(i)=0.
	3954	ds01(i)=0.
	3955	ds02(i)=0.
	3956	ds03(i)=0.
	3957	ds04(i)=0.
	3958	ds05(i)=0.
	3959	dh01(i)=0.
	3960	dh02(i)=0.
	3961	dh03(i)=0.
	3962	dh04(i)=0.
	3963	dh05(i)=0.
	3964	dtrconv01(i)=0.
	3965	dtrconv02(i)=0.
	3966	dtrconv03(i)=0.
	3967	dtrconv04(i)=0.
	3968	dtrconv05(i)=0.
	3969	dtherm01(i)=0.
	3970	dtherm02(i)=0.
	3971	dtherm03(i)=0.
	3972	dtherm04(i)=0.
	3973	dtherm05(i)=0.
	3974	dhkecv01(i)=0.
	3975	dhkecv02(i)=0.
	3976	dhkecv03(i)=0.
	3977	dhkecv04(i)=0.
	3978	dhkecv05(i)=0.
	3979	d_tr_ds01(i)=0.
	3980	d_tr_ds02(i)=0.
	3981	d_tr_ds03(i)=0.
	3982	d_tr_ds04(i)=0.
	3983	d_tr_ds05(i)=0.
	3984	dhkelsc01(i)=0.
	3985	dhkelsc02(i)=0.
	3986	dhkelsc03(i)=0.
	3987	dhkelsc04(i)=0.
	3988	dhkelsc05(i)=0.
	3989	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	3990
	3991	if(id_prec>0) trm01(i)=trm(i,id_prec)
	3992	if(id_fine>0) trm02(i)=trm(i,id_fine)
	3993	if(id_coss>0) trm03(i)=trm(i,id_coss)
	3994	if(id_codu>0) trm04(i)=trm(i,id_codu)
	3995	if(id_scdu>0) trm05(i)=trm(i,id_scdu)
	3996	if(id_prec>0) sconc01(i)=sconc_seri(i,id_prec)
	3997	if(id_fine>0) sconc02(i)=sconc_seri(i,id_fine)
	3998	if(id_coss>0) sconc03(i)=sconc_seri(i,id_coss)
	3999	if(id_codu>0) sconc04(i)=sconc_seri(i,id_codu)
	4000	if(id_scdu>0) sconc05(i)=sconc_seri(i,id_scdu)
	4001	if(id_prec>0) flux01(i)=flux_tr(i,id_prec)
	4002	if(id_fine>0) flux02(i)=flux_tr(i,id_fine)
	4003	if(id_coss>0) flux03(i)=flux_tr(i,id_coss)
	4004	if(id_codu>0) flux04(i)=flux_tr(i,id_codu)
	4005	if(id_scdu>0) flux05(i)=flux_tr(i,id_scdu)
	4006	if(id_prec>0) ds01(i)=his_ds(i,id_prec)
	4007	if(id_fine>0) ds02(i)=his_ds(i,id_fine)
	4008	if(id_coss>0) ds03(i)=his_ds(i,id_coss)
	4009	if(id_codu>0) ds04(i)=his_ds(i,id_codu)
	4010	if(id_scdu>0) ds05(i)=his_ds(i,id_scdu)
	4011	if(id_prec>0) dh01(i)=his_dh(i,id_prec)
	4012	if(id_fine>0) dh02(i)=his_dh(i,id_fine)
	4013	if(id_coss>0) dh03(i)=his_dh(i,id_coss)
	4014	if(id_codu>0) dh04(i)=his_dh(i,id_codu)
	4015	if(id_scdu>0) dh05(i)=his_dh(i,id_scdu)
	4016	if(id_prec>0) dtrconv01(i)=dtrconv(i,id_prec)
	4017	if(id_fine>0) dtrconv02(i)=dtrconv(i,id_fine)
	4018	if(id_coss>0) dtrconv03(i)=dtrconv(i,id_coss)
	4019	if(id_codu>0) dtrconv04(i)=dtrconv(i,id_codu)
	4020	if(id_scdu>0) dtrconv05(i)=dtrconv(i,id_scdu)
	4021	if(id_prec>0) dtherm01(i)=his_th(i,id_prec)
	4022	if(id_fine>0) dtherm02(i)=his_th(i,id_fine)
	4023	if(id_coss>0) dtherm03(i)=his_th(i,id_coss)
	4024	if(id_codu>0) dtherm04(i)=his_th(i,id_codu)
	4025	if(id_scdu>0) dtherm05(i)=his_th(i,id_scdu)
	4026	if(id_prec>0) dhkecv01(i)=his_dhkecv(i,id_prec)
	4027	if(id_fine>0) dhkecv02(i)=his_dhkecv(i,id_fine)
	4028	if(id_coss>0) dhkecv03(i)=his_dhkecv(i,id_coss)
	4029	if(id_codu>0) dhkecv04(i)=his_dhkecv(i,id_codu)
	4030	if(id_scdu>0) dhkecv05(i)=his_dhkecv(i,id_scdu)
	4031	if(id_prec>0) d_tr_ds01(i)=his_ds(i,id_prec)
	4032	if(id_fine>0) d_tr_ds02(i)=his_ds(i,id_fine)
	4033	if(id_coss>0) d_tr_ds03(i)=his_ds(i,id_coss)
	4034	if(id_codu>0) d_tr_ds04(i)=his_ds(i,id_codu)
	4035	if(id_scdu>0) d_tr_ds05(i)=his_ds(i,id_scdu)
	4036	if(id_prec>0) dhkelsc01(i)=his_dhkelsc(i,id_prec)
	4037	if(id_fine>0) dhkelsc02(i)=his_dhkelsc(i,id_fine)
	4038	if(id_coss>0) dhkelsc03(i)=his_dhkelsc(i,id_coss)
	4039	if(id_codu>0) dhkelsc04(i)=his_dhkelsc(i,id_codu)
	4040	if(id_scdu>0) dhkelsc05(i)=his_dhkelsc(i,id_scdu)
	4041	u10m_ss(i)=u10m_ec(i)
	4042	v10m_ss(i)=v10m_ec(i)
	4043	ENDDO
	4044	! 3D outs
	4045	DO i=1, klon
	4046	DO k=1,klev
	4047	d_tr_cv01(i,k) =0.
	4048	d_tr_cv02(i,k) =0.
	4049	d_tr_cv03(i,k) =0.
	4050	d_tr_cv04(i,k) =0.
	4051	d_tr_cv05(i,k) =0.
	4052	d_tr_trsp01(i,k) =0.
	4053	d_tr_trsp02(i,k) =0.
	4054	d_tr_trsp03(i,k) =0.
	4055	d_tr_trsp04(i,k) =0.
	4056	d_tr_trsp05(i,k) =0.
	4057	d_tr_sscav01(i,k)=0.
	4058	d_tr_sscav02(i,k)=0.
	4059	d_tr_sscav03(i,k)=0.
	4060	d_tr_sscav04(i,k)=0.
	4061	d_tr_sscav05(i,k)=0.
	4062	d_tr_sat01(i,k) =0.
	4063	d_tr_sat02(i,k) =0.
	4064	d_tr_sat03(i,k) =0.
	4065	d_tr_sat04(i,k) =0.
	4066	d_tr_sat05(i,k) =0.
	4067	d_tr_uscav01(i,k)=0.
	4068	d_tr_uscav02(i,k)=0.
	4069	d_tr_uscav03(i,k)=0.
	4070	d_tr_uscav04(i,k)=0.
	4071	d_tr_uscav05(i,k)=0.
	4072	d_tr_insc01(i,k)=0.!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	4073	d_tr_insc02(i,k)=0.
	4074	d_tr_insc03(i,k)=0.
	4075	d_tr_insc04(i,k)=0.
	4076	d_tr_insc05(i,k)=0.
	4077	d_tr_bcscav01(i,k)=0.
	4078	d_tr_bcscav02(i,k)=0.
	4079	d_tr_bcscav03(i,k)=0.
	4080	d_tr_bcscav04(i,k)=0.
	4081	d_tr_bcscav05(i,k)=0.
	4082	d_tr_evapls01(i,k)=0.
	4083	d_tr_evapls02(i,k)=0.
	4084	d_tr_evapls03(i,k)=0.
	4085	d_tr_evapls04(i,k)=0.
	4086	d_tr_evapls05(i,k)=0.
	4087	d_tr_ls01(i,k)=0.
	4088	d_tr_ls02(i,k)=0.
	4089	d_tr_ls03(i,k)=0.
	4090	d_tr_ls04(i,k)=0.
	4091	d_tr_ls05(i,k)=0.!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	4092	d_tr_dyn01(i,k)=0.
	4093	d_tr_dyn02(i,k)=0.
	4094	d_tr_dyn03(i,k)=0.
	4095	d_tr_dyn04(i,k)=0.
	4096	d_tr_dyn05(i,k)=0.!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	4097	d_tr_cl01(i,k)=0.
	4098	d_tr_cl02(i,k)=0.
	4099	d_tr_cl03(i,k)=0.
	4100	d_tr_cl04(i,k)=0.
	4101	d_tr_cl05(i,k)=0.!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	4102	d_tr_th01(i,k)=0.
	4103	d_tr_th02(i,k)=0.
	4104	d_tr_th03(i,k)=0.
	4105	d_tr_th04(i,k)=0.
	4106	d_tr_th05(i,k)=0.!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[3786]	4107	ENDDO
	4108	ENDDO
	4109
	4110	IF(1==0) THEN
	4111	! calcul in original trunk version; problem: budget not closed. Corrected in "ELSE"
	4112	DO i=1, klon
	4113	DO k=1,klev
	4114
[2630]	4115	if(id_prec>0) d_tr_cv01(i,k) =d_tr_cv_o(i,k,id_prec)
	4116	if(id_fine>0) d_tr_cv02(i,k) =d_tr_cv_o(i,k,id_fine)
	4117	if(id_coss>0) d_tr_cv03(i,k) =d_tr_cv_o(i,k,id_coss)
	4118	if(id_codu>0) d_tr_cv04(i,k) =d_tr_cv_o(i,k,id_codu)
	4119	if(id_scdu>0) d_tr_cv05(i,k) =d_tr_cv_o(i,k,id_scdu)
	4120	if(id_prec>0) d_tr_trsp01(i,k) =d_tr_trsp_o(i,k,id_prec)
	4121	if(id_fine>0) d_tr_trsp02(i,k) =d_tr_trsp_o(i,k,id_fine)
	4122	if(id_coss>0) d_tr_trsp03(i,k) =d_tr_trsp_o(i,k,id_coss)
	4123	if(id_codu>0) d_tr_trsp04(i,k) =d_tr_trsp_o(i,k,id_codu)
	4124	if(id_scdu>0) d_tr_trsp05(i,k) =d_tr_trsp_o(i,k,id_scdu)
	4125	if(id_prec>0) d_tr_sscav01(i,k)=d_tr_sscav_o(i,k,id_prec)
	4126	if(id_fine>0) d_tr_sscav02(i,k)=d_tr_sscav_o(i,k,id_fine)
	4127	if(id_coss>0) d_tr_sscav03(i,k)=d_tr_sscav_o(i,k,id_coss)
	4128	if(id_codu>0) d_tr_sscav04(i,k)=d_tr_sscav_o(i,k,id_codu)
	4129	if(id_scdu>0) d_tr_sscav05(i,k)=d_tr_sscav_o(i,k,id_scdu)
	4130	if(id_prec>0) d_tr_sat01(i,k) =d_tr_sat_o(i,k,id_prec)
	4131	if(id_fine>0) d_tr_sat02(i,k) =d_tr_sat_o(i,k,id_fine)
	4132	if(id_coss>0) d_tr_sat03(i,k) =d_tr_sat_o(i,k,id_coss)
	4133	if(id_codu>0) d_tr_sat04(i,k) =d_tr_sat_o(i,k,id_codu)
	4134	if(id_scdu>0) d_tr_sat05(i,k) =d_tr_sat_o(i,k,id_scdu)
	4135	if(id_prec>0) d_tr_uscav01(i,k)=d_tr_uscav_o(i,k,id_prec)
	4136	if(id_fine>0) d_tr_uscav02(i,k)=d_tr_uscav_o(i,k,id_fine)
	4137	if(id_coss>0) d_tr_uscav03(i,k)=d_tr_uscav_o(i,k,id_coss)
	4138	if(id_codu>0) d_tr_uscav04(i,k)=d_tr_uscav_o(i,k,id_codu)
	4139	if(id_scdu>0) d_tr_uscav05(i,k)=d_tr_uscav_o(i,k,id_scdu)
	4140	if(id_prec>0) d_tr_insc01(i,k)=d_tr_insc_o(i,k,id_prec)
	4141	if(id_fine>0) d_tr_insc02(i,k)=d_tr_insc_o(i,k,id_fine)
	4142	if(id_coss>0) d_tr_insc03(i,k)=d_tr_insc_o(i,k,id_coss)
	4143	if(id_codu>0) d_tr_insc04(i,k)=d_tr_insc_o(i,k,id_codu)
	4144	if(id_scdu>0) d_tr_insc05(i,k)=d_tr_insc_o(i,k,id_scdu)
	4145	if(id_prec>0) d_tr_bcscav01(i,k)=d_tr_bcscav_o(i,k,id_prec)
	4146	if(id_fine>0) d_tr_bcscav02(i,k)=d_tr_bcscav_o(i,k,id_fine)
	4147	if(id_coss>0) d_tr_bcscav03(i,k)=d_tr_bcscav_o(i,k,id_coss)
	4148	if(id_codu>0) d_tr_bcscav04(i,k)=d_tr_bcscav_o(i,k,id_codu)
	4149	if(id_scdu>0) d_tr_bcscav05(i,k)=d_tr_bcscav_o(i,k,id_scdu)
	4150	if(id_prec>0) d_tr_evapls01(i,k)=d_tr_evapls_o(i,k,id_prec)
	4151	if(id_fine>0) d_tr_evapls02(i,k)=d_tr_evapls_o(i,k,id_fine)
	4152	if(id_coss>0) d_tr_evapls03(i,k)=d_tr_evapls_o(i,k,id_coss)
	4153	if(id_codu>0) d_tr_evapls04(i,k)=d_tr_evapls_o(i,k,id_codu)
	4154	if(id_scdu>0) d_tr_evapls05(i,k)=d_tr_evapls_o(i,k,id_scdu)
	4155	ENDDO
	4156	ENDDO
[3786]	4157	ELSE ! correction pour fermeture de bilan, par FH dans les simus de Binta pour Habib
[2630]	4158	DO i=1, klon
	4159	DO k=1,klev
[3786]	4160	if(id_prec>0) d_tr_cv01(i,k) =d_tr_cv(i,k,id_prec)/pdtphys
	4161	if(id_fine>0) d_tr_cv02(i,k) =d_tr_cv(i,k,id_fine)/pdtphys
	4162	if(id_coss>0) d_tr_cv03(i,k) =d_tr_cv(i,k,id_coss)/pdtphys
	4163	if(id_codu>0) d_tr_cv04(i,k) =d_tr_cv(i,k,id_codu)/pdtphys
	4164	if(id_scdu>0) d_tr_cv05(i,k) =d_tr_cv(i,k,id_scdu)/pdtphys
	4165	if(id_prec>0) d_tr_trsp01(i,k) =d_tr_trsp(i,k,id_prec)/pdtphys
	4166	if(id_fine>0) d_tr_trsp02(i,k) =d_tr_trsp(i,k,id_fine)/pdtphys
	4167	if(id_coss>0) d_tr_trsp03(i,k) =d_tr_trsp(i,k,id_coss)/pdtphys
	4168	if(id_codu>0) d_tr_trsp04(i,k) =d_tr_trsp(i,k,id_codu)/pdtphys
	4169	if(id_scdu>0) d_tr_trsp05(i,k) =d_tr_trsp(i,k,id_scdu)/pdtphys
	4170	if(id_prec>0) d_tr_sscav01(i,k)=d_tr_sscav(i,k,id_prec)/pdtphys
	4171	if(id_fine>0) d_tr_sscav02(i,k)=d_tr_sscav(i,k,id_fine)/pdtphys
	4172	if(id_coss>0) d_tr_sscav03(i,k)=d_tr_sscav(i,k,id_coss)/pdtphys
	4173	if(id_codu>0) d_tr_sscav04(i,k)=d_tr_sscav(i,k,id_codu)/pdtphys
	4174	if(id_scdu>0) d_tr_sscav05(i,k)=d_tr_sscav(i,k,id_scdu)/pdtphys
	4175	if(id_prec>0) d_tr_sat01(i,k) =d_tr_sat(i,k,id_prec)/pdtphys
	4176	if(id_fine>0) d_tr_sat02(i,k) =d_tr_sat(i,k,id_fine)/pdtphys
	4177	if(id_coss>0) d_tr_sat03(i,k) =d_tr_sat(i,k,id_coss)/pdtphys
	4178	if(id_codu>0) d_tr_sat04(i,k) =d_tr_sat(i,k,id_codu)/pdtphys
	4179	if(id_scdu>0) d_tr_sat05(i,k) =d_tr_sat(i,k,id_scdu)/pdtphys
	4180	if(id_prec>0) d_tr_uscav01(i,k)=d_tr_uscav(i,k,id_prec)/pdtphys
	4181	if(id_fine>0) d_tr_uscav02(i,k)=d_tr_uscav(i,k,id_fine)/pdtphys
	4182	if(id_coss>0) d_tr_uscav03(i,k)=d_tr_uscav(i,k,id_coss)/pdtphys
	4183	if(id_codu>0) d_tr_uscav04(i,k)=d_tr_uscav(i,k,id_codu)/pdtphys
	4184	if(id_scdu>0) d_tr_uscav05(i,k)=d_tr_uscav(i,k,id_scdu)/pdtphys
	4185	if(id_prec>0) d_tr_insc01(i,k)=d_tr_insc(i,k,id_prec)/pdtphys
	4186	if(id_fine>0) d_tr_insc02(i,k)=d_tr_insc(i,k,id_fine)/pdtphys
	4187	if(id_coss>0) d_tr_insc03(i,k)=d_tr_insc(i,k,id_coss)/pdtphys
	4188	if(id_codu>0) d_tr_insc04(i,k)=d_tr_insc(i,k,id_codu)/pdtphys
	4189	if(id_scdu>0) d_tr_insc05(i,k)=d_tr_insc(i,k,id_scdu)/pdtphys
	4190	if(id_prec>0) d_tr_bcscav01(i,k)=d_tr_bcscav(i,k,id_prec)/pdtphys
	4191	if(id_fine>0) d_tr_bcscav02(i,k)=d_tr_bcscav(i,k,id_fine)/pdtphys
	4192	if(id_coss>0) d_tr_bcscav03(i,k)=d_tr_bcscav(i,k,id_coss)/pdtphys
	4193	if(id_codu>0) d_tr_bcscav04(i,k)=d_tr_bcscav(i,k,id_codu)/pdtphys
	4194	if(id_scdu>0) d_tr_bcscav05(i,k)=d_tr_bcscav(i,k,id_scdu)/pdtphys
	4195	if(id_prec>0) d_tr_evapls01(i,k)=d_tr_evapls(i,k,id_prec)/pdtphys
	4196	if(id_fine>0) d_tr_evapls02(i,k)=d_tr_evapls(i,k,id_fine)/pdtphys
	4197	if(id_coss>0) d_tr_evapls03(i,k)=d_tr_evapls(i,k,id_coss)/pdtphys
	4198	if(id_codu>0) d_tr_evapls04(i,k)=d_tr_evapls(i,k,id_codu)/pdtphys
	4199	if(id_scdu>0) d_tr_evapls05(i,k)=d_tr_evapls(i,k,id_scdu)/pdtphys
	4200	ENDDO
	4201	ENDDO
	4202	ENDIF
	4203
	4204	IF(1==0) THEN ! This "if" is as in original trunk
	4205	DO i=1, klon
	4206	DO k=1,klev
[2630]	4207	if(id_prec>0) d_tr_ls01(i,k)=d_tr_ls_o(i,k,id_prec)
	4208	if(id_fine>0) d_tr_ls02(i,k)=d_tr_ls_o(i,k,id_fine)
	4209	if(id_coss>0) d_tr_ls03(i,k)=d_tr_ls_o(i,k,id_coss)
	4210	if(id_codu>0) d_tr_ls04(i,k)=d_tr_ls_o(i,k,id_codu)
	4211	if(id_scdu>0) d_tr_ls05(i,k)=d_tr_ls_o(i,k,id_scdu)
	4212	if(id_prec>0) d_tr_dyn01(i,k)=d_tr_dyn_o(i,k,id_prec)
	4213	if(id_fine>0) d_tr_dyn02(i,k)=d_tr_dyn_o(i,k,id_fine)
	4214	if(id_coss>0) d_tr_dyn03(i,k)=d_tr_dyn_o(i,k,id_coss)
	4215	if(id_codu>0) d_tr_dyn04(i,k)=d_tr_dyn_o(i,k,id_codu)
	4216	if(id_scdu>0) d_tr_dyn05(i,k)=d_tr_dyn_o(i,k,id_scdu)
	4217	if(id_prec>0) d_tr_cl01(i,k)=d_tr_cl_o(i,k,id_prec)
	4218	if(id_fine>0) d_tr_cl02(i,k)=d_tr_cl_o(i,k,id_fine)
	4219	if(id_coss>0) d_tr_cl03(i,k)=d_tr_cl_o(i,k,id_coss)
	4220	if(id_codu>0) d_tr_cl04(i,k)=d_tr_cl_o(i,k,id_codu)
	4221	if(id_scdu>0) d_tr_cl05(i,k)=d_tr_cl_o(i,k,id_scdu)
	4222	if(id_prec>0) d_tr_th01(i,k)=d_tr_th_o(i,k,id_prec)
	4223	if(id_fine>0) d_tr_th02(i,k)=d_tr_th_o(i,k,id_fine)
	4224	if(id_coss>0) d_tr_th03(i,k)=d_tr_th_o(i,k,id_coss)
	4225	if(id_codu>0) d_tr_th04(i,k)=d_tr_th_o(i,k,id_codu)
	4226	if(id_scdu>0) d_tr_th05(i,k)=d_tr_th_o(i,k,id_scdu)
	4227	ENDDO
	4228	ENDDO
	4229	ELSE
	4230	DO i=1, klon
	4231	DO k=1,klev
	4232	if(id_prec>0) d_tr_ls01(i,k)=d_tr_ls(i,k,id_prec)/pdtphys
	4233	if(id_fine>0) d_tr_ls02(i,k)=d_tr_ls(i,k,id_fine)/pdtphys
	4234	if(id_coss>0) d_tr_ls03(i,k)=d_tr_ls(i,k,id_coss)/pdtphys
	4235	if(id_codu>0) d_tr_ls04(i,k)=d_tr_ls(i,k,id_codu)/pdtphys
	4236	if(id_scdu>0) d_tr_ls05(i,k)=d_tr_ls(i,k,id_scdu)/pdtphys
	4237	if(id_prec>0) d_tr_dyn01(i,k)=d_tr_dyn(i,k,id_prec)/pdtphys
	4238	if(id_fine>0) d_tr_dyn02(i,k)=d_tr_dyn(i,k,id_fine)/pdtphys
	4239	if(id_coss>0) d_tr_dyn03(i,k)=d_tr_dyn(i,k,id_coss)/pdtphys
	4240	if(id_codu>0) d_tr_dyn04(i,k)=d_tr_dyn(i,k,id_codu)/pdtphys
	4241	if(id_scdu>0) d_tr_dyn05(i,k)=d_tr_dyn(i,k,id_scdu)/pdtphys
	4242	if(id_prec>0) d_tr_cl01(i,k)=d_tr_cl(i,k,id_prec)/pdtphys
	4243	if(id_fine>0) d_tr_cl02(i,k)=d_tr_cl(i,k,id_fine)/pdtphys
	4244	if(id_coss>0) d_tr_cl03(i,k)=d_tr_cl(i,k,id_coss)/pdtphys
	4245	if(id_codu>0) d_tr_cl04(i,k)=d_tr_cl(i,k,id_codu)/pdtphys
	4246	if(id_scdu>0) d_tr_cl05(i,k)=d_tr_cl(i,k,id_scdu)/pdtphys
	4247	if(id_prec>0) d_tr_th01(i,k)=d_tr_th(i,k,id_prec)/pdtphys
	4248	if(id_fine>0) d_tr_th02(i,k)=d_tr_th(i,k,id_fine)/pdtphys
	4249	if(id_coss>0) d_tr_th03(i,k)=d_tr_th(i,k,id_coss)/pdtphys
	4250	if(id_codu>0) d_tr_th04(i,k)=d_tr_th(i,k,id_codu)/pdtphys
	4251	if(id_scdu>0) d_tr_th05(i,k)=d_tr_th(i,k,id_scdu)/pdtphys
	4252	ENDDO
	4253	ENDDO
	4254	ENDIF
	4255
	4256
	4257	IF (logitime) THEN
	4258	CALL SYSTEM_CLOCK(COUNT=clock_end)
	4259
	4260	dife=clock_end-clock_start
	4261	ti_outs=dife*MAX(0,SIGN(1,dife)) &
	4262	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	4263	tia_outs=tia_outs+REAL(ti_outs)/REAL(clock_rate)
	4264	ENDIF
	4265
	4266	IF (logitime) THEN
	4267	CALL SYSTEM_CLOCK(COUNT=clock_end)
	4268
	4269	dife=clock_end-clock_start_spla
	4270	ti_spla=dife*MAX(0,SIGN(1,dife)) &
	4271	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	4272
	4273
	4274	tia_spla=tia_spla+REAL(ti_spla)/REAL(clock_rate)
	4275	print *,'times for this timestep:timeproc,timeproc/time_pytracr_spl-'
	4276	print *,'time spla',REAL(ti_spla)/REAL(clock_rate) &
	4277	,REAL(ti_spla)/REAL(ti_spla)
	4278	print *,'time init',REAL(ti_init)/REAL(clock_rate) &
	4279	,REAL(ti_init)/REAL(ti_spla)
	4280	print *,'time inittype',REAL(ti_inittype)/REAL(clock_rate) &
	4281	,REAL(ti_inittype)/REAL(ti_spla)
	4282	print *,'time inittwrite',REAL(ti_inittwrite)/REAL(clock_rate) &
	4283	,REAL(ti_inittwrite)/REAL(ti_spla)
	4284	print *,'time emis',REAL(ti_emis)/REAL(clock_rate) &
	4285	,REAL(ti_emis)/REAL(ti_spla)
	4286	print *,'time depo ',REAL(ti_depo)/REAL(clock_rate) &
	4287	,REAL(ti_depo)/REAL(ti_spla)
	4288	print *,'time cltr',REAL(ti_cltr)/REAL(clock_rate) &
	4289	,REAL(ti_cltr)/REAL(ti_spla)
	4290	print *,'time ther',REAL(ti_ther)/REAL(clock_rate) &
	4291	,REAL(ti_ther)/REAL(ti_spla)
	4292	print *,'time sedi',REAL(ti_sedi)/REAL(clock_rate) &
	4293	,REAL(ti_sedi)/REAL(ti_spla)
	4294	print *,'time gas to part',REAL(ti_gasp)/REAL(clock_rate) &
	4295	,REAL(ti_gasp)/REAL(ti_spla)
	4296	print *,'time AP wet',REAL(ti_wetap)/REAL(clock_rate) &
	4297	,REAL(ti_wetap)/REAL(ti_spla)
	4298	print *,'time convective',REAL(ti_cvltr)/REAL(clock_rate) &
	4299	,REAL(ti_cvltr)/REAL(ti_spla)
	4300	print *,'time NP lsc scav',REAL(ti_lscs)/REAL(clock_rate) &
	4301	,REAL(ti_lscs)/REAL(ti_spla)
	4302	print *,'time opt,brdn,etc',REAL(ti_brop)/REAL(clock_rate) &
	4303	,REAL(ti_brop)/REAL(ti_spla)
	4304	print *,'time outputs',REAL(ti_outs)/REAL(clock_rate) &
	4305	,REAL(ti_outs)/REAL(ti_spla)
	4306
	4307
	4308	print *,'--time accumulated: time proc, time proc/time phytracr_spl--'
	4309	print *,'time spla',tia_spla
	4310	print *,'time init',tia_init,tia_init/tia_spla
	4311	print *,'time inittype',tia_inittype,tia_inittype/tia_spla
	4312	print *,'time inittwrite',tia_inittwrite,tia_inittwrite/tia_spla
	4313	print *,'time emis',tia_emis,tia_emis/tia_spla
	4314	print *,'time depo',tia_depo,tia_depo/tia_spla
	4315	print *,'time cltr',tia_cltr,tia_cltr/tia_spla
	4316	print *,'time ther',tia_ther,tia_ther/tia_spla
	4317	print *,'time sedi',tia_sedi,tia_sedi/tia_spla
	4318	print *,'time gas to part',tia_gasp,tia_gasp/tia_spla
	4319	print *,'time AP wet',tia_wetap,tia_wetap/tia_spla
	4320	print *,'time convective',tia_cvltr,tia_cvltr/tia_spla
	4321	print *,'time NP lsc scav',tia_lscs,tia_lscs/tia_spla
	4322	print *,'time opt,brdn,etc',tia_brop,tia_brop/tia_spla
	4323	print *,'time outputs',tia_outs,tia_outs/tia_spla
	4324
	4325
	4326
	4327	dife=clock_end_outphytracr-clock_start_outphytracr
	4328	ti_nophytracr=dife*MAX(0,SIGN(1,dife)) &
	4329	+(dife+clock_per_max)*MAX(0,SIGN(1,-dife))
	4330	tia_nophytracr=tia_nophytracr+REAL(ti_nophytracr)/REAL(clock_rate)
	4331	print *,'Time outside phytracr; Time accum outside phytracr'
	4332	print*,REAL(ti_nophytracr)/REAL(clock_rate),tia_nophytracr
	4333
	4334	clock_start_outphytracr=clock_end
	4335
	4336	ENDIF
	4337	print *,'END PHYTRACR_SPL '
	4338	print *,'lmt_so2ff_l FIN' , MINVAL(lmt_so2ff_l), MAXVAL(lmt_so2ff_l)
	4339
	4340	! CALL abort_gcm('TEST1', 'OK1', 1)
	4341
	4342	RETURN
	4343	END SUBROUTINE phytracr_spl
	4344
	4345	SUBROUTINE readregionsdims2_spl(nbreg,fileregions)
	4346
	4347	USE mod_grid_phy_lmdz
	4348	USE mod_phys_lmdz_para
	4349
	4350	IMPLICIT NONE
	4351	CHARACTER*800 fileregions
	4352	CHARACTER*800 auxstr
	4353	INTEGER nbreg
	4354
	4355	IF (is_mpi_root .AND. is_omp_root) THEN
	4356
	4357	OPEN (UNIT=1,FILE=trim(adjustl(fileregions)))
	4358	READ(1,'(a)') auxstr
	4359	READ(1,'(i10)') nbreg
	4360	CLOSE(UNIT=1)
	4361	ENDIF
	4362	CALL bcast(nbreg)
	4363
	4364	END SUBROUTINE readregionsdims2_spl
	4365
	4366	SUBROUTINE readregionsdims_spl(nbreg_ind,fileregionsdimsind, &
	4367	nbreg_dust,fileregionsdimsdust, &
	4368	nbreg_bb,fileregionsdimsbb)
	4369	USE mod_grid_phy_lmdz
	4370	USE mod_phys_lmdz_para
	4371
	4372	IMPLICIT NONE
	4373	CHARACTER*800 fileregionsdimsind
	4374	CHARACTER*800 fileregionsdimsdust
	4375	CHARACTER*800 fileregionsdimsbb
	4376	CHARACTER*800 auxstr
	4377	INTEGER nbreg_ind,nbreg_dust,nbreg_bb
	4378
	4379	IF (is_mpi_root .AND. is_omp_root) THEN
	4380
	4381	OPEN (UNIT=1,FILE=trim(adjustl(fileregionsdimsind)))
	4382	READ(1,'(a)') auxstr
	4383	READ(1,'(i10)') nbreg_ind
	4384	CLOSE(UNIT=1)
	4385
	4386	OPEN (UNIT=1,FILE=trim(adjustl(fileregionsdimsdust)))
	4387	READ(1,'(a)') auxstr
	4388	READ(1,'(i10)') nbreg_dust
	4389	CLOSE(UNIT=1)
	4390
	4391	OPEN (UNIT=1,FILE=trim(adjustl(fileregionsdimsbb)))
	4392	READ(1,'(a)') auxstr
	4393	READ(1,'(i10)') nbreg_bb
	4394	CLOSE(UNIT=1)
	4395
	4396
	4397	ENDIF
	4398	CALL bcast(nbreg_ind)
	4399	CALL bcast(nbreg_dust)
	4400	CALL bcast(nbreg_bb)
	4401
	4402	END SUBROUTINE readregionsdims_spl
	4403
	4404	SUBROUTINE readregions_spl(iregion,filenameregion)
	4405	USE dimphy
	4406	USE mod_grid_phy_lmdz
	4407	USE mod_phys_lmdz_para
	4408
	4409	IMPLICIT NONE
	4410	CHARACTER() filenameregion
	4411	INTEGER iregion(klon)
	4412	INTEGER iregion_glo(klon_glo)
	4413	INTEGER k
	4414
	4415	IF (is_mpi_root .AND. is_omp_root) THEN
	4416
	4417	print *,trim(adjustl(filenameregion))
	4418	OPEN(1,file=trim(adjustl(filenameregion)))
	4419	DO k=1,klon_glo
	4420	READ(1,'(i10)') iregion_glo(k)
	4421	ENDDO
	4422	CLOSE(UNIT=1)
	4423	ENDIF
	4424	CALL scatter(iregion_glo,iregion)
	4425
	4426	END SUBROUTINE readregions_spl
	4427
[3806]	4428	!! AS: SUBROUTINE readscaleparams_spl pas appellee
[2630]	4429	SUBROUTINE readscaleparams_spl(scale_param, nbreg, &
	4430	filescaleparams)
	4431	USE mod_grid_phy_lmdz
	4432	USE mod_phys_lmdz_para
	4433	IMPLICIT NONE
	4434
	4435	CHARACTER*800 filescaleparams
	4436	INTEGER nbreg
	4437	REAL scale_param(nbreg)
	4438	INTEGER k
	4439
	4440	IF (is_mpi_root .AND. is_omp_root) THEN
	4441	OPEN(1,file=trim(adjustl(filescaleparams)),form='unformatted')
	4442	do k=1,nbreg
	4443	read(1) scale_param(k)
	4444	enddo
	4445	CLOSE(1)
	4446	ENDIF
	4447	CALL bcast(scale_param)
	4448	! print *,'holaaaaaaaaaaaa'
	4449	! print *,scale_param
	4450
	4451	END SUBROUTINE readscaleparams_spl
	4452
	4453	SUBROUTINE readscaleparamsnc_spl(scale_param_ind, &
	4454	nbreg_ind, paramname_ind, &
	4455	scale_param_ff, nbreg_ff,paramname_ff, &
	4456	scale_param_bb, nbreg_bb,paramname_bb, &
	4457	scale_param_dustacc, nbreg_dustacc,paramname_dustacc, &
	4458	scale_param_dustcoa, nbreg_dustcoa,paramname_dustcoa, &
	4459	scale_param_dustsco, nbreg_dustsco,paramname_dustsco, &
	4460	param_wstarBLperregion, nbreg_wstardustBL, paramname_wstarBL, &
	4461	param_wstarWAKEperregion, nbreg_wstardustWAKE, paramname_wstarWAKE, &
	4462	scale_param_ssacc , paramname_ssacc, &
	4463	scale_param_sscoa , paramname_sscoa, &
	4464	filescaleparams,julien,jH_phys, pdtphys,debutphy)
	4465	! SUBROUTINE readscaleparamsnc_spl(scale_param, nbreg, &
	4466	! filescaleparams,paramname,&
	4467	! julien,jH_phys, pdtphys,debutphy)
	4468	USE mod_grid_phy_lmdz
	4469	USE mod_phys_lmdz_para
	4470	IMPLICIT NONE
	4471
	4472	CHARACTER*800 filescaleparams
	4473	CHARACTER*100 paramname_ind,paramname_ff,paramname_bb
	4474	CHARACTER*100 paramname_dustacc, paramname_dustcoa
	4475	CHARACTER*100 paramname_dustsco
	4476	CHARACTER*100 paramname_ssacc
	4477	CHARACTER*100 paramname_sscoa
	4478	CHARACTER*100 paramname_wstarBL
	4479	CHARACTER*100 paramname_wstarWAKE
	4480
	4481	INTEGER nbreg,iday
	4482	INTEGER nbreg_ind, nbreg_ff, nbreg_bb , nbreg_dustacc
	4483	INTEGER nbreg_dustcoa , nbreg_dustsco, nbreg_wstardustBL
	4484	INTEGER nbreg_wstardustWAKE
	4485	INTEGER,PARAMETER :: nbreg_ssacc=1
	4486	INTEGER,PARAMETER :: nbreg_sscoa=1
	4487	REAL,PARAMETER :: sca_resol = 24. ! resolution of scalig params in hours
	4488	REAL scale_param_ind(nbreg_ind)
	4489	REAL scale_param_bb(nbreg_bb)
	4490	REAL scale_param_ff(nbreg_ff)
	4491	REAL scale_param_dustacc(nbreg_dustacc)
	4492	REAL scale_param_dustcoa(nbreg_dustcoa)
	4493	REAL scale_param_dustsco(nbreg_dustsco)
	4494	REAL param_wstarBLperregion(nbreg_wstardustBL)
	4495	REAL param_wstarWAKEperregion(nbreg_wstardustWAKE)
	4496	REAL scale_param_ssacc
	4497	REAL scale_param_ssacc_tmp(nbreg_ssacc)
	4498	REAL scale_param_sscoa
	4499	REAL scale_param_sscoa_tmp(nbreg_sscoa)
	4500
	4501	INTEGER k,step_sca,test_sca
	4502	REAL :: jH_phys, pdtphys
	4503	REAL,SAVE :: jH_sca, jH_ini
	4504	INTEGER julien
	4505	LOGICAL debutphy
	4506	SAVE step_sca,test_sca,iday
	4507	!$OMP THREADPRIVATE(step_sca,test_sca,iday)
	4508	!$OMP THREADPRIVATE(jH_sca,jH_ini)
	4509
	4510	IF (debutphy) THEN
	4511	iday=julien
	4512	step_sca=1
	4513	test_sca=0
	4514	jH_ini=jH_phys
	4515	jH_sca=jH_phys
	4516	ENDIF
	4517
	4518	IF (test_sca .EQ. 0 ) THEN
	4519	! READ file!!
	4520	call read_scalenc(filescaleparams,paramname_ind, &
	4521	nbreg_ind,step_sca, &
	4522	scale_param_ind)
	4523	call read_scalenc(filescaleparams,paramname_bb, &
	4524	nbreg_bb,step_sca, &
	4525	scale_param_bb)
	4526	call read_scalenc(filescaleparams,paramname_ff, &
	4527	nbreg_ff,step_sca, &
	4528	scale_param_ff)
	4529	call read_scalenc(filescaleparams,paramname_dustacc, &
	4530	nbreg_dustacc,step_sca, &
	4531	scale_param_dustacc)
	4532	call read_scalenc(filescaleparams,paramname_dustcoa, &
	4533	nbreg_dustcoa,step_sca, &
	4534	scale_param_dustcoa)
	4535	call read_scalenc(filescaleparams,paramname_dustsco, &
	4536	nbreg_dustsco,step_sca, &
	4537	scale_param_dustsco)
	4538	call read_scalenc(filescaleparams,paramname_wstarBL, &
	4539	nbreg_wstardustBL,step_sca, &
	4540	param_wstarBLperregion)
	4541	call read_scalenc(filescaleparams,paramname_wstarWAKE, &
	4542	nbreg_wstardustWAKE,step_sca, &
	4543	param_wstarWAKEperregion)
	4544	call read_scalenc(filescaleparams,paramname_ssacc, &
	4545	nbreg_ssacc,step_sca, &
	4546	scale_param_ssacc_tmp)
	4547	call read_scalenc(filescaleparams,paramname_sscoa, &
	4548	nbreg_sscoa,step_sca, &
	4549	scale_param_sscoa_tmp)
	4550	scale_param_ssacc=scale_param_ssacc_tmp(1)
	4551	scale_param_sscoa=scale_param_sscoa_tmp(1)
	4552
	4553	!print *,'JEREADFILE',julien,jH_phys
	4554	step_sca= step_sca + 1
	4555	test_sca=1
	4556	ENDIF
	4557
	4558	jH_sca=jH_sca+pdtphys/(24.*3600.)
	4559	IF (jH_sca.GT.(sca_resol)/24.) THEN
	4560	test_sca=0
	4561	jH_sca=jH_ini
	4562	ENDIF
	4563
	4564	END SUBROUTINE readscaleparamsnc_spl
	4565
	4566	SUBROUTINE read_scalenc(filescaleparams,paramname,nbreg,step_sca, &
	4567	scale_param)
	4568
	4569	USE mod_grid_phy_lmdz
	4570	USE mod_phys_lmdz_para
	4571	IMPLICIT NONE
	4572
	4573	include "netcdf.inc"
	4574
	4575	CHARACTER*800 filescaleparams
	4576	CHARACTER*100 paramname
	4577	INTEGER nbreg, step_sca
	4578	REAL scale_param(nbreg)
	4579	!local vars
	4580	integer nid,ierr,nvarid
	4581	real rcode,auxreal
	4582	integer start(4),count(4), status
	4583	! local
	4584	integer debutread,countread
	4585	CHARACTER*104 varname
	4586	CHARACTER*2 aux_2s
	4587	integer i, j, ig
	4588	!$OMP MASTER
	4589	IF (is_mpi_root .AND. is_omp_root) THEN
	4590	!nci=NCOPN(trim(adjustl(filescaleparams)),NCNOWRIT,rcode)
	4591	ierr = NF_OPEN (trim(adjustl(filescaleparams)),NF_NOWRITE, nid)
	4592	if (ierr .EQ. NF_NOERR) THEN
	4593	debutread=step_sca
	4594	countread=1
	4595
	4596	do i=1,nbreg
	4597	WRITE(aux_2s,'(i2.2)') i
	4598	varname= trim(adjustl(paramname))//aux_2s
	4599	print *,varname
	4600	ierr = NF_INQ_VARID (nid,trim(adjustl(varname)), nvarid)
	4601	ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debutread, &
	4602	countread, auxreal)
	4603	IF (ierr .NE. NF_NOERR) THEN
	4604	PRINT*, 'Pb de lecture pour modvalues'
	4605	print *,'JE scale_var, step_sca',trim(adjustl(varname)),step_sca
	4606	CALL HANDLE_ERR(ierr)
	4607	print *,'error ierr= ',ierr
	4608	CALL exit(1)
	4609	call abort_gcm('read_scalenc','error reading variable',1)
	4610	ENDIF
	4611
	4612	print *,auxreal
	4613	scale_param(i)=auxreal
	4614	enddo
	4615
	4616	ierr = NF_CLOSE(nid)
	4617	else
	4618	print *,'File '//trim(adjustl(filescaleparams))//' not found'
	4619	print *,'doing nothing...'
	4620	endif
	4621
	4622	ENDIF ! mpi_root
	4623	!$OMP END MASTER
	4624	!$OMP BARRIER
	4625	! CALL scatter(var local _glo,var local) o algo asi
	4626	call bcast(scale_param)
	4627	END SUBROUTINE read_scalenc
	4628
	4629
	4630
	4631	END MODULE

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

Download in other formats:

Original Format