<fix> default startmag_type should be 1, for element magnetic instead of atomic magnetic.

<range> only for atomic_rho

Author: dyzheng

source/module_cell/read_atoms.cpp

@@ -618,7 +618,7 @@ bool UnitCell_pseudo::read_atom_positions(std::ifstream &ifpos, std::ofstream &o
                                         mv.z = true ;
                                         atoms[it].vel[ia].set(0,0,0);
 #ifndef __CMD
-										atoms[it].mag[ia]=magnet.start_magnetization[it];
+										//atoms[it].mag[ia]=magnet.start_magnetization[it];//if this line is used, default startmag_type would be 2
 #endif										
 										atoms[it].angle1[ia]=0;
 										atoms[it].angle2[ia]=0;

source/src_pw/charge.cpp

@@ -331,7 +331,7 @@ void Charge::atomic_rho(const int spin_number_need, double** rho_in)const		// Pe
 				else if(spin_number_need==2)
 				{
 					if(startmag_type==1)
-					{ cout<<"this should not be called";
+					{
 						for (int ig = 0; ig < GlobalC::pw.ngmc ; ig++)
 						{
 							const std::complex<double> swap = GlobalC::pw.strucFac(it, ig)* rho_lgl[GlobalC::pw.ig2ngg[ig]];
@@ -375,7 +375,7 @@ void Charge::atomic_rho(const int spin_number_need, double** rho_in)const		// Pe
 				{
 					//noncolinear case
 					if(startmag_type == 1)
-					{cout<<"this would not be call";
+					{
 						for (int ig = 0; ig < GlobalC::pw.ngmc ; ig++)
 						{
 							const std::complex<double> swap = GlobalC::pw.strucFac(it, ig)* rho_lgl[GlobalC::pw.ig2ngg[ig]];