Changeset 2078

Timestamp:

Jul 4, 2014, 12:44:45 PM (10 years ago)

Author:

idelkadi

Message:

Rajout des directives OPENMP pour la lecture des cles et parametres de convection et wakes dans les fichier conv_param.data, ep_param.data et wake_param.data

Location:

LMDZ5/trunk/libf/phylmd

Files:

• cv3_routines.F90 (modified) (4 diffs)
• wake.F90 (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/cv3_routines.F90

@@ -5 +5 @@
 
 SUBROUTINE cv3_param(nd, delt)
+
+  use mod_phys_lmdz_para
   IMPLICIT NONE
 
@@ -91 +93 @@
     betad = 10.0 ! original value (from convect 4.3)
 
+   !$OMP MASTER
     OPEN (99, FILE='conv_param.data', STATUS='old', FORM='formatted', ERR=9999)
     READ (99, *, END=9998) dpbase
@@ -124 +127 @@
     WRITE (*, *) 'tlcrit=', tlcrit
 ! IM end: ajout fis. reglage ep
+  !$OMP END MASTER
+
+   CALL bcast(dpbase)
+   CALL bcast(pbcrit)
+   CALL bcast(ptcrit)
+   CALL bcast(sigdz)
+   CALL bcast(spfac)
+   CALL bcast(tau)
+   CALL bcast(flag_wb)
+   CALL bcast(wbmax)
+
+   CALL bcast(flag_epkeorig)
+   CALL bcast(elcrit)
+   CALL bcast(tlcrit)
 
     first = .FALSE.
@@ -129 +146 @@
   END IF ! (first)
 
-! print*,'tau=',tau
   beta = 1.0 - delt/tau
   alpha1 = 1.5E-3

LMDZ5/trunk/libf/phylmd/wake.F90

@@ -19 +19 @@
 
   USE dimphy
+  use mod_phys_lmdz_para
   IMPLICIT NONE
   ! ============================================================================
@@ -151 +152 @@
   ! Variables à fixer
   REAL alon
-  REAL coefgw
-  REAL :: wdens_ref
-  REAL stark
-  REAL alpk
+  LOGICAL, SAVE :: first = .TRUE.
+  !$OMP THREADPRIVATE(first)
+  REAL, SAVE ::  stark, wdens_ref, coefgw, alpk, crep_upper, crep_sol  
+  !$OMP THREADPRIVATE(stark, wdens_ref, coefgw, alpk, crep_upper, crep_sol)
   REAL delta_t_min
   INTEGER nsub
@@ -232 +233 @@
 
   REAL, DIMENSION (klon, klev) :: crep
-  REAL crep_upper, crep_sol
 
   REAL, DIMENSION (klon, klev) :: ppi
@@ -280 +280 @@
   ! alpk = 0.05
 
+ if (first) then
   stark = 0.33
   alpk = 0.25
@@ -288 +289 @@
 
   ! cc nrlmd Lecture du fichier wake_param.data
+ !$OMP MASTER
   OPEN (99, FILE='wake_param.data', STATUS='old', FORM='formatted', ERR=9999)
   READ (99, *, END=9998) stark
@@ -296 +298 @@
   CLOSE (99)
 9999 CONTINUE
+ !$OMP END MASTER
+  CALL bcast(stark)
+  CALL bcast(alpk)
+  CALL bcast(wdens_ref)
+  CALL bcast(coefgw)
+
+  first=.false.
+ endif
 
   ! Initialisation de toutes des densites a wdens_ref.