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

Quentin 26/09/2024: add kernels ice4_nucleation

parent c4040f2e
No related branches found
No related tags found
No related merge requests found
......@@ -76,8 +76,10 @@ ELSEWHERE
ENDWHERE
!$mnh_end_expand_where(JI=1:KSIZE)
!$acc kernels
ZUSW(:)=0.
ZZW(:)=0.
!$acc end kernels
!$mnh_expand_where(JI=1:KSIZE)
WHERE(GNEGT(:))
ZZW(:)=ALOG(PT(:))
......@@ -85,8 +87,9 @@ WHERE(GNEGT(:))
ZZW(:)=EXP(CST%XALPI - CST%XBETAI/PT(:) - CST%XGAMI*ZZW(:)) ! es_i
END WHERE
!$mnh_end_expand_where(JI=1:KSIZE)
!$acc kernels
ZSSI(:)=0.
!$acc end kernels
!$mnh_expand_where(JI=1:KSIZE)
WHERE(GNEGT(:))
ZZW(:)=MIN(PPABST(:)/2., ZZW(:)) ! safety limitation
......@@ -104,7 +107,9 @@ WHERE(GNEGT(:))
END WHERE
!$mnh_end_expand_where(JI=1:KSIZE)
!$acc kernels
ZZW(:)=0.
!$acc loop independent
DO JI=1,KSIZE
IF(GNEGT(JI)) THEN
IF(PT(JI)<CST%XTT-5.0 .AND. ZSSI(JI)>0.0) THEN
......@@ -115,14 +120,16 @@ DO JI=1,KSIZE
ENDIF
ENDIF
ENDDO
!$acc end kernels
!$mnh_expand_where(JI=1:KSIZE)
WHERE(GNEGT(:))
ZZW(:)=ZZW(:)-PCIT(:)
ZZW(:)=MIN(ZZW(:), 50.E3) ! limitation provisoire a 50 l^-1
END WHERE
!$mnh_end_expand_where(JI=1:KSIZE)
!$acc kernels
PRVHENI_MR(:)=0.
!$acc end kernels
!$mnh_expand_where(JI=1:KSIZE)
WHERE(GNEGT(:))
!
......@@ -134,7 +141,9 @@ END WHERE
!$mnh_end_expand_where(JI=1:KSIZE)
!Limitation due to 0 crossing of temperature
IF(PARAMI%LFEEDBACKT) THEN
!$acc kernels
ZW(:)=0.
!$acc end kernels
!$mnh_expand_where(JI=1:KSIZE)
WHERE(GNEGT(:))
ZW(:)=MIN(PRVHENI_MR(:), &
......
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