Changeset 662

Timestamp:

May 21, 2012, 10:49:27 AM (13 years ago)

Author:

acolaitis

Message:

minor modif to pbl_param. Now that we compute wstar taking into account diffusion, we dont need to recompute it in pbl_param.

File:

• trunk/LMDZ.MARS/libf/phymars/pbl_parameters.F (modified) (5 diffs)

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

@@ -25 +25 @@
 !     pu(ngrid,nlay)   u component of the wind
 !     pv(ngrid,nlay)   v component of the wind
-!     wstar_in(ngrid)  free convection velocity in thermals
+!     wstar_in(ngrid)  free convection velocity in PBL
 !     hfmax(ngrid)     maximum vertical turbulent heat flux in thermals
 !     zmax(ngrid)      height reached by the thermals (pbl height)
@@ -40 +40 @@
 !     ustar(ngrid)     friction velocity
 !     tstar(ngrid)     friction temperature
-!     wstar(ngrid)     free convection velocity
 !     L_mo(ngrid)      monin_obukhov length
 !
@@ -70 +69 @@
       REAL Teta_out(ngrid,n_out)
       REAL, INTENT(OUT) :: ustar(ngrid), tstar(ngrid)
-      REAL wstar(ngrid)
       REAL, INTENT(OUT) :: L_mo(ngrid)
 
@@ -434 +432 @@
       ENDDO
 
-! Recompute wstar
-! We follow Spiga et. al 2010 (QJRMS)
-! ------------
-
-      DO ig=1, ngrid
-         IF (zmax(ig) .gt. 0.) THEN
-            wstar(ig)=(pg*zmax(ig)*hfmax(ig)/pbl_teta(ig))**(1./3.)
-         ELSE
-            wstar(ig)=0.
-         ENDIF
-      ENDDO
-
 !------------------------------------------------------------------------
 !------------------------------------------------------------------------
@@ -482 +468 @@
 
       vhf(:) = dvhf(:)*hfmax(:)
-      vvv(:) = dvvv(:)*(wstar(:))**2
+      vvv(:) = dvvv(:)*(wstar_in(:))**2
 
       ENDIF ! of if calltherm