Changeset 852 for LMDZ4/trunk/libf/phytherm/thermcell_plume.F90

Timestamp:

Oct 11, 2007, 3:43:42 PM (17 years ago)

Author:

Laurent Fairhead

Message:

Mise a jour de la physique avec thermiques avec la version de FH d'aout 2007
LF

File:

• LMDZ4/trunk/libf/phytherm/thermcell_plume.F90 (modified) (5 diffs)

LMDZ4/trunk/libf/phytherm/thermcell_plume.F90

@@ -1 +1 @@
       SUBROUTINE thermcell_plume(ngrid,klev,ztv,zthl,po,zl,rhobarz,  &
      &           zlev,pplev,pphi,zpspsk,l_mix,r_aspect,alim_star,  &
-     &           lentr,zmax_sec,f0,detr_star,entr_star,f_star,ztva,  &
+     &           lalim,zmax_sec,f0,detr_star,entr_star,f_star,ztva,  &
      &           ztla,zqla,zqta,zha,zw2,zqsatth,lmix,linter,lev_out)
 
@@ -29 +29 @@
       REAL l_mix
       REAL r_aspect
-      INTEGER lentr(ngrid)
+      INTEGER lalim(ngrid)
       integer lev_out                           ! niveau pour les print
 
@@ -133 +133 @@
 !tests sur la definition du detr
 !calcul de detr_star et entr_star
-                w_est(ig,3)=zw2(ig,2)*  &
-     &                   ((f_star(ig,2))**2)  &
-     &                   /(f_star(ig,2)+alim_star(ig,2))**2+  &
-     &                   2.*RG*(ztva(ig,1)-ztv(ig,2))/ztv(ig,2)  &
-     &                   *(zlev(ig,3)-zlev(ig,2))
+
+
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! FH le test miraculeux de Catherine ? Le bout du tunel ?
+!               w_est(ig,3)=zw2(ig,2)*  &
+!    &                   ((f_star(ig,2))**2)  &
+!    &                   /(f_star(ig,2)+alim_star(ig,2))**2+  &
+!    &                   2.*RG*(ztva(ig,1)-ztv(ig,2))/ztv(ig,2)  &
+!    &                   *(zlev(ig,3)-zlev(ig,2))
+!     if (l.gt.2) then
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
       if (l.gt.2) then
+          w_est(ig,3)=zw2(ig,2)* &
+     &      ((f_star(ig,2))**2) &
+     &      /(f_star(ig,2)+alim_star(ig,2))**2+ &
+     &      2.*RG*(ztva(ig,2)-ztv(ig,2))/ztv(ig,2) &
+     &      *(zlev(ig,3)-zlev(ig,2))
+
+
 !---------------------------------------------------------------------------
 !calcul de l entrainement et du detrainement lateral
@@ -217 +232 @@
 !calcul de entr_star
 !
-        if (detr_star(ig,l).gt.f_star(ig,l)) then
-           detr_star(ig,l)=f_star(ig,l)
-!a decommenter ou pas?
-!           entr_star(ig,l)=0.
-        endif
-
 ! Deplacement du calcul de entr_star pour eviter d'avoir aussi
 ! entr_star > fstar.
+! Redeplacer suite a la transformation du cas detr>f
 ! FH
-          entr_star(ig,l)=0.4*detr_star(ig,l)
-!
+        if(l.gt.lalim(ig)) then
+         entr_star(ig,l)=0.4*detr_star(ig,l)
+        else
+
+! FH :
+! Cette ligne doit permettre de garantir qu'on a toujours un flux = 1
+! en haut de la couche d'alimentation.
+! A remettre en questoin a la premiere occasion mais ca peut aider a 
+! ecrire un code robuste.
+! Que ce soit avec ca ou avec l'ancienne facon de faire (e* = 0 mais
+! d* non nul) on a une discontinuité de e* ou d* en haut de la couche
+! d'alimentation, ce qui n'est pas forcement heureux.
+         entr_star(ig,l)=max(detr_star(ig,l)-alim_star(ig,l),0.)
+         detr_star(ig,l)=entr_star(ig,l)
+        endif
+
+!
+        if (detr_star(ig,l).gt.f_star(ig,l)) then
+
+!  Ce test est là entre autres parce qu'on passe par des valeurs
+!  delirantes de detr_star.
+!  ca vaut sans doute le coup de verifier pourquoi.
+
+           detr_star(ig,l)=f_star(ig,l)
+           if (l.gt.lalim(ig)+1) then
+               entr_star(ig,l)=0.
+               alim_star(ig,l)=0.
+! FH ajout pour forcer a stoper le thermique juste sous le sommet
+! de la couche (voir calcul de finter)
+               zw2(ig,l+1)=-1.e-10
+               linter(ig)=l+1
+            else
+               entr_star(ig,l)=detr_star(ig,l)
+            endif
+        endif
+
       else
          detr_star(ig,l)=0.
          entr_star(ig,l)=0.
       endif
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! FH inutile si on conserve comme on l'a fait plus haut entr=detr
+! dans la couche d'alimentation
 !pas d entrainement dans la couche alim
-      if ((l.lt.lentr(ig))) then
-           entr_star(ig,l)=0.
-      endif
+!      if ((l.le.lalim(ig))) then
+!           entr_star(ig,l)=0.
+!      endif
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !
 !prise en compte du detrainement et de l entrainement dans le calcul du flux
@@ -308 +357 @@
 !
 !initialisations pour le calcul de la hauteur du thermique, de l'inversion et de la vitesse verticale max 
+
+            if (zw2(ig,l+1)>0. .and. zw2(ig,l+1).lt.1.e-10) then
+!               stop'On tombe sur le cas particulier de thermcell_dry'
+                print*,'On tombe sur le cas particulier de thermcell_plume'
+                zw2(ig,l+1)=0.
+                linter(ig)=l+1
+            endif
+
+
         if (zw2(ig,l+1).lt.0.) then 
            linter(ig)=(l*(zw2(ig,l+1)-zw2(ig,l))  &
      &               -zw2(ig,l))/(zw2(ig,l+1)-zw2(ig,l))
            zw2(ig,l+1)=0.
-        else
+        endif
+
            wa_moy(ig,l+1)=sqrt(zw2(ig,l+1)) 
-        endif
+
         if (wa_moy(ig,l+1).gt.wmaxa(ig)) then
 !   lmix est le niveau de la couche ou w (wa_moy) est maximum