***********************************************************************
*      max512 - calculate mean and search maximum                     *
*        menu : input,integer  clean menu                             *
*        map  : input,real(512,512)  dirty map                        *
*        amx  : output,real    maximum of map                         *
*        imx  : output,integer x-coordinate of the maximum point      *
*        lmx  : output,integer y-coordinate of the maximum point      *
*                                         Feb. 1999  K. Fujiki        *
***********************************************************************
c
      subroutine max512(menu,map,amx,imx,lmx)
c
      real map(512,512)
c
      amx=0.
      if (menu.gt.0) then
       do 100 l=1,512
       do 100 i=1,512
        if (amx.lt.map(i,l)) then
         amx=map(i,l)
         imx=i
         lmx=l
        end if
  100  continue
      else
       do 110 l=1,512
       do 110 i=1,512
        if (abs(amx).lt.abs(map(i,l))) then
         amx=map(i,l)
         imx=i
         lmx=l
        end if
  110  continue
      end if
c
      return
      end
c
***********************************************************************
*      max256 - calculate mean and search maximum                     *
*        menu : input,integer  clean menu                             *
*        map  : input,real(256,256)  dirty map                        *
*        amx  : output,real    maximum of map                         *
*        imx  : output,integer x-coordinate of the maximum point      *
*        lmx  : output,integer y-coordinate of the maximum point      *
*                                         Feb. 1999  K. Fujiki        *
***********************************************************************
c
      subroutine max256(menu,map,amx,imx,lmx)
c
      real map(256,256)
c
      amx=0.
      if (menu.gt.0) then
       do 100 l=1,256
       do 100 i=1,256
        if (amx.lt.map(i,l)) then
         amx=map(i,l)
         imx=i
         lmx=l
        end if
  100  continue
      else
       do 110 l=1,256
       do 110 i=1,256
        if (abs(amx).lt.abs(map(i,l))) then
         amx=map(i,l)
         imx=i
         lmx=l
        end if
  110  continue
      end if
c
      return
      end
c
***********************************************************************
*      max128 - calculate mean and search maximum                     *
*        menu : input,integer  clean menu                             *
*        map  : input,real(128,128)  dirty map                        *
*        amx  : output,real    maximum of map                         *
*        imx  : output,integer x-coordinate of the maximum point      *
*        lmx  : output,integer y-coordinate of the maximum point      *
*                                         Feb. 1999  K. Fujiki        *
***********************************************************************
c
      subroutine max128(menu,map,amx,imx,lmx)
c
      real map(128,128)
c
      amx=0.
      if (menu.gt.0) then
       do 100 l=1,128
       do 100 i=1,128
        if (amx.lt.map(i,l)) then
         amx=map(i,l)
         imx=i
         lmx=l
        end if
  100  continue
      else
       do 110 l=1,128
       do 110 i=1,128
        if (abs(amx).lt.abs(map(i,l))) then
         amx=map(i,l)
         imx=i
         lmx=l
        end if
  110  continue
      end if
c
      return
      end
c