Changeset 472 for codes/icosagcm/trunk/src/physics_dcmip2016.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_dcmip2016.f90 (modified) (1 diff)

codes/icosagcm/trunk/src/physics_dcmip2016.f90

@@ -6 +6 @@
 !$OMP THREADPRIVATE(testcase)
 
-  TYPE(t_field),POINTER :: f_out_i(:)
-  REAL(rstd),POINTER :: out_i(:,:)
-
-  TYPE(t_field),POINTER  :: f_precl(:)
-  REAL(rstd),ALLOCATABLE :: precl_packed(:)
-
-  TYPE(t_field),POINTER  :: f_Q1(:)
-  TYPE(t_field),POINTER  :: f_Q2(:)
-  TYPE(t_field),POINTER  :: f_PS(:)
-  TYPE(t_field),POINTER  :: f_rhodz(:)
-  TYPE(t_field),POINTER  :: f_Q1_col_int(:)
-  TYPE(t_field),POINTER  :: f_Q2_col_int(:)
+  TYPE(t_field),SAVE,POINTER :: f_out_i(:)
+  REAL(rstd),SAVE,POINTER :: out_i(:,:)
+
+  TYPE(t_field),SAVE,POINTER  :: f_precl(:)
+  REAL(rstd),SAVE,ALLOCATABLE :: precl_packed(:)
+!$OMP THREADPRIVATE(precl_packed)  
+
+  TYPE(t_field),SAVE,POINTER  :: f_Q1(:)
+  TYPE(t_field),SAVE,POINTER  :: f_Q2(:)
+  TYPE(t_field),SAVE,POINTER  :: f_PS(:)
+  TYPE(t_field),SAVE,POINTER  :: f_rhodz(:)
+  TYPE(t_field),SAVE,POINTER  :: f_Q1_col_int(:)
+  TYPE(t_field),SAVE,POINTER  :: f_Q2_col_int(:)
   PUBLIC :: init_physics, full_physics, write_physics