geometry={h1;o,h1,r;h2,o,r,h1,theta} r=1.9 theta=102 runmrci={if(.not.scfdone) hf;multi;mrci} runmrci extrapolate,procedure=runmrci,bases=vtz:vqz