Changeset 1834 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Dec 4, 2017, 6:35:57 PM (7 years ago)

Author:

mlefevre

Message:

MESOSCALE GENERIC. added turb_mod and changed turbdiff to send info on ustar to dynamical core in LES mode. by the way ustar is now calculated.

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• turb_mod.F90 (added)
• turbdiff.F90 (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phystd/turbdiff.F90

@@ -13 +13 @@
       use comcstfi_mod, only: rcp, g, r, cpp
       use callkeys_mod, only: water,tracer,nosurf
-
+#ifdef MESOSCALE
+      use turb_mod, only : ustar
+      use comm_wrf, only : comm_LATENT_HF
+#endif
       implicit none
 
@@ -259 +262 @@
       ENDDO
 !we treat only winds, energy and tracers coefficients will be computed with upadted winds
- 
 !JL12 calculate the flux coefficients (tables multiplied elements by elements)
       zfluxv(1:ngrid,1:nlay)=zkv(1:ngrid,1:nlay)*zb0(1:ngrid,1:nlay)
@@ -359 +361 @@
          zcdh(ig) = zcdh_true(ig)*sqrt(zu2)
          zkh(ig,1)= zcdh(ig)
-      ENDDO
+         ustar(ig)=sqrt(zcdv_true(ig))*sqrt(zu2)
+      ENDDO
+
 
 !     JL12 calculate the flux coefficients (tables multiplied elements by elements)
@@ -719 +723 @@
 
       if(water)then
+#ifndef MESOSCALE
          call writediagfi(ngrid,'beta','Dryness coefficient',' ',2,dryness)
+#endif
          if (tracer) then
+#ifndef MESOSCALE
             call writediagfi(ngrid,'dqevap','evaporated water vapor specific concentration','s-1',3,pdqevap)
+#else
+            comm_LATENT_HF(:)=0.0
+            comm_LATENT_HF(1:ngrid)=RLVTT*dqsdif_total(1:ngrid)/ptimestep
+#endif
          endif
       endif