--- trunk/dyn3d/ismax.f 2013/11/15 18:45:49 76 +++ trunk/dyn3d/ismax.f90 2014/03/05 14:38:41 81 @@ -1,24 +1,25 @@ -! -! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/ismax.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $ -! - function ismax(n,sx,incx) -c - IMPLICIT NONE -c - INTEGER n,i,incx,ismax,ix - real sx((n-1)*incx+1),sxmax -c - ix=1 - ismax=1 - sxmax=sx(1) - do 10 i=1,n-1 - ix=ix+incx - if(sx(ix).gt.sxmax) then - sxmax=sx(ix) - ismax=i+1 - endif -10 continue -c - return - end + +! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/ismax.F,v 1.1.1.1 2004/05/19 +! 12:53:07 lmdzadmin Exp $ + +FUNCTION ismax(n, sx, incx) + + IMPLICIT NONE + + INTEGER n, i, incx, ismax, ix + REAL sx((n-1)*incx+1), sxmax + + ix = 1 + ismax = 1 + sxmax = sx(1) + DO i = 1, n - 1 + ix = ix + incx + IF (sx(ix)>sxmax) THEN + sxmax = sx(ix) + ismax = i + 1 + END IF + END DO + + RETURN +END FUNCTION ismax