Changeset 319 for trunk/LMDZ.MARS/libf

Timestamp:

Oct 14, 2011, 11:59:58 AM (13 years ago)

Author:

aslmd

Message:

LMDZ.MARS : minor changes to thermals [output L_mo + default settings]. transparent to users.

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• calltherm_interface.F90 (modified) (3 diffs)
• physiq.F (modified) (3 diffs)
• surflayer_interpol.F (modified) (4 diffs)

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

@@ -8 +8 @@
      & pdu_th,pdv_th,pdt_th,pdq_th,lmax_th,zmax_th,pbl_dtke,hfmax,wmax)
 
-       USE ioipsl_getincom
+       !USE ioipsl_getincom
 
       implicit none
@@ -107 +107 @@
 
             pq_therm(:,:,:)=0.
-            qtransport_thermals=.true.
-            call getin("qtransport_thermals",qtransport_thermals)
+            qtransport_thermals=.true. !! default setting
+            !call getin("qtransport_thermals",qtransport_thermals)
             if(qtransport_thermals) then
             if(tracer) then
@@ -124 +124 @@
             buoyancyEst(:,:)=0.
 
-       dtke_thermals=.false.
-       call getin("dtke_thermals",dtke_thermals)
+       dtke_thermals=.false. !! default setting
+       !call getin("dtke_thermals",dtke_thermals)
          if(dtke_thermals) then
 

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

@@ -327 +327 @@
       REAL zcdv(ngridmx), zcdh(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 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 zu2(ngridmx)
 c=======================================================================
@@ -1658 +1659 @@
 
 
-         z_out=0.
+         z_out=1. !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)
+     &              ,tsurf,zh,z_out,Teta_out,u_out,ustar,tstar,L_mo)
 
          zu2(:)=sqrt(zu(:,1)*zu(:,1)+zv(:,1)*zv(:,1))
@@ -1680 +1681 @@
      &              'friction potential temperature','K',
      &                         2,tstar)
+         call WRITEDIAGFI(ngrid,'L',
+     &              'Monin Obukhov length','m',
+     &                         2,L_mo)
          else
            if((.not. calltherm).and.(z_out .gt. 0.)) then

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

@@ -1 +1 @@
       SUBROUTINE surflayer_interpol(ngrid,nlay,pz0, 
-     & pg,pz,pu,pv,wmax,pts,ph,z_out,Teta_out,u_out,ustar,tstar)
+     & pg,pz,pu,pv,wmax,pts,ph,z_out,Teta_out,u_out,ustar,tstar,L_mo)
       IMPLICIT NONE
 !=======================================================================
@@ -45 +45 @@
       REAL, INTENT(OUT) :: Teta_out(ngrid),u_out(ngrid)! interpolated fields at z_out : potential temperature and norm(uv)
       REAL, INTENT(OUT) :: ustar(ngrid), tstar(ngrid) ! u* and teta*
+      REAL, INTENT(OUT) :: L_mo(ngrid)                ! Monin-Obukhov length
 
 !   Local:
@@ -81 +82 @@
 
       REAL zu2(ngrid)                  ! Large-scale wind at first layer
-      REAL L_mo(ngrid)                ! Monin-Obukhov length
 !-----------------------------------------------------------------------
 !   couche de surface:
@@ -331 +331 @@
 ! Usefull diagnostics for the interpolation routine :
 
-!         call WRITEDIAGFI(ngrid,'L',
-!     &              'Monin Obukhov length','m',
-!     &                         2,L_mo)
 !         call WRITEDIAGFI(ngrid,'z0T',
 !     &              'thermal roughness length','m',