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

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/pbl_surface_mod.F90

    r1227 r1235  
    12871287    INCLUDE "dimsoil.h"
    12881288    INCLUDE "clesphys.h"
     1289    INCLUDE "compbl.h"
    12891290
    12901291! Input variables
     
    13451346                u10m(i,nsrf)  = u10m(i,nsrf_comp1)
    13461347                v10m(i,nsrf)  = v10m(i,nsrf_comp1)
    1347                 tke(i,:,nsrf) = tke(i,:,nsrf_comp1)
     1348                if (iflag_pbl > 1) then
     1349                 tke(i,:,nsrf) = tke(i,:,nsrf_comp1)
     1350                endif
    13481351                mfois(nsrf) = mfois(nsrf) + 1
    13491352             ELSE
     
    13571360                u10m(i,nsrf)  = u10m(i,nsrf_comp2) *pctsrf_old(i,nsrf_comp2) + u10m(i,nsrf_comp3) *pctsrf_old(i,nsrf_comp3)
    13581361                v10m(i,nsrf)  = v10m(i,nsrf_comp2) *pctsrf_old(i,nsrf_comp2) + v10m(i,nsrf_comp3) *pctsrf_old(i,nsrf_comp3)
    1359                 tke(i,:,nsrf) = tke(i,:,nsrf_comp2)*pctsrf_old(i,nsrf_comp2) + tke(i,:,nsrf_comp3)*pctsrf_old(i,nsrf_comp3)
     1362                if (iflag_pbl > 1) then
     1363                 tke(i,:,nsrf) = tke(i,:,nsrf_comp2)*pctsrf_old(i,nsrf_comp2) + tke(i,:,nsrf_comp3)*pctsrf_old(i,nsrf_comp3)
     1364                endif
    13601365           
    13611366                ! Security abort. This option has never been tested. To test, comment the following line.
Note: See TracChangeset for help on using the changeset viewer.