Changeset 499 for trunk/LMDZ.MARS/libf/phymars/vdif_cd.F

Timestamp:

Jan 16, 2012, 10:16:46 AM (13 years ago)

Author:

acolaitis

Message:

M 498 libf/phymars/calltherm_interface.F90
^{------------------> Added computation of wstar. Wmax is no longer an output.}

M 498 libf/phymars/physiq.F
^{------------------> Minor changes.}

M 498 libf/phymars/vdif_cd.F
^{------------------> Changed computation of bulk Richardson number to allow more unstable configurations.}

M 498 libf/phymars/vdifc.F
^{------------------> Changed computation of surface conductance coefficients following a new approach. Comparisons with LES are much better.}

File:

• trunk/LMDZ.MARS/libf/phymars/vdif_cd.F (modified) (3 diffs)

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

@@ -1 +1 @@
       SUBROUTINE vdif_cd(ngrid,nlay,pz0,
-     &           pg,pz,pu,pv,wmax,pts,ph,pcdv,pcdh)
+     &           pg,pz,pu,pv,wstar,pts,ph,pcdv,pcdh)
       IMPLICIT NONE
 c=======================================================================
@@ -46 +46 @@
       REAL, INTENT(IN) :: pu(ngrid,nlay),pv(ngrid,nlay)
       REAL, INTENT(IN) :: pts(ngrid),ph(ngrid,nlay)
-      REAL, INTENT(IN) :: wmax(ngrid)
+      REAL, INTENT(IN) :: wstar(ngrid)
       REAL, INTENT(OUT) :: pcdv(ngrid),pcdh(ngrid) ! momentum and heat drag coefficient
 
@@ -146 +146 @@
 c Richardson number formulation proposed by D.E. England et al. (1995)
 
-!         zu2=MAX(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1),0.25*wmax(ig)**2)
+!         zu2=MAX(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1),0.25*wstar(ig)**2)
 !         zu2=pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1)
-         zu2(ig)=MAX(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1),wmax(ig)**2)
-!       zu2(ig)=pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) + (0.5*wmax(ig))**2
-
-               ! we add the wmax to simulate
+         zu2(ig)=MAX(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1),             &
+     &      (0.3*wstar(ig))**2)
+!       zu2(ig)=pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) + (0.5*wstar(ig))**2
+
+               ! we add the wstar to simulate
                ! bulk Ri changes due to subgrid wind feeding the thermals