Changeset 5569

Timestamp:

Mar 7, 2025, 1:00:54 PM (5 months ago)

Author:

Laurent Fairhead

Message:

The THREADPRIVATE statements that were added in r5483 are erroneous as all variables concerned are
initialized under a OMP MASTER flag and thus not initialised for non master threads

Location:

LMDZ6/trunk/libf/phylmd

Files:

• iophy.F90 (modified) (2 diffs)
• phys_local_var_mod.F90 (modified) (1 diff)
• phys_state_var_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/iophy.F90

@@ -13 +13 @@
   INTEGER, SAVE :: itau_iophy
   LOGICAL :: check_dim = .false.
-!$OMP THREADPRIVATE(io_lat,io_lon,phys_domain_id,npstn,nptabij,itau_iophy)
+!$OMP THREADPRIVATE(itau_iophy)
 
   INTERFACE histwrite_phy
@@ -971 +971 @@
   REAL, ALLOCATABLE, DIMENSION(:) :: fieldok
   logical, save :: is_active = .true.
-!$OMP THREADPRIVATE(is_active)
 
   IF (check_dim .AND. is_master) WRITE(lunout,*)'histwrite2d_phy for ',trim(var%name)

LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -1105 +1105 @@
       ALLOCATE(proba_notrig(klon), random_notrig(klon))
       ALLOCATE(coef_clos(klon), coef_clos_eff(klon))
+      coef_clos(:)=0.
+      coef_clos_eff(:)=0.
 
       ALLOCATE(dnwd0(klon, klev))

LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -657 +657 @@
 !
       ALLOCATE(Mipsh(klon,klev))
+      Mipsh(:,:)=0.0
       ALLOCATE(Ma(klon,klev))
       ALLOCATE(qcondc(klon,klev))