diff --git a/lib/potential/mkernel_tricone.cpp b/lib/potential/mkernel_tricone.cpp
index 1e00e48..556429b 100644
--- a/lib/potential/mkernel_tricone.cpp
+++ b/lib/potential/mkernel_tricone.cpp
@@ -212,8 +212,8 @@ void gctl::magkernel(matrix<double> &kernel, const array<mag_tricone> &top_ele,
 		for (j = 0; j < e_size; j++)
 		{
 			mag_b = magkernel_single(top_ele[j], obsp[i], Rs.get(i));
-			if (comp_type == Bx) kernel[i][j] = mag_b.z;
-			if (comp_type == By) kernel[i][j] = mag_b.y;
+			if (comp_type == Bx) kernel[i][j] = mag_b.y;
+			if (comp_type == By) kernel[i][j] = mag_b.z;
 			if (comp_type == Bz) kernel[i][j] = mag_b.x;
 		}
 	}
@@ -227,8 +227,8 @@ void gctl::magkernel(matrix<double> &kernel, const array<mag_tricone> &top_ele,
 		for (j = 0; j < e_size; j++)
 		{
 			mag_b = magkernel_single(btm_ele[j], obsp[i], Rs.get(i));
-			if (comp_type == Bx) kernel[i][j] -= mag_b.z;
-			if (comp_type == By) kernel[i][j] -= mag_b.y;
+			if (comp_type == Bx) kernel[i][j] -= mag_b.y;
+			if (comp_type == By) kernel[i][j] -= mag_b.z;
 			if (comp_type == Bz) kernel[i][j] -= mag_b.x;
 		}
 	}
@@ -300,8 +300,8 @@ void gctl::magkernel(spmat<double> &kernel, const array<mag_tricone> &top_ele, c
 			if (geometry3d::angle(obsp[i].s2c(), cen) < cut_angle*GCTL_Pi/180.0)
 			{
 				mag_b = magkernel_single(top_ele[j], obsp[i], Rs.get(i));
-				if (comp_type == Bx) triplts[tri_idx[j + i*e_size]].val = mag_b.z;
-				if (comp_type == By) triplts[tri_idx[j + i*e_size]].val = mag_b.y;
+				if (comp_type == Bx) triplts[tri_idx[j + i*e_size]].val = mag_b.y;
+				if (comp_type == By) triplts[tri_idx[j + i*e_size]].val = mag_b.z;
 				if (comp_type == Bz) triplts[tri_idx[j + i*e_size]].val = mag_b.x;
 			}
 		}
@@ -319,8 +319,8 @@ void gctl::magkernel(spmat<double> &kernel, const array<mag_tricone> &top_ele, c
 			if (geometry3d::angle(obsp[i].s2c(), cen) < cut_angle*GCTL_Pi/180.0)
 			{
 				mag_b = magkernel_single(btm_ele[j], obsp[i], Rs.get(i));
-				if (comp_type == Bx) triplts[tri_idx[j + i*e_size]].val -= mag_b.z;
-				if (comp_type == By) triplts[tri_idx[j + i*e_size]].val -= mag_b.y;
+				if (comp_type == Bx) triplts[tri_idx[j + i*e_size]].val -= mag_b.y;
+				if (comp_type == By) triplts[tri_idx[j + i*e_size]].val -= mag_b.z;
 				if (comp_type == Bz) triplts[tri_idx[j + i*e_size]].val -= mag_b.x;
 			}
 		}