Changeset 267 for trunk

Timestamp:

Aug 8, 2011, 2:35:47 PM (14 years ago)

Author:

acolaitis

Message:

--- AC 03/08/2011 ---
M 265 libf/phymars/physiq.F
<> Added a PBL section for outputs, with a call to SL outputs via surflayer_interpol and thermals outputs

A 0 libf/phymars/surflayer_interpol.F
<> New subroutine to interpolate horizontal velocity norm and potential temperature in the surface layer.

THIS ROUTINE IS NOT VALIDATED YET. IT IS TURNED OFF BY DEFAULT AND IS HERE FOR DEVELOPMENT PURPOSES ONLY FOR NOW.

M 265 libf/phymars/vdif_cd.F
<> Important modification to the Reynolds formula : due to a confusion of symbols used in the litterature, a wrong simplification of

numerical values had been implicitely done... It is now corrected and yields better results for the mixed layer temperature.

M 265 libf/phymars/vdifc.F
<> Cosmetic modif (added comment for more clarity)

---

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/physiq.F (modified) (4 diffs)
• libf/phymars/surflayer_interpol.F (added)
• libf/phymars/vdif_cd.F (modified) (2 diffs)
• libf/phymars/vdifc.F (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -847 +847 @@
 - vdifc.F Now takes into account sub-grid gustiness, evaluated from thermals activity (it's proxy being the maximum vertical velocity)
 - Minor modification in meso_inc_les.F: u* is now taken from the new vdifc and not recomputed from a simple law
+== 08/08/2011 == AC
+- Formula correction for the Reynolds in vdif_cd : will affect the thermal roughness length. Results obtained with the precedent formula (previous commit) probably underestimate by about 1K the temperature in the mixing layer.
+- Added a new subroutine to libf/phymars called surflayer_interpol.F. This routine can be called as a diagnostic for interpolation of horizontal velocity norm and temperature in the surface layer, in 1D or 3D. The output height is controlled by a parameter z_out, to be set in physiq.F. This might be later added to a .def file. z_out is set to 0. by default, resulting in no call to this subroutine.

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

@@ -318 +318 @@
       REAL time_phys
 
-c Variables from thermal
+c Variables for PBL
 
       REAL lmax_th_out(ngridmx),zmax_th(ngridmx)
       REAL, SAVE :: wmax_th(ngridmx)
-      REAL, SAVE :: hfmax_th(ngridmx)  !! AS: pourquoi en save ?
+      REAL hfmax_th(ngridmx)
       REAL pdu_th(ngridmx,nlayermx),pdv_th(ngridmx,nlayermx)
       REAL pdt_th(ngridmx,nlayermx),pdq_th(ngridmx,nlayermx,nqmx)
@@ -328 +328 @@
       REAL dtke_th(ngridmx,nlayermx+1)
       REAL zcdv(ngridmx)
+      REAL Teta_out(ngridmx),u_out(ngridmx)  ! Interpolated teta and u at z_out
+      REAL z_out                          ! height of interpolation between z0 and z1
+      REAL ustar(ngridmx),tstar(ngridmx)  ! friction velocity and friction potential temp
+      REAL zu2(ngridmx)
 c=======================================================================
 
@@ -1553 +1557 @@
 
 c        ----------------------------------------------------------
+c        ----------------------------------------------------------
+c        PBL OUTPUS
+c        ----------------------------------------------------------
+c        ----------------------------------------------------------
+
+
+c        ----------------------------------------------------------
+c        Outputs of surface layer
+c        ----------------------------------------------------------
+
+
+         z_out=0.
+         if (calltherm .and. (z_out .gt. 0.)) then
+         call surflayer_interpol(ngrid,nlayer,z0,g,zzlay,zu,zv,wmax_th
+     &              ,tsurf,zt(:,:)*(zplay(:,:)/zplev(:,:))**rcp
+     &              ,z_out,Teta_out,u_out,ustar,tstar)
+
+         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)
+         else
+           if((.not. calltherm).and.(z_out .gt. 0.)) then
+            print*, 'WARNING : no interpolation in surface-layer :'
+            print*, 'Outputing surface-layer quantities without thermals
+     & does not make sense'
+           endif
+         endif
+
+c        ----------------------------------------------------------
 c        Outputs of thermals
 c        ----------------------------------------------------------
@@ -1589 +1636 @@
 
 c        ----------------------------------------------------------
+c        ----------------------------------------------------------
+c        END OF PBL OUTPUS
+c        ----------------------------------------------------------
+c        ----------------------------------------------------------
+
+
+c        ----------------------------------------------------------
 c        Output in netcdf file "diagsoil.nc" for subterranean
 c          variables (output every "ecritphy", as for writediagfi)

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

@@ -52 +52 @@
       INTEGER ig
 
-      REAL karman
+      REAL karman,nu    ! Von Karman constant and fluid kinematic viscosity
       LOGICAL firstcal
-      DATA karman/.41/
+      DATA karman,nu/.41,0.001/
       DATA firstcal/.true./
-      SAVE karman
+      SAVE karman,nu
 
 c    Local(2):
@@ -179 +179 @@
       endif
 
-       reynolds(ig)=sqrt(fm(ig))*sqrt(pu(ig,1)**2 + pv(ig,1)**2)*pz0(ig)
-     &                   /(log(z1z0))
+       reynolds(ig)=karman*sqrt(fm(ig))*sqrt(pu(ig,1)**2 + pv(ig,1)**2)
+     &       *pz0(ig)/(log(z1z0)*nu)
        pz0tcomp(ig)=pz0(ig)*exp(-karman*7.3*
      &              (reynolds(ig)**0.25)*(prandtl(ig)**0.5))

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

@@ -238 +238 @@
 
         zu2=pu(ig,1)*pu(ig,1)+pv(ig,1)*pv(ig,1)
-     &             +(1.2*wmax(ig))**2
-
+     &             +(1.2*wmax(ig))**2        !subgrid gustiness is used to enhance surface flux only, and not u*,t* computations
+        
         zcdv(ig)=zcdv_true(ig)*sqrt(zu2)
         zcdh(ig)=zcdh_true(ig)*sqrt(zu2)