implicit  real*8 (a-h,o-z)
 include  'parameter'
 dimension xlat(mxy), xlon(mxy), ylat(mxy), ylon(mxy)
 include  'constant'

 do i=1,464 {
     xlat(i)= -35. + i/8. - 1/16.
     write(11,10) xlat(i)
     }
 do j=1,464 {
     xlon(j)= -80. + j/8. - 1/16.
     write(12,10) xlon(j)
     }

 read(5,2)
 do i=1,200000 {
    read(5,*,end=100) ylat(i), ylon(i)
    }
 100 continue
 ix= i - 1

 do i=1,464 {
    xlat1= xlat(i)
    do j=1,464 {
       xlon1= xlon(j)
       call inside( xlat1, xlon1, ylat, ylon, ix, ins1 )
       if (ins1 == 0) { }
       else {
         write(6,20) xlat1, xlon1
          }
       }
    }

  2 format( )
 10 format( f9.3, f9.3)
 20 format( f12.4, f12.4, '    10  0  Helvetica-Bold MC  X')
 
 end