Changeset 1232 for LMDZ4/trunk

Timestamp:

Aug 27, 2009, 3:23:40 PM (15 years ago)

Author:

musat

Message:

Correction Bug originel sur lb, lt : lb(klon), lt(klon) ==> lb(klon,klev), lt(klon,klev)
IM

File:

• LMDZ4/trunk/libf/phylmd/plevel_new.F (modified) (6 diffs)

LMDZ4/trunk/libf/phylmd/plevel_new.F

@@ -46 +46 @@
 
 cym      save lt,lb,ptop,pbot,aist,aisb
-      INTEGER,ALLOCATABLE,SAVE,DIMENSION(:) :: lt,lb
+      INTEGER,ALLOCATABLE,SAVE,DIMENSION(:,:) :: lt,lb
       REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: aist,aisb
 c$OMP THREADPRIVATE(lt,lb,aist,aisb)      
@@ -57 +57 @@
 c
       if (first) then
-         allocate(lt(klon),lb(klon))
+         allocate(lt(klon,klevSTD),lb(klon,klevSTD))
          allocate(aist(klon,klevSTD),aisb(klon, klevSTD))
          first=.false.
@@ -78 +78 @@
                IF ( ABS(pres(nlev)-pgcm(i,ilev) ) .LT.
      &              ABS(pres(nlev)-pgcm(i,1)) ) THEN
-                  lt(i) = ilev  ! 2
-                  lb(i) = ilev-1 ! 1
+                  lt(i,nlev) = ilev  ! 2
+                  lb(i,nlev) = ilev-1 ! 1
                ELSE
-                  lt(i) = 2
-                  lb(i) = 1
+                  lt(i,nlev) = 2
+                  lb(i,nlev) = 1
                ENDIF
             ENDDO
@@ -90 +90 @@
                   ptop = pgcm(i,k+1)
                   IF (ptop.LE.pres(nlev) .AND. pbot.GE.pres(nlev)) THEN
-                     lt(i) = k+1
-                     lb(i) = k
+                     lt(i,nlev) = k+1
+                     lb(i,nlev) = k
                   ENDIF
                ENDDO
@@ -103 +103 @@
 c     Modif Frederic Hourdin (3/01/02)
                
-               aist(i,nlev) = LOG( pgcm(i,lb(i))/ pres(nlev) )
-     &              / LOG( pgcm(i,lb(i))/ pgcm(i,lt(i)) )
-               aisb(i,nlev) = LOG( pres(nlev) / pgcm(i,lt(i)) )
-     &              / LOG( pgcm(i,lb(i))/ pgcm(i,lt(i)))
+               aist(i,nlev) = LOG( pgcm(i,lb(i,nlev))/ pres(nlev) )
+     &              / LOG( pgcm(i,lb(i,nlev))/ pgcm(i,lt(i,nlev)) )
+               aisb(i,nlev) = LOG( pres(nlev) / pgcm(i,lt(i,nlev)) )
+     &              / LOG( pgcm(i,lb(i,nlev))/ pgcm(i,lt(i,nlev)))
             ENDDO
          ENDDO
@@ -124 +124 @@
             ELSE
                Qpres(i,nlev) = 
-     &              Qgcm(i,lb(i))*aisb(i,nlev) +
-     &              Qgcm(i,lt(i))*aist(i,nlev)
+     &              Qgcm(i,lb(i,nlev))*aisb(i,nlev) +
+     &              Qgcm(i,lt(i,nlev))*aist(i,nlev)
             ENDIF
          ENDDO