Fix for bug with improper magnitude returned of earth's field

Modules/LemmaCore/src/EarthModel.cpp

 		BMag = bfield.norm(  );
 	    BInc = std::acos (bfield.dot(Vector3r(0,0,1)) / BMag) ;
 	    BDec = std::acos (bfield.dot(Vector3r(1,0,0)) / BMag) ;
-		BFieldUnit = BField.array() / BMag;
+		BFieldUnit = BField / BField.norm();
  	}
 
 	void EarthModel::SetMagneticFieldIncDecMag(const Real &inc,
 		BField(0) = BMag * cos(BInc*(PI/180.)) * cos(BDec*(PI/180.));
 		BField(1) = BMag * cos(BInc*(PI/180.)) * sin(BDec*(PI/180.));
 		BField(2) = BMag * sin(BInc*(PI/180.))  ;
-		BFieldUnit = BField.array() / BMag;
+		BFieldUnit = BField / BField.norm();
 	}
 
 	// ====================  INQUIRY       =======================
 	}
 
 	Vector3r EarthModel::GetMagneticFieldUnitVector() {
-		return this->BFieldUnit;
+		return BField/BField.norm();
 	}
 
 	Real  EarthModel::GetMagneticFieldMagnitude() {