c
***********************************************************************
*      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      *
*                                         March 18, 1995 K. Fujiki    *
***********************************************************************

      subroutine max512(menu,map,amx,imx,lmx)
c
      real map(512,512)
      integer iloc(2)
c
      if (menu.gt.0) then
        amx=maxval(map)
        iloc=maxloc(map)
        imx=iloc(1)
        lmx=iloc(2)
      else
        amx=maxval(abs(map))
        iloc=maxloc(abs(map))
        imx=iloc(1)
        lmx=iloc(2)
        amx=sign(1.,map(imx,lmx))*amx
      end if
c
      return
      end

***********************************************************************
*      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      *
*                                         March 18, 1995 K. Fujiki    *
***********************************************************************


      subroutine max256(menu,map,amx,imx,lmx)
c
      real map(256,256)
      integer iloc(2)
c
      if (menu.gt.0) then
        amx=maxval(map)
        iloc=maxloc(map)
        imx=iloc(1)
        lmx=iloc(2)
      else
        amx=maxval(abs(map))
        iloc=maxloc(abs(map))
        imx=iloc(1)
        lmx=iloc(2)
        amx=sign(1.,map(imx,lmx))*amx
      end if
c
      return
      end
c
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      *
*                                         March 18, 1995 K. Fujiki    *
***********************************************************************


      subroutine max128(menu,map,amx,imx,lmx)
c
      real map(128,128)
      integer iloc(2)
c
      if (menu.gt.0) then
        amx=maxval(map)
        iloc=maxloc(map)
        imx=iloc(1)
        lmx=iloc(2)
      else
        amx=maxval(abs(map))
        iloc=maxloc(abs(map))
        imx=iloc(1)
        lmx=iloc(2)
        amx=sign(1.,map(imx,lmx))*amx
      end if
c
      return
      end