orbs = {"3s","3p","3d","4d"} f = ReadFPLOBasisFunctions(orbs,"+fval.001.1") orbs = {"3d3/2","3d5/2"} G, F = ReadFPLOBasisFunctions(orbs,"+fval.001.1","+fval.001.2")