Changeset 292 for trunk/LMDZ.MARS/libf

Timestamp:

Sep 9, 2011, 5:31:42 PM (13 years ago)

Author:

acolaitis

Message:

--- AC 07/09/2011 ---

libf/phymars/surflayer_interpol.F has been corrected. u* and t* are now computed at "gcm level", i.e. we consider that surface flux enhancement by subgrid gustiness has already been taken into account in vdifc. Hence, we do not use beta coefficients for u* and teta* in the interpolation. However, we still use it for the Richardson computation, to avoid problems with low winds conditions.

---

File:

• trunk/LMDZ.MARS/libf/phymars/surflayer_interpol.F (modified) (4 diffs)

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

@@ -134 +134 @@
           zu2(ig)=
 !     &  (MAX(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1),wmax(ig)**2))
-     &  ( pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) + wmax(ig)**2)
+     &  ( pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1))
 !     &  pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1)
 !        ELSE
@@ -144 +144 @@
      &      *sqrt(pz(ig,1)*pz0(ig))
      &      *(((log(pz(ig,1)/pz0(ig)))**2)/(log(pz(ig,1)/pz0t)))
-     &      *(ph(ig,1)-pts(ig))/zu2(ig)
+     &      *(ph(ig,1)-pts(ig))/(zu2(ig) + (0.5*wmax(ig))**2)
 
 !     &  /(MAX(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1),wmax(ig)**2))
@@ -186 +186 @@
 
        reynolds(ig)=karman*sqrt(fm(ig))
-     &              *sqrt(zu2(ig))
+     &              *sqrt(zu2(ig)+(0.5*wmax(ig))**2)
 !    &               *sqrt(pu(ig,1)**2 + pv(ig,1)**2)
      &              *pz0(ig)/(log(z1z0)*nu)
@@ -226 +226 @@
 !simpler definition of u* and teta*, equivalent to the formula above :
 
-            ustar(ig)=sqrt(pcdv(ig))*sqrt(zu2(ig))                   
-            tstar(ig)=-pcdh(ig)*(pts(ig)-ph(ig,1))/sqrt(pcdv(ig))
+            ustar(ig)=sqrt(pcdv(ig))*sqrt(zu2(ig))
+            tstar(ig)=-pcdh(ig)*(pts(ig)-ph(ig,1))
+     &     /sqrt(pcdv(ig))
 
            if (tstar(ig) .lt. -50) then