Molar weight calculator

//EDIT: 2/23/16; optimised function "ElementCount"
// The function molwt(str) returns the molar weight of a chemical formula specified by "str", e.g. molwt("H2O") returns 18.0148 (in grams per mole). 
// The function is case-sensitive, e.g. molwt("h2o") will not return what you expect.
// The function ElementCount(input) returns a list of elements and their abundance, e.g. ElementCount("H2O") returns "H;2;O;1;"
// The element abundance in the input string MUST be a positive integer

function molwt(str)
	string str
	
	variable mw
	variable AN
	variable weight
	variable n
	variable i
	
	string List = ElementCount(str)
	string Element =""
	
	//initialise elements and atomic weights: list index equals atomic number	
	string el=""
	el += ";H;He;Li;Be;B;C;N;O;F;Ne;Na;Mg;Al;Si;P;S;Cl;Ar;"
	el += "K;Ca;Sc;Ti;V;Cr;Mn;Fe;Co;Ni;Cu;Zn;Ga;Ge;As;Se;Br;Kr;"
	el += "Rb;Sr;Y;Zr;Nb;Mo;Tc;Ru;Rh;Pd;Ag;Cd;In;Sn;Sb;Te;I;Xe;"
	el += "Cs;Ba;La;Ce;Pr;Nd;Pm;Sm;Eu;Gd;Tb;Dy;Ho;Er;Tm;Yb;Lu;"
	el += "Hf;Ta;W;Re;Os;Ir;Pt;Au;Hg;Tl;Pb;Bi;Po;At;Rn;"
	el += "Fr;Ra;Ac;Th;Pa;U;Np;Pu;"
	
	string aw=""
	aw+=";1.0079;4.0026;6.941;9.0122;10.811;12.011;14.007;15.999;18.998;20.18;22.99;24.305;26.982;28.085;30.974;32.065;35.453;39.948;"
	aw+="39.098;40.078;44.956;47.867;50.941;51.996;54.938;55.845;58.933;58.693;63.546;65.409;69.723;72.64;74.922;78.96;79.904;83.798;"
	aw+="85.468;87.62;88.906;91.224;92.906;95.94;98.906;101.07;102.91;106.42;107.87;112.41;114.82;118.71;121.76;127.6;126.9;131.29;"
	aw+="132.91;137.33;138.91;140.12;140.91;144.24;146.92;150.36;151.96;157.25;158.93;162.5;164.93;167.26;168.93;173.04;174.97;"
	aw+="178.49;180.95;183.84;186.21;190.23;192.22;195.08;196.97;200.59;204.38;207.2;208.98;208.98;209.99;222.02;223.02;"
	aw+="226.03;227.03;232.04;231.04;238.03;237.05;244.06;"
	

	do
		Element = StringFromList(i, List) 
		
		if ( strlen(Element) == 0) 
			break
		endif
		
		AN = WhichListItem(Element, el)
		weight = str2num(StringFromList(AN, aw))
		n =  str2num(StringFromList(i+1, List))
		
		mw += (n * weight) 
				
		i+=2
	while(1)
	
	return mw
end

function/S ElementCount(input)  // analyses input string for elements and number of elements
	String input
 
	string list = ""
	variable i
	variable TheEnd = StrLen(input)
	
	for(i=0;i<TheEnd; i+=1)
 
		// consider cases where string item is a letter and upper case		
		if (GrepString( input[i], "[[:alpha:]]")  &&  GrepString(input[i], "[[:upper:]]") )						
				list += input[i]
 
				// is it a two character element, e.g. Fe instead of F?
				if (GrepString(input[i+1], "[[:alpha:]]")  && GrepString(input[i+1], "[[:lower:]]") )	
					list +=  input[i+1]
					i+=1
				endif
 
 				//this element is covered
 				list +=";"
 				 				
				// test if next list item is numeric or string; if it is string, then the present element has a count of 1, add this to list ....
				// ...or check if i+1 is the last entry
				if (GrepString(input[i+1], "[[:alpha:]]") || i+1 == TheEnd)
					list += "1;"
				endif
				
				// we are done
				continue
 
		endif
 
		// consider cases where string item is a number			
		if (GrepString(input[i], "[[:digit:]]") )
 
				list +=input[i]
 
				// is there more than one numeric character?
				do
					if (GrepString(input[i+1] "[[:digit:]]"))
						list += input[i+1]
						i+=1						
					else
						break
					endif					
 
				while(1)
 
				List +=";"		
		endif			
	endfor

	return list
end