casa6

*=======================================================================

*     Copyright (C) 1999,2001,2002

*     Associated Universities, Inc. Washington DC, USA.

*

*     This library is free software; you can redistribute it and/or

*     modify it under the terms of the GNU Library General Public

*     License as published by the Free Software Foundation; either

*     version 2 of the License, or (at your option) any later version.

*

*     This library is distributed in the hope that it will be useful,

*     but WITHOUT ANY WARRANTY; without even the implied warranty of

*     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

*     GNU Library General Public License for more details.

*

*     You should have received a copy of the GNU Library General Public

*     License along with this library; if not, write to the Free

*     Software Foundation, Inc., 675 Massachusetts Ave, Cambridge,

*     MA 02139, USA.

*

*     Correspondence concerning AIPS++ should be addressed as follows:

*            Internet email: aips2-request@nrao.edu.

*            Postal address: AIPS++ Project Office

*                            National Radio Astronomy Observatory

*                            520 Edgemont Road

*                            Charlottesville, VA 22903-2475 USA

*

*     $Id$

*-----------------------------------------------------------------------

C

C Calculate gridded coordinates

C

      subroutine sgridsdclip (xy, sampling, pos, loc, off)

      implicit none

      integer sampling

      integer loc(2), off(2)

      double precision xy(2)

      real pos(2)

      integer idim

​

      do idim=1,2

         pos(idim)=xy(idim)+1.0

         loc(idim)=nint(pos(idim))

         off(idim)=nint((loc(idim)-pos(idim))*sampling)

      end do

      return 

      end

C

C Is this on the grid?

C

      logical function ogridsdclip (nx, ny, loc, support)

      implicit none

      integer nx, ny, loc(2), support

      ogridsdclip=(loc(1)-support.ge.1).and.(loc(1)+support.le.nx).and.

     $     (loc(2)-support.ge.1).and.(loc(2)+support.le.ny)

      return

      end

C

C Grid a number of visibility records: single dish gridding

C but with complex images including additional process for

C min/max clipping 

C

C This subroutine assumes specific clipping process after

C all the data are accumulated