Changeset 1759 for LMDZ5/trunk/libf/phydev

Timestamp:

May 24, 2013, 11:52:12 AM (11 years ago)

Author:

Ehouarn Millour

Message:

IOIPSL routine getin is not threadsafe, so when running in OpenMP, it should be called by only one thread (and the result copied to other threads in the case of threadprivate variables).
EM

File:

• LMDZ5/trunk/libf/phydev/physiq.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phydev/physiq.F90

@@ -75 +75 @@
 integer,save :: iwrite_phys ! output every iwrite_phys physics step
 !$OMP THREADPRIVATE(iwrite_phys)
+integer :: iwrite_phys_omp ! intermediate variable to read iwrite_phys
 real :: t_ops ! frequency of the IOIPSL operations (eg average over...)
 real :: t_wrt ! frequency of the IOIPSL outputs
@@ -86 +87 @@
 ! Initialize outputs:
   itau0=0
-  iwrite_phys=1 !default: output every physics timestep
-  call getin("iwrite_phys",iwrite_phys)
+!$OMP MASTER
+  iwrite_phys_omp=1 !default: output every physics timestep
+  ! NB: getin() is not threadsafe; only one thread should call it.
+  call getin("iwrite_phys",iwrite_phys_omp)
+!$OMP END MASTER
+!$OMP BARRIER
+  iwrite_phys=iwrite_phys_omp
   t_ops=pdtphys*iwrite_phys ! frequency of the IOIPSL operation
   t_wrt=pdtphys*iwrite_phys ! frequency of the outputs in the file