We apply real space full multiple scattering theory to interpret the Fe K-edge XANES spectra of TPP[Fe(Pc)L 2 ] 2 (L=CN, Cl and Br) systems, which show the giant magnetoresistance (GMR) at the low temperatures. In the previous paper, we have reported the absorption edge shift of the XANES spectra, whose origin remains unclear, for TPP[Fe(Pc)L 2 ] 2 systems. In order to clarify the relation between the charge of the Fe atom, the local structure of the axial ligand and the XANES spectra, we improve the calculation of the XANES spectra by taking into account the wider region including the neighboring Fe(Pc)L 2 and TPP molecules. Our multiple scattering analyses suggest that the spectral shape is strongly influenced by the distance between a central Fe and axial ligands L. The number of Fe 3d electrons obtained by density functional theory calculations show weak dependence on the axial ligands L. The EXAFS spectra, the polarization dependence and the temperature dependence of the XANES spectra are also discussed.