Changeset 3960 for LMDZ6/branches/LMDZ-tracers/libf

Timestamp:

Jul 13, 2021, 12:41:02 PM (3 years ago)

Author:

acozic

Message:

remove !$OMP directives to allow parallelisation in openmp
infotrac is call outside of a parallel omp zone, all variables in SAVE are shared by all process when we enter in OMP zone. If we put these variables in THREADPRIVATE we need to call bcast to share them.

File:

• LMDZ6/branches/LMDZ-tracers/libf/dyn3d_common/infotrac.F90 (modified) (3 diffs)

LMDZ6/branches/LMDZ-tracers/libf/dyn3d_common/infotrac.F90

@@ -112 +112 @@
                               nbIso                        !--- Number of available isotopes family
   CHARACTER(LEN=256), SAVE :: type_trac                    !--- Keyword for tracers type
-!$OMP THREADPRIVATE(nqtot, nbtr, nqo, nbIso, type_trac)
 
   !=== DERIVED TYPES EMBEDDING MOST INFORMATIONS ABOUT TRACERS AND ISOTOPES FAMILIES
   TYPE(tra), TARGET,  SAVE, ALLOCATABLE ::  tracers(:)     !=== TRACERS DESCRIPTORS VECTOR
   TYPE(iso), TARGET,  SAVE, ALLOCATABLE :: isotopes(:)     !=== ISOTOPES PARAMETERS VECTOR
-!$OMP THREADPRIVATE(tracers, isotopes)
 
   !=== ALIASES FOR CURRENTLY SELECTED ISOTOPES FAMILY OF VARIABLES EMBEDDED IN THE VECTOR "isotopes"
@@ -131 +129 @@
   INTEGER,            SAVE, POINTER     :: iZonIso(:,:)    !--- INDEX IN "isoTrac" AS f(tagging zone, isotope)
   INTEGER,            SAVE, POINTER     :: iTraPha(:,:)    !--- INDEX IN "isoTrac" AS f(isotopic tracer, phase)
-!$OMP THREADPRIVATE(isotope, ixIso,iH2O, isoCheck, isoKeys, isoName,isoZone,isoPhas, niso,nzon,npha,nitr, iZonIso,iTraPha)
 
   !=== VARIABLES EMBEDDED IN "tracers", BUT DUPLICATED, AS THEY ARE RATHER FREQUENTLY USED + VARIABLES FOR INCA
@@ -141 +138 @@
                                               niadv(:)     !--- Indexes of true tracers  (<=nqtot, such that iadv(idx)>0)
   CHARACTER(LEN=8),   SAVE, ALLOCATABLE ::   solsym(:)     !--- Names from INCA                                    (nbtr)
-!OMP THREADPRIVATE(tnat, alpha_ideal, conv_flg, pbl_flg, itr_indice, niadv, solsym)
 
 #ifdef CPP_StratAer
   !=== SPECIFIC TO STRATOSPHERIC AEROSOLS (CK/OB)
   INTEGER, SAVE :: nbtr_bin, nbtr_sulgas, id_OCS_strat, id_H2SO4_strat, id_SO2_strat, id_BIN01_strat, id_TEST_strat
-!OMP THREADPRIVATE(nbtr_bin, nbtr_sulgas, id_OCS_strat, id_H2SO4_strat, id_SO2_strat, id_BIN01_strat, id_TEST_strat)
 #endif