Skip to content
Snippets Groups Projects
Commit 816b30e3 authored by WAUTELET Philippe's avatar WAUTELET Philippe
Browse files

Philippe 08/06/2018: SURFEX correction: AINT(...,KIND=16) crashes PGI and does not seem useful

parent 2f3d9fd8
No related branches found
No related tags found
No related merge requests found
...@@ -99,7 +99,7 @@ ENDDO ...@@ -99,7 +99,7 @@ ENDDO
DO JP = 1,SIZE(PW,2) DO JP = 1,SIZE(PW,2)
DO JI = 1,SIZE(PW,1) DO JI = 1,SIZE(PW,1)
IF (PW(JI,JP)/=XUNDEF) THEN IF (PW(JI,JP)/=XUNDEF) THEN
ZINT = AINT(PW(JI,JP),KIND=16) ZINT = AINT(PW(JI,JP))
IF (PW(JI,JP)/=ZINT) THEN IF (PW(JI,JP)/=ZINT) THEN
PW(JI,JP) = ZINT + NINT((PW(JI,JP)-ZINT)*ZPREC)/ZPREC PW(JI,JP) = ZINT + NINT((PW(JI,JP)-ZINT)*ZPREC)/ZPREC
ENDIF ENDIF
......
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