Changeset 6120 for LMDZ6/trunk/libf

Timestamp:

Mar 23, 2026, 11:29:20 PM (9 days ago)

Author:

lguez

Message:

Set ncid_climoz to OpenMP threadprivate

Set ncid_climoz to OpenMP threadprivate in procedure
physiq. ncid_climoz is only used by the master thread of root MPI
process so we can set it to threadprivate and define it only for the
master thread of root MPI process.

Location:

LMDZ6/trunk/libf/phylmd

Files:

• open_climoz_m.f90 (modified) (1 diff)
• physiq_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/open_climoz_m.f90

@@ -25 +25 @@
 ! and in strictly ascending order.
 !-------------------------------------------------------------------------------
-! Arguments (OpenMP shared):
   INTEGER, INTENT(OUT):: ncid      !--- "climoz_LMDZ.nc" identifier
+  ! Following arguments are OpenMP shared:
   REAL, allocatable, intent(out):: press_in_cen(:) !--- at cells centers
   REAL, allocatable, INTENT(OUT):: press_in_edg(:) !--- at the interfaces (pressure intervals)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -1069 +1069 @@
     !     2: read two ozone climatologies, the average day and night
     !     climatology and the daylight climatology
-    !$OMP THREADPRIVATE(read_climoz)
+    INTEGER,SAVE :: ncid_climoz                ! NetCDF file identifier
+    !$OMP THREADPRIVATE(read_climoz, ncid_climoz)
 
     !  Following variables are OpenMP shared
-    INTEGER,SAVE :: ncid_climoz                ! NetCDF file identifier
     REAL, ALLOCATABLE, SAVE :: press_cen_climoz(:) ! Pressure levels
     REAL, ALLOCATABLE, SAVE :: press_edg_climoz(:) ! Edges of pressure intervals
@@ -2124 +2124 @@
        ENDIF
 
-       !$omp single
+       !$omp master
        IF (read_climoz >= 1) CALL open_climoz(ncid_climoz, press_cen_climoz,   &
             press_edg_climoz, time_climoz, ok_daily_climoz, adjust_tropopause)
-       !$omp end single
+       !$omp end master
        !
        !IM betaCRF