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

Timestamp:

Jan 11, 2012, 2:21:08 PM (13 years ago)

Author:

acolaitis

Message:

M 489 libf/phymars/thermcell_main_mars.F90
^{--------------------> New parameters for thermals, following latest version of the relevant article}

M 489 libf/phymars/physiq.F
^{--------------------> Minor modifications}

D 489 libf/phymars/surflayer_interpol.F
A 0 libf/phymars/pbl_parameters.F
^{--------------------> Replaced surflayer_interpol.F by a cleaner and richer version called pbl_parameters.F}

This routine is the base of what will later be implemented in the MCD and as a tool to
compute surface properties from output fields (so that it will ultimately disappear from libf
and might become an utility)

M 489 libf/phymars/vdif_cd.F
^{--------------------> Minor modification to the Richardson computation}

M 489 libf/phymars/vdifc.F
^{--------------------> Minor modification to the aerodynamic conductances computation, replacing wmax by hfmax, which}

is a better proxy for convective activity. Might be replaced by w_star later.

File:

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

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

@@ -111 +111 @@
       REAL zdz,zbuoy(ngridmx,nlayermx),zw2m
       LOGICAL active(ngridmx),activetmp(ngridmx)
-      REAL a1,b1,ae,be,ad,bd
+      REAL a1,b1,ae,be,ad,bd,fdfu
       INTEGER tic
 
@@ -304 +304 @@
       linter(:)=1.
 
+! --------------------------------------------------------------------------
+! -------------- MAIN PARAMETERS FOR THERMALS MODEL ------------------------
+! --------------  see thermiques.pro and getfit.py -------------------------
+
 !      a1=2.5 ; b1=0.0015 ; ae=0.045 ; be = 0.6  ! svn baseline
 
-      a1=1.60226 ; b1=0.0006 ; ae=0.0454 ; be = 0.57 !improved fits
-      ad = 0.0005114  ; bd = -0.662
+! Using broad downdraft selection
+!      a1=1.60226 ; b1=0.0006 ; ae=0.0454 ; be = 0.57
+!      ad = 0.0005114  ; bd = -0.662
+!      fdfu = -1.9
+
+! Using conditional sampling downdraft selection
+      a1=1.4716 ; b1=0.0005698 ; ae=0.03683 ; be = 0.57421
+      ad = 0.00048088  ; bd = -0.6697
+      fdfu = -1.3
+
+! --------------------------------------------------------------------------
+! --------------------------------------------------------------------------
+! --------------------------------------------------------------------------
 
 ! Initialisation des variables entieres
@@ -513 +528 @@
 !---------------------------------------------------------------------------
 
-      DO tic=0,6  ! internal convergence loop
+      DO tic=0,1  ! internal convergence loop
       activetmp(:)=active(:) .and. f_star(:,l+1)>1.e-10
       do ig=1,ngridmx
@@ -640 +655 @@
 ! ================= FIN PLUME ===============================================
 ! ===========================================================================
-
 
 ! ===========================================================================
@@ -1089 +1103 @@
          do l=1,lmax(ig)
               if(zlay(ig,l) .le. 0.7*zmax(ig)) then
-              fm_down(ig,l) =-1.9*fm(ig,l)
+              fm_down(ig,l) =fm(ig,l)* &
+     &      max(fdfu,-3*max(0.,(zlay(ig,l)/zmax(ig)))-0.9)
               endif
 
@@ -1102 +1117 @@
             endif
 
-!             if(zlay(ig,l) .le. 0.06*zmax(ig)) then
-!          ztvd(ig,l)=ztv(ig,l)*max(0.,(zlay(ig,l)/zmax(ig))/19.231 + 0.9938)
-!             elseif(zlay(ig,l) .le. 0.6*zmax(ig)) then
-!          ztvd(ig,l)=ztv(ig,l)*max(0.,(zlay(ig,l)/zmax(ig)-0.075)/187.931 + 0.9982)
-!             else
-!!             elseif(zlay(ig,l) .le. 0.7*zmax(ig)) then
-!          ztvd(ig,l)=ztv(ig,l)*max(0.,(zlay(ig,l)/zmax(ig) -0.60)/(-1333) + 1.00025)
-!!             else
-!!          ztvd(ig,l)=ztv(ig,l)
-!            endif
 
          enddo
@@ -1126 +1131 @@
          if(lmax(ig) .gt. 1) then
 ! No downdraft in the very-near surface layer, we begin at k=3
-         do k=2,lmax(ig)
+ 
+         do k=3,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))
@@ -1332 +1338 @@
          do ig=1,ngridmx
            pdtadj(ig,l)=(zdthladj(ig,l)+zdthladj_down(ig,l))*zpopsk(ig,l)
+!           pdtadj(ig,l)=zdthladj(ig,l)*zpopsk(ig,l)
          enddo
       enddo