Changeset 766 for LMDZ4/trunk/libf/phylmd/iniphysiq.F

Timestamp:

Jun 4, 2007, 4:34:47 PM (17 years ago)

Author:

Laurent Fairhead

Message:

Merge entre la version V3_conv et le HEAD
YM, JG, LF

File:

• LMDZ4/trunk/libf/phylmd/iniphysiq.F (modified) (4 diffs)

LMDZ4/trunk/libf/phylmd/iniphysiq.F

@@ -9 +9 @@
      $           plat,plon,parea,pcu,pcv,
      $           prad,pg,pr,pcpp)
+      USE dimphy,klon2=>klon,klon=>klon2
+      USE comgeomphy
       IMPLICIT NONE
 c
@@ -43 +45 @@
 c   -------------
  
-#include "dimensions.h"
-#include "dimphy.h"
-#include "comgeomphy.h"
-
+cym#include "dimensions.h"
+cym#include "dimphy.h"
+cym#include "comgeomphy.h"
+#include "YOMCST.h"
       REAL prad,pg,pr,pcpp,punjours
  
@@ -52 +54 @@
       REAL plat(ngrid),plon(ngrid),parea(klon),pcu(klon),pcv(klon)
       INTEGER pdayref
+      INTEGER :: ibegin,iend
  
       REAL ptimestep
@@ -70 +73 @@
          STOP
       ENDIF
+c$OMP PARALLEL PRIVATE(ibegin,iend) 
+c$OMP+         SHARED(parea,pcu,pcv,plon,plat)
+      ibegin=klon_begin+klon_omp_begin(omp_rank)-1
+      iend=ibegin+klon_omp-1
 
-      airephy=parea
-      cuphy=pcu
-      cvphy=pcv
-      rlond = plon
-      rlatd = plat
+      airephy(1:klon_omp)=parea(ibegin:iend)
+      cuphy(1:klon_omp)=pcu(ibegin:iend)
+      cvphy(1:klon_omp)=pcv(ibegin:iend)
+      rlond(1:klon_omp) = plon(ibegin:iend)
+      rlatd(1:klon_omp) = plat(ibegin:iend)
 
       call suphec
+
+c$OMP END PARALLEL
+
       print*,'ATTENTION !!! TRAVAILLER SUR INIPHYSIQ'
       print*,'CONTROLE DES LATITUDES, LONGITUDES, PARAMETRES ...'