Changeset 1146 for LMDZ4/trunk/libf/phylmd/cv3_cine.F

Timestamp:

Apr 9, 2009, 12:11:35 PM (15 years ago)

Author:

Laurent Fairhead

Message:

Réintegration dans le tronc des modifications issues de la branche LMDZ-dev
comprises entre la révision 1074 et 1145
Validation: une simulation de 1 jour en séquentiel sur PC donne les mêmes
résultats entre la trunk et la dev
LF

Location:

LMDZ4/trunk

Files:

• . (modified) (1 prop)
• libf/phylmd/cv3_cine.F (modified) (12 diffs)

LMDZ4/trunk/libf/phylmd/cv3_cine.F

@@ -33 +33 @@
       integer ifst(nloc),isublcl(nloc)
       logical lswitch(nloc),lswitch1(nloc),lswitch2(nloc)
+      logical exist_lfc(nloc)
+      real plfc(nloc)
       real dpmax
       real deltap,dcin
       real buoylcl(nloc),tvplcl(nloc),tvlcl(nloc)
-      real plfc(nloc),p0(nloc)
+      real p0(nloc)
       real buoyz(nloc), buoy(nloc,nd)
 c
@@ -50 +52 @@
 c      Recompute buoyancies
 c--------------------------------------------------------------
-      DO k = 1,nl
+      DO k = 1,nd
         DO il = 1,ncum
+!      print*,'tvp tv=',tvp(il,k),tv(il,k)
           buoy(il,k) = tvp(il,k) - tv(il,k)
         ENDDO
       ENDDO
-c
-c---------------------------------------------------------------
-c premiere couche contenant un  niveau de flotabilite positive
-c et premiere couche contenant un  niveau de flotabilite negative
-c  au dessus du niveau de condensation
-c---------------------------------------------------------------
-      do il = 1,ncum
-        itop(il) =nl-1
-        ineg(il) = nl-1
-      enddo
-      do 100 k=nl,1,-1
-       do 110 il=1,ncum
-        if (k .ge. icb(il)) then
-         if (buoy(il,k) .gt. 0.) then
-          itop(il)=k
-         else
-          ineg(il)=k
-         endif
-        endif
-110    continue
-100   continue
-c      print *,' itop, ineg, icb ',itop(1),ineg(1), icb(1)
-c
 c---------------------------------------------------------------
 c
@@ -109 +89 @@
 c
 c---------------------------------------------------------------
+c premiere couche contenant un  niveau de flotabilite positive
+c et premiere couche contenant un  niveau de flotabilite negative
+c  au dessus du niveau de condensation
+c---------------------------------------------------------------
+      do il = 1,ncum
+        itop(il) =nl-1
+        ineg(il) = nl-1
+        exist_lfc(il) = .FALSE.
+      enddo
+      do 100 k=nl-1,1,-1
+       do 110 il=1,ncum
+        if (k .ge. ifst(il)) then
+         if (buoy(il,k) .gt. 0.) then
+          itop(il)=k
+          exist_lfc(il) = .TRUE.
+         else
+          ineg(il)=k
+         endif
+        endif
+110    continue
+100   continue
+c
+c---------------------------------------------------------------
+c When there is no positive buoyancy level, set Plfc, Cina and Cinb
+c to arbitrary extreme values.
+c---------------------------------------------------------------
+      DO il = 1,ncum
+       IF (.NOT.exist_lfc(il)) THEN
+         Plfc(il) = 1.111
+         Cinb(il) = -1111.
+         Cina(il) = -1112.
+       ENDIF
+      ENDDO
+c
+c
+c---------------------------------------------------------------
 c -- Two cases : BUOYlcl >= 0 and BUOYlcl < 0.
 c---------------------------------------------------------------
@@ -118 +134 @@
       DPMAX = 50.
       DO il = 1,ncum
-        lswitch1(il)=BUOYlcl(il) .GE. 0.
+        lswitch1(il)=BUOYlcl(il) .GE. 0. .AND. exist_lfc(il)
         lswitch(il) = lswitch1(il)
       ENDDO
@@ -233 +249 @@
 C
       DO il = 1,ncum
-        lswitch1(il)=BUOYlcl(il) .LT. 0.
+        lswitch1(il)=BUOYlcl(il) .LT. 0. .AND. exist_lfc(il)
         lswitch(il) = lswitch1(il)
       ENDDO
@@ -239 +255 @@
 c 2.0.1 Premiere  couche ou la flotabilite est negative au dessus du sol
 c ----------------------------------------------------
-c    au cas ou il existe  sinon ilow=1 (nk apres)
+c    au cas ou elle existe  sinon ilow=1 (nk apres)
 c      on suppose que la parcelle part de la premiere couche
 c
@@ -248 +264 @@
       ENDDO
 c
-      do 200 i=nl,1,-1
+      do 200 k=nl,1,-1
         DO il = 1,ncum
         IF (lswitch(il) .AND. k .LE.icb(il)-1) THEN
@@ -292 +308 @@
         dcin = RD*(BUOYz(il)+BUOYlcl(il))*deltap/(P0(il)+Plcl(il))
         CINB(il) = min(0.,dcin)
-cc        print *,'buoyz(il),buoylcl(il),deltap,p0(il),plcl(il),dcin ',
-cc     $           buoyz(il),buoylcl(il),deltap,p0(il),plcl(il),dcin
-      ENDIF
-      ENDDO
-c        print*, 'CINB ',CINB(1),'DCIN ',DCIN,I,BUOYz(1),BUOYlcl(1)
+      ENDIF
+      ENDDO
 c
       DO il = 1,ncum
@@ -316 +329 @@
       ENDDO
 c
-      IF (lswitch(1)) THEN
-c        print*,'ilow= ',ilow(1),'DCIN0 ',DCIN,P0(1),P(1,ilow(1))
-c        print*,'buoy',(BUOY(1,k),k=1,itop(1))
-      ENDIF
 c
 C  Middle part of CINB : integral from P(ilow) to P(isublcl)
@@ -332 +341 @@
         ENDIF
         ENDDO
-c      print*, 'CINB ', CINB(1), 'DCIN',DCIN,k,BUOY(1,k),BUOY(1,k+1)
       ENDDO
 c
@@ -345 +353 @@
       ENDDO
 C
-c        print*, ' CINB ', CINB(1), 'Dcin ',dcin
 c
 cc      ENDIF
@@ -439 +446 @@
       ENDDO
 cc      ENDIF
-c       Print *,' Plcl,P(itop-1),P(itop),PLFC,BUOYlcl'
-c     $      ,Plcl(1),P(1,itop(1)-1),P(1,itop(1)),PLFC(1),BUOYlcl(1)
-C
-c        print*, 'CIN above', CINA(1), 'CIN below',CINB(1)
 c