function gravity = dial2gravity (dial_reading)

%gravimeter data
datashit = [ 	[2900 3000 3100 3200 ]
		[3058.5 3164.06 3269.63 3375.22] 
		[1.05562 1.05574 1.05587 1.005598]];

%loop over readings
for j=1:length(dial_reading)
	dial = dial_reading(j);
	
	if ((dial < 2900) | (dial > 3299))
		printf ("ERROR!!!!!!!!!haha! value out of range")
		continue;
	%find Milligal and interval values
	end
	
	for i=1:4
		if (floor( dial_reading(j)/100) = datashit(1, i)/100 )
			intervalvalue= datashit(3, i);
			millgalvalue = datashit(2, i);
		end
		tableline =i;
	end
	
	%calculate offset (reading - table)
	offset = dial - datashit(1,tableline);

	%multiply by interval value
	offset = offset * intervalvalue;
	
	%add to reading
	gravity(j) = millgalvalue + offset;
end