Changeset 1235 for LMDZ4/branches/LMDZ4-dev

Timestamp:

Aug 31, 2009, 2:48:59 PM (15 years ago)

Author:

musat

Message:

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

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/pbl_surface_mod.F90 (modified) (3 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/pbl_surface_mod.F90

@@ -1287 +1287 @@
     INCLUDE "dimsoil.h"
     INCLUDE "clesphys.h"
+    INCLUDE "compbl.h"
 
 ! Input variables
@@ -1345 +1346 @@
                 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
@@ -1357 +1360 @@
                 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.