source: trunk/LMDZ.GENERIC/libf/misc/isamax.F @ 1974

Last change on this file since 1974 was 135, checked in by aslmd, 14 years ago

CHANGEMENT ARBORESCENCE ETAPE 2 -- NON COMPLET

File size: 385 bytes
Line 
1      function isamax(n,sx,incx)
2c
3      IMPLICIT NONE
4c
5      INTEGER n,i,incx,isamax,ix
6      real sx((n-1)*incx+1),sxmax,abssx
7c
8      ix=1
9      isamax=1
10      sxmax=ABS(sx(1))
11      do 10 i=1,n-1
12         ix=ix+incx
13         abssx=ABS(sx(ix))
14         if(abssx.gt.sxmax) then
15             sxmax=abssx
16             isamax=i+1
17         endif
1810    continue
19c
20      return
21      end
Note: See TracBrowser for help on using the repository browser.