diff --git a/src/MNH/convect_updraft.f90 b/src/MNH/convect_updraft.f90 index 13b65f82a29a5f071f864d5d4153785357541f50..d860b62401954899c0df857b3178df86e8ea38c3 100644 --- a/src/MNH/convect_updraft.f90 +++ b/src/MNH/convect_updraft.f90 @@ -144,7 +144,6 @@ END MODULE MODI_CONVECT_UPDRAFT !! Original 07/11/95 !! Last modified 10/12/97 !! V.Masson, C.Lac, Sept. 2010 : Correction of a loop for reproducibility -! P. Wautelet 12/02/2019: bugfix: PURR/PURS were sometimes =-0. -> problems later !------------------------------------------------------------------------------- ! !* 0. DECLARATIONS @@ -380,9 +379,7 @@ DO JK = MAX( IKB + 1, JKMIN ), IKE - 1 PUTPR(:) = PUTPR(:) + PUPR(:,JKP) ! total precipitation rate ZWORK2(:) = PURR(:,JKP) / MAX( 1.E-8, PURC(:,JKP) + PURI(:,JKP) ) PURR(:,JKP) = ZWORK2(:) * PURC(:,JKP) ! liquid precipitation - PURR(:,JKP)=MAX(TINY(0.),PURR(:,JKP)) !To prevent problems when =-0. (occur sometimes) PURS(:,JKP) = ZWORK2(:) * PURI(:,JKP) ! solid precipitation - PURS(:,JKP)=MAX(TINY(0.),PURS(:,JKP)) !To prevent problems when =-0. (occur sometimes) ! ! !* 7. Update r_c, r_i, enthalpy, r_w for precipitation