diff --git a/src/MNH/lima_ice_deposition.f90 b/src/MNH/lima_ice_deposition.f90 index 8c7c57e4091305b31d8906cec050cd18427f63cb..752feb93be82551252500320e25b21cf6c02217b 100644 --- a/src/MNH/lima_ice_deposition.f90 +++ b/src/MNH/lima_ice_deposition.f90 @@ -65,6 +65,7 @@ SUBROUTINE LIMA_ICE_DEPOSITION (PTSTEP, LDCOMPUTE, & !! MODIFICATIONS !! ------------- !! Original 15/03/2018 +!! B. Vié 30/08/2021 Disable CNVS if LSNOW=F !! !------------------------------------------------------------------------------- ! @@ -171,5 +172,9 @@ P_CI_CNVS(:) = - ZZW2(:) ! END WHERE ! +IF (.NOT.LSNOW) THEN + P_RI_CNVS(:) = 0. + P_CI_CNVS(:) = 0. +END IF ! END SUBROUTINE LIMA_ICE_DEPOSITION diff --git a/src/MNH/lima_snow_deposition.f90 b/src/MNH/lima_snow_deposition.f90 index 697f9ee74f5f9101579f9724421e76ee6f93d614..d39f714d83639730bbff492e14ae8001b06cdd64 100644 --- a/src/MNH/lima_snow_deposition.f90 +++ b/src/MNH/lima_snow_deposition.f90 @@ -149,7 +149,7 @@ WHERE( GMICRO ) ! ZZW(:) = 0.0 WHERE ( (PRST(:)>XRTMIN(5)) ) - ZZW(:) = ( PSSI(:)/(PAI(:)) ) * & + ZZW(:) = ( PSSI(:)/(PAI(:)*PRHODREF(:)) ) * & ( X0DEPS*PLBDS(:)**XEX0DEPS + X1DEPS*PCJ(:)*PLBDS(:)**XEX1DEPS ) ZZW(:) = ZZW(:)*(0.5+SIGN(0.5,ZZW(:))) - ABS(ZZW(:))*(0.5-SIGN(0.5,ZZW(:))) END WHERE