function addCRC(w)
	wave w
	// w should be an 8 bit wave
	Redimension/B/U w // treat data as unsigned
	variable LSB, MSB
	variable CRC=0xFFFF
	variable carryflag, i, j, byte
	for (i=0;i<numpnts(w); i+=1)
		byte=w[i]
		CRC=CRC%^byte  // XOR
		for (j=0; j<8; j+=1)		
			carryflag=CRC & 1 // if bit1 is set, carryflag=1 
			CRC=floor(CRC/2) // shift bits right
			if (carryflag) // XOR with A001
				CRC = CRC %^ 0xA001
			endif
		endfor		
	endfor
	LSB=floor(CRC/256)
	MSB=mod(CRC,256)
	w[numpnts(w)] = {MSB, LSB}
end