Changeset 2647 for trunk/LMDZ.GENERIC/libf/phystd/thermcell_flux.F90

Timestamp:

Mar 24, 2022, 4:02:39 PM (3 years ago)

Author:

dbardet

Message:

Thermal plume model: Fix an index error for calculating the maximal authorized interlayer mass flux of the plume.

File:

• trunk/LMDZ.GENERIC/libf/phystd/thermcell_flux.F90 (modified) (1 diff)

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

@@ -278 +278 @@
          DO ig=1,ngrid
             fff0 = fm(ig,l+1)
-            fmax = rhobarz(ig,l+1) * zw2(ig,l+1) * alpha_max
+            fmax = rhobarz(ig,l) * zw2(ig,l+1) * alpha_max
+            ! rhobarz is an arithmetic mean of the mid-layer rho, and thus evaluated
+            ! at the interlayer
+            ! so rhobarz(ig,l) can be taken to evaluate the inter-layer maximal
+            ! authorized mass flux fmax (with the plume vertical velocity zw2).
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 ! AB: The plume mass flux can be reduced.