Changeset 1026 for LMDZ4/trunk/libf/phylmd/thermcell_height.F90

Timestamp:

Oct 24, 2008, 5:09:39 PM (16 years ago)

Author:

lmdzadmin

Message:

Modifs thermiques
FH/IM

File:

• LMDZ4/trunk/libf/phylmd/thermcell_height.F90 (modified) (4 diffs)

LMDZ4/trunk/libf/phylmd/thermcell_height.F90

@@ -7 +7 @@
       IMPLICIT NONE
 #include "iniprint.h"
+#include "thermcell.h"
 
       INTEGER ig,l
@@ -23 +24 @@
       REAL zmax0(ngrid)
       REAL zmix(ngrid)
+      REAL num(ngrid)
+      REAL denom(ngrid)
 
       REAL zlevinter(ngrid)
@@ -70 +73 @@
          zlevinter(ig)=zlev(ig,1)
       enddo
+
+      if (iflag_thermals_ed.ge.1) then
+
+         num(:)=0.
+         denom(:)=0.
+         do ig=1,ngrid
+          do l=1,nlay
+             num(ig)=num(ig)+zw2(ig,l)*zlev(ig,l)*(zlev(ig,l+1)-zlev(ig,l))
+             denom(ig)=denom(ig)+zw2(ig,l)*(zlev(ig,l+1)-zlev(ig,l))
+          enddo
+       enddo
+       do ig=1,ngrid
+       if (denom(ig).gt.1.e-10) then
+          zmax(ig)=2.*num(ig)/denom(ig)
+          zmax0(ig)=zmax(ig)
+       endif 
+       enddo
+
+       else
+
       do  ig=1,ngrid
 ! calcul de zlevinter
@@ -80 +103 @@
        zmax0(ig)=zmax(ig)
       enddo
+
+
+      endif
+!endif iflag_thermals_ed
 !
 ! def de  zmix continu (profil parabolique des vitesses)