smeanToZero.F

      subroutine smeanToZero
     &     (A, m, n, meanA, maxm, maxn,VALEX)

c---------------------------------------------------------------------c
c     "iterativeEOF", Version 2.0, December, 2002.                    c
c      A Krylov-based code for computing EOF                          c
c      Vincent TOUMAZOU                                               c
c      Vincent.Toumazou@cnes.fr                                       c
c      CNES/MERCATOR-Ocean                                            c
c      18, Avenue Edouard Belin                                       c
c      F-31401 TOULOUSE Cedex 4                                       c
c      URL: http://www.mercator.com.fr                                c
c      URL: http://www.cnes.fr                                        c
c---------------------------------------------------------------------c
c      This subroutine transforms the matrix A such that the means    c
c      of the raws of the output matrix are equal to 0.               c
c                                                                     c
c      INPUT  :                                                       c
c      ------                                                         c
c         - A        : SINGLE PRECISION                               c
c                      m by n matrix                                  c
c         - m, n     : INTEGER                                        c
c                      Dimensions of A                                c
c         - maxm, maxn are the leading dimension of A                 c
c         - meanA    : SINGLE PRECISION                               c
c                      vector of length m whose components are the    c
c                      means of the raws of A.                        c
c                                                                     c
c      OUTPUT :                                                       c
c      ------                                                         c
c         - A        : SINGLE PRECISION                               c
c                      Transformation of A such that the raws have    c
c                      zero means.                                    c
c---------------------------------------------------------------------c

      integer          m, n, maxm, maxn
      real             A(maxm, maxn), meanA(maxm),VALEX

      integer          i,j

      do i=1,m
         do j=1,n
            if(A(i,j).ne.VALEX) then
            A(i,j)=A(i,j)-meanA(i)
            endif
         enddo
      enddo

      end