Changeset 300

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.

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/calltherm_mars.F90 (modified) (1 diff)
• libf/phymars/thermcell_main_mars.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -1011 +1011 @@
      which needed to be updated every time...
 --> This is completely transparent to the casual GCM users and only appears in the WRITEDIAGFI section of physiq.F
+
+== 21/09/11 == AC
+
+Revision on several settings for the thermals model. This version relies on fits done for an article to be published, and is more precise.
+
+ M             299   libf/phymars/thermcell_main_mars.F90
+ ^-----------------> Changed downdraft to updraft mass flux ratio from -1.8 to -1.9
+                     Changed first level for downdraft mixing from k=3 to k=2. Only level 1 is non-mixed now.
+                     Changed coefficients for downdraft to updraft thermal buoyancy ratios.
+
+ M             299   libf/phymars/calltherm_mars.F90
+ ^-----------------> Changed r_aspect_thermals from 2. to 1.5 for the GCM version to better match buoyancy profiles.

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

@@ -100 +100 @@
 #else
          nsplit_thermals=50
-         r_aspect_thermals=2.
+         r_aspect_thermals=1.5
 #endif
 

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))