source: LMDZ4/trunk/libf/dyn3dpar/cray.F @ 5427

Last change on this file since 5427 was 1279, checked in by Laurent Fairhead, 15 years ago

Merged LMDZ4-dev branch changes r1241:1278 into the trunk
Running trunk and LMDZ4-dev in LMDZOR configuration on local
machine (sequential) and SX8 (4-proc) yields identical results
(restart and restartphy are identical binarily)
Log history from r1241 to r1278 is available by switching to
source:LMDZ4/branches/LMDZ4-dev-20091210

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 836 bytes
RevLine 
[630]1!
2! $Header$
3!
4#ifdef CRAY
5      SUBROUTINE riencray
6      END
7#else
8      subroutine scopy(n,sx,incx,sy,incy)
9c
10      IMPLICIT NONE
11c
12      integer n,incx,incy,ix,iy,i
13      real sx((n-1)*incx+1),sy((n-1)*incy+1)
14c
[1279]15      if (incx.eq.1.and.incy.eq.1) then
16      do 10 i=1,n
17         sy(i)=sx(i)
1810    continue
19      else
[630]20      iy=1
21      ix=1
[1279]22      do 11 i=1,n
[630]23         sy(iy)=sx(ix)
24         ix=ix+incx
25         iy=iy+incy
[1279]2611    continue
27      endif
[630]28c
29      return
30      end
31
32      function ssum(n,sx,incx)
33c
34      IMPLICIT NONE
35c
36      integer n,incx,i,ix
37      real ssum,sx((n-1)*incx+1)
38c
39      ssum=0.
[1279]40      if (incx.eq.1) then
41      do 10 i=1,n
42         ssum=ssum+sx(i)
4310    continue
44      else
[630]45      ix=1
[1279]46      do 11 i=1,n
[630]47         ssum=ssum+sx(ix)
48         ix=ix+incx
[1279]4911    continue
50      endif
[630]51c
52      return
53      end
54#endif
Note: See TracBrowser for help on using the repository browser.