Skip to content
Snippets Groups Projects
Commit 3817d727 authored by VIE Benoît's avatar VIE Benoît
Browse files

bugfix low values of ri everywhere

parent be045b6f
No related branches found
No related tags found
No related merge requests found
......@@ -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 ! Not used
XCTMIN(6) = 1.0E-3 ! Not used
XCTMIN(7) = 1.0E-3 ! Not used
......
......@@ -275,7 +275,7 @@ CLOUD: select case ( hcloud )
! Correction where ri<0 or Ni<0
if ( lcold_lima ) 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(:, :, :) / &
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment