Changeset 1372

Timestamp:

Apr 27, 2010, 12:56:22 PM (15 years ago)

Author:

musat

Message:

Bug fix: overflow for pplay field
IM

Location:

LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd

Files:

• thermcellV0_main.F90 (modified) (2 diffs)
• thermcell_main.F90 (modified) (2 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcellV0_main.F90

@@ -593 +593 @@
       pcon(ig)=pplay(ig,1)*(zo(ig,1)/zqsat(ig,1))**CHI
       enddo
-      do k=1,nlay
+!IM   do k=1,nlay
+      do k=1,nlay-1
          do ig=1,ngrid
          if ((pcon(ig).le.pplay(ig,k))  &
@@ -600 +601 @@
          endif
          enddo
+      enddo
+!IM
+      do ig=1,ngrid
+        if (pcon(ig).le.pplay(ig,nlay)) then 
+           zcon2(ig)=zlay(ig,nlay)-(pcon(ig)-pplay(ig,nlay))/(RG*rho(ig,nlay))/100.
+           abort_message = 'thermcellV0_main: les thermiques vont trop haut '
+           CALL abort_gcm (modname,abort_message,1)
+        endif
       enddo
       if (prt_level.ge.1) print*,'14b OK convect8'

LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_main.F90

@@ -610 +610 @@
       pcon(ig)=pplay(ig,1)*(zo(ig,1)/zqsat(ig,1))**CHI
       enddo
-      do k=1,nlay
+!IM   do k=1,nlay
+      do k=1,nlay-1
          do ig=1,ngrid
          if ((pcon(ig).le.pplay(ig,k))  &
@@ -617 +618 @@
          endif
          enddo
+      enddo
+!IM
+      do ig=1,ngrid
+        if (pcon(ig).le.pplay(ig,nlay)) then 
+           zcon2(ig)=zlay(ig,nlay)-(pcon(ig)-pplay(ig,nlay))/(RG*rho(ig,nlay))/100.
+           abort_message = 'thermcellV0_main: les thermiques vont trop haut '
+           CALL abort_gcm (modname,abort_message,1)
+        endif
       enddo
       if (prt_level.ge.1) print*,'14b OK convect8'