Changeset 1023 for LMDZ4/trunk/libf/phylmd/surf_land_orchidee_mod.F90

Timestamp:

Oct 23, 2008, 11:26:49 AM (16 years ago)

Author:

lsce

Message:

Bug parallélisme en MPI/OPENMP lorsque pour un grand nombre de CPU, le processus maître et/ou la tâche maitre se retrouve sans points de terre, et sont donc exlus de la liste des domaines dans orchidee.
YM

File:

• LMDZ4/trunk/libf/phylmd/surf_land_orchidee_mod.F90 (modified) (3 diffs)

LMDZ4/trunk/libf/phylmd/surf_land_orchidee_mod.F90

@@ -279 +279 @@
        ENDIF
    
-
+       CALL bcast(lon_scat)
+       CALL bcast(lat_scat)
 !
 ! Allouer et initialiser le tableau des voisins et des fraction de continents
@@ -366 +367 @@
        CALL Init_synchro_omp
        
-       IF (knon_mpi > 0) THEN
+       IF (knon > 0) THEN
          CALL Init_intersurf(nbp_lon,nbp_lat,knon,ktindex,offset,orch_omp_size,orch_omp_rank,orch_comm)
        ENDIF
@@ -505 +506 @@
       CALL MPI_COMM_SPLIT(COMM_LMDZ_PHY,color,mpi_rank,orch_comm,ierr)
 #endif
-
-   ENDIF
+    
+    ENDIF
+    CALL bcast_omp(orch_comm)
+    
+    IF (knon_mpi /= 0) THEN
+      orch_omp_size=0
+      DO i=0,omp_size-1
+        IF (knon_omp_para(i) /=0) THEN
+          orch_omp_size=orch_omp_size+1
+          IF (i==omp_rank) orch_omp_rank=orch_omp_size-1
+        ENDIF
+      ENDDO
+    ENDIF
    
-   IF (knon_mpi /= 0) THEN
-     orch_omp_size=0
-     DO i=0,omp_size-1
-       IF (knon_omp_para(i) /=0) THEN
-         orch_omp_size=orch_omp_size+1
-         IF (i==omp_rank) orch_omp_rank=orch_omp_size-1
-       ENDIF
-     ENDDO
-   ENDIF
-       
     
   END SUBROUTINE Get_orchidee_communicator