Changeset 836 for LMDZ4/trunk/libf

Timestamp:

Aug 30, 2007, 10:55:25 AM (17 years ago)

Author:

lsce

Message:

AC + ACo + JG : correction d'un bug suite a la parallelisation de la routine

gath2cpl qui entrainait un mauvais masque pour un decoupage
de l'atmosphere sur plus de deux processeurs.

File:

• LMDZ4/trunk/libf/phylmd/cpl_mod.F90 (modified) (6 diffs)

LMDZ4/trunk/libf/phylmd/cpl_mod.F90

@@ -17 +17 @@
 ! Use statements
 !*************************************************************************************
-  USE dimphy, ONLY : klon, zmasq
+  USE dimphy, ONLY : klon
   USE mod_phys_lmdz_para
   USE ioipsl
@@ -71 +71 @@
   REAL, ALLOCATABLE, DIMENSION(:,:), SAVE   :: pctsrf_sav   
   !$OMP THREADPRIVATE(pctsrf_sav)
-  REAL, ALLOCATABLE, DIMENSION(:,:), SAVE   :: zmasq2D
-  !$OMP THREADPRIVATE(zmasq2D)
   INTEGER, ALLOCATABLE, DIMENSION(:), SAVE  :: unity
   !$OMP THREADPRIVATE(unity)
@@ -181 +179 @@
     ALLOCATE(read_alb_sic(iim, jj_nb), stat = error)
     sum_error = sum_error + error
-    ALLOCATE(zmasq2D(iim, jj_nb), stat = error)
-    sum_error = sum_error + error    
-   
 
     IF (sum_error /= 0) THEN
@@ -202 +197 @@
     cpl_taux = 0.   ; cpl_tauy = 0.  ; cpl_rriv2D = 0. ; cpl_rcoa2D = 0.
     cpl_rlic2D = 0. ; cpl_windsp = 0.
-
-!*************************************************************************************
-! Transform the land-ocean mask into 2D grid. 
-! Colorize zmasq2D with 99 so that after gath2cpl points not valid can be recognized. 
-!
-!*************************************************************************************
-    zmasq2D(:,:) = 99.
-    CALL gath2cpl(zmasq, zmasq2D, klon, unity)
 
 !*************************************************************************************
@@ -936 +923 @@
     INTEGER, DIMENSION(iim*(jjm+1))                      :: ndexct
     REAL                                                 :: Up, Down
-    REAL, DIMENSION(iim, jj_nb)                       :: tmp_lon, tmp_lat
-    REAL, DIMENSION(iim, jj_nb, 4)                    :: pctsrf2D
-    REAL, DIMENSION(iim, jj_nb)                       :: deno
+    REAL, DIMENSION(iim, jj_nb)                          :: tmp_lon, tmp_lat
+    REAL, DIMENSION(iim, jj_nb, 4)                       :: pctsrf2D
+    REAL, DIMENSION(iim, jj_nb)                          :: deno
     CHARACTER(len = 20)                                  :: modname = 'cpl_send_all'
     CHARACTER(len = 80)                                  :: abort_message
    
 ! Variables with fields to coupler
-    REAL, DIMENSION(iim, jj_nb)                       :: tmp_taux
-    REAL, DIMENSION(iim, jj_nb)                       :: tmp_tauy
-    REAL, DIMENSION(iim, jj_nb)                       :: tmp_calv
+    REAL, DIMENSION(iim, jj_nb)                          :: tmp_taux
+    REAL, DIMENSION(iim, jj_nb)                          :: tmp_tauy
+    REAL, DIMENSION(iim, jj_nb)                          :: tmp_calv
 ! Table with all fields to send to coupler
-    REAL, DIMENSION(iim, jj_nb, jpflda2o1+jpflda2o2)  :: tab_flds
+    REAL, DIMENSION(iim, jj_nb, jpflda2o1+jpflda2o2)     :: tab_flds
 #ifdef CPP_PARA
     INCLUDE 'mpif.h'
@@ -1036 +1023 @@
     tmp_tauy(:,:)    = 0.0
     
-    ! For all valid grid cells not entier land
-    WHERE (zmasq2D /= 1. .AND. zmasq2D /=99. )
-       deno =  pctsrf2D(:,:,is_oce) + pctsrf2D(:,:,is_sic) ! fraction oce+seaice
+
+    ! fraction oce+seaice
+    deno =  pctsrf2D(:,:,is_oce) + pctsrf2D(:,:,is_sic) 
+    ! For all valid grid cells containing some fraction of ocean or sea-ice
+    WHERE ( deno(:,:) /= 0 )
+       tab_flds(:,:,15) = cpl_rain2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) +    &
+            cpl_rain2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
+       tab_flds(:,:,16) = cpl_snow2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) +    &
+            cpl_snow2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
        
-       tab_flds(:,:,15) = cpl_rain2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno +    &
-            cpl_rain2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno
-       tab_flds(:,:,16) = cpl_snow2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno +    &
-            cpl_snow2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno
-       
-       tmp_taux = cpl_taux2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno +    &
-            cpl_taux2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno
-       tmp_tauy = cpl_tauy2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno +    &
-            cpl_tauy2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno
+       tmp_taux = cpl_taux2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) +    &
+            cpl_taux2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
+       tmp_tauy = cpl_tauy2D(:,:,1) * pctsrf2D(:,:,is_oce) / deno(:,:) +    &
+            cpl_tauy2D(:,:,2) * pctsrf2D(:,:,is_sic) / deno(:,:)
     ENDWHERE