c
c "read.f" reads Dr. Scambos' Antarctic Grounding and Coast Line data files as 
c		posted at the National Snow and Ice Center, e.g.,
c		http://nsidc.org/data/docs/agdc/nsdic0280_scambos/
c	  	and extracts data for a region of interest.
c
c Jan.-17, 2012 Andreas Muenchow, University of Delaware
c
	integer n,ix,iy,idev
	real x,y,xlat,xlon,xlat0,xlat1,xlon0,xlon1
	character*11 dummy
c
c Pine Island region, Antarctica
c
	xlat0 = -74
	xlat1 = -75.5
	xlon0 = -98.0
	xlon1 = -104.0
c	
	write(6,*)'Antarctic MOA: 1-Grounding, 2-Coast Line'
	read(5,*)idev
	if (idev .eq. 1) then
	    open (unit=7,file='MOA_groundingline.txt')
	    open (unit=8,file='Grounding.dat')
	else
	    open (unit=7,file='MOA_coast_file.txt')
	    open (unit=8,file='Coast.dat')
	end if
c
	read(7,101)dummy,n
	do 1 i=1,n
	   read(7,*)line,ix,iy,x,y,xlat,xlon,b1
	   if (xlat.lt.xlat0 .and. xlat.gt.xlat1) then
	   if (xlon.lt.xlon0 .and. xlon.gt.xlon1) then
		write(8,102)xlon,xlat
	   end if
	   end if
1	continue
	close (7)
102	format(2f10.5)
101	format(//////,a11,i7/)
	stop
	end