- Timestamp:
- Nov 5, 2018, 3:24:59 PM (6 years ago)
- Location:
- LMDZ6/branches/DYNAMICO-conv
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/DYNAMICO-conv
- Property svn:mergeinfo changed
/LMDZ6/trunk removed
- Property svn:mergeinfo changed
-
LMDZ6/branches/DYNAMICO-conv/libf/phylmd/ener_conserv.F90
r3356 r3411 222 222 DO k=1,klev 223 223 bils_ec(:)=bils_ec(:)-d_t_ec(:,k)*masse(:,k) 224 bils_tke(:)=bils_tke(:)+0.5*(dtke(:,k)+dtke(:,k+1))*masse(:,k) 224 225 bils_diss(:)=bils_diss(:)-d_t_diss(:,k)*masse(:,k) 225 226 bils_kinetic(:)=bils_kinetic(:)+masse(:,k)* & … … 234 235 ENDDO 235 236 bils_ec(:)=rcpd*bils_ec(:)/pdtphys 237 bils_tke(:)=bils_tke(:)/pdtphys 236 238 bils_diss(:)=rcpd*bils_diss(:)/pdtphys 237 239 bils_kinetic(:)= 0.5*bils_kinetic(:)/pdtphys 238 240 bils_enthalp(:)=rcpd*bils_enthalp(:)/pdtphys 239 241 bils_latent(:)=rlvtt*bils_latent(:)/pdtphys 240 !jyg<241 IF (iflag_pbl > 1) THEN242 DO k=1,klev243 bils_tke(:)=bils_tke(:)+0.5*(dtke(:,k)+dtke(:,k+1))*masse(:,k)244 ENDDO245 bils_tke(:)=bils_tke(:)/pdtphys246 ENDIF ! (iflag_pbl > 1)247 !>jyg248 242 249 243 IF (iflag_ener_conserv>=1) THEN
Note: See TracChangeset
for help on using the changeset viewer.