/*
** (c) 1996-2000 The Regents of the University of California (through
** E.O. Lawrence Berkeley National Laboratory), subject to approval by
** the U.S. Department of Energy.  Your use of this software is under
** license -- the license agreement is attached and included in the
** directory as license.txt or you may contact Berkeley Lab's Technology
** Transfer Department at TTD@lbl.gov.  NOTICE OF U.S. GOVERNMENT RIGHTS.
** The Software was developed under funding from the U.S. Government
** which consequently retains certain rights as follows: the
** U.S. Government has been granted for itself and others acting on its
** behalf a paid-up, nonexclusive, irrevocable, worldwide license in the
** Software to reproduce, prepare derivative works, and perform publicly
** and display publicly.  Beginning five (5) years after the date
** permission to assert copyright is obtained from the U.S. Department of
** Energy, and subject to any subsequent five (5) year renewals, the
** U.S. Government is granted for itself and others acting on its behalf
** a paid-up, nonexclusive, irrevocable, worldwide license in the
** Software to reproduce, prepare derivative works, distribute copies to
** the public, perform publicly and display publicly, and to permit
** others to do so.
*/

#include "REAL.H"
#include "ArrayLim.H"

#define SDIM 3

      subroutine AVGPLANE(v,DIMS(v),avg,DIMS(avg),jlo,jhi)
      implicit none
      integer DIMDEC(v)
      integer DIMDEC(avg)
      integer jlo, jhi
      REAL_T   v(DIMV(v))
      REAL_T avg(DIMV(avg))
      integer i,j,k,jj,N

      jj = ARG_L2(avg)
      do k=ARG_L3(v),ARG_H3(v)
         do i=ARG_L1(v),ARG_H1(v)
            avg(i,jj,k) = 0.d0
            do j=jlo,jhi
               avg(i,jj,k) = avg(i,jj,k) + v(i,j,k)
            enddo
            N = jhi - jlo + 1
            avg(i,jj,k) = avg(i,jj,k) / N
         enddo
      enddo
      end

