Changeset 472 for codes/icosagcm/trunk/src/physics_interface.f90

Timestamp:

09/05/16 14:59:46 (8 years ago)

Author:

ymipsl

Message:

Repairing openMP :
physic_column and physic_dcmip2016 seems OK. Only for thread assigned to domain.

YM

File:

• codes/icosagcm/trunk/src/physics_interface.f90 (modified) (4 diffs)

codes/icosagcm/trunk/src/physics_interface.f90

@@ -26 +26 @@
 
   TYPE(t_physics_inout), SAVE :: physics_inout
+!$OMP THREADPRIVATE(physics_inout) 
  
 ! pack_info contains indices used by pack/unpack routines
@@ -39 +40 @@
 
   TYPE(t_pack_info), ALLOCATABLE, SAVE :: pack_info(:)
-  
+!$OMP THREADPRIVATE(pack_info)
+
+
   INTERFACE pack_field
      MODULE PROCEDURE pack_2D
@@ -76 +79 @@
     IMPLICIT NONE
     INTEGER :: ind, offset
-    !$OMP MASTER
+
     offset=0
     ALLOCATE(pack_info(ndomain))
     DO ind=1,ndomain
+       IF (.NOT. assigned_domain(ind)) CYCLE
        CALL swap_dimensions(ind)
        CALL swap_geometry(ind)
@@ -87 +91 @@
     END DO
     physics_inout%ngrid = offset
-    !$OMP END MASTER
-    !$OMP BARRIER
+
   END SUBROUTINE init_pack_before