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

Timestamp:

Apr 6, 2010, 2:49:00 PM (14 years ago)

Author:

idelkadi

Message:

Nouvelle version des thermiques

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_main.F90 (modified) (6 diffs)

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

@@ -130 +130 @@
 
       real wmax(klon)
+      real wmax_tmp(klon)
       real wmax_sec(klon)
       real fm0(klon,klev+1),entr0(klon,klev),detr0(klon,klev)
@@ -191 +192 @@
 
 
-! #define wrgrads_thermcell
-#undef wrgrads_thermcell
+#define wrgrads_thermcell
 #ifdef wrgrads_thermcell
 ! Initialisation des sorties grads pour les thermiques.
@@ -208 +208 @@
 
       fm=0. ; entr=0. ; detr=0.
+
+      print*,'THERMCELL MAIN OPT7'
 
       icount=icount+1
@@ -395 +397 @@
 
       elseif (iflag_thermals_ed<=19) then
-! Version d'Arnaud Jam
-!         print*,'THERM RIO et al 2010'
+!        print*,'THERM RIO et al 2010, version d Arnaud'
          CALL thermcellV1_plume(itap,ngrid,nlay,ptimestep,ztv,zthl,po,zl,rhobarz,&
      &    zlev,pplev,pphi,zpspsk,alim_star,alim_star_tot,  &
@@ -426 +427 @@
       CALL thermcell_height(ngrid,nlay,lalim,lmin,linter,lmix,zw2,  &
      &           zlev,lmax,zmax,zmax0,zmix,wmax,lev_out)
+! Attention, w2 est transforme en sa racine carree dans cette routine
+! Le probleme vient du fait que linter et lmix sont souvent égaux à 1.
+      wmax_tmp=0.
+      do  l=1,nlay
+         wmax_tmp(:)=max(wmax_tmp(:),zw2(:,l))
+      enddo
+!     print*,"ZMAX ",lalim,lmin,linter,lmix,lmax,zmax,zmax0,zmix,wmax
 
 
@@ -665 +673 @@
          enddo
       enddo
-!calcul de ale_bl et alp_bl
-!pour le calcul d'une valeur integree entre la surface et lmax
-      do ig=1,ngrid
-      alp_int(ig)=0.
-      ale_int(ig)=0.
-      n_int(ig)=0
-      enddo
-!
-      do l=1,nlay
-      do ig=1,ngrid
-       if(l.LE.lmax(ig)) THEN
-        alp_int(ig)=alp_int(ig)+0.5*rhobarz(ig,l)*wth3(ig,l)
-        ale_int(ig)=ale_int(ig)+0.5*zw2(ig,l)**2
-        n_int(ig)=n_int(ig)+1
-       endif
-      enddo
-      enddo
+!
 !      print*,'avant calcul ale et alp' 
 !calcul de ALE et ALP pour la convection
+      Alp_bl(:)=0.
+      Ale_bl(:)=0.
+!          print*,'ALE,ALP ,l,zw2(ig,l),Ale_bl(ig),Alp_bl(ig)'
+      do l=1,nlay
       do ig=1,ngrid
-!      Alp_bl(ig)=0.5*rhobarz(ig,lmix_bis(ig))*wth3(ig,lmix(ig))
-!          Alp_bl(ig)=0.5*rhobarz(ig,nivcon(ig))*wth3(ig,nivcon(ig))
-!      Alp_bl(ig)=0.5*rhobarz(ig,lmix(ig))*wth3(ig,lmix(ig)) 
-!     &           *0.1
-!valeur integree de alp_bl * 0.5:
-       if (n_int(ig).gt.0) then
-       Alp_bl(ig)=0.5*alp_int(ig)/n_int(ig)
-!       if (Alp_bl(ig).lt.0.) then
-!       Alp_bl(ig)=0.
-       endif
-!       endif
-!         write(18,*),'rhobarz,wth3,Alp',rhobarz(ig,nivcon(ig)),
-!     s               wth3(ig,nivcon(ig)),Alp_bl(ig)
-!            write(18,*),'ALP_BL',Alp_bl(ig),lmix(ig)
-!      Ale_bl(ig)=0.5*zw2(ig,lmix_bis(ig))**2
-!      if (nivcon(ig).eq.1) then
-!       Ale_bl(ig)=0.
-!       else
-!valeur max de ale_bl:
-       Ale_bl(ig)=0.5*zw2(ig,lmix(ig))**2 
-!     & /2.
-!     & *0.1
-!        Ale_bl(ig)=0.5*zw2(ig,lmix_bis(ig))**2 
-!       if (n_int(ig).gt.0) then
-!       Ale_bl(ig)=ale_int(ig)/n_int(ig)
-!        Ale_bl(ig)=4.
-!       endif
-!       endif
-!            Ale_bl(ig)=0.5*wth2(ig,lmix_bis(ig))
-!          Ale_bl(ig)=wth2(ig,nivcon(ig))
-!          write(19,*),'wth2,ALE_BL',wth2(ig,nivcon(ig)),Ale_bl(ig)
-      enddo
+           Alp_bl(ig)=max(Alp_bl(ig),0.5*rhobarz(ig,l)*wth3(ig,l) )
+           Ale_bl(ig)=max(Ale_bl(ig),0.5*zw2(ig,l)**2)
+!          print*,'ALE,ALP',l,zw2(ig,l),Ale_bl(ig),Alp_bl(ig)
+      enddo
+      enddo
+
+!     print*,'AAAAAAA ',Alp_bl,Ale_bl,lmix
+
+
+! TEST. IL FAUT REECRIRE LES ALE et ALP
+!     Ale_bl(:)=0.5*wmax(:)*wmax(:)
+!     Alp_bl(:)=0.1*wmax(:)*wmax(:)*wmax(:)
+
 !test:calcul de la ponderation des couches pour KE
 !initialisations