Index: LMDZ6/trunk/libf/phylmd/Dust/dustemission_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/Dust/dustemission_mod.F90	(revision 3805)
+++ LMDZ6/trunk/libf/phylmd/Dust/dustemission_mod.F90	(revision 3806)
@@ -636,8 +636,4 @@
     ALLOCATE( srel(nats,nclass) )
     ALLOCATE( srel2(nats,nclass) )
-    !!AS: moved to subroutine dustemission_output_init
-    !ALLOCATE( m1dflux(klon) )
-    !ALLOCATE( m2dflux(klon) )
-    !ALLOCATE( m3dflux(klon) )
 
 
@@ -684,5 +680,4 @@
       !&      .and.pctsrf(i)>0.5.and.Pini(i,nts)>0.)THEN
   ! JE20150605<< easier to read
-  !    IF(pctsrf(i)>0.5.and.Pini(i,nts)>0.)THEN
       IF(pctsrf(i,is_ter)>0.5.and.Pini(i,nts)>0.)THEN
   ! JE20150605>>
@@ -1275,17 +1270,5 @@
                       modwm=sqrt((wind10ms(kwb)**2)+(1.2*zwstar(i))**2)
                       ustarns=cdnms*modwm*100.
-!JE20150202 <<
-! Do not have too much sense.. and is not anymore in the chimere14b version.
-!
-!                      utmin=umin/(cdnms*ceff)
-!                   IF(wind10cm(kwb).ge.utmin)THEN
-!                      ustarsalt=ustarns+  &
-!                    (0.3*(wind10cm(kwb)/100.-utmin/100.)**2.)
-!                   ELSE
-!                      ustarsalt=ustarns
-!                   ENDIF
-! ustarsalt should be :
                     ustarsalt=ustarns
-!JE20150202 >>
 
 
Index: LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90	(revision 3805)
+++ LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90	(revision 3806)
@@ -34,7 +34,4 @@
   CHARACTER*800 fileregionsdimsbb
   CHARACTER*800 fileregionsdimswstar
-!  CHARACTER*800 filescaleparamsind
-!  CHARACTER*800 filescaleparamsdust
-!  CHARACTER*800 filescaleparamsbb
   CHARACTER*100 paramname_ind
   CHARACTER*100 paramname_bb
@@ -423,21 +420,9 @@
 !$OMP THREADPRIVATE(flux_sparam_sscoa,u10m_ss,v10m_ss)
 
-! Select dust emission scheme ver the Sahara:
+! Select dust emission scheme for the Sahara:
 !      LOGICAL,PARAMETER,SAVE ::  ok_chimeredust=.FALSE.
       LOGICAL,PARAMETER ::  ok_chimeredust=.TRUE.
 !!!!!! !$OMP THREADPRIVATE(ok_chimeredust)
 
-!OH   REAL,SAVE :: scale_param_ssacc  !Scaling parameter for Fine Sea Salt
-!OH   REAL,SAVE :: scale_param_sscoa  !Scaling parameter for Coarse Sea Salt
-!OH   REAL,ALLOCATABLE,SAVE :: scale_param_ind(nbreg_ind) !Scaling parameter for industrial emissionsi of SO2
-!OH   REAL,ALLOCATABLE,SAVE :: scale_param_bb(nbreg_bb)  !Scaling parameter for biomas burning (SO2, BC & OM)
-!OH   REAL,ALLOCATABLE,SAVE :: scale_param_ff(nbreg_ind)  !Scaling parameter for industrial emissions (fossil fuel)
-!OH   REAL,ALLOCATABLE,SAVE :: scale_param_dustacc(nbreg_dust)  !Scaling parameter for Fine Dust
-!OH   REAL,ALLOCATABLE,SAVE :: scale_param_dustcoa(nbreg_dust)  !Scaling parameter for Coarse Dust
-!OH   REAL,ALLOCATABLE,SAVE :: scale_param_dustsco(nbreg_dust)  !Scaling parameter for SCoarse Dust
-!OH   REAL,ALLOCATABLE,SAVE :: param_wstarBLperregion(nbreg_wstardust)
-!OH   REAL,ALLOCATABLE,SAVE :: param_wstarWAKEperregion(nbreg_wstardust)
-!!!! !$OMP THREADPRIVATE( scale_param_ssacc, scale_param_sscoa, scale_param_ind, scale_param_bb, scale_param_ff, scale_param_dustacc, scale_param_dustcoa, scale_param_dustsco, param_wstarBLperregion, param_wstarWAKEperregion)
-
 
 CONTAINS
@@ -446,6 +431,6 @@
 SUBROUTINE phytracr_spl_out_init()
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-!AS : This new subroutine centralises les ALLOCATE dont on a besoin au premier appel 
-!     de phys_output_write_spl dans physiq
+!AS : This subroutine centralises the ALLOCATE needed for the 1st call of
+!     phys_output_write_spl in physiq
 
     USE dimphy
@@ -455,7 +440,4 @@
     ! pour les variables m[1-3]dflux
     CALL dustemis_out_init() 
-
-    !AS: moved here from  section INITIALISATIONS of subroutine phytracr_spl_ini
-    !    ALLOCATE variables from spla_output_write.h, section "traceur"
 
 !traceur
@@ -686,5 +668,5 @@
       ALLOCATE( v10m_ss(klon))
 
-!AS: dans phys_output_write_spl, hors spla_output_write.h
+!AS: in phys_output_write_spl, but not in spla_output_write.h
 !------------------------------------------------------
       ALLOCATE(d_tr_cl(klon,klev,nbtr))
@@ -876,6 +858,4 @@
       INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: masque_terra !
 !$OMP THREADPRIVATE(masque_aqua,masque_terra)
-!!$OMP THREADPRIVATE(aod550_aqua,aod550_terra,aod670_aqua,aod670_terra)
-!!$OMP THREADPRIVATE(aod865_aqua,aod865_terra)
 
   INTEGER, SAVE :: nbreg_dust, nbreg_ind, nbreg_bb, nbreg_ss,nbreg_wstardust
@@ -1060,6 +1040,4 @@
       REAL, intent(in) ::  pmflxr(klon,klev+1), pmflxs(klon,klev+1)   !--convection
       REAL, intent(in) ::  prfl(klon,klev+1),   psfl(klon,klev+1)     !--large-scale
-! JE      REAL pmflxr(klon,klev), pmflxs(klon,klev)   !--convection       ! Titane
-! JE      REAL prfl(klon,klev),   psfl(klon,klev)     !--large-scale      ! Titane
       REAL :: ql_incl ! contenu en eau liquide nuageuse dans le nuage ! ql_incl=oliq/rneb
       REAL  :: ql_incloud_ref    ! ref value of in-cloud condensed water content
@@ -1131,7 +1109,4 @@
       REAL flux_tr(klon,nbtr)
       REAL m_conc(klon,klev)
-!      REAL sed_ss(klon)    ! corresponds to tracer 3
-!      REAL sed_dust(klon)  ! corresponds to tracer 4
-!      REAL sed_dustsco(klon)  ! corresponds to tracer 4
       REAL henry(nbtr)  !--cste de Henry  mol/l/atm
       REAL kk(nbtr)     !--coefficient de var avec T (K)
@@ -1157,9 +1132,7 @@
    
       CHARACTER*2 str2
-      LOGICAL ok_histrac
-!JE2014124      PARAMETER (ok_histrac=.true.)
-      PARAMETER (ok_histrac=.false.)
-!      PARAMETER (ok_chimeredust=.false.) 
-!      PARAMETER (ok_chimeredust=.true.) 
+!!AS:      LOGICAL ok_histrac
+!!!JE2014124      PARAMETER (ok_histrac=.true.)
+!!      PARAMETER (ok_histrac=.false.)
       INTEGER ndex2d(iim*(jjm+1)), ndex3d(iim*(jjm+1)*klev)
       INTEGER nhori1, nhori2, nhori3, nhori4, nhori5, nvert
@@ -1176,5 +1149,4 @@
       REAL zx_tmp_2d(iim,jjm+1), zx_tmp_3d(iim,jjm+1,klev)
       REAL zx_tmp_fi2d(klon), zx_tmp_fi3d(klon, klev)
-!      REAL zx_lon(iim,jjm+1), zx_lat(iim,jjm+1)
       REAL zx_lon_glo(nbp_lon,nbp_lat), zx_lat_glo(nbp_lon,nbp_lat)
       REAL zsto, zout, zout_h, zout_m, zjulian
@@ -1208,7 +1180,4 @@
 ! JE for updating in  cltrac
       REAL,DIMENSION(klon,klev)             :: delp     ! epaisseur de couche (Pa)
-!JE20140507      REAL,DIMENSION(klon,nbtr)       :: d_tr_dry ! Td depot sec/traceur (1st layer),ALLOCATABLE,SAVE  jyg
-!JE20140507      REAL,DIMENSION(klon,nbtr)        ::  flux_tr_dry
-!      SAVE  d_tr_dry
 !! JE for include gas to particle conversion in output
 !      REAL his_g2pgas(klon)      ! gastoparticle in gas units (check!)
@@ -1617,512 +1586,12 @@
          print *,'ecrit_tra=', pdtphys, ecrit_tra
 
-         IF (ok_histrac) THEN
-           IF (is_mpi_root .AND. is_omp_root) THEN
-  
-           itra=0
-!
-           CALL ymds2ju(1900, 1, 1, 0.0, zjulian)
-!
-!           print *, 'klon,iim,jjm+1 = ',klon,iim,jjm+1
-           print *, 'glo klon,iim,jjm+1 = ',klon_glo,nbp_lon,nbp_lat
-           CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,rlon,zx_lon_glo)
-!
-!           DO i = 1, iim
-           DO i = 1, nbp_lon
-             zx_lon_glo(i,1) = rlon(i+1)
-             zx_lon_glo(i,nbp_lat) = rlon(i+1)
-           ENDDO
-!
-      CALL histbeg("histrac_spl", nbp_lon,zx_lon_glo,            &
-                       nbp_lat,zx_lat_glo,                       &
-                       1,nbp_lon,1,nbp_lat, 0, zjulian, pdtphys, &
-                       nhori1, nid_tra1)
-!
-      CALL histbeg("lessivage_spl", nbp_lon,zx_lon_glo,            &
-                       nbp_lat,zx_lat_glo,                         &
-                       1,nbp_lon,1,nbp_lat, 0, zjulian, pdtphys,   &
-                       nhori2, nid_tra2)
-! 
-      CALL histbeg("traceur_spl", nbp_lon,zx_lon_glo,               &
-                       nbp_lat,zx_lat_glo,                         &
-                      1,nbp_lon,1,nbp_lat, 0, zjulian, pdtphys,    &
-                       nhori3, nid_tra3)
-!
-      CALL histvert(nid_tra1, "presnivs", "Vertical levels", "mb",  &
-                      nbp_lev, presnivs, nvert)
-!
-      CALL histvert(nid_tra2, "presnivs", "Vertical levels", "mb",  &
-                      nbp_lev, presnivs, nvert)
-!
-      CALL histvert(nid_tra3, "presnivs", "Vertical levels", "mb",  &
-                      nbp_lev, presnivs, nvert)
-!
-           zsto = pdtphys
-           zout = pdtphys * FLOAT(ecrit_tra)
-           zout_h = pdtphys * FLOAT(ecrit_tra_h)
-           zout_m = pdtphys * FLOAT(ecrit_tra_m)
-           print *,'zsto zout=', zsto, zout
-
-!
-!----------------- HISTORY FILES OF TRACER EMISSIONS -------------------
-!
-! HISTRAC
-!
-       CALL histdef(nid_tra1, "fluxbb", "Flux BB", "mg/m2/s",       &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,     &
-                       "ave(X)", zsto,zout)                         
-!                                                                   
-      CALL histdef(nid_tra1, "fluxff", "Flux FF", "mg/m2/s",        & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,      &
-                       "ave(X)", zsto,zout)                          
-!                                                                    
-      CALL histdef(nid_tra1, "fluxbcbb", "Flux BC-BB", "mg/m2/s",    &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,       &
-                       "ave(X)", zsto,zout)                           
-!                                                                     
-      CALL histdef(nid_tra1, "fluxbcff", "Flux BC-FF", "mg/m2/s",     &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,        &
-                       "ave(X)", zsto,zout)                            
-!                                                                      
-      CALL histdef(nid_tra1, "fluxbcnff", "Flux BC-NFF", "mg/m2/s",    &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                             
-!                                                                       
-      CALL histdef(nid_tra1, "fluxbcba", "Flux BC-BA", "mg/m2/s",       &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,          &
-                       "ave(X)", zsto,zout)                              
-!                                                                        
-      CALL histdef(nid_tra1, "fluxbc", "Flux BC", "mg/m2/s",    &         
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,  &
-                       "ave(X)", zsto,zout)                      
-!                                                                
-      CALL histdef(nid_tra1, "fluxombb", "Flux OM-BB", "mg/m2/s" ,  &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,      &
-                       "ave(X)", zsto,zout)                          
-!                                                                    
-      CALL histdef(nid_tra1, "fluxomff", "Flux OM-FF", "mg/m2/s",    &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,       &
-                       "ave(X)", zsto,zout)                           
-!                                                                     
-      CALL histdef(nid_tra1, "fluxomnff", "Flux OM-NFF", "mg/m2/s",  & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,       &
-                       "ave(X)", zsto,zout)                           
-!                                                                     
-      CALL histdef(nid_tra1, "fluxomba", "Flux OM-BA", "mg/m2/s",    & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,       &
-                       "ave(X)", zsto,zout)                           
-!                                                                     
-      CALL histdef(nid_tra1, "fluxomnat", "Flux OM-NT", "mg/m2/s",   & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,       &
-                       "ave(X)", zsto,zout)                           
-!                                                                     
-      CALL histdef(nid_tra1, "fluxom", "Flux OM", "mg/m2/s",         & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,       &
-                       "ave(X)", zsto,zout)                           
-!                                                                     
-      CALL histdef(nid_tra1,"fluxh2sff","Flux H2S FF","mgS/m2/s",    & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,       &
-                       "ave(X)", zsto,zout)                           
-!                                                                     
-      CALL histdef(nid_tra1,"fluxh2snff","Flux H2S non-FF",          & 
-                       "mgS/m2/s",nbp_lon,nbp_lat,nhori1, 1,1,1,     &
-                        -99, 32,                                     & 
-                       "ave(X)", zsto,zout)                           
-!                                                                     
-      CALL histdef(nid_tra1,"fluxso2ff","Flux SO2 FF","mgS/m2/s",    &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,       &
-                       "ave(X)", zsto,zout)                           
-!                                                                     
-      CALL histdef(nid_tra1,"fluxso2nff","Flux SO2 non-FF",          & 
-                       "mgS/m2/s",nbp_lon,nbp_lat,nhori1, 1,1,1,     &
-                        -99, 32,                                     &
-                       "ave(X)", zsto,zout)                           
-!                                                                      
-      CALL histdef(nid_tra1, "fluxso2bb", "Flux SO2 BB","mgS/m2/s",   & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,        &  
-                       "ave(X)", zsto,zout)                           
-!                                                                      
-      CALL histdef(nid_tra1,"fluxso2vol","Flux SO2 Vol","mgS/m2/s",    &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1, "fluxso2ba", "Flux SO2 Ba","mgS/m2/s",    & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1, "fluxso2", "Flux SO2","mgS/m2/s",         & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1,"fluxso4ff","Flux SO4 FF","mgS/m2/s",      & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1,"fluxso4nff","Flux SO4 non-FF",            & 
-                   "mgS/m2/s", nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32, &
-                   "ave(X)", zsto,zout)                                
-!                                                                       
-      CALL histdef(nid_tra1, "fluxso4bb", "Flux SO4 BB","mgS/m2/s",    & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1, "fluxso4ba", "Flux SO4 Ba","mgS/m2/s",    & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1, "fluxso4", "Flux SO4","mgS/m2/s",         & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1, "fluxdms", "Flux DMS", "mgS/m2/s",        & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1,"fluxh2sbio","Flux H2S Bio","mgS/m2/s",    & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1, "fluxdustec",                             & 
-                                      "Flux Dust EC", "mg/m2/s",       &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                             
-!                                                                       
-      CALL histdef(nid_tra1,"fluxddfine","DD Fine Mode","mg/m2/s",     &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1,"fluxddcoa","DD Coarse Mode","mg/m2/s",    & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1,"fluxddsco","DD SCoarse Mode","mg/m2/s",   & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1,"fluxdd","Flux DD","mg/m2/s",              & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1,"fluxssfine","SS Fine Mode","mg/m2/s",     & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1,"fluxsscoa","SS Coarse Mode","mg/m2/s",    & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-      CALL histdef(nid_tra1,"fluxss","Flux SS","mg/m2/s",              & 
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                            
-!                                                                       
-!nhl          CALL histdef(nid_tra1,"fluxso4chem","SO4 chem prod",      
-!nhl    .                  "gAer/kgAir",
-!nhl    .                  nbp_lon,nbp_lat,nhori1, nbp_lev,1,nbp_lev,nvert, 32,
-!nhl    .                  "ave(X)", zsto,zout)
-!
-          CALL histdef(nid_tra1,"flux_sparam_ind","Ind emiss",      &
-                       "mg/m2/s",                                   &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,      & 
-                       "ave(X)", zsto,zout)                           
-!                                                                    
-          CALL histdef(nid_tra1,"flux_sparam_bb","BB emiss",        &
-                       "mg/m2/s",                                   &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,      &
-                       "ave(X)", zsto,zout)                          
-!                                                                    
-          CALL histdef(nid_tra1,"flux_sparam_ff","FF emiss",        &
-                       "mg/m2/s",                                   &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,      &
-                       "ave(X)", zsto,zout)                          
-!                                                                    
-          CALL histdef(nid_tra1,"flux_sparam_ddfine","DD fine emiss",  &
-                       "mg/m2/s",                                      &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                             
-!                                                                       
-          CALL histdef(nid_tra1,"flux_sparam_ddcoa","DD coarse emiss",  & 
-                       "mg/m2/s",                                       &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,          &
-                       "ave(X)", zsto,zout)                              
-!                                                                        
-          CALL histdef(nid_tra1,"flux_sparam_ddsco","DD Scoarse emiss", &
-                       "mg/m2/s",                                       &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,          &
-                       "ave(X)", zsto,zout)                              
-!                                                                        
-          CALL histdef(nid_tra1,"flux_sparam_ssfine","SS fine emiss",   &
-                       "mg/m2/s",                                       &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,          &
-                       "ave(X)", zsto,zout)                              
-!                                                                        
-          CALL histdef(nid_tra1,"flux_sparam_sscoa","SS coarse emiss",  &
-                       "mg/m2/s",                                       &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,          &
-                       "ave(X)", zsto,zout)                              
-!                                                                        
-          CALL histdef(nid_tra1,"u10m","Zonal wind at 10 m",            &
-                       "m/s",                                           &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,          &
-                       "ave(X)", zsto,zout)                              
-!                                                                        
-          CALL histdef(nid_tra1,"v10m","Meridional wind at 10 m",       &
-                       "m/s",                                           &
-                       nbp_lon,nbp_lat,nhori1, 1,1,1, -99, 32,          &
-                       "ave(X)", zsto,zout)                              
-!                                                                        
-!nhl          CALL histdef(nid_tra1,"flux_sparam_sulf","SO4 chem prod", 
-!nhl    .                  "gAer/kgAir",
-!nhl    .                  nbp_lon,nbp_lat,nhori1, nbp_lev,1,nbp_lev,nvert, 32,
-!nhl    .                  "ave(X)", zsto,zout)
-!
-! TRACEUR
-!
-          CALL histdef(nid_tra3, "taue550", "Tau ext 550", " ",           &
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,            &
-                       "ave(X)", zsto,zout)                                  
-!                                                                           
-          CALL histdef(nid_tra3, "taue670", "Tau ext 670", " ",            &  
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,             &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-          CALL histdef(nid_tra3, "taue865", "Tau ext 865", " ",            & 
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,             &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-          CALL histdef(nid_tra3, "taue550_tr2", "Tau ext 550tr2", " ",     & 
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,             &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-          CALL histdef(nid_tra3, "taue670_tr2", "Tau ext 670tr2", " ",     & 
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,             &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-          CALL histdef(nid_tra3, "taue865_tr2", "Tau ext 865tr2", " ",     & 
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,             &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-          CALL histdef(nid_tra3, "taue550_ss", "Tau ext 550ss", " ",       & 
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,             &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-          CALL histdef(nid_tra3, "taue670_ss", "Tau ext 670ss", " ",       & 
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,             &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-          CALL histdef(nid_tra3, "taue865_ss", "Tau ext 865ss", " ",       & 
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,             &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-          CALL histdef(nid_tra3, "taue550_dust", "Tau ext 550dust", " "    & 
-                       ,nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,            &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-          CALL histdef(nid_tra3, "taue670_dust", "Tau ext 670dust", " "    & 
-                       ,nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,            &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-          CALL histdef(nid_tra3, "taue865_dust", "Tau ext 865dust", " "    & 
-                       ,nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,            &
-                       "ave(X)", zsto,zout)                                 
-                                                                            
-          CALL histdef(nid_tra3, "taue550_dustsco",                     &   
-                       "Tau ext 550dustsco", " "                        &
-                       ,nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                              
-!                                                                        
-           CALL histdef(nid_tra3, "taue670_dustsco",                    &
-                       "Tau ext 670dustsco", " "                        &
-                       ,nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                              
-!                                                                        
-           CALL histdef(nid_tra3, "taue865_dustsco",                    &
-                       "Tau ext 865dustsco", " "                        &
-                       ,nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,         &
-                       "ave(X)", zsto,zout)                              
-                                                                         
-                                                                        
-        CALL histdef(nid_tra3, "taue550_aqua", "Tau ext 550 aqua", " ",   & 
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,            &
-                       "inst(X)", zout,zout)                               
-      CALL histdef(nid_tra3, "taue550_terra", "Tau ext 550 terra", " ",   & 
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,            & 
-                       "inst(X)", zout,zout)                               
-        CALL histdef(nid_tra3, "taue670_aqua", "Tau ext 670 aqua", " ",   & 
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,            & 
-                       "inst(X)", zout,zout)                               
-      CALL histdef(nid_tra3, "taue670_terra", "Tau ext 670 terra", " ",   & 
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,            & 
-                       "inst(X)", zout,zout)                               
-        CALL histdef(nid_tra3, "taue865_aqua", "Tau ext 865 aqua", " ",   & 
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,            & 
-                       "inst(X)", zout,zout)                               
-      CALL histdef(nid_tra3, "taue865_terra", "Tau ext 865 terra", " ",   & 
-                       nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,            & 
-                       "inst(X)", zout,zout)                               
-                                                                           
-                                                                           
-          DO it=1, nbtr
-!
-          WRITE(str2,'(i2.2)') it
-!
-          CALL histdef(nid_tra3, "trm"//str2, "Burden No."//str2,         & 
-                     "mgS/m2", nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,    &
-                       "ave(X)", zsto,zout)                                
-!                                                                          
-          CALL histdef(nid_tra3, "sconc"//str2, "Surf Conc. No."//str2,   & 
-                       "mg/m3", nbp_lon,nbp_lat,nhori3, 1,1,1, -99, 32,   &
-                       "ave(X)", zsto,zout)                                
-!                                                                          
-! LESSIVAGE                                                                 
-!
-          CALL histdef(nid_tra2, "flux"//str2, "emission"//str2,           &  
-                   "mgS/m2/s", nbp_lon,nbp_lat,nhori2, 1,1,1, -99, 32,     &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-          CALL histdef(nid_tra2, "ds"//str2, "Depot sec No."//str2,        & 
-                   "mgS/m2/s", nbp_lon,nbp_lat,nhori2, 1,1,1, -99, 32,     &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-          CALL histdef(nid_tra2,"dh"//str2,                                 &
-                    "Depot hum No."//str2,                                 &
-                   "mgS/m2/s", nbp_lon,nbp_lat,nhori2, 1,1,1, -99, 32,     &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-          CALL histdef(nid_tra2,"dtrconv"//str2,                           &
-                     "Tiedke convective"//str2,                            &
-                  "mgS/m2/s", nbp_lon,nbp_lat,nhori2, 1,1,1, -99, 32,      &
-                       "ave(X)", zsto,zout)                                 
-                                                                            
-          CALL histdef(nid_tra2,"dtherm"//str2,                            &
-                       "Thermals dtracer"//str2,                           &
-                   "mgS/m2/s", nbp_lon,nbp_lat,nhori2, 1,1,1, -99, 32,     &
-                       "ave(X)", zsto,zout)                                 
-                                                                            
-          CALL histdef(nid_tra2,"dhkecv"//str2,                            &
-                       "KE dep hum convective"//str2,                      &
-                   "mgS/m2/s", nbp_lon,nbp_lat,nhori2, 1,1,1, -99, 32,     &
-                       "ave(X)", zsto,zout)                                 
-          CALL histdef(nid_tra2,"dhkelsc"//str2,                            &
-                       "KE dep hum large scale"//str2,                      &
-                   "mgS/m2/s", nbp_lon,nbp_lat,nhori2, 1,1,1, -99, 32,      &
-                       "ave(X)", zsto,zout)                                  
-                                                                             
-                                
-          CALL histdef(nid_tra2,"d_tr_ds"//str2,                            &
-                       " Tendance dep sec"//str2,                      &
-                   "mgS/m2/s", nbp_lon,nbp_lat,nhori2, 1,1,1, -99, 32,     &
-                       "ave(X)", zsto,zout)                                 
-
-                                            
-          CALL histdef(nid_tra2,"d_tr_cv"//str2,                          & 
-                       "cvltr d_tr_cv"//str2,                             &
-                       "mgS/m2/s",                                        &
-                   nbp_lon,nbp_lat,nhori2, nbp_lev,1,nbp_lev,nvert, 32,   &
-                       "ave(X)", zsto,zout)                                 
-          CALL histdef(nid_tra2,"d_tr_trsp"//str2                         & 
-                       ,"cvltr d_tr_trsp"//str2,                          &
-                       "mgS/m2/s",                                        &
-                   nbp_lon,nbp_lat,nhori2, nbp_lev,1,nbp_lev,nvert, 32,   &
-                       "ave(X)", zsto,zout)                                
-          CALL histdef(nid_tra2,"d_tr_sscav"//str2                        & 
-                       ,"cvltr d_tr_sscav"//str2,"mgS/m2/s",                 &
-                   nbp_lon,nbp_lat,nhori2, nbp_lev,1,nbp_lev,nvert, 32,      &
-                       "ave(X)", zsto,zout)                                  
-          CALL histdef(nid_tra2,"d_tr_sat"//str2                            &  
-                       ,"cvltr d_tr_sat"//str2,                             &  
-                       "mgS/m2/s",                                          & 
-                   nbp_lon,nbp_lat,nhori2, nbp_lev,1,nbp_lev,nvert, 32,     &
-                       "ave(X)", zsto,zout)                                  
-        CALL histdef(nid_tra2,"d_tr_uscav"//str2,                           & 
-                    "cvltr d_tr_uscav"//str2,                               &
-                       "mgS/m2/s",                                          &
-                   nbp_lon,nbp_lat,nhori2, nbp_lev,1,nbp_lev,nvert, 32,     &
-                       "ave(X)", zsto,zout)                                  
-        CALL histdef(nid_tra2,"d_tr_insc"//str2,                           &   !!!
-                    "cvltr d_tr_insc"//str2,                               &
-                       "mgS/m2/s",                                          &
-                   nbp_lon,nbp_lat,nhori2, nbp_lev,1,nbp_lev,nvert, 32,     &
-                       "ave(X)", zsto,zout)                                  
-        CALL histdef(nid_tra2,"d_tr_bcscav"//str2,                           & 
-                    "cvltr d_tr_bcscav"//str2,                               &
-                       "mgS/m2/s",                                          &
-                   nbp_lon,nbp_lat,nhori2, nbp_lev,1,nbp_lev,nvert, 32,     &
-                       "ave(X)", zsto,zout)                                  
-        CALL histdef(nid_tra2,"d_tr_evapls"//str2,                           & 
-                    "cvltr d_tr_evapls"//str2,                               &
-                       "mgS/m2/s",                                          &
-                   nbp_lon,nbp_lat,nhori2, nbp_lev,1,nbp_lev,nvert, 32,     &
-                       "ave(X)", zsto,zout)                                  
-        CALL histdef(nid_tra2,"d_tr_ls"//str2,                           & 
-                    "cvltr d_tr_ls"//str2,                               &
-                       "mgS/m2/s",                                          &
-                   nbp_lon,nbp_lat,nhori2, nbp_lev,1,nbp_lev,nvert, 32,     &
-                       "ave(X)", zsto,zout)                                     !!
-        CALL histdef(nid_tra2,"d_tr_dyn"//str2,                           & 
-                    "large-scale d_tr_dyn"//str2,                               &
-                       "mgS/m2/s",                                          &
-                   nbp_lon,nbp_lat,nhori2, nbp_lev,1,nbp_lev,nvert, 32,     &
-                       "ave(X)", zsto,zout)                                     !!
-        CALL histdef(nid_tra2,"d_tr_cl"//str2,                           & 
-                    "cvltr d_tr_cl"//str2,                               &
-                       "mgS/m2/s",                                          &
-                   nbp_lon,nbp_lat,nhori2, nbp_lev,1,nbp_lev,nvert, 32,     &
-                       "ave(X)", zsto,zout)                                 !!
-        CALL histdef(nid_tra2,"d_tr_th"//str2,                           & 
-                    "cvltr d_tr_th"//str2,                               &
-                       "mgS/m2/s",                                          &
-                   nbp_lon,nbp_lat,nhori2, nbp_lev,1,nbp_lev,nvert, 32,     &
-                       "ave(X)", zsto,zout)                                 !!
-                                                                             
-
-
-!
-          ENDDO
-!
-          CALL histdef(nid_tra2, "sed_ss", "Sedmet. Tr3",                   & 
-                       "mg/m2/s", nbp_lon,nbp_lat,nhori2, 1,1,1, -99,       &
-                         32,                                                &
-                       "ave(X)", zsto,zout)                                  
-!                                                                            
-          CALL histdef(nid_tra2, "sed_dust", "Sedmet. Tr4",                 &
-                       "mg/m2/s", nbp_lon,nbp_lat,nhori2, 1,1,1,            &
-                        -99, 32,                                            &
-                       "ave(X)", zsto,zout)                                  
-!                                                                            
-          CALL histdef(nid_tra2, "sed_dustsco", "Sedmet. Tr5",              &
-                       "mg/m2/s", nbp_lon,nbp_lat,nhori2, 1,1,1,            &
-                        -99, 32,                                            &
-                       "ave(X)", zsto,zout)                                  
-!                                                                            
-          CALL histdef(nid_tra2, "g2p_gas", "Gas2particle gas sink",       & 
-                   "mg-S/m2/s", nbp_lon,nbp_lat,nhori2, 1,1,1,-99, 32,     &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-          CALL histdef(nid_tra2, "g2p_aer", "Gas2particle tr2 src",        & 
-                       "mg/m2/s", nbp_lon,nbp_lat,nhori2, 1,1,1,-99,32,    &
-                       "ave(X)", zsto,zout)                                 
-!                                                                           
-!-------------------------------------------------------------------        
-!
-          CALL histend(nid_tra1)
-!
-          CALL histend(nid_tra2)
-!
-          CALL histend(nid_tra3)
-!
-!-------------------------------------------------------------------
-
-!       nbjour=1
-         ENDIF ! mpi root
-         ENDIF !--ok_histrac
-
-!
-!        IF (.NOT.edgar.AND.bateau) THEN
-!        PRINT *,'ATTENTION risque de compter double les bateaux'
-!        STOP
-!        ENDIF
-!
-!
-!
+!!AS deleting lines
+!!         IF (ok_histrac) THEN
+!!           IF (is_mpi_root .AND. is_omp_root) THEN
+!!-----many deleted lines----
+!!!       nbjour=1
+!!         ENDIF ! mpi root
+!!         ENDIF !--ok_histrac
+
 !$OMP END MASTER
 !$OMP BARRIER
@@ -2138,215 +1607,13 @@
 !$OMP MASTER
 
-      !ALLOCATE(d_tr_cl(klon,klev,nbtr),d_tr_dry(klon,nbtr))
       ALLOCATE(d_tr_dry(klon,nbtr))
       ALLOCATE(flux_tr_dry(klon,nbtr),d_tr_dec(klon,klev,nbtr))
-      !ALLOCATE(d_tr_cv(klon,klev,nbtr))
-      !ALLOCATE(d_tr_insc(klon,klev,nbtr),d_tr_bcscav(klon,klev,nbtr))
-      !ALLOCATE(d_tr_evapls(klon,klev,nbtr),d_tr_ls(klon,klev,nbtr))
-      !ALLOCATE(qPrls(klon,nbtr),d_tr_trsp(klon,klev,nbtr))
       ALLOCATE(qPrls(klon,nbtr),qPr(klon,klev,nbtr))
-      !ALLOCATE(d_tr_sscav(klon,klev,nbtr),d_tr_sat(klon,klev,nbtr))
-      !ALLOCATE(d_tr_uscav(klon,klev,nbtr),qPr(klon,klev,nbtr))
       ALLOCATE(qDi(klon,klev,nbtr))
       ALLOCATE(qPa(klon,klev,nbtr),qMel(klon,klev,nbtr))
       ALLOCATE(qTrdi(klon,klev,nbtr),dtrcvMA(klon,klev,nbtr))
-      !ALLOCATE(d_tr_th(klon,klev,nbtr))
-      !ALLOCATE(d_tr_lessi_impa(klon,klev,nbtr))
-      !ALLOCATE(d_tr_lessi_nucl(klon,klev,nbtr))
-
-      !ALLOCATE( diff_aod550_tot(klon)     )
-      !ALLOCATE( diag_aod670_tot(klon)     )
-      !ALLOCATE( diag_aod865_tot(klon)     )
-      !ALLOCATE( diff_aod550_tr2(klon)     )
-      !ALLOCATE( diag_aod670_tr2(klon)     )
-      !ALLOCATE( diag_aod865_tr2(klon)     )
-      !ALLOCATE( diag_aod550_ss(klon)      )
-      !ALLOCATE( diag_aod670_ss(klon)      )
-      !ALLOCATE( diag_aod865_ss(klon)      )
-      !ALLOCATE( diag_aod550_dust(klon)    )
-      !ALLOCATE( diag_aod670_dust(klon)    )
-      !ALLOCATE( diag_aod865_dust(klon)    )
-      !ALLOCATE( diag_aod550_dustsco(klon)  )
-      !ALLOCATE( diag_aod670_dustsco(klon)  )
-      !ALLOCATE( diag_aod865_dustsco(klon)  )
-
-
-      !ALLOCATE(  sconc01(klon)     )
-      !ALLOCATE(  trm01(klon)     )
-      !ALLOCATE(  sconc02(klon)     )
-      !ALLOCATE(  trm02(klon)     )
-      !ALLOCATE(  sconc03(klon)     )
-      !ALLOCATE(  trm03(klon)     )
-      !ALLOCATE(  sconc04(klon)     )
-      !ALLOCATE(  trm04(klon)     )
-      !ALLOCATE(  sconc05(klon)     )
-      !ALLOCATE(  trm05(klon)     )
-
-
-      !ALLOCATE(  flux01(klon)     )
-      !ALLOCATE(  flux02(klon)     )
-      !ALLOCATE(  flux03(klon)     )
-      !ALLOCATE(  flux04(klon)     )
-      !ALLOCATE(  flux05(klon)     )
-      !ALLOCATE(  ds01(klon)     )
-      !ALLOCATE(  ds02(klon)     )
-      !ALLOCATE(  ds03(klon)     )
-      !ALLOCATE(  ds04(klon)     )
-      !ALLOCATE(  ds05(klon)     )
-      !ALLOCATE(  dh01(klon)     )
-      !ALLOCATE(  dh02(klon)     )
-      !ALLOCATE(  dh03(klon)     )
-      !ALLOCATE(  dh04(klon)     )
-      !ALLOCATE(  dh05(klon)     )
-      !ALLOCATE(  dtrconv01(klon)     )
-      !ALLOCATE(  dtrconv02(klon)     )
-      !ALLOCATE(  dtrconv03(klon)     )
-      !ALLOCATE(  dtrconv04(klon)     )
-      !ALLOCATE(  dtrconv05(klon)     )
-      !ALLOCATE(  dtherm01(klon)     )
-      !ALLOCATE(  dtherm02(klon)     )
-      !ALLOCATE(  dtherm03(klon)     )
-      !ALLOCATE(  dtherm04(klon)     )
-      !ALLOCATE(  dtherm05(klon)     )
-      !ALLOCATE(  dhkecv01(klon)     )
-      !ALLOCATE(  dhkecv02(klon)     )
-      !ALLOCATE(  dhkecv03(klon)     )
-      !ALLOCATE(  dhkecv04(klon)     )
-      !ALLOCATE(  dhkecv05(klon)     )
-      !ALLOCATE(  d_tr_ds01(klon)     )
-      !ALLOCATE(  d_tr_ds02(klon)     )
-      !ALLOCATE(  d_tr_ds03(klon)     )
-      !ALLOCATE(  d_tr_ds04(klon)     )
-      !ALLOCATE(  d_tr_ds05(klon)     )
-      !ALLOCATE(  dhkelsc01(klon)     )
-      !ALLOCATE(  dhkelsc02(klon)     )
-      !ALLOCATE(  dhkelsc03(klon)     )
-      !ALLOCATE(  dhkelsc04(klon)     )
-      !ALLOCATE(  dhkelsc05(klon)     )
-      !ALLOCATE(  d_tr_cv01(klon,klev))
-      !ALLOCATE(  d_tr_cv02(klon,klev))
-      !ALLOCATE(  d_tr_cv03(klon,klev))
-      !ALLOCATE(  d_tr_cv04(klon,klev))
-      !ALLOCATE(  d_tr_cv05(klon,klev))
-      !ALLOCATE(  d_tr_trsp01(klon,klev))
-      !ALLOCATE(  d_tr_trsp02(klon,klev))
-      !ALLOCATE(  d_tr_trsp03(klon,klev))
-      !ALLOCATE(  d_tr_trsp04(klon,klev))
-      !ALLOCATE(  d_tr_trsp05(klon,klev))
-      !ALLOCATE(  d_tr_sscav01(klon,klev))
-      !ALLOCATE(  d_tr_sscav02(klon,klev))
-      !ALLOCATE(  d_tr_sscav03(klon,klev))
-      !ALLOCATE(  d_tr_sscav04(klon,klev))
-      !ALLOCATE(  d_tr_sscav05(klon,klev))
-      !ALLOCATE(  d_tr_sat01(klon,klev))
-      !ALLOCATE(  d_tr_sat02(klon,klev))
-      !ALLOCATE(  d_tr_sat03(klon,klev))
-      !ALLOCATE(  d_tr_sat04(klon,klev))
-      !ALLOCATE(  d_tr_sat05(klon,klev))
-      !ALLOCATE(  d_tr_uscav01(klon,klev))
-      !ALLOCATE(  d_tr_uscav02(klon,klev))
-      !ALLOCATE(  d_tr_uscav03(klon,klev))
-      !ALLOCATE(  d_tr_uscav04(klon,klev))
-      !ALLOCATE(  d_tr_uscav05(klon,klev))
-      
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-      !ALLOCATE(  d_tr_insc01(klon,klev))
-      !ALLOCATE(  d_tr_insc02(klon,klev))
-      !ALLOCATE(  d_tr_insc03(klon,klev))
-      !ALLOCATE(  d_tr_insc04(klon,klev))
-      !ALLOCATE(  d_tr_insc05(klon,klev))
-      !ALLOCATE(  d_tr_bcscav01(klon,klev))
-      !ALLOCATE(  d_tr_bcscav02(klon,klev))
-      !ALLOCATE(  d_tr_bcscav03(klon,klev))
-      !ALLOCATE(  d_tr_bcscav04(klon,klev))
-      !ALLOCATE(  d_tr_bcscav05(klon,klev))
-      !ALLOCATE(  d_tr_evapls01(klon,klev))
-      !ALLOCATE(  d_tr_evapls02(klon,klev))
-      !ALLOCATE(  d_tr_evapls03(klon,klev))
-      !ALLOCATE(  d_tr_evapls04(klon,klev))
-      !ALLOCATE(  d_tr_evapls05(klon,klev))
-      !ALLOCATE(  d_tr_ls01(klon,klev))
-      !ALLOCATE(  d_tr_ls02(klon,klev))
-      !ALLOCATE(  d_tr_ls03(klon,klev))
-      !ALLOCATE(  d_tr_ls04(klon,klev))
-      !ALLOCATE(  d_tr_ls05(klon,klev))
-      !ALLOCATE(  d_tr_dyn01(klon,klev))
-      !ALLOCATE(  d_tr_dyn02(klon,klev))
-      !ALLOCATE(  d_tr_dyn03(klon,klev))
-      !ALLOCATE(  d_tr_dyn04(klon,klev))
-      !ALLOCATE(  d_tr_dyn05(klon,klev))
-      !ALLOCATE(  d_tr_cl01(klon,klev))
-      !ALLOCATE(  d_tr_cl02(klon,klev))
-      !ALLOCATE(  d_tr_cl03(klon,klev))
-      !ALLOCATE(  d_tr_cl04(klon,klev))
-      !ALLOCATE(  d_tr_cl05(klon,klev))
-      !ALLOCATE(  d_tr_th01(klon,klev))
-      !ALLOCATE(  d_tr_th02(klon,klev))
-      !ALLOCATE(  d_tr_th03(klon,klev))
-      !ALLOCATE(  d_tr_th04(klon,klev))
-      !ALLOCATE(  d_tr_th05(klon,klev))
-
-      !ALLOCATE( sed_ss3D(klon,klev))
-      !ALLOCATE( sed_dust3D(klon,klev))
-      !ALLOCATE( sed_dustsco3D(klon,klev))
-
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-
-      !ALLOCATE( sed_ss(klon))
-      !ALLOCATE( sed_dust(klon))
-      !ALLOCATE( sed_dustsco(klon))
-      !ALLOCATE( his_g2pgas(klon))
-      !ALLOCATE( his_g2paer(klon))
-
-      !ALLOCATE( fluxbb(klon))
-      !ALLOCATE( fluxff(klon))
-      !ALLOCATE( fluxbcbb(klon))
-      !ALLOCATE( fluxbcff(klon))
-      !ALLOCATE( fluxbcnff(klon))
-      !ALLOCATE( fluxbcba(klon))
-      !ALLOCATE( fluxbc(klon))
-      !ALLOCATE( fluxombb(klon))
-      !ALLOCATE( fluxomff(klon))
-      !ALLOCATE( fluxomnff(klon))
-      !ALLOCATE( fluxomba(klon))
-      !ALLOCATE( fluxomnat(klon))
-      !ALLOCATE( fluxom(klon))
-      !ALLOCATE( fluxh2sff(klon))
-      !ALLOCATE( fluxh2snff(klon))
-      !ALLOCATE( fluxso2ff(klon))
-      !ALLOCATE( fluxso2nff(klon))
-      !ALLOCATE( fluxso2bb(klon))
-      !ALLOCATE( fluxso2vol(klon))
-      !ALLOCATE( fluxso2ba(klon))
-      !ALLOCATE( fluxso2(klon))
-      !ALLOCATE( fluxso4ff(klon))
-      !ALLOCATE( fluxso4nff(klon))
-      !ALLOCATE( fluxso4bb(klon))
-      !ALLOCATE( fluxso4ba(klon))
-      !ALLOCATE( fluxso4(klon))
-      !ALLOCATE( fluxdms(klon))
-      !ALLOCATE( fluxh2sbio(klon))
-      !ALLOCATE( fluxdustec(klon))
-      !ALLOCATE( fluxddfine(klon))
-      !ALLOCATE( fluxddcoa(klon))
-      !ALLOCATE( fluxddsco(klon))
-      !ALLOCATE( fluxdd(klon))
-      !ALLOCATE( fluxssfine(klon))
-      !ALLOCATE( fluxsscoa(klon))
-      !ALLOCATE( fluxss(klon))
-      !ALLOCATE( flux_sparam_ind(klon))
-      !ALLOCATE( flux_sparam_bb(klon))
-      !ALLOCATE( flux_sparam_ff(klon))
-      !ALLOCATE( flux_sparam_ddfine(klon))
-      !ALLOCATE( flux_sparam_ddcoa(klon))
-      !ALLOCATE( flux_sparam_ddsco(klon))
-      !ALLOCATE( flux_sparam_ssfine(klon))
-      !ALLOCATE( flux_sparam_sscoa(klon))
-      !ALLOCATE( u10m_ss(klon))
-      !ALLOCATE( v10m_ss(klon))
-
-
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
        ALLOCATE(d_tr_cv_o(klon,klev,nbtr))
        ALLOCATE(d_tr_trsp_o(klon,klev,nbtr))
@@ -2375,52 +1642,8 @@
        ALLOCATE(masque_aqua(klon))  
        ALLOCATE(masque_terra(klon))  
-!       ALLOCATE(aod550_aqua(klon))  
-!       ALLOCATE(aod550_terra(klon))  
-!       ALLOCATE(aod670_aqua(klon))  
-!       ALLOCATE(aod670_terra(klon))  
-!       ALLOCATE(aod865_aqua(klon))  
-!       ALLOCATE(aod865_terra(klon)) 
-
-      !ALLOCATE( aod550_terra(klon))  
-      !ALLOCATE( aod550_tr2_terra(klon))  
-      !ALLOCATE( aod550_ss_terra(klon))   
-      !ALLOCATE( aod550_dust_terra(klon))   
-      !ALLOCATE( aod550_dustsco_terra(klon))   
-      !ALLOCATE( aod670_terra(klon))   
-      !ALLOCATE( aod670_tr2_terra(klon))  
-      !ALLOCATE( aod670_ss_terra(klon))  
-      !ALLOCATE( aod670_dust_terra(klon))  
-      !ALLOCATE( aod670_dustsco_terra(klon))  
-      !ALLOCATE( aod865_terra(klon))   
-      !ALLOCATE( aod865_tr2_terra(klon))  
-      !ALLOCATE( aod865_ss_terra(klon))  
-      !ALLOCATE( aod865_dust_terra(klon))  
-      !ALLOCATE( aod865_dustsco_terra(klon))  
-
-      !ALLOCATE( aod550_aqua(klon))  
-      !ALLOCATE( aod550_tr2_aqua(klon))  
-      !ALLOCATE( aod550_ss_aqua(klon))   
-      !ALLOCATE( aod550_dust_aqua(klon))   
-      !ALLOCATE( aod550_dustsco_aqua(klon))   
-      !ALLOCATE( aod670_aqua(klon))   
-      !ALLOCATE( aod670_tr2_aqua(klon))  
-      !ALLOCATE( aod670_ss_aqua(klon))  
-      !ALLOCATE( aod670_dust_aqua(klon))  
-      !ALLOCATE( aod670_dustsco_aqua(klon))  
-      !ALLOCATE( aod865_aqua(klon))   
-      !ALLOCATE( aod865_tr2_aqua(klon))  
-      !ALLOCATE( aod865_ss_aqua(klon))  
-      !ALLOCATE( aod865_dust_aqua(klon))  
-      !ALLOCATE( aod865_dustsco_aqua(klon))  
  
 
        masque_aqua(:)=0
        masque_terra(:)=0
-!       aod550_aqua(:)=0.
-!       aod550_terra(:)=0.
-!       aod670_aqua(:)=0.
-!       aod670_terra(:)=0.
-!       aod865_aqua(:)=0.
-!       aod865_terra(:)=0.
 
       aod550_terra(:)=0.  
@@ -2715,4 +1938,5 @@
       CLOSE(76)
       CLOSE(77)
+      CLOSE(95)
 
       ENDIF ! mpi root
@@ -4116,10 +3340,4 @@
        !new utc day: put in 0 everything
 !JE20150518<<
-!       aod550_aqua(:) =0.
-!       aod550_terra(:) =0.
-!       aod670_aqua(:) =0.
-!       aod670_terra(:) =0.
-!       aod865_aqua(:) =0.
-!       aod865_terra(:) =0.
        masque_aqua(:) =0
        masque_terra(:) =0
@@ -4158,10 +3376,4 @@
 
       DO i=1,klon
-!         aod550_aqua(i)=aod550_aqua(i)+   &
-!                       masque_aqua_cur(i)*diff_aod550_tot(i)
-!         aod670_aqua(i)=aod670_aqua(i)+   &
-!                        masque_aqua_cur(i)*diag_aod670_tot(i)
-!         aod865_aqua(i)=aod865_aqua(i)+   &
-!                       masque_aqua_cur(i)*diag_aod865_tot(i)
 
        aod550_terra(i)=aod550_terra(i)+   &
@@ -4228,11 +3440,6 @@
        aod865_dustsco_aqua(i)= aod865_dustsco_aqua(i) + &
                        masque_aqua_cur(i)*diag_aod865_dustsco(i)
-!         aod550_terra(i)=aod550_terra(i)+  &
-!                       masque_terra_cur(i)*diff_aod550_tot(i)
-!         aod670_terra(i)=aod670_terra(i)+  &
-!                       masque_terra_cur(i)*diag_aod670_tot(i)
-!         aod865_terra(i)=aod865_terra(i)+   &
-!                       masque_terra_cur(i)*diag_aod865_tot(i)
-         masque_aqua(i)=masque_aqua(i)+masque_aqua_cur(i)
+
+               masque_aqua(i)=masque_aqua(i)+masque_aqua_cur(i)
          masque_terra(i)=masque_terra(i)+masque_terra_cur(i)
       ENDDO
@@ -4307,68 +3514,16 @@
                    aod865_dustsco_terra(i)= -999.
                ENDIF
-!              IF (masque_terra(i).GT. 0) THEN
-!                   aod550_terra(i) = aod550_terra(i)/masque_terra(i)
-!                   aod670_terra(i)=aod670_terra(i)/masque_terra(i)
-!                   aod865_terra(i)=aod865_terra(i)/masque_terra(i)
-!
-!               ELSE
-!                   aod550_terra(i) = -999.
-!                   aod670_terra(i) = -999.
-!                   aod865_terra(i) = -999.
-!               ENDIF
           ENDDO          
-!      !write  dbg
-!       CALL writefield_phy("aod550_aqua",aod550_aqua,1) 
-!       CALL writefield_phy("aod550_terra",aod550_terra,1) 
-!       CALL writefield_phy("masque_aqua",float(masque_aqua),1)
-!       CALL writefield_phy("masque_terra",float(masque_terra),1)
-
-
-      IF (ok_histrac) THEN
-!      write in output file
-      call gather(aod550_aqua,aod550_aqua_glo)
-      call gather(aod550_terra,aod550_terra_glo)
-      call gather(aod670_aqua,aod670_aqua_glo)
-      call gather(aod670_terra,aod670_terra_glo)
-      call gather(aod865_aqua,aod865_aqua_glo)
-      call gather(aod865_terra,aod865_terra_glo)
-
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-
-      CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, aod550_aqua_glo ,zx_tmp_2d)
-      CALL histwrite(nid_tra3,"taue550_aqua",itra,zx_tmp_2d, &
-                                      nbp_lon*(nbp_lat),ndex2d)
-
-      CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, aod550_terra_glo ,zx_tmp_2d)
-      CALL histwrite(nid_tra3,"taue550_terra",itra,zx_tmp_2d, &
-                                      nbp_lon*(nbp_lat),ndex2d)
-      CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, aod670_aqua_glo ,zx_tmp_2d)
-      CALL histwrite(nid_tra3,"taue670_aqua",itra,zx_tmp_2d, &
-                                      nbp_lon*(nbp_lat),ndex2d)
-      CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, aod670_terra_glo ,zx_tmp_2d)
-      CALL histwrite(nid_tra3,"taue670_terra",itra,zx_tmp_2d, &
-                                      nbp_lon*(nbp_lat),ndex2d)
-
-      CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, aod865_aqua_glo ,zx_tmp_2d)
-      CALL histwrite(nid_tra3,"taue865_aqua",itra,zx_tmp_2d, &
-                                      nbp_lon*(nbp_lat),ndex2d)
-      CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, aod865_terra_glo ,zx_tmp_2d)
-      CALL histwrite(nid_tra3,"taue865_terra",itra,zx_tmp_2d, &
-                                      nbp_lon*(nbp_lat),ndex2d)
-      ENDIF
-!$OMP END MASTER
-!$OMP BARRIER
-      ENDIF
-!       !put in 0 everything
-!       aod550_aqua(:) =0.
-!       aod550_terra(:) =0.
-!       aod670_aqua(:) =0.
-!       aod670_terra(:) =0.
-!       aod865_aqua(:) =0.
-!       aod865_terra(:) =0.
-!       masque_aqua(:) =0 
-!       masque_terra(:) =0 
-      ENDIF 
+
+!!AS deleting lines
+!!      IF (ok_histrac) THEN
+!!!!      write in output file
+!!----many deleted lines
+!!      ENDIF  !mpi_root
+!!!$OMP END MASTER
+!!!$OMP BARRIER
+!!      ENDIF  !--ok_histrac
+
+      ENDIF ! jH_cur...
 
 
@@ -4743,954 +3898,10 @@
       ENDDO
 
-      IF (ok_histrac) THEN
-!
-! SAVING VARIABLES IN TRACEUR
-!
-     call gather(diff_aod550_tot  ,auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo ,zx_tmp_2d_glo)
-     CALL histwrite(nid_tra3,"taue550",itra,zx_tmp_2d_glo,                 &
-                                      nbp_lon*(nbp_lat),ndex2d)              
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( diag_aod670_tot  , auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo,zx_tmp_2d_glo)     
-     CALL histwrite(nid_tra3,"taue670",itra,zx_tmp_2d_glo,                 &   
-                                      nbp_lon*(nbp_lat),ndex2d)              
-!                                                                       
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( diag_aod865_tot  , auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo,zx_tmp_2d_glo)     
-     CALL histwrite(nid_tra3,"taue865",itra,zx_tmp_2d_glo,                 &  
-                                      nbp_lon*(nbp_lat),ndex2d)              
-!                                                                       
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather(  diff_aod550_tr2 , auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo,zx_tmp_2d_glo)     
-     CALL histwrite(nid_tra3,"taue550_tr2",itra,zx_tmp_2d_glo,             &  
-                                      nbp_lon*(nbp_lat),ndex2d)              
-!                                                                       
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather(  diag_aod670_tr2 , auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo,zx_tmp_2d_glo)     
-     CALL histwrite(nid_tra3,"taue670_tr2",itra,zx_tmp_2d_glo,             &  
-                                      nbp_lon*(nbp_lat),ndex2d)              
-!                                                                       
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( diag_aod865_tr2  , auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo,zx_tmp_2d_glo)     
-     CALL histwrite(nid_tra3,"taue865_tr2",itra,zx_tmp_2d_glo,             &  
-                                      nbp_lon*(nbp_lat),ndex2d)              
-!                                                                       
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather(  diag_aod550_ss, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo ,zx_tmp_2d_glo)      
-     CALL histwrite(nid_tra3,"taue550_ss",itra,zx_tmp_2d_glo,              &  
-                                      nbp_lon*(nbp_lat),ndex2d)              
-!                                                                       
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( diag_aod670_ss , auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo ,zx_tmp_2d_glo)      
-     CALL histwrite(nid_tra3,"taue670_ss",itra,zx_tmp_2d_glo,              &  
-                                      nbp_lon*(nbp_lat),ndex2d)              
-!                                                                       
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( diag_aod865_ss, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo ,zx_tmp_2d_glo)      
-     CALL histwrite(nid_tra3,"taue865_ss",itra,zx_tmp_2d_glo,              &  
-                                      nbp_lon*(nbp_lat),ndex2d)              
-!                                                                       
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( diag_aod550_dust, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo ,zx_tmp_2d_glo)    
-     CALL histwrite(nid_tra3,"taue550_dust",itra,zx_tmp_2d_glo,             & 
-                                      nbp_lon*(nbp_lat),ndex2d)               
-!                                                                        
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( diag_aod670_dust, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo ,zx_tmp_2d_glo)     
-     CALL histwrite(nid_tra3,"taue670_dust",itra,zx_tmp_2d_glo,             &  
-                                      nbp_lon*(nbp_lat),ndex2d)               
-!                                                                        
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( diag_aod865_dust, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo ,zx_tmp_2d_glo)     
-     CALL histwrite(nid_tra3,"taue865_dust",itra,zx_tmp_2d_glo,             &  
-                                      nbp_lon*(nbp_lat),ndex2d)               
-!                                                                        
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( diag_aod550_dustsco, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo ,zx_tmp_2d_glo)  
-     CALL histwrite(nid_tra3,"taue550_dustsco",itra,zx_tmp_2d_glo,          &  
-                                      nbp_lon*(nbp_lat),ndex2d)               
-!                                                                        
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( diag_aod670_dustsco, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo ,zx_tmp_2d_glo)  
-     CALL histwrite(nid_tra3,"taue670_dustsco",itra,zx_tmp_2d_glo,          &  
-                                      nbp_lon*(nbp_lat),ndex2d)               
-!                                                                        
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( diag_aod865_dustsco, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1, klon_glo,nbp_lon,nbp_lat, auxklon_glo ,zx_tmp_2d_glo)  
-     CALL histwrite(nid_tra3,"taue865_dustsco",itra,zx_tmp_2d_glo,          &  
-                                      nbp_lon*(nbp_lat),ndex2d)               
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-                                                                         
-!$OMP MASTER
-     DO it=1,nbtr                                                        
-!                                                                        
-     WRITE(str2,'(i2.2)') it
-! 
-     call gather( trm, auxklonnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklonnbtr_glo(1,it) , zx_tmp_2d_glo)
-     CALL histwrite(nid_tra3,"trm"//str2,itra,zx_tmp_2d_glo,              & 
-                                         nbp_lon*(nbp_lat),ndex2d)          
-!                                                                      
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-     call gather( sconc_seri, auxklonnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklonnbtr_glo(1,it) ,zx_tmp_2d_glo)     
-     CALL histwrite(nid_tra3,"sconc"//str2,itra,zx_tmp_2d_glo,            &  
-                                         nbp_lon*(nbp_lat),ndex2d)          
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-!                                                                      
-! SAVING VARIABLES IN LESSIVAGE                                         
-!                                                                       
-     call gather( flux_tr, auxklonnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklonnbtr_glo(1,it) ,zx_tmp_2d_glo)
-     CALL histwrite(nid_tra2,"flux"//str2,itra,zx_tmp_2d_glo,               & 
-                    nbp_lon*(nbp_lat),ndex2d)                                 
-!                                                                        
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-     call gather( his_ds, auxklonnbtr_glo )
-!! $OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklonnbtr_glo(1,it) ,zx_tmp_2d_glo)           
-     CALL histwrite(nid_tra2,"ds"//str2,itra,zx_tmp_2d_glo,                 &  
-                    nbp_lon*(nbp_lat),ndex2d)                                 
-!                                                                        
-      ENDIF
-! !$OMP END MASTER
-! !$OMP BARRIER
-      ENDDO
-
-     DO it=1,nbtr
-     WRITE(str2,'(i2.2)') it
-      DO i=1, klon                                                        
-       zx_tmp_fi2d(i) = his_dhlsc(i,it)+his_dhcon(i,it)+               &  
-                        his_dhbclsc(i,it)+his_dhbccon(i,it)
-       his_dh(i,it)= his_dhlsc(i,it)+his_dhcon(i,it)+               &  
-                   his_dhbclsc(i,it)+his_dhbccon(i,it)
-
-      ENDDO
-!
-     call gather( zx_tmp_fi2d, auxklon_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)
-     CALL histwrite(nid_tra2,"dh"//str2,itra,zx_tmp_2d_glo,                  &
-                    nbp_lon*(nbp_lat),ndex2d)                                  
-!                                                                         
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-     call gather( his_dhkecv, auxklonnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklonnbtr_glo(1,it) ,zx_tmp_2d_glo)        
-     CALL histwrite(nid_tra2,"dhkecv"//str2,itra,zx_tmp_2d_glo,              &  
-                    nbp_lon*(nbp_lat),ndex2d)       
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!                           
-!                                                                         
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-     call gather( his_dhkelsc, auxklonnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklonnbtr_glo(1,it) ,zx_tmp_2d_glo)       
-     CALL histwrite(nid_tra2,"dhkelsc"//str2,itra,zx_tmp_2d_glo,             &  
-                    nbp_lon*(nbp_lat),ndex2d)                                  
-!                                                                         
-                                                                          
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-!    call gather( d_tr_cv_o,  auxklonklevnbtr_glo )
-     call gather( d_tr_cv,  auxklonklevnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(nbp_lev,klon_glo,nbp_lon,nbp_lat,auxklonklevnbtr_glo(1,1,it) ,             &  
-                      zx_tmp_3d_glo)                                         
-     CALL histwrite(nid_tra2,"d_tr_cv"//str2,itra,zx_tmp_3d_glo,             &  
-                                  nbp_lon*(nbp_lat)*nbp_lev,ndex3d)                  
-!
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-     call gather( d_tr_trsp_o, auxklonklevnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(nbp_lev,klon_glo,nbp_lon,nbp_lat,auxklonklevnbtr_glo(1,1,it) ,           &    
-                      zx_tmp_3d_glo)                                            
-     CALL histwrite(nid_tra2,"d_tr_trsp"//str2,itra,zx_tmp_3d_glo,           &    
-                                  nbp_lon*(nbp_lat)*nbp_lev,ndex3d)                  
-!
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-     call gather( d_tr_sscav_o, auxklonklevnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(nbp_lev,klon_glo,nbp_lon,nbp_lat,auxklonklevnbtr_glo(1,1,it) ,          &    
-                      zx_tmp_3d_glo)                                            
-     CALL histwrite(nid_tra2,"d_tr_sscav"//str2,itra,zx_tmp_3d_glo,          &    
-                                  nbp_lon*(nbp_lat)*nbp_lev,ndex3d)                  
-!
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-     call gather( d_tr_sat_o, auxklonklevnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(nbp_lev,klon_glo,nbp_lon,nbp_lat,auxklonklevnbtr_glo(1,1,it) ,            &    
-                      zx_tmp_3d_glo)                                            
-     CALL histwrite(nid_tra2,"d_tr_sat"//str2,itra,zx_tmp_3d_glo,            &    
-                                  nbp_lon*(nbp_lat)*nbp_lev,ndex3d)                  
-!
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-     call gather( d_tr_uscav_o, auxklonklevnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(nbp_lev,klon_glo,nbp_lon,nbp_lat,auxklonklevnbtr_glo(1,1,it) ,          &    
-                       zx_tmp_3d_glo)                                           
-     CALL histwrite(nid_tra2,"d_tr_uscav"//str2,itra,zx_tmp_3d_glo,          &    
-                                  nbp_lon*(nbp_lat)*nbp_lev,ndex3d)                  
-!                                                                            
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-        !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-     call gather( d_tr_insc_o, auxklonklevnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(nbp_lev,klon_glo,nbp_lon,nbp_lat,auxklonklevnbtr_glo(1,1,it) ,          &    
-                       zx_tmp_3d_glo)                                           
-     CALL histwrite(nid_tra2,"d_tr_insc"//str2,itra,zx_tmp_3d_glo,          &    
-                                  nbp_lon*(nbp_lat)*nbp_lev,ndex3d)                  
-!                                                                            
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-      
-      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-     call gather( d_tr_bcscav_o, auxklonklevnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(nbp_lev,klon_glo,nbp_lon,nbp_lat,auxklonklevnbtr_glo(1,1,it) ,          &    
-                       zx_tmp_3d_glo)                                           
-     CALL histwrite(nid_tra2,"d_tr_bcscav"//str2,itra,zx_tmp_3d_glo,          &    
-                                  nbp_lon*(nbp_lat)*nbp_lev,ndex3d)                  
-!                                                                            
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-     call gather( d_tr_evapls_o, auxklonklevnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(nbp_lev,klon_glo,nbp_lon,nbp_lat,auxklonklevnbtr_glo(1,1,it) ,          &    
-                       zx_tmp_3d_glo)                                           
-     CALL histwrite(nid_tra2,"d_tr_evapls"//str2,itra,zx_tmp_3d_glo,          &    
-                                  nbp_lon*(nbp_lat)*nbp_lev,ndex3d)                  
-!                                                                            
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-!    call gather( d_tr_ls_o, auxklonklevnbtr_glo )
-     call gather( d_tr_ls, auxklonklevnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(nbp_lev,klon_glo,nbp_lon,nbp_lat,auxklonklevnbtr_glo(1,1,it) ,          &    
-                       zx_tmp_3d_glo)                                           
-     CALL histwrite(nid_tra2,"d_tr_ls"//str2,itra,zx_tmp_3d_glo,          &    
-                                  nbp_lon*(nbp_lat)*nbp_lev,ndex3d)                  
-                                                                            
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-
-
-      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-!    call gather( d_tr_dyn_o, auxklonklevnbtr_glo )
-     call gather( d_tr_dyn, auxklonklevnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(nbp_lev,klon_glo,nbp_lon,nbp_lat,auxklonklevnbtr_glo(1,1,it) ,          &    
-                       zx_tmp_3d_glo)                                           
-     CALL histwrite(nid_tra2,"d_tr_dyn"//str2,itra,zx_tmp_3d_glo,          &    
-                                  nbp_lon*(nbp_lat)*nbp_lev,ndex3d)                  
-                                                                            
-      print*,'ECRTIURES TENDANCES MODIFIEES NON MAIS'
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-!    call gather( d_tr_cl_o, auxklonklevnbtr_glo )
-     call gather( d_tr_cl, auxklonklevnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(nbp_lev,klon_glo,nbp_lon,nbp_lat,auxklonklevnbtr_glo(1,1,it) ,          &    
-                       zx_tmp_3d_glo)                                           
-     CALL histwrite(nid_tra2,"d_tr_cl"//str2,itra,zx_tmp_3d_glo,          &    
-                                  nbp_lon*(nbp_lat)*nbp_lev,ndex3d)                  
-                                                                            
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-!    call gather( d_tr_th_o, auxklonklevnbtr_glo )
-     call gather( d_tr_th, auxklonklevnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(nbp_lev,klon_glo,nbp_lon,nbp_lat,auxklonklevnbtr_glo(1,1,it) ,          &    
-                       zx_tmp_3d_glo)                                           
-     CALL histwrite(nid_tra2,"d_tr_th"//str2,itra,zx_tmp_3d_glo,          &    
-                                  nbp_lon*(nbp_lat)*nbp_lev,ndex3d)                  
-                                                                            
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-     call gather( dtrconv,auxklonnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklonnbtr_glo(1,it) ,zx_tmp_2d_glo)           
-     CALL histwrite(nid_tra2,"dtrconv"//str2,itra,zx_tmp_2d_glo,            & 
-                    nbp_lon*(nbp_lat),ndex2d)                                 
-!                                                                        
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-     call gather( his_th, auxklonnbtr_glo )
-! !$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklonnbtr_glo(1,it) ,zx_tmp_2d_glo)           
-     CALL histwrite(nid_tra2,"dtherm"//str2,itra,zx_tmp_2d_glo,             &  
-                    nbp_lon*(nbp_lat),ndex2d)                                 
-      ENDIF ! mpi root
-! !$OMP END MASTER
-! !$OMP BARRIER
-!                                                                        
-                                                                         
-     ENDDO                                                               
-!
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( sed_ss, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)
-     CALL histwrite(nid_tra2,"sed_ss",itra,zx_tmp_2d_glo,                & 
-                    nbp_lon*(nbp_lat),ndex2d)                              
-!                                                                     
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( sed_dust, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)            
-     CALL histwrite(nid_tra2,"sed_dust",itra,zx_tmp_2d_glo,               & 
-                    nbp_lon*(nbp_lat),ndex2d)                               
-!                                                                      
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( sed_dustsco, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)          
-     CALL histwrite(nid_tra2,"sed_dustsco",itra,zx_tmp_2d_glo,              &
-                    nbp_lon*(nbp_lat),ndex2d)                                 
-!                                                                        
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( his_g2pgas, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)             
-     CALL histwrite(nid_tra2,"g2p_gas",itra,zx_tmp_2d_glo,                   & 
-                    nbp_lon*(nbp_lat),ndex2d)                                  
-!                                                                         
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( his_g2paer, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)              
-     CALL histwrite(nid_tra2,"g2p_aer",itra,zx_tmp_2d_glo,                   &  
-                    nbp_lon*(nbp_lat),ndex2d)                                  
-! SAVING VARIABLES IN HISTRAC                                              
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxbb, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-      CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)                  
-      CALL histwrite(nid_tra1,"fluxbb",itra,zx_tmp_2d_glo,                   & 
-                                    nbp_lon*(nbp_lat),ndex2d)                 
-!                                                                        
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxff, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)                 
-     CALL histwrite(nid_tra1,"fluxff",itra,zx_tmp_2d_glo,                   &  
-                                    nbp_lon*(nbp_lat),ndex2d)                 
-!                                                                        
-! ======================== BC =============================              
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxbcbb, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxbcbb",itra,zx_tmp_2d_glo,                 & 
-                                    nbp_lon*(nbp_lat),ndex2d)                 
-!                                                                        
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxbcff, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxbcff",itra,zx_tmp_2d_glo,                  & 
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-!                                                                         
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxbcnff, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxbcnff",itra,zx_tmp_2d_glo,                 &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-!                                                                         
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxbcba, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)                
-     CALL histwrite(nid_tra1,"fluxbcba",itra,zx_tmp_2d_glo,                  &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-!                                                                         
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxbc, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)                  
-     CALL histwrite(nid_tra1,"fluxbc",itra,zx_tmp_2d_glo,                    &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-! ======================== OM =============================               
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxombb, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)                
-     CALL histwrite(nid_tra1,"fluxombb",itra,zx_tmp_2d_glo,                  &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-!                                                                         
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxomff, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)                
-     CALL histwrite(nid_tra1,"fluxomff",itra,zx_tmp_2d_glo,                  &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-!                                                                         
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxomnff, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxomnff",itra,zx_tmp_2d_glo,                 &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-!                                                                         
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxomba, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)                
-     CALL histwrite(nid_tra1,"fluxomba",itra,zx_tmp_2d_glo,                  &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-!                                                                         
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxomnat, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxomnat",itra,zx_tmp_2d_glo,                 &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-!                                                                         
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxom, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)                  
-     CALL histwrite(nid_tra1,"fluxom",itra,zx_tmp_2d_glo,                    &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-! ======================== SO4 =============================              
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxso4ff, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxso4ff",itra,zx_tmp_2d_glo,                 &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-!                                                                         
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxso4nff, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)              
-     CALL histwrite(nid_tra1,"fluxso4nff",itra,zx_tmp_2d_glo,                &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-!                                                                         
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxso4bb, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxso4bb",itra,zx_tmp_2d_glo,                 &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-!                                                                         
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxso4ba, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxso4ba",itra,zx_tmp_2d_glo,                 &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-!                                                                         
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxso4, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo ,zx_tmp_2d_glo)                 
-     CALL histwrite(nid_tra1,"fluxso4",itra,zx_tmp_2d_glo,                   &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-! ======================== H2S =============================              
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxh2sff, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxh2sff",itra,zx_tmp_2d_glo,                 &  
-                                    nbp_lon*(nbp_lat),ndex2d)                  
-!                                                                         
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxh2snff, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)              
-     CALL histwrite(nid_tra1,"fluxh2snff",itra,zx_tmp_2d_glo,                 & 
-                                    nbp_lon*(nbp_lat),ndex2d)                   
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxh2sbio, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxh2sbio",itra,zx_tmp_2d_glo,                 &  
-                                    nbp_lon*(nbp_lat),ndex2d)                   
-! ======================== SO2 =============================               
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxso2ff, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)                
-     CALL histwrite(nid_tra1,"fluxso2ff",itra,zx_tmp_2d_glo,                  &  
-                                    nbp_lon*(nbp_lat),ndex2d)                   
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxso2nff, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxso2nff",itra,zx_tmp_2d_glo,                 &  
-                                    nbp_lon*(nbp_lat),ndex2d)                   
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxso2bb, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)                
-     CALL histwrite(nid_tra1,"fluxso2bb",itra,zx_tmp_2d_glo,                  &  
-                                    nbp_lon*(nbp_lat),ndex2d)                   
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxso2vol, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxso2vol",itra,zx_tmp_2d_glo,                 &  
-                                    nbp_lon*(nbp_lat),ndex2d)                   
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxso2ba, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)                
-     CALL histwrite(nid_tra1,"fluxso2ba",itra,zx_tmp_2d_glo,                  &  
-                                    nbp_lon*(nbp_lat),ndex2d)                   
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxso2, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)                  
-     CALL histwrite(nid_tra1,"fluxso2",itra,zx_tmp_2d_glo,                    &  
-                                    nbp_lon*(nbp_lat),ndex2d)                   
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxdms, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)                  
-     CALL histwrite(nid_tra1,"fluxdms",itra,zx_tmp_2d_glo,                    &  
-                                    nbp_lon*(nbp_lat),ndex2d)                   
-! ======================== DD =============================                
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxdustec, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxdustec",itra,zx_tmp_2d_glo,                 &  
-                                    nbp_lon*(nbp_lat),ndex2d)                   
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxddfine, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxddfine",itra,zx_tmp_2d_glo,                 &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxddcoa, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)                
-     CALL histwrite(nid_tra1,"fluxddcoa",itra,zx_tmp_2d_glo,                  &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxddsco, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)                
-     CALL histwrite(nid_tra1,"fluxddsco",itra,zx_tmp_2d_glo,                  &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxdd, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)                   
-     CALL histwrite(nid_tra1,"fluxdd",itra,zx_tmp_2d_glo,                     &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-! ======================== SS =============================                
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxssfine, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)               
-     CALL histwrite(nid_tra1,"fluxssfine",itra,zx_tmp_2d_glo,                 &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxsscoa, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)                
-     CALL histwrite(nid_tra1,"fluxsscoa",itra,zx_tmp_2d_glo,                  &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( fluxss, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)                   
-     CALL histwrite(nid_tra1,"fluxss",itra,zx_tmp_2d_glo,                     &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-
-!     call gather( , auxklon_glo )
-!!!!      IF (is_mpi_root .AND. is_omp_root) THEN
-!nhl     CALL gr_fi_ecrit(nbp_lev,klon_glo,nbp_lon,nbp_lat,fluxso4chem,zx_tmp_3d_glo)       
-!nhl     CALL histwrite(nid_tra1,"fluxso4chem",itra,zx_tmp_3d_glo,            &  
-!nhl    .                             nbp_lon*(nbp_lat)*nbp_lev,ndex3d)            
-!                                                                          
-     call gather( flux_sparam_ind, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)          
-     CALL histwrite(nid_tra1,"flux_sparam_ind",itra,zx_tmp_2d_glo,            &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( flux_sparam_bb, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)           
-     CALL histwrite(nid_tra1,"flux_sparam_bb",itra,zx_tmp_2d_glo,             &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( flux_sparam_ff, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)           
-     CALL histwrite(nid_tra1,"flux_sparam_ff",itra,zx_tmp_2d_glo,             &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( flux_sparam_ddfine, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)       
-     CALL histwrite(nid_tra1,"flux_sparam_ddfine",itra,zx_tmp_2d_glo,         &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( flux_sparam_ddcoa, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)        
-     CALL histwrite(nid_tra1,"flux_sparam_ddcoa",itra,zx_tmp_2d_glo,          &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( flux_sparam_ddsco, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)        
-     CALL histwrite(nid_tra1,"flux_sparam_ddsco",itra,zx_tmp_2d_glo,          &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( flux_sparam_ssfine, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)       
-     CALL histwrite(nid_tra1,"flux_sparam_ssfine",itra,zx_tmp_2d_glo,         &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( flux_sparam_sscoa, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)        
-     CALL histwrite(nid_tra1,"flux_sparam_sscoa",itra,zx_tmp_2d_glo,          &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( u10m_ec, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)                  
-     CALL histwrite(nid_tra1,"u10m",itra,zx_tmp_2d_glo,                       &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-     call gather( v10m_ec, auxklon_glo )
-!$OMP MASTER
-      IF (is_mpi_root .AND. is_omp_root) THEN
-     CALL gr_fi_ecrit(1,klon_glo,nbp_lon,nbp_lat,auxklon_glo,zx_tmp_2d_glo)                  
-     CALL histwrite(nid_tra1,"v10m",itra,zx_tmp_2d_glo,                       &  
-                                  nbp_lon*(nbp_lat),ndex2d)                     
-!                                                                          
-!     call gather( , auxklon_glo )
-!!!   !$OMP MASTER
-!      IF (is_mpi_root .AND. is_omp_root) THEN
-!nhl     CALL gr_fi_ecrit(nbp_lev,klon_glo,nbp_lon,nbp_lat,flux_sparam_sulf,zx_tmp_3d_glo)  
-!nhl     CALL histwrite(nid_tra1,"flux_sparam_sulf",itra,zx_tmp_3d_glo,       &  
-!nhl    .                             nbp_lon*(nbp_lat)*nbp_lev,ndex3d)            
-!                                                                          
-      ENDIF ! mpi root
-!$OMP END MASTER
-!$OMP BARRIER
-
-      ENDIF ! ok_histrac                                                    
+!AS: commenting out and deleting lines 
+!!      IF (ok_histrac) THEN
+!!
+!! SAVING VARIABLES IN TRACEUR
+!!----- many lines deleted----
+!!      ENDIF ! ok_histrac                                                    
                                                                             
 
@@ -5751,6 +3962,4 @@
        dhkelsc04(i)=0. 
        dhkelsc05(i)=0.
-!       u10m_ss(i)=u10m_ec(i)
-!       v10m_ss(i)=v10m_ec(i)
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
@@ -6192,4 +4401,5 @@
       END SUBROUTINE readregions_spl
 
+!! AS: SUBROUTINE readscaleparams_spl pas appellee
       SUBROUTINE readscaleparams_spl(scale_param, nbreg, &
                                              filescaleparams)
Index: LMDZ6/trunk/libf/phylmd/Dust/read_vent.F
===================================================================
--- LMDZ6/trunk/libf/phylmd/Dust/read_vent.F	(revision 3805)
+++ LMDZ6/trunk/libf/phylmd/Dust/read_vent.F	(revision 3806)
@@ -33,7 +33,7 @@
 c
          ncidu1=NCOPN('u10m.nc',NCNOWRIT,rcode)
-         varidu1=NCVID(ncidu1,'UWND',rcode)
+         varidu1=NCVID(ncidu1,'U10M',rcode)
          ncidv1=NCOPN('v10m.nc',NCNOWRIT,rcode)
-         varidv1=NCVID(ncidv1,'VWND',rcode)
+         varidv1=NCVID(ncidv1,'V10M',rcode)
 c
       endif
