Skip to content
Snippets Groups Projects
Commit 9714ce05 authored by RIETTE Sébastien's avatar RIETTE Sébastien
Browse files

S. Riette 13/6/2022 Bugfix for unphysical points

parent bead35b8
No related branches found
No related tags found
No related merge requests found
......@@ -69,7 +69,11 @@ INTEGER :: JI
IF (LHOOK) CALL DR_HOOK('ICE4_NUCLEATION', 0, ZHOOK_HANDLE)!
!
!$mnh_expand_where(JI=1:KSIZE)
GNEGT(:)=PT(:)<CST%XTT .AND. PRVT(:)>ICED%XRTMIN(1) .AND. ODCOMPUTE(:)
WHERE(ODCOMPUTE(:))
GNEGT(:)=PT(:)<CST%XTT .AND. PRVT(:)>ICED%XRTMIN(1)
ELSEWHERE
GNEGT(:)=.FALSE.
ENDWHERE
!$mnh_end_expand_where(JI=1:KSIZE)
ZUSW(:)=0.
......
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