Index: LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90	(revision 3813)
+++ LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90	(revision 3814)
@@ -695,4 +695,11 @@
 
   ALLOCATE(  tsol(klon)              )
+
+!AS: IF permettant le debranchage des coefs de Jeronimo Escribano: fichiers *_meta
+! nbreg_* sont initialisés à 1 dans phytracr_spl, if debutphy, 
+! avant d'appeler la subroutine presente, phytracr_spl_ini 
+! (phytracr_spl_ini appele readregionsdims2_spl, 
+! qui lit et fait "bcast" de nbreg_ind,_bb,_dust,_wstardust dans fichiers regions_*_meta)
+IF("ASSIM"=="YES") THEN
   fileregionsdimsind='regions_ind_meta'
   fileregionsdimsdust='regions_dustacc_meta'
@@ -704,4 +711,6 @@
   call  readregionsdims2_spl(nbreg_bb,fileregionsdimsbb)
   call  readregionsdims2_spl(nbreg_wstardust,fileregionsdimswstar)
+  ENDIF ! ASSIM
+! fin debranchage
 
 !readregions_spl()
@@ -748,4 +757,6 @@
 
   !temporal hardcoded null inicialization of assimilation emmision factors
+!AS: scale_param sont ensuite lus dans modvalues.nc 
+! par la subroutine read_scalenc, appelee par readscaleparamsnc_spl
   scale_param_ssacc=1.
   scale_param_sscoa=1.
@@ -758,5 +769,4 @@
   param_wstarBLperregion(:)=0.
   param_wstarWAKEperregion(:)=0.
-
 
 
@@ -926,5 +936,5 @@
       INTEGER :: aux_mask1
       INTEGER :: aux_mask2
-      INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: iregion_so4 !Defines regions for SO4
+      INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: iregion_so4 !Defines regions for SO4 ; AS: PAS UTILISE!
       INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: iregion_ind  !Defines regions for SO2, BC & OM
       INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: iregion_bb   !Defines regions for SO2, BC & OM
@@ -1277,5 +1287,7 @@
 
   filescaleparams='modvalues.nc'
-  CALL readscaleparamsnc_spl(scale_param_ind,                        &
+!AS: debranchage de lecture des coefs d'assmilation de Jeronimo Escribano
+  IF("ASSIM"=="YES") THEN
+    CALL readscaleparamsnc_spl(scale_param_ind,                        &
         nbreg_ind, paramname_ind,                                    &
         scale_param_ff, nbreg_ind,paramname_ff,                      &
@@ -1289,5 +1301,9 @@
         scale_param_sscoa  ,  paramname_sscoa,                    &
            filescaleparams,ijulday,jH_cur, pdtphys,debutphy)
-! add seasalt
+  ENDIF ! ASSIM
+!AS: le commentaire suivant "add seasalt" ne semble pas avoir ete mis en pratique.
+! Des fichiers regions_ssacc et _sscoa existent mais ne semblent pas lus.
+! Ca reste donc aux valeurs initialisées: nbreg_ss=1, scale_param_ss*=1, cf fichiers ss et modvalues
+!! add seasalt
 
   print *,'JE : check scale_params'
@@ -1853,7 +1869,15 @@
 
 
-
-       IF (debutphy) then
-
+      IF (debutphy) then
+
+! AS: initialisation des indices par point de grille physique iregion_*
+! (variables tenant de l'assimilation, a eliminer dans un 2eme temps) 
+       iregion_dust(:)=1
+       iregion_ind(:)=1
+       iregion_bb(:)=1
+       iregion_wstardust(:)=1
+
+!AS: lecture des indices dans fichiers "regions_*" eliminee par IF("ASSIM"="YES") (faux donc)
+       IF("ASSIM"=="YES") THEN
       c_FullName1='regions_dustacc'
       !c_FullName1='regions_dust'
@@ -1943,4 +1967,6 @@
 !$OMP END MASTER
 !$OMP BARRIER
+
+      ENDIF  ! ASSIM 
 
       ENDIF  ! debutphy
@@ -3565,6 +3591,4 @@
 !  SAVING AEROSOL RELATED VARIABLES INTO FILE
 !======================================================================
-!
-!JE20141224      IF (ok_histrac) THEN
 !
       ndex2d = 0
@@ -3702,4 +3726,5 @@
          fluxss(i)=fluxssfine(i)+fluxsscoa(i)
       ENDDO
+
 !      prepare outputs cvltr