Program n : grgg, grll, grps

      program progn
      implicit none

      real xlat(2,3),xlon(2,3),roots(3)
      integer i,j

      print *,'This example define three grids:'
# define a global gaussian grid
      print *,'The gaussian grid defined:'
      print *,''
      call dgauss(roots,3,1)
      call grgg(xlat,xlon,2,3,roots,0)
      do i=1,2
         do j=1,3
            print *,'(',xlat(i,j),xlon(i,j),')'
         enddo
         print *,''
      enddo   
      print *,''

# define a latidude-logitude grid
      print *,'The latitude-longitude grid defined:'
      print *,''
      call grll(xlat,xlon,2,3,45.5,-73.0,1.0,1.0)
      
      do i=1,2
         do j=1,3
            print *,'(',xlat(i,j),xlon(i,j),')'
         enddo
         print *,''
      enddo   
      print *,''

# define a polar stereographic grid
      print *,'The polar stereographic grid defined:'
      print *,''
      call grps(xlat,xlon,2,3,5.0,6.0,20000.0,20.0,1)
      do i=1,2
         do j=1,3
            print *,'(',xlat(i,j),xlon(i,j),')'
         enddo
         print *,''
      enddo   
      
      stop
      end