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

Timestamp:

Aug 8, 2011, 2:35:47 PM (13 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)

---

File:

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

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))