bugfix unrealistic regions of low ri

src/MNH/ini_lima.f90

@@ -147,7 +147,7 @@ ALLOCATE( XCTMIN(7) )
 XCTMIN(1) = 1.0       ! Not used
 XCTMIN(2) = 1.0E-3    ! Nc
 XCTMIN(3) = 1.0E-3    ! Nr
-XCTMIN(4) = 1.0E-3    ! Ni
+XCTMIN(4) = 1.0E-0    ! Ni
 XCTMIN(5) = 1.0E-3    ! Ns
 XCTMIN(6) = 1.0E-3    ! Ng
 XCTMIN(7) = 1.0E-3    ! Nh

src/MNH/sources_neg_correct.f90

@@ -276,7 +276,7 @@ CLOUD: select case ( hcloud )
 ! Correction where ri<0 or Ni<0
      if ( krr.GE.4 ) then
         zmask(:,:,:)=(prrs(:, :, :, 4) < xrtmin_lima(4) / ptstep)
-        if (nmom_i.ge.2) zmask(:,:,:)=(zmask(:,:,:) .or. prsvs(:, :, :, nsv_lima_ni) < 0. )
+        if (nmom_i.ge.2) zmask(:,:,:)=(zmask(:,:,:) .or. prsvs(:, :, :, nsv_lima_ni) < xctmin_lima(4) / ptstep)
         where ( zmask(:,:,:) )
            prrs(:, :, :, 1) = prrs(:, :, :, 1) + prrs(:, :, :, 4)
            prths(:, :, :) = prths(:, :, :) - prrs(:, :, :, 4) * zls(:, :, :) /  &