pro pinacate, infile,outfile
	openr, un, /get, infile
	openw, oun, /get, outfile
		
	s=''
	while s ne 'Vent' do begin
		readf, un, s
		s=strsplit(s,' ',/extract)
		s=s(where(s ne ' '))
		s=s[0]
	endwhile
			
	northing=0d
	easting=0d
	i=0
	j=0
	while not eof(un) do begin
		i=i+1
		s=''
		readf, un, s
		s=strsplit(s,' ',/extract)
		s=s(where(s ne ' '))
		s=s[0:2]
		
		tmp=s(0)
		reads, tmp, j		
		
		if j ne i then begin 
			print, j
			i=j
		endif
		
		tmp=s(1)
		reads, tmp, easting
		tmp=s(2)
		reads, tmp, northing
		
		easting=1000*(easting+200.)
		northing=1000*(northing+3500.)

		printf, oun, j, convertUTMtoLL([easting,northing])
		if j gt 417 then while not eof(un) do readf, un, tmp
		
	endwhile

	free_lun, oun, un
end
		

function convertUTMtoLL, locations

	units=envi_translate_projection_units('Degrees')

;; changed datum from WGS84 to NAD27 10/26/2005 after dad emailed
;; corrected info
;	oproj=envi_proj_create(/geographic, units=units, datum='WGS-84')
	oproj=envi_proj_create(/geographic, units=units, datum='NAD27')
	units=envi_translate_projection_units('Meters')
;	iproj=envi_proj_create(/utm, units=units, datum='WGS-84', zone=12)
	iproj=envi_proj_create(/utm, units=units, datum='NAD27', zone=12)

	envi_convert_projection_coordinates, locations(0), locations(1), iproj, $
						newXmap, newYmap, oproj
	map=dblarr(2)
	map(0)=newXmap
	map(1)=newYmap

	return, map
end