Changeset 704 for LMDZ4/branches/V3_test/libf/phylmd/plevel.F

Timestamp:

Aug 17, 2006, 5:41:51 PM (18 years ago)

Author:

Laurent Fairhead

Message:

Inclusion des modifs de Y. Meurdesoif pour la version V3
LF

File:

• LMDZ4/branches/V3_test/libf/phylmd/plevel.F (modified) (3 diffs)

LMDZ4/branches/V3_test/libf/phylmd/plevel.F

@@ -7 +7 @@
 c================================================================
 c================================================================
-
+      USE dimphy
       IMPLICIT none
 
-#include "dimensions.h"
-#include "dimphy.h"
+cym#include "dimensions.h"
+cy#include "dimphy.h"
 
 c================================================================
@@ -42 +42 @@
 c   -------
 
-      INTEGER lt(klon), lb(klon)
-      REAL ptop, pbot, aist(klon), aisb(klon)
+cym      INTEGER lt(klon), lb(klon)
+cym      REAL ptop, pbot, aist(klon), aisb(klon)
 
-      save lt,lb,ptop,pbot,aist,aisb
-
+cym      save lt,lb,ptop,pbot,aist,aisb
+      INTEGER,ALLOCATABLE,SAVE,DIMENSION(:) :: lt,lb,aist,aisb
+c$OMP THREADPRIVATE(lt,lb,aist,aisb)      
+      REAL,SAVE :: ptop, pbot
+c$OMP THREADPRIVATE(ptop, pbot)      
+      LOGICAL,SAVE :: first = .true.
+c$OMP THREADPRIVATE(first)
       INTEGER i, k
 c
-
+      if (first) then
+        allocate(lt(klon),lb(klon),aist(klon),aisb(klon))
+        first=.false.
+      endif
+      
 c=====================================================================
       if (lnew) then
-c   on réinitialise les réindicages et les poids
+c   on r�nitialise les r�ndicages et les poids
 c=====================================================================
 
@@ -87 +96 @@
 c
 c ...   Modif . P. Le Van    ( 20/01/98) ....
-c       Modif Frédéric Hourdin (3/01/02)
+c       Modif Fr��ic Hourdin (3/01/02)
 
         aist(i) = LOG( pgcm(i,lb(i))/ pres )