Changeset 3814 for LMDZ6

Timestamp:

Jan 29, 2021, 1:20:22 PM (3 years ago)

Author:

asima

Message:

SPLA : unplugging Jeronimo Escribano’s assimilation emission coefficient files, in a minimalist, non-distructive way.
Each file reading is put under an IF("ASSIM"=="YES"), always false.
Thorough removal of assimilation coefficients will follow.

File:

• LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90 (modified) (11 diffs)

LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90

@@ -695 +695 @@
 
   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 +711 @@
   call  readregionsdims2_spl(nbreg_bb,fileregionsdimsbb)
   call  readregionsdims2_spl(nbreg_wstardust,fileregionsdimswstar)
+  ENDIF ! ASSIM
+! fin debranchage
 
 !readregions_spl()
@@ -748 +757 @@
 
   !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 +769 @@
   param_wstarBLperregion(:)=0.
   param_wstarWAKEperregion(:)=0.
-
 
 
@@ -926 +936 @@
       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 +1287 @@
 
   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 +1301 @@
         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 +1869 @@
 
 
-
-       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 +1967 @@
 !$OMP END MASTER
 !$OMP BARRIER
+
+      ENDIF  ! ASSIM 
 
       ENDIF  ! debutphy
@@ -3565 +3591 @@
 !  SAVING AEROSOL RELATED VARIABLES INTO FILE
 !======================================================================
-!
-!JE20141224      IF (ok_histrac) THEN
 !
       ndex2d = 0
@@ -3702 +3726 @@
          fluxss(i)=fluxssfine(i)+fluxsscoa(i)
       ENDDO
+
 !      prepare outputs cvltr