Changeset 1858 for LMDZ5/trunk/libf

Timestamp:

Sep 3, 2013, 11:32:25 AM (11 years ago)

Author:

yann meurdesoif

Message:

omp_chunk has now a default value of (llm+1)/omp_size instead of 1.

YM

Location:

LMDZ5/trunk/libf/dyn3dmem

Files:

• conf_gcm.F (modified) (3 diffs)
• parallel_lmdz.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/dyn3dmem/conf_gcm.F

@@ -17 +17 @@
       use mod_filtre_fft_loc, ONLY : use_filtre_fft_loc=>use_filtre_fft
       use mod_hallo, ONLY : use_mpi_alloc
-      USE parallel_lmdz, ONLY : omp_chunk
       USE infotrac, ONLY : type_trac
       use assert_m, only: assert
@@ -908 +907 @@
       use_mpi_alloc=.FALSE.
       CALL getin('use_mpi_alloc',use_mpi_alloc)
-
-!Config  Key  = omp_chunk
-!Config  Desc = taille des blocs openmp
-!Config  Def  = 1
-!Config  Help = defini la taille des packets d'it�ration openmp
-!Config         distribu�e � chaque t�che lors de l'entr�e dans une
-!Config         boucle parall�lis�e
-  
-      omp_chunk=1
-      CALL getin('omp_chunk',omp_chunk)
 
 !Config key = ok_strato
@@ -1006 +995 @@
       write(lunout,*)' use_filtre_fft = ', use_filtre_fft
       write(lunout,*)' use_mpi_alloc = ', use_mpi_alloc
-      write(lunout,*)' omp_chunk = ', omp_chunk
       write(lunout,*)' ok_strato = ', ok_strato
       write(lunout,*)' ok_gradsfile = ', ok_gradsfile

LMDZ5/trunk/libf/dyn3dmem/parallel_lmdz.F90

@@ -4 +4 @@
   MODULE parallel_lmdz
   USE mod_const_mpi
-    
+#ifdef CPP_IOIPSL
+      use IOIPSL
+#else
+! if not using IOIPSL, we still need to use (a local version of) getin
+      use ioipsl_getincom
+#endif    
     INTEGER,PARAMETER :: halo_max=3
     
@@ -221 +226 @@
         omp_size=OMP_GET_NUM_THREADS()
 !$OMP END MASTER
-        omp_rank=OMP_GET_THREAD_NUM()    
+        omp_rank=OMP_GET_THREAD_NUM()  
+
+!Config  Key  = omp_chunk
+!Config  Desc = taille des blocs openmp
+!Config  Def  = 1
+!Config  Help = defini la taille des packets d'it�ration openmp
+!Config         distribue a chaque tache lors de l'entree dans une
+!Config         boucle parallelisee
+  
+      omp_chunk=(llm+1)/omp_size
+      IF (MOD(llm+1,omp_size)/=0) omp_chunk=omp_chunk+1
+      CALL getin('omp_chunk',omp_chunk)
+        
 #else    
         omp_size=1