- Timestamp:
- Jun 29, 2018, 12:31:11 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/branches/DYNAMICO-conv/libf/phylmd/ener_conserv.F90
r2903 r3356 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)225 224 bils_diss(:)=bils_diss(:)-d_t_diss(:,k)*masse(:,k) 226 225 bils_kinetic(:)=bils_kinetic(:)+masse(:,k)* & … … 235 234 ENDDO 236 235 bils_ec(:)=rcpd*bils_ec(:)/pdtphys 237 bils_tke(:)=bils_tke(:)/pdtphys238 236 bils_diss(:)=rcpd*bils_diss(:)/pdtphys 239 237 bils_kinetic(:)= 0.5*bils_kinetic(:)/pdtphys 240 238 bils_enthalp(:)=rcpd*bils_enthalp(:)/pdtphys 241 239 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 242 248 243 249 IF (iflag_ener_conserv>=1) THEN
Note: See TracChangeset
for help on using the changeset viewer.