Changeset 1794 for LMDZ5/trunk/libf/phy1d/1DUTILS.h_with_writelim
- Timestamp:
- Jul 18, 2013, 9:52:31 AM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phy1d/1DUTILS.h_with_writelim
r1763 r1794 710 710 RETURN 711 711 END 712 subroutine scopy(n,sx,incx,sy,incy)713 !714 IMPLICIT NONE715 !716 integer n,incx,incy,ix,iy,i717 real sx((n-1)*incx+1),sy((n-1)*incy+1)718 !719 iy=1720 ix=1721 do 10 i=1,n722 sy(iy)=sx(ix)723 ix=ix+incx724 iy=iy+incy725 10 continue726 !727 return728 end729 712 subroutine wrgradsfi(if,nl,field,name,titlevar) 730 713 implicit none … … 2669 2652 return 2670 2653 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.