30   int q,q0,al0,al1,pmax;  in Initial_CntCoes()  local
1007 	  al0 = tmp_index2[L0][Mul0][M0];   in Initial_CntCoes()
1011 	  for (p=0; p<Spe_Specified_Num[wan][al0]; p++){  in Initial_CntCoes()
1019 	    for (p=0; p<Spe_Specified_Num[wan][al0]; p++){  in Initial_CntCoes()
1020 	      sum = sum + CntCoes[Mc_AN][al0][p]*CntCoes[Mc_AN][al1][p];  in Initial_CntCoes()
1023 	    for (p=0; p<Spe_Specified_Num[wan][al0]; p++){  in Initial_CntCoes()
1028 	  for (p=0; p<Spe_Specified_Num[wan][al0]; p++){  in Initial_CntCoes()
1029 	    CntCoes[Mc_AN][al0][p] = CntCoes[Mc_AN][al0][p] - Tmp_CntCoes[p];  in Initial_CntCoes()
1037 	  for (p=0; p<Spe_Specified_Num[wan][al0]; p++){  in Initial_CntCoes()
1038 	    sum = sum + CntCoes[Mc_AN][al0][p]*CntCoes[Mc_AN][al0][p];  in Initial_CntCoes()
1039 	    if (Max0<fabs(CntCoes[Mc_AN][al0][p])){  in Initial_CntCoes()
1040 	      Max0 = fabs(CntCoes[Mc_AN][al0][p]);  in Initial_CntCoes()
1050 	  tmp1 = sgn(CntCoes[Mc_AN][al0][pmax]);  in Initial_CntCoes()
1052 	  for (p=0; p<Spe_Specified_Num[wan][al0]; p++){  in Initial_CntCoes()
1053 	    CntCoes[Mc_AN][al0][p] = tmp0*tmp1*CntCoes[Mc_AN][al0][p];  in Initial_CntCoes()
1083           al0 = al;   in Initial_CntCoes()
1096           al = al0;  in Initial_CntCoes()
1126 	    al0 = tmp_index2[L0][Mul0][M0];   in Initial_CntCoes()
1130 	    for (p=0; p<Spe_Specified_Num[wan][al0]; p++){  in Initial_CntCoes()
1138 	      for (p=0; p<Spe_Specified_Num[wan][al0]; p++){  in Initial_CntCoes()
1139 		sum = sum + CntCoes[Mc_AN][al0][p]*CntCoes[Mc_AN][al1][p];  in Initial_CntCoes()
1142 	      for (p=0; p<Spe_Specified_Num[wan][al0]; p++){  in Initial_CntCoes()
1147 	    for (p=0; p<Spe_Specified_Num[wan][al0]; p++){  in Initial_CntCoes()
1148 	      CntCoes[Mc_AN][al0][p] = CntCoes[Mc_AN][al0][p] - Tmp_CntCoes[p];  in Initial_CntCoes()
1155 	    for (p=0; p<Spe_Specified_Num[wan][al0]; p++){  in Initial_CntCoes()
1156 	      sum = sum + CntCoes[Mc_AN][al0][p]*CntCoes[Mc_AN][al0][p];  in Initial_CntCoes()
1157 	      if (Max0<fabs(CntCoes[Mc_AN][al0][p])){  in Initial_CntCoes()
1158 		Max0 = fabs(CntCoes[Mc_AN][al0][p]);  in Initial_CntCoes()
1163 	    tmp1 = sgn(CntCoes[Mc_AN][al0][pmax]);  in Initial_CntCoes()
1165 	    for (p=0; p<Spe_Specified_Num[wan][al0]; p++){  in Initial_CntCoes()
1166 	      CntCoes[Mc_AN][al0][p] = tmp0*tmp1*CntCoes[Mc_AN][al0][p];  in Initial_CntCoes()