Ignore:
Timestamp:
Jul 18, 2013, 9:52:31 AM (12 years ago)
Author:
Ehouarn Millour
Message:

Some cleanup around ismin/ismax/scopy/ssum routines which are defined in multiple places and should definitely only be in 'bibio'.
EM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phy1d/1DUTILS.h_with_writelim

    r1763 r1794  
    710710      RETURN
    711711      END
    712       subroutine scopy(n,sx,incx,sy,incy)
    713 !
    714       IMPLICIT NONE
    715 !
    716       integer n,incx,incy,ix,iy,i
    717       real sx((n-1)*incx+1),sy((n-1)*incy+1)
    718 !
    719       iy=1
    720       ix=1
    721       do 10 i=1,n
    722          sy(iy)=sx(ix)
    723          ix=ix+incx
    724          iy=iy+incy
    725 10    continue
    726 !
    727       return
    728       end
    729712      subroutine wrgradsfi(if,nl,field,name,titlevar)
    730713      implicit none
     
    26692652        return
    26702653        end
    2671 !===============================================================
    2672       function ismin(n,sx,incx)
    2673 
    2674       implicit none
    2675       integer n,i,incx,ismin,ix
    2676       real sx((n-1)*incx+1),sxmin
    2677 
    2678       ix=1
    2679       ismin=1
    2680       sxmin=sx(1)
    2681       do i=1,n-1
    2682          ix=ix+incx
    2683          if(sx(ix).lt.sxmin) then
    2684              sxmin=sx(ix)
    2685              ismin=i+1
    2686          endif
    2687       enddo
    2688 
    2689       return
    2690       end
    2691 
    2692 !===============================================================
    2693       function ismax(n,sx,incx)
    2694 
    2695       implicit none
    2696       integer n,i,incx,ismax,ix
    2697       real sx((n-1)*incx+1),sxmax
    2698 
    2699       ix=1
    2700       ismax=1
    2701       sxmax=sx(1)
    2702       do i=1,n-1
    2703          ix=ix+incx
    2704          if(sx(ix).gt.sxmax) then
    2705              sxmax=sx(ix)
    2706              ismax=i+1
    2707          endif
    2708       enddo
    2709 
    2710       return
    2711       end
    2712 
     2654
Note: See TracChangeset for help on using the changeset viewer.