Changeset 3071


Ignore:
Timestamp:
Nov 14, 2017, 4:17:16 PM (7 years ago)
Author:
jyg
Message:

Correction in ener_conserv.F90: bils_tke is
set to zero if iflag_pbl <= 1.

Cosmetic correction in physiq_mod.F90: pbl_tke is
initialized to zero if iflag_pbl <= 1 and
klon_glo=1.

Location:
LMDZ6/trunk/libf/phylmd
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/ener_conserv.F90

    r2903 r3071  
    222222      DO k=1,klev
    223223        bils_ec(:)=bils_ec(:)-d_t_ec(:,k)*masse(:,k)
    224         bils_tke(:)=bils_tke(:)+0.5*(dtke(:,k)+dtke(:,k+1))*masse(:,k)
    225224        bils_diss(:)=bils_diss(:)-d_t_diss(:,k)*masse(:,k)
    226225        bils_kinetic(:)=bils_kinetic(:)+masse(:,k)* &
    … …  
    235234      ENDDO
    236235      bils_ec(:)=rcpd*bils_ec(:)/pdtphys
    237       bils_tke(:)=bils_tke(:)/pdtphys
    238236      bils_diss(:)=rcpd*bils_diss(:)/pdtphys
    239237      bils_kinetic(:)= 0.5*bils_kinetic(:)/pdtphys
    240238      bils_enthalp(:)=rcpd*bils_enthalp(:)/pdtphys
    241239      bils_latent(:)=rlvtt*bils_latent(:)/pdtphys
     240!jyg<
     241      IF (iflag_pbl > 1) THEN
     242        DO k=1,klev
     243          bils_tke(:)=bils_tke(:)+0.5*(dtke(:,k)+dtke(:,k+1))*masse(:,k)
     244        ENDDO
     245        bils_tke(:)=bils_tke(:)/pdtphys
     246      ENDIF  ! (iflag_pbl > 1)
     247!>jyg
    242248
    243249IF (iflag_ener_conserv>=1) THEN

  • LMDZ6/trunk/libf/phylmd/physiq_mod.F90

    r3048 r3071  
    14001400!jyg<
    14011401       IF (klon_glo==1) THEN
    1402           pbl_tke(:,:,is_ave) = 0.
    1403           DO nsrf=1,nbsrf
    1404             DO k = 1,klev+1
    1405                  pbl_tke(:,k,is_ave) = pbl_tke(:,k,is_ave) &
    1406                      +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
    1407             ENDDO
    1408           ENDDO
     1402          IF (iflag_pbl > 1) THEN         
     1403              pbl_tke(:,:,is_ave) = 0.
     1404              DO nsrf=1,nbsrf
     1405                DO k = 1,klev+1
     1406                     pbl_tke(:,k,is_ave) = pbl_tke(:,k,is_ave) &
     1407                         +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
     1408                ENDDO
     1409              ENDDO
     1410          ELSE   ! (iflag_pbl > 1)
     1411              pbl_tke(:,:,:) = 0.
     1412          ENDIF  ! (iflag_pbl > 1)
    14091413!>jyg
    14101414       ENDIF
Note: See TracChangeset for help on using the changeset viewer.