Changeset 3924 for LMDZ6/trunk/libf/phylmd

Timestamp:

Jun 8, 2021, 6:33:45 PM (4 years ago)

Author:

Laurent Fairhead

Message:

Version modifiee par Camille pour version isotopique

File:

• LMDZ6/trunk/libf/phylmd/infotrac_phy.F90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/infotrac_phy.F90

@@ -20 +20 @@
   INTEGER, SAVE :: nbtr
 !$OMP THREADPRIVATE(nbtr)
+
+  INTEGER, SAVE :: nqtottr
+!$OMP THREADPRIVATE(nqtottr)
 
 ! ThL : number of CO2 tracers                   ModThL
@@ -97 +100 @@
     INTEGER,SAVE :: niso,ntraceurs_zone,ntraciso
 !$OMP THREADPRIVATE(niso,ntraceurs_zone,ntraciso)
+
+    INTEGER, ALLOCATABLE, DIMENSION(:), SAVE ::  itr_indice ! numéro iq entre 1 et nqtot qui correspond au traceur itr entre 1 et nqtottr
+!$OMP THREADPRIVATE(itr_indice)
  
 CONTAINS
 
-  SUBROUTINE init_infotrac_phy(nqtot_,nqo_,nbtr_,nqCO2_,tname_,ttext_,type_trac_,&
+  SUBROUTINE init_infotrac_phy(nqtot_,nqo_,nbtr_,nqtottr_,nqCO2_,tname_,ttext_,type_trac_,&
                                niadv_,conv_flg_,pbl_flg_,solsym_,&
                                nqfils_,nqdesc_,nqdesc_tot_,iqfils_,iqpere_,&
@@ -108 +114 @@
                                iso_indnum_,zone_num_,phase_num_,&
                                indnum_fn_num_,index_trac_,&
-                               niso_,ntraceurs_zone_,ntraciso_&
+                               niso_,ntraceurs_zone_,ntraciso_,itr_indice_&
 #ifdef CPP_StratAer
                                ,nbtr_bin_,nbtr_sulgas_&
@@ -122 +128 @@
     INTEGER,INTENT(IN) :: nqo_
     INTEGER,INTENT(IN) :: nbtr_
+    INTEGER,INTENT(IN) :: nqtottr_
     INTEGER,INTENT(IN) :: nqCO2_
 #ifdef CPP_StratAer
@@ -162 +169 @@
     INTEGER,INTENT(IN) :: ntraceurs_zone_
     INTEGER,INTENT(IN) :: ntraciso_
+    INTEGER,INTENT(IN) :: itr_indice_(nqtottr_)
 
     CHARACTER(LEN=30) :: modname="init_infotrac_phy"
@@ -169 +177 @@
     nbtr=nbtr_
     nqCO2=nqCO2_
+    nqtottr=nqtottr_
 #ifdef CPP_StratAer
     nbtr_bin=nbtr_bin_
@@ -242 +251 @@
       ALLOCATE(index_trac(ntraceurs_zone,niso))
       index_trac(:,:)=index_trac_(:,:)
+
+      ALLOCATE(itr_indice(nqtottr))
+      itr_indice(:)=itr_indice_(:)
     ENDIF ! of IF(ok_isotopes)