Ignore:
Timestamp:
May 22, 2014, 5:30:38 PM (10 years ago)
Author:
musat
Message:

Calcul de bils_ech seulement si iflag_ener_conserv >=1
IM

  • This line, and those below, will be ignored--

M ener_conserv.F90
M phys_output_write_mod.F90

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

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/ener_conserv.F90

    r2042 r2051  
    153153
    154154      bils_ec(:)=0.
    155       bils_ech(:)=0.
    156155      bils_tke(:)=0.
    157156      bils_diss(:)=0.
     
    161160      DO k=1,klev
    162161        bils_ec(:)=bils_ec(:)-d_t_ec(:,k)*masse(:,k)
    163         bils_ech(:)=bils_ech(:)-d_t_ech(:,k)*masse(:,k)
    164162        bils_tke(:)=bils_tke(:)+0.5*(dtke(:,k)+dtke(:,k+1))*masse(:,k)
    165163        bils_diss(:)=bils_diss(:)-d_t_diss(:,k)*masse(:,k)
     
    175173      ENDDO
    176174      bils_ec(:)=rcpd*bils_ec(:)/pdtphys
    177       bils_ech(:)=rcpd*bils_ech(:)/pdtphys
    178175      bils_tke(:)=bils_tke(:)/pdtphys
    179176      bils_diss(:)=rcpd*bils_diss(:)/pdtphys
     
    181178      bils_enthalp(:)=rcpd*bils_enthalp(:)/pdtphys
    182179      bils_latent(:)=rlvtt*bils_latent(:)/pdtphys
     180
     181IF (iflag_ener_conserv>=1) THEN
     182      bils_ech(:)=0.
     183      DO k=1,klev
     184        bils_ech(:)=bils_ech(:)-d_t_ech(:,k)*masse(:,k)
     185      ENDDO
     186      bils_ech(:)=rcpd*bils_ech(:)/pdtphys
     187ENDIF
     188
    183189RETURN
    184190
  • LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

    r2042 r2051  
    467467       CALL histwrite_phy(o_bils_diss, bils_diss)
    468468       CALL histwrite_phy(o_bils_ec, bils_ec)
    469        CALL histwrite_phy(o_bils_ech, bils_ech)
     469       IF (iflag_ener_conserv>=1) THEN
     470         CALL histwrite_phy(o_bils_ech, bils_ech)
     471       ENDIF
    470472       CALL histwrite_phy(o_bils_tke, bils_tke)
    471473       CALL histwrite_phy(o_bils_kinetic, bils_kinetic)
Note: See TracChangeset for help on using the changeset viewer.