Changeset 496 for trunk/LMDZ.MARS/libf/phymars/physiq.F

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/physiq.F (modified) (9 diffs)

trunk/LMDZ.MARS/libf/phymars/physiq.F

@@ -333 +333 @@
       REAL lmax_th_out(ngridmx),zmax_th(ngridmx)
       REAL, SAVE :: wmax_th(ngridmx)
-      REAL hfmax_th(ngridmx)
+      REAL, SAVE :: hfmax_th(ngridmx)
       REAL pdu_th(ngridmx,nlayermx),pdv_th(ngridmx,nlayermx)
       REAL pdt_th(ngridmx,nlayermx),pdq_th(ngridmx,nlayermx,nqmx)
@@ -342 +342 @@
       REAL z_out                          ! height of interpolation between z0 and z1 [meters]
       REAL ustar(ngridmx),tstar(ngridmx)  ! friction velocity and friction potential temp
-      REAL L_mo(ngridmx)
+      REAL L_mo(ngridmx),wstar(ngridmx),vhf(ngridmx),vvv(ngridmx)
       REAL zu2(ngridmx)
 c=======================================================================
@@ -569 +569 @@
                  enddo
               endif
+           else
+              zdtnlte(:,:)=0.
            endif
 
@@ -746 +748 @@
              IF (zh(ig,1) .lt. tsurf(ig)) THEN
                wmax_th(ig)=1.
-             ENDIF        
+               hfmax_th(ig)=0.2
+             ELSE
+               wmax_th(ig)=0.
+               hfmax_th(ig)=0.
+             ENDIF      
           ENDDO
         ENDIF
@@ -762 +768 @@
      $        zdum1,zdum2,zdh,pdq,zflubid,
      $        zdudif,zdvdif,zdhdif,zdtsdif,q2,
-     &        zdqdif,zdqsdif,wmax_th,zcdv,zcdh)
+     &        zdqdif,zdqsdif,wmax_th,zcdv,zcdh,hfmax_th)
 
 #ifdef MESOSCALE
@@ -853 +859 @@
         lmax_th(:)=0
         wmax_th(:)=0.
+        hfmax_th(:)=0.
         lmax_th_out(:)=0.
         end if
@@ -1803 +1810 @@
          z_out=0.
          if (calltherm .and. (z_out .gt. 0.)) then
-         call surflayer_interpol(ngrid,nlayer,z0,g,zzlay,zu,zv,wmax_th
-     &              ,tsurf,zh,z_out,Teta_out,u_out,ustar,tstar,L_mo)
-
-         zu2(:)=sqrt(zu(:,1)*zu(:,1)+zv(:,1)*zv(:,1))
-         call WRITEDIAGFI(ngridmx,'sqrt(zu2)',
-     &              'horizontal velocity norm','m/s',
-     &                         2,zu2)
-
-         call WRITEDIAGFI(ngridmx,'Teta_out',
-     &              'potential temperature at z_out','K',
-     &                         2,Teta_out)
-         call WRITEDIAGFI(ngridmx,'u_out',
-     &              'horizontal velocity norm at z_out','m/s',
-     &                         2,u_out)
-         call WRITEDIAGFI(ngridmx,'u*',
-     &              'friction velocity','m/s',
-     &                         2,ustar)
-         call WRITEDIAGFI(ngridmx,'teta*',
-     &              'friction potential temperature','K',
-     &                         2,tstar)
-         call WRITEDIAGFI(ngrid,'L',
-     &              'Monin Obukhov length','m',
-     &                         2,L_mo)
+
+         call pbl_parameters(ngrid,nlayer,z0,
+     & g,zzlay,zu,zv,wmax_th,hfmax_th,zmax_th,tsurf,zh,z_out,
+     & Teta_out,u_out,ustar,tstar,wstar,L_mo,vhf,vvv)
+
          else
            if((.not. calltherm).and.(z_out .gt. 0.)) then
@@ -1914 +1903 @@
          z_out=0.
          if (calltherm .and. (z_out .gt. 0.)) then
-         call surflayer_interpol(ngrid,nlayer,z0,g,zzlay,zu,zv,wmax_th
-     &              ,tsurf,zh,z_out,Teta_out,u_out,ustar,tstar,L_mo)
-
-         zu2(:)=sqrt(zu(:,1)*zu(:,1)+zv(:,1)*zv(:,1))
-         call WRITEDIAGFI(ngridmx,'sqrt(zu2)',
-     &              'horizontal velocity norm','m/s',
-     &                         0,zu2)
-
-         call WRITEDIAGFI(ngridmx,'Teta_out',
-     &              'potential temperature at z_out','K',
-     &                         0,Teta_out)
-         call WRITEDIAGFI(ngridmx,'u_out',
-     &              'horizontal velocity norm at z_out','m/s',
-     &                         0,u_out)
-         call WRITEDIAGFI(ngridmx,'u*',
-     &              'friction velocity','m/s',
-     &                         0,ustar)
-         call WRITEDIAGFI(ngridmx,'teta*',
-     &              'friction potential temperature','K',
-     &                         0,tstar)
+
+         call pbl_parameters(ngrid,nlayer,z0,
+     & g,zzlay,zu,zv,wmax_th,hfmax_th,zmax_th,tsurf,zh,z_out,
+     & Teta_out,u_out,ustar,tstar,wstar,L_mo,vhf,vvv)
+
          else
            if((.not. calltherm).and.(z_out .gt. 0.)) then
@@ -1947 +1921 @@
      &              'hauteur du thermique','point',
      &                         0,lmax_th_out)
+        call WRITEDIAGFI(ngridmx,'zmax_th',
+     &              'hauteur du thermique','m',
+     &                         0,zmax_th)
         call WRITEDIAGFI(ngridmx,'hfmax_th',
      &              'maximum TH heat flux','K.m/s',