if(.not.geom_tag_to_element(aname, symb, elem, atn)) + call errquit('cosmo def radii:cannot resolve tag' + //aname, 0,0)