Changeset 4022

Timestamp:

Jan 20, 2026, 1:07:53 AM (5 days ago)

Author:

aslmd

Message:

MESOSCALE: sensible flux computation is broken, has to fix it by putting the right zcdh back in physiq_mod again

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (1 diff)

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

@@ -3211 +3211 @@
      .         - capcal(1:ngrid,1)*zdtsdif(1:ngrid,1)
       else
-        sensibFlux(1:ngrid) =
-     &   (pplay(1:ngrid,1)/(r*pt(1:ngrid,1)))*cpp
-     &   *sqrt(pu(1:ngrid,1)*pu(1:ngrid,1)+pv(1:ngrid,1)*pv(1:ngrid,1)
-     &         +(log(1.+0.7*wstar(1:ngrid) + 2.3*wstar(1:ngrid)**2))**2)
-     &   *zcdh(1:ngrid,1)*(tsurf(1:ngrid,1)-zh(1:ngrid,1))
+        print*,"sensible flux computations broken! zcdh no longer here"
+!        sensibFlux(1:ngrid) =
+!     &   (pplay(1:ngrid,1)/(r*pt(1:ngrid,1)))*cpp
+!     &   *sqrt(pu(1:ngrid,1)*pu(1:ngrid,1)+pv(1:ngrid,1)*pv(1:ngrid,1)
+!     &         +(log(1.+0.7*wstar(1:ngrid) + 2.3*wstar(1:ngrid)**2))**2)
+!     &   *zcdh(1:ngrid,1)*(tsurf(1:ngrid,1)-zh(1:ngrid,1))
       endif