Changeset 3743 for trunk/LMDZ.MARS

Timestamp:

Apr 28, 2025, 6:30:07 PM (8 weeks ago)

Author:

emillour

Message:

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

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/aeronomars/moldiff_MPF.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4823 +4823 @@
 == 28/04/2025 == JBC
 Simplification and clarification for the prescription of the atmospheric water profile in 1D: the boolean 'prescribed_h2ovap' activates the option, the file "profile_prescribed_h2ovap" defines the prescribed profile and the value 'h2ovap_relax_time', if positive, activates the relaxation and gives the constant of time.
+
+== 28/04/2025 == EM
+OpenMP bug fix: add some missing threadprivate statements in moldiff_MPF

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: