function [calcconcs calcabs]=normalise(calcconcs,obsdata,totconc)
	disp("normalise")	

	#concs
	# normalise concentrations
	for n=1:rows(calcconcs)
		scaling=sum(calcconcs(n,:));
		calcconcs(n,:)=(totconc/scaling).*calcconcs(n,:);
	endfor
	calcabs=pinv(calcconcs)*obsdata;

endfunction
