Changeset 1318 for LMDZ4/branches/LMDZ4V5.0-dev/libf/cosp/cosp.F90

Timestamp:

Feb 23, 2010, 9:12:08 PM (14 years ago)

Author:

yann meurdesoif

Message:

YM : Parallelisation COSP MPI+OpenMP

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/cosp/cosp.F90 (modified) (6 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/cosp/cosp.F90

@@ -26 +26 @@
   USE MOD_COSP_TYPES
   USE MOD_COSP_SIMULATOR
+  USE mod_phys_lmdz_para
+  USE mod_grid_phy_lmdz
   IMPLICIT NONE
 
@@ -61 +63 @@
   real :: minv,maxv
   real :: maxp,minp
-  integer,dimension(:),allocatable :: & ! Dimensions nPoints
+  integer,dimension(:),save,  allocatable :: & ! Dimensions nPoints
                   seed    !  It is recommended that the seed is set to a different value for each model
                           !  gridbox it is called on, as it is possible that the choice of the same 
                           !  seed value every time may introduce some statistical bias in the results, 
                           !  particularly for low values of NCOL.
+!$OMP THREADPRIVATE(seed)
+  real,dimension(:),allocatable :: rseed    !  It is recommended that the seed is set to a different value for each model
   ! Types used in one iteration
   type(cosp_gridbox) :: gbx_it
@@ -76 +80 @@
   type(cosp_radarstats) :: stradar_it
   type(cosp_lidarstats) :: stlidar_it
+  
+  logical,save :: first_cosp=.TRUE.
+!$OMP THREADPRIVATE(first_cosp)
   
   !++++++++++ Dimensions ++++++++++++
@@ -162 +169 @@
   endif
 
-   
+  if (first_cosp) then   
    ! We base the seed in the decimal part of the surface pressure.
-   allocate(seed(Npoints))
-   seed = int(gbx%psfc) ! This is to avoid division by zero when Npoints = 1   
+     allocate(seed(Npoints))
+
+     allocate(rseed(klon_glo))
+     CALL gather(gbx%psfc,rseed)
+     call bcast(rseed)
+!   seed = int(gbx%psfc) ! This is to avoid division by zero when Npoints = 1   
       ! Roj Oct/2008 ... Note: seed value of 0 caused me some problems + I want to 
       ! randomize for each call to COSP even when Npoints ==1
-   minp = minval(gbx%psfc)
-   maxp = maxval(gbx%psfc)
-   if (Npoints .gt. 1) seed=int((gbx%psfc-minp)/(maxp-minp)*100000) + 1
-
-  
+     minp = minval(rseed)
+     maxp = maxval(rseed)
+   
+     if (Npoints .gt. 1) THEN
+       seed=int((gbx%psfc-minp)/(maxp-minp)*100000) + 1
+     else
+       seed=int(gbx%psfc-minp)
+     endif
+
+     deallocate(rseed)
+     first_cosp=.false.
+   endif
+    
    if (gbx%Npoints_it >= gbx%Npoints) then ! One iteration gbx%Npoints
         call cosp_iter(overlap,seed,cfg,vgrid,gbx,sgx,sgradar,sglidar,isccp,misr,stradar,stlidar)
@@ -273 +292 @@
         call free_cosp_lidarstats(stlidar_it)
    endif
-   deallocate(seed)
 
     
@@ -305 +323 @@
   integer :: i,j,k
   real,dimension(:,:),pointer :: column_frac_out ! Array with one column of frac_out
-  integer :: scops_debug=0    !  set to non-zero value to print out inputs for debugging in SCOPS
+  integer,parameter :: scops_debug=0    !  set to non-zero value to print out inputs for debugging in SCOPS
+  
   real,dimension(:, :),allocatable :: cca_scops,ls_p_rate,cv_p_rate, &
                      tca_scops ! Cloud cover in each model level (HORIZONTAL gridbox fraction) of total cloud.