Skip to content
Snippets Groups Projects
Commit 582bf442 authored by RODIER Quentin's avatar RODIER Quentin
Browse files

Quentin 11/12/2020: temporary force the possible uncorrected negative values...

Quentin 11/12/2020: temporary force the possible uncorrected negative values to be > 0. This commit should be reverted soon and the missing right correction of negative values of rv must be found
parent 7d5539c7
No related branches found
No related tags found
No related merge requests found
...@@ -90,6 +90,7 @@ CLOUD: select case ( hcloud ) ...@@ -90,6 +90,7 @@ CLOUD: select case ( hcloud )
( zcph(:, :, :) * zexn(:, :, :) ) ( zcph(:, :, :) * zexn(:, :, :) )
prrs(:, :, :, 2) = 0. prrs(:, :, :, 2) = 0.
end where end where
prrs(:,:,:,1) = max(0.0,prrs(:,:,:,1))
case( 'ICE3', 'ICE4' ) case( 'ICE3', 'ICE4' )
...@@ -141,6 +142,8 @@ CLOUD: select case ( hcloud ) ...@@ -141,6 +142,8 @@ CLOUD: select case ( hcloud )
prrs(:, :, :, 4) = prrs(:, :, :, 4) - zcor(:, :, :) prrs(:, :, :, 4) = prrs(:, :, :, 4) - zcor(:, :, :)
end where end where
end if end if
prrs(:,:,:,1) = max(0.0,prrs(:,:,:,1))
! !
! !
case( 'C2R2', 'KHKO' ) case( 'C2R2', 'KHKO' )
...@@ -163,6 +166,8 @@ CLOUD: select case ( hcloud ) ...@@ -163,6 +166,8 @@ CLOUD: select case ( hcloud )
prrs(:, :, :, 2) = 0. prrs(:, :, :, 2) = 0.
prsvs(:, :, :, nsv_c2r2beg + 1) = 0. prsvs(:, :, :, nsv_c2r2beg + 1) = 0.
end where end where
prrs(:,:,:,1) = max(0.0,prrs(:,:,:,1))
! !
! !
case( 'LIMA' ) case( 'LIMA' )
...@@ -224,7 +229,7 @@ CLOUD: select case ( hcloud ) ...@@ -224,7 +229,7 @@ CLOUD: select case ( hcloud )
deallocate( zcor ) deallocate( zcor )
end if end if
end if end if
prrs(:,:,:,1) = max(0.0,prrs(:,:,:,1))
prsvs(:, :, :, nsv_lima_beg : nsv_lima_end) = Max( 0.0, prsvs(:, :, :, nsv_lima_beg : nsv_lima_end) ) prsvs(:, :, :, nsv_lima_beg : nsv_lima_end) = Max( 0.0, prsvs(:, :, :, nsv_lima_beg : nsv_lima_end) )
end select CLOUD end select CLOUD
......
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