Changeset 3743 for trunk/LMDZ.MARS/libf

Timestamp:

Apr 28, 2025, 6:30:07 PM (9 months ago)

Author:

emillour

Message:

Mars PCM:
OpenMP bug fix: add some missing threadprivate statements in moldiff_MPF
EM

File:

• trunk/LMDZ.MARS/libf/aeronomars/moldiff_MPF.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/aeronomars/moldiff_MPF.F90

@@ -63 +63 @@
       real*8 zz(nlayer)
       real*8,dimension(:,:),allocatable,save :: qq,qnew,qint,FacMass
+!$OMP THREADPRIVATE(qq,qnew,qint,FacMass)
       real*8,dimension(:,:),allocatable,save :: rhoK,rhokinit
+!$OMP THREADPRIVATE(rhoK,rhokinit)
       real*8 rhoT(nlayer)
       real*8 dmmeandz(nlayer)
@@ -76 +78 @@
       real*8,dimension(:),allocatable :: alpha,beta,delta,ksi,eps,zeta,prod,loss
       real*8,dimension(:),allocatable,save :: wi,Wad,Uthermal,Lambdaexo,Hspecie,alphaT
+!$OMP THREADPRIVATE(wi,Wad,Uthermal,Lambdaexo,Hspecie,alphaT)
       real*8,dimension(:),allocatable,save :: Mtot1,Mtot2,Mraf1,Mraf2
+!$OMP THREADPRIVATE(Mtot1,Mtot2,Mraf1,Mraf2)
       integer,parameter :: ListeDiffNb=16
       character(len=20),dimension(ListeDiffNb) :: ListeDiff
@@ -92 +96 @@
 ! vertical index limit for the molecular diffusion
       integer,save :: il0  
-
-!$OMP THREADPRIVATE(qq,qnew,qint,FacMass,rhoK,rhokinit,wi,Wad,Uthermal,Lambdaexo,Hspecie,Mtot1,Mtot2,Mraf1,Mraf2,i_h2,i_h,i_d,il0)
+!$OMP THREADPRIVATE(i_h2,i_h,i_d,i_hd,il0)
 
 ! Tracer indexes in the GCM: