0845	Arrive. Safety meeting at 9. First, going to try and figure out why
		the extract.g script was failing for Sophia yesterday...?
	Looking at crop100 and crop200 in the SFDregrid folder, they are
		BAD!!! 
0857	Toney arrives... Followed by Johnny... safety training!
0925	Feeling safe - all done.
	crop100 and crop200 have PROBLEMS with reference pixels... when
		cropping, they get coordinates thrown way off. bad! delete!
	deleted.
0950	dropped off new travel form, for trip home in Aug.
	DJ says that triangles come from undefined pixels. try to mask/avoid/
		remove them, and things should work better.
	DJ also said to re-download the GASS data cubes to my machine and try
		running extract.g on them to see if it works or not... they're
		going, but will take a while to finish.
1030	after re-downloading, extract.g works fine on my machine.
1035	also tried cropping the SFD image first, and then re-gridding, but
		coordinate system still doesn't carry over right, and results
		are totally garbage.
1040	heading up to talk to DJ before tour of GBT at 1045
1245	back from tour, time for lunch
1315	back from lunch. going to keep working on undefined/masked pixels
		and try to eliminate triangle problem.
1340	command to mask image based on itself, with logic:
	- im.calcmask(mask='$im>0')
		...though this only works for certain types of images/
		arrays...
1350	There might be another problem since when you load the SFD (sgp or
		ngp) data into aips++, it won't display it unless you zoom
		in slightly to avoid looking at the edges...
1410	Better way to edit image masks:
		im is just some image
	    - local p, m
	    - im.getregion(pixels=p, pixelmask=m)
	    - m[1:10,1:10] := F
	    - im.putregion (pixelmask=m)
	this will mask out a 10x10 region in bottom left of image
1435	SFD data actually goes 0.5 degrees further than it's supposed to, in
		some places... at 0,90,180,270, it stops right at 0deg, but
		at 45, etc, it goes 0.5deg further past 0!
1445	Sophia's going to try logging into Ninkasi here to see if extract.g
		will work then. Seems to be a tricky problem...
1505	When I'm on makar, I don't have file permissions and can't run the
		scripts anyway. When logged into makar, but SSHed into ninkasi
		it works. When she's logged into ninkasi, but SSHED into makar
		it works. Therefore, we have learned nothing new...
	Also have been reading the Rupen article more, and some of the one
		by Westpfahl, too.
1515	tiny fragment of code to mask all 0 points in re-gridded SFD data:
		im := image('SFD_i100_sgp_regrid.img')
		local p,m
		im.getregion(pixels=p,pixelmask=m)
		pl := len(p)
		ml := len(m)
		for (i in 1:pl) {
			if (p[i]==0) {
				m[i] := F
			}
		}
		im.putregion(pixelmask=m)
		im.tofits('SFD_i100_sgp_regrid_mask.fits')
1550	and again for the ngp!
	to combine these two, use:
		im1 := image('SFD_i100_ngp_regrid_mask.fits')
		im2 := image('SFD_i100_sgp_regrid_mask.fits')
		...
		...?   continued below >>>>

1635	Talked to Toney, and then to Sue Ann about borrowing an optical
		telescope sometime for fun. She said come back Monday and
		she'd get me, Ron Maddalena, and Cara Ross (?) and teach me
		how to use it. Then sometime we could pull it out of the
		visitor center, and have the janitor turn the lights off
		and use it up here. Sounds good! Just have to go back Mon...
1650	Talked with DJ about progress... he said he can't find the "undefined"
		pixels either (in the SFD IR data), but thinks the weird
		display/grid-line issue	is related. Even if we could crop
		the pixels around the edges, we'd lose the galactic disk,
		and then there'd be a gap when we put the ngp and sgp back
		together. Going to talk to Jay tomorrow about it, and maybe
		just email the SFD group and ask for data gridded in the
		same way our GASS data is.
	continued attempts to seam together ngp and sgp masked/regridded...
		im1 := image('SFD_i100_ngp_regrid_mask.fits')
		im2 := image('SFD_i100_sgp_regrid_mask.fits')
		im1a := im1.getchunk()
		im2a := im2.getchunk()
		cs := im1.coordsys()
		for (i in 1:2497) {
			for (j in 1:2496) {
				if (!(im1a[i,j] > 0)) {
					im1a[i,j]:=0
				}
				if (!(im2a[i,j] > 0)) {
					im2a[i,j]:=0
				}
			}
		}
		im3a:=im1a+im2a
		im3 := imagefromarray(pixels=im3a,outfile='SFD_i100_all_regrid_mask.img',csys=cs)
	other problem: when adding these images together, the overlap matters!
		if doing simple addition, there will be double-sampled areas,
		PROBLEMS!!
1740	Heading out... going with Jonathan and Courtney to Trent's, then home.