Changeset 594 for trunk/LMDZ.GENERIC/libf/phystd/vdifc.F

Timestamp:

Mar 22, 2012, 8:25:03 PM (13 years ago)

Author:

jleconte

Message:

• New turbulent diffusion scheme solving "most" energy conservation problems:
  • Turbulent energy created by buoyancy effects is now dissipated back into enthalpy
  • the scheme is now written in an enthalpy conservative way
• Turbulent diffusion now treated in routine turbdiff (in F90).
• Temporarily, for comparison, the old vdifc can be used

by setting UseTurbDiff?=.false. in physiq.F90

• The sensible heat flux is now an output
• Corrected evaporation at the surface when all the surface water is evaporated.

File:

• trunk/LMDZ.GENERIC/libf/phystd/vdifc.F (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/vdifc.F

@@ -4 +4 @@
      &     pu,pv,ph,pq,ptsrf,pemis,pqsurf,
      &     pdufi,pdvfi,pdhfi,pdqfi,pfluxsrf,
-     &     pdudif,pdvdif,pdhdif,pdtsrf,pq2,
+     &     pdudif,pdvdif,pdhdif,pdtsrf,sensibFlux,pq2,
      &     pdqdif,pdqsdif,lastcall)
 
@@ -55 +55 @@
       REAL pfluxsrf(ngrid)
       REAL pdudif(ngrid,nlay),pdvdif(ngrid,nlay),pdhdif(ngrid,nlay)
-      REAL pdtsrf(ngrid),pcapcal(ngrid)
+      REAL pdtsrf(ngrid),sensibFlux(ngrid),pcapcal(ngrid)
       REAL pq2(ngrid,nlay+1)
       
@@ -674 +674 @@
          enddo
       enddo
-      
+
+      DO ig=1,ngrid  ! computing sensible heat flux (atm => surface)
+         sensibFlux(ig)=cpp*zb(ig,1)/ptimestep*(zh(ig,1)-ztsrf2(ig))
+      ENDDO      
+
       if (tracer) then
          do iq = 1, nq