Changeset 1236 for LMDZ4/trunk/libf

Timestamp:

Aug 31, 2009, 2:55:42 PM (15 years ago)

Author:

musat

Message:

Petit bug tke: evolution tke uniquement si iflag_pbl>1
IM

File:

• LMDZ4/trunk/libf/phylmd/pbl_surface_mod.F90 (modified) (3 diffs)

LMDZ4/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -1278 +1278 @@
     INCLUDE "dimsoil.h"
     INCLUDE "clesphys.h"
+    INCLUDE "compbl.h"
 
 ! Input variables
@@ -1336 +1337 @@
                 u10m(i,nsrf)  = u10m(i,nsrf_comp1)
                 v10m(i,nsrf)  = v10m(i,nsrf_comp1)
-                tke(i,:,nsrf) = tke(i,:,nsrf_comp1)
+                if (iflag_pbl > 1) then
+                 tke(i,:,nsrf) = tke(i,:,nsrf_comp1)
+                endif
                 mfois(nsrf) = mfois(nsrf) + 1
              ELSE
@@ -1348 +1351 @@
                 u10m(i,nsrf)  = u10m(i,nsrf_comp2) *pctsrf_old(i,nsrf_comp2) + u10m(i,nsrf_comp3) *pctsrf_old(i,nsrf_comp3)
                 v10m(i,nsrf)  = v10m(i,nsrf_comp2) *pctsrf_old(i,nsrf_comp2) + v10m(i,nsrf_comp3) *pctsrf_old(i,nsrf_comp3)
-                tke(i,:,nsrf) = tke(i,:,nsrf_comp2)*pctsrf_old(i,nsrf_comp2) + tke(i,:,nsrf_comp3)*pctsrf_old(i,nsrf_comp3)
+                if (iflag_pbl > 1) then
+                 tke(i,:,nsrf) = tke(i,:,nsrf_comp2)*pctsrf_old(i,nsrf_comp2) + tke(i,:,nsrf_comp3)*pctsrf_old(i,nsrf_comp3)
+                endif
             
                 ! Security abort. This option has never been tested. To test, comment the following line.