Changeset 1236 for LMDZ4/trunk/libf


Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phylmd/pbl_surface_mod.F90

    r1146 r1236  
    12781278    INCLUDE "dimsoil.h"
    12791279    INCLUDE "clesphys.h"
     1280    INCLUDE "compbl.h"
    12801281
    12811282! Input variables
     
    13361337                u10m(i,nsrf)  = u10m(i,nsrf_comp1)
    13371338                v10m(i,nsrf)  = v10m(i,nsrf_comp1)
    1338                 tke(i,:,nsrf) = tke(i,:,nsrf_comp1)
     1339                if (iflag_pbl > 1) then
     1340                 tke(i,:,nsrf) = tke(i,:,nsrf_comp1)
     1341                endif
    13391342                mfois(nsrf) = mfois(nsrf) + 1
    13401343             ELSE
     
    13481351                u10m(i,nsrf)  = u10m(i,nsrf_comp2) *pctsrf_old(i,nsrf_comp2) + u10m(i,nsrf_comp3) *pctsrf_old(i,nsrf_comp3)
    13491352                v10m(i,nsrf)  = v10m(i,nsrf_comp2) *pctsrf_old(i,nsrf_comp2) + v10m(i,nsrf_comp3) *pctsrf_old(i,nsrf_comp3)
    1350                 tke(i,:,nsrf) = tke(i,:,nsrf_comp2)*pctsrf_old(i,nsrf_comp2) + tke(i,:,nsrf_comp3)*pctsrf_old(i,nsrf_comp3)
     1353                if (iflag_pbl > 1) then
     1354                 tke(i,:,nsrf) = tke(i,:,nsrf_comp2)*pctsrf_old(i,nsrf_comp2) + tke(i,:,nsrf_comp3)*pctsrf_old(i,nsrf_comp3)
     1355                endif
    13511356           
    13521357                ! Security abort. This option has never been tested. To test, comment the following line.
Note: See TracChangeset for help on using the changeset viewer.