Changeset 300 for trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90

Timestamp:

Sep 21, 2011, 2:54:49 PM (13 years ago)

Author:

acolaitis

Message:

Revision on several settings for the thermals model. This version relies on fits done for an article to be published, and is more precise. See README for more infos.

File:

• trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90

@@ -1139 +1139 @@
          do l=1,lmax(ig)
               if(zlay(ig,l) .le. 0.7*zmax(ig)) then
-              fm_down(ig,l) =-1.8*fm(ig,l)
+              fm_down(ig,l) =-1.9*fm(ig,l)
               endif
 
              if(zlay(ig,l) .le. 0.06*zmax(ig)) then
           ztvd(ig,l)=ztv(ig,l)*max(0.,(1.+(sqrt((zlay(ig,l)/zmax(ig))/0.122449) - 1.)*(ztva(ig,l)/ztv(ig,l) - 1.)))
-             elseif(zlay(ig,l) .le. 0.6*zmax(ig)) then
-          ztvd(ig,l)=ztv(ig,l)*max(0.,1.-0.3*(ztva(ig,l)/ztv(ig,l) - 1.))
+             elseif(zlay(ig,l) .le. 0.4*zmax(ig)) then
+          ztvd(ig,l)=ztv(ig,l)*max(0.,1.-0.25*(ztva(ig,l)/ztv(ig,l) - 1.))
              elseif(zlay(ig,l) .le. 0.7*zmax(ig)) then
-          ztvd(ig,l)=ztv(ig,l)*max(0.,(1.+(((zlay(ig,l)/zmax(ig))-0.7)/0.333333)*(ztva(ig,l)/ztv(ig,l) - 1.)))
+          ztvd(ig,l)=ztv(ig,l)*max(0.,(1.+(((zlay(ig,l)/zmax(ig))-0.7)/1.)*(ztva(ig,l)/ztv(ig,l) - 1.)))
              else
           ztvd(ig,l)=ztv(ig,l)
@@ -1176 +1176 @@
          if(lmax(ig) .gt. 1) then
 ! No downdraft in the very-near surface layer, we begin at k=3
-         do k=3,lmax(ig)
+         do k=2,lmax(ig)
             zdthladj_down(ig,k)=(1./masse(ig,k))*(fm_down(ig,k+1)*ztv(ig,k+1)- &
      & fm_down(ig,k)*ztv(ig,k)+fm_down(ig,k)*ztvd(ig,k)-fm_down(ig,k+1)*ztvd(ig,k+1))