diff --git a/src/CMakeLists.txt b/src/CMakeLists.txt
index e67b23e..6e584f2 100644
--- a/src/CMakeLists.txt
+++ b/src/CMakeLists.txt
@@ -13,6 +13,8 @@ add_executable(gm3d_gui WIN32 MACOSX_BUNDLE ${DIR_SRC})
 
 include_directories(/usr/local/include)
 
-find_library(FLTK_LIBRARY fltk HINTS /usr/local/lib)
+find_library(HINTS /usr/local/lib)
+#find_library(FLTK_LIBRARY fltk HINTS /usr/local/lib)
+
 
 target_link_libraries(gm3d_gui PUBLIC ${FLTK_LIBRARY})
\ No newline at end of file
diff --git a/src/add_interface_block.cpp b/src/add_interface_block.cpp
index 51f3550..80fe1b2 100644
--- a/src/add_interface_block.cpp
+++ b/src/add_interface_block.cpp
@@ -1,5 +1,5 @@
 #include "gm3d.h"
-//我们读入一个界面数据 插值计算每个块体中心位置的值 然后按情况赋值
+//���Ƕ���һ���������� ��ֵ����ÿ����������λ�õ�ֵ Ȼ�������ֵ
 int GM3D::AddInterfaceBlock(modelist para_list){
 	int m, n, xnum, ynum;
 	double xs,xe,xmin,xmax,dx;
@@ -18,7 +18,7 @@ int GM3D::AddInterfaceBlock(modelist para_list){
 		if (open_infile(infile,filename)) return -1;
 
 		while (getline(infile,temp_str)){
-			//#range必须出现在数据之前
+			//#range�������������֮ǰ
 			if (*(temp_str.begin()) == '#'){
 				if (6 == sscanf(temp_str.c_str(),"# range=%lf/%lf/%lf/%lf/%lf/%lf",&xs,&dx,&xe,&ys,&dy,&ye)){
 					xmin = MIN(xs,xe); xmax = MAX(xs,xe);
@@ -75,7 +75,7 @@ int GM3D::AddInterfaceBlock(modelist para_list){
 
 				if (model_cube_[i].cen.z >= temp_topo.z){
 					if (model_block_val_[i] == BDL_MAX)
-						model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
+						model_block_val_[i] = para_list.mod_value; //ע���ظ���ֵ�Ŀ���Ḳ��
 					else
 						model_block_val_[i] += para_list.mod_value;
 					model_added = true;
@@ -92,7 +92,7 @@ int GM3D::AddInterfaceBlock(modelist para_list){
 
 				if (model_cube_[i].cen.z <= temp_topo.z){
 					if (model_block_val_[i] == BDL_MAX)
-						model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
+						model_block_val_[i] = para_list.mod_value; //ע���ظ���ֵ�Ŀ���Ḳ��
 					else
 						model_block_val_[i] += para_list.mod_value;
 					model_added = true;
diff --git a/src/add_model_gui.cpp b/src/add_model_gui.cpp
index fed6592..2e33ce9 100644
--- a/src/add_model_gui.cpp
+++ b/src/add_model_gui.cpp
@@ -18,7 +18,7 @@ Fl_Return_Button *can_add_btn=(Fl_Return_Button *)0;
 Fl_Button *sig_add_btn=(Fl_Button *)0;
 Fl_Input *mod_val_input=(Fl_Input *)0;
 
-// 声明三个字符串用于保存新建的模型类型和赋值类型
+// ���������ַ������ڱ����½���ģ�����ͺ͸�ֵ����
 char add_mod_type[1024];
 char add_val_type[1024];
 char agn_val_type[1024];
@@ -81,8 +81,8 @@ void cb_bot_val_rbtn(Fl_Button*, void*){
   return;
 }
 
-// 声明一个字符串用于保存新建的模型参数并添加至列表中
-// 模型参数的排列为<模型类型> <赋值类型> <物理参数字符串> <几何参数字符串>
+// ����һ���ַ������ڱ����½���ģ�Ͳ������������б���
+// ģ�Ͳ���������Ϊ<ģ������> <��ֵ����> <���������ַ���> <���β����ַ���>
 char add_mod_para[1024];
 
 void cb_sig_add_btn(Fl_Widget*, void*){
@@ -101,7 +101,7 @@ void cb_sig_add_btn(Fl_Widget*, void*){
   mod_para_brw->add(add_mod_para);
   return;
 }
-//取消添加模型
+//ȡ������ģ��
 void cb_can_add_btn(Fl_Widget*, void*){
 	add_mod_win->hide();
 	return;
diff --git a/src/add_models.cpp b/src/add_models.cpp
index e594adf..ec61292 100644
--- a/src/add_models.cpp
+++ b/src/add_models.cpp
@@ -10,7 +10,7 @@ int GM3D::AddModels(char* filename){
 	while(getline(infile,temp_str)){
 		if (*(temp_str.begin()) == '#') continue;
 		else{
-			//按每行5个数据解析 初始化为含观测值与不确定度的观测点
+			//��ÿ��5�����ݽ��� ��ʼ��Ϊ���۲�ֵ�벻ȷ���ȵĹ۲��
 			if (4 == sscanf(temp_str.c_str(),"%s %s %lf %s",
 				temp_list.mod_type,temp_list.val_type,&temp_list.mod_value,temp_list.mod_para)){
 				model_list_.push_back(temp_list);
diff --git a/src/add_regular_block.cpp b/src/add_regular_block.cpp
index 4834976..abb80a3 100644
--- a/src/add_regular_block.cpp
+++ b/src/add_regular_block.cpp
@@ -16,7 +16,7 @@ int GM3D::AddRegularBlock(modelist para_list){
 				if (model_cube_[i].cen.x >= xmin && model_cube_[i].cen.x <= xmax &&
 					model_cube_[i].cen.y >= ymin && model_cube_[i].cen.y <= ymax &&
 					model_cube_[i].cen.z >= zmin && model_cube_[i].cen.z <= zmax){
-					model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
+					model_block_val_[i] = para_list.mod_value; //ע���ظ���ֵ�Ŀ���Ḳ��
 					model_added = true;
 				}
 			}
@@ -27,7 +27,7 @@ int GM3D::AddRegularBlock(modelist para_list){
 					model_cube_[i].cen.y >= ymin && model_cube_[i].cen.y <= ymax &&
 					model_cube_[i].cen.z >= zmin && model_cube_[i].cen.z <= zmax){
 					if (model_block_val_[i] == BDL_MAX)
-						model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
+						model_block_val_[i] = para_list.mod_value; //ע���ظ���ֵ�Ŀ���Ḳ��
 					else
 						model_block_val_[i] += para_list.mod_value;
 					model_added = true;
@@ -39,7 +39,7 @@ int GM3D::AddRegularBlock(modelist para_list){
 				if (model_cube_[i].cen.x >= xmin && model_cube_[i].cen.x <= xmax &&
 					model_cube_[i].cen.y >= ymin && model_cube_[i].cen.y <= ymax &&
 					model_cube_[i].cen.z >= zmin && model_cube_[i].cen.z <= zmax){
-					model_block_val_[i] = BDL_MAX; //注意重复赋值的块体会覆盖
+					model_block_val_[i] = BDL_MAX; //ע���ظ���ֵ�Ŀ���Ḳ��
 					model_added = true;
 				}
 			}
diff --git a/src/add_sphere_block.cpp b/src/add_sphere_block.cpp
index 9bba9bf..521526b 100644
--- a/src/add_sphere_block.cpp
+++ b/src/add_sphere_block.cpp
@@ -20,7 +20,7 @@ int GM3D::AddSphereBlock(modelist para_list){
 				rad_limit = rad_x*rad_y*rad_z/sqrt(pow(rad_y*rad_z*sin(theta)*cos(phi),2) + pow(rad_x*rad_z*sin(theta)*sin(phi),2) + pow(rad_x*rad_y*cos(theta),2));
 
 				if (dist <= rad_limit){
-					model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
+					model_block_val_[i] = para_list.mod_value; //ע���ظ���ֵ�Ŀ���Ḳ��
 					model_added = true;
 				}
 			}
@@ -39,7 +39,7 @@ int GM3D::AddSphereBlock(modelist para_list){
 
 				if (dist <= rad_limit){
 					if (model_block_val_[i] == BDL_MAX)
-						model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
+						model_block_val_[i] = para_list.mod_value; //ע���ظ���ֵ�Ŀ���Ḳ��
 					else
 						model_block_val_[i] += para_list.mod_value;
 					model_added = true;
@@ -59,7 +59,7 @@ int GM3D::AddSphereBlock(modelist para_list){
 				rad_limit = rad_x*rad_y*rad_z/sqrt(pow(rad_y*rad_z*sin(theta)*cos(phi),2) + pow(rad_x*rad_z*sin(theta)*sin(phi),2) + pow(rad_x*rad_y*cos(theta),2));
 
 				if (dist <= rad_limit){
-					model_block_val_[i] = BDL_MAX; //注意重复赋值的块体会覆盖
+					model_block_val_[i] = BDL_MAX; //ע���ظ���ֵ�Ŀ���Ḳ��
 					model_added = true;
 				}
 			}
diff --git a/src/add_tilted_block.cpp b/src/add_tilted_block.cpp
index 904548c..444b4d0 100644
--- a/src/add_tilted_block.cpp
+++ b/src/add_tilted_block.cpp
@@ -19,7 +19,7 @@ int GM3D::AddTiltedBlock(modelist para_list){
 
 		if (!strcmp(para_list.val_type,"replace")){
 			for (int i = 0; i < model_num_; i++){
-				//计算当前层的x与y范围
+				//���㵱ǰ���x��y��Χ
 				layer_xmin = (model_cube_[i].cen.z - zmin)*(xmin_2 - xmin_1)/(zmax - zmin) + xmin_1;
 				layer_xmax = (model_cube_[i].cen.z - zmin)*(xmax_2 - xmax_1)/(zmax - zmin) + xmax_1;
 				layer_ymin = (model_cube_[i].cen.z - zmin)*(ymin_2 - ymin_1)/(zmax - zmin) + ymin_1;
@@ -28,14 +28,14 @@ int GM3D::AddTiltedBlock(modelist para_list){
 				if (model_cube_[i].cen.x >= layer_xmin && model_cube_[i].cen.x <= layer_xmax &&
 					model_cube_[i].cen.y >= layer_ymin && model_cube_[i].cen.y <= layer_ymax &&
 					model_cube_[i].cen.z >= zmin && model_cube_[i].cen.z <= zmax){
-					model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
+					model_block_val_[i] = para_list.mod_value; //ע���ظ���ֵ�Ŀ���Ḳ��
 					model_added = true;
 				}
 			}
 		}
 		else if (!strcmp(para_list.val_type,"add")){
 			for (int i = 0; i < model_num_; i++){
-				//计算当前层的x与y范围
+				//���㵱ǰ���x��y��Χ
 				layer_xmin = (model_cube_[i].cen.z - zmin)*(xmin_2 - xmin_1)/(zmax - zmin) + xmin_1;
 				layer_xmax = (model_cube_[i].cen.z - zmin)*(xmax_2 - xmax_1)/(zmax - zmin) + xmax_1;
 				layer_ymin = (model_cube_[i].cen.z - zmin)*(ymin_2 - ymin_1)/(zmax - zmin) + ymin_1;
@@ -45,7 +45,7 @@ int GM3D::AddTiltedBlock(modelist para_list){
 					model_cube_[i].cen.y >= layer_ymin && model_cube_[i].cen.y <= layer_ymax &&
 					model_cube_[i].cen.z >= zmin && model_cube_[i].cen.z <= zmax){
 					if (model_block_val_[i] == BDL_MAX)
-						model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
+						model_block_val_[i] = para_list.mod_value; //ע���ظ���ֵ�Ŀ���Ḳ��
 					else
 						model_block_val_[i] += para_list.mod_value;
 					model_added = true;
@@ -54,7 +54,7 @@ int GM3D::AddTiltedBlock(modelist para_list){
 		}
 		else if (!strcmp(para_list.val_type,"erase")){
 			for (int i = 0; i < model_num_; i++){
-				//计算当前层的x与y范围
+				//���㵱ǰ���x��y��Χ
 				layer_xmin = (model_cube_[i].cen.z - zmin)*(xmin_2 - xmin_1)/(zmax - zmin) + xmin_1;
 				layer_xmax = (model_cube_[i].cen.z - zmin)*(xmax_2 - xmax_1)/(zmax - zmin) + xmax_1;
 				layer_ymin = (model_cube_[i].cen.z - zmin)*(ymin_2 - ymin_1)/(zmax - zmin) + ymin_1;
@@ -63,7 +63,7 @@ int GM3D::AddTiltedBlock(modelist para_list){
 				if (model_cube_[i].cen.x >= layer_xmin && model_cube_[i].cen.x <= layer_xmax &&
 					model_cube_[i].cen.y >= layer_ymin && model_cube_[i].cen.y <= layer_ymax &&
 					model_cube_[i].cen.z >= zmin && model_cube_[i].cen.z <= zmax){
-					model_block_val_[i] = BDL_MAX; //注意重复赋值的块体会覆盖
+					model_block_val_[i] = BDL_MAX; //ע���ظ���ֵ�Ŀ���Ḳ��
 					model_added = true;
 				}
 			}
diff --git a/src/build_regular_grid.cpp b/src/build_regular_grid.cpp
index 060610c..1027563 100644
--- a/src/build_regular_grid.cpp
+++ b/src/build_regular_grid.cpp
@@ -46,26 +46,26 @@ int GM3D::BuildRegularGrid(char* space_para){
 		while(x >= xmin && x <= xmax){
 			z = zs;
 			while(z >= zmin && z <= zmax){
-				//添加denMod
+				//����denMod
 				temp_cu.cen.id = model_cube_.size();
 				temp_cu.cen.x = x; temp_cu.cen.y = y; temp_cu.cen.z = z;
-				//添加mshVert
+				//����mshVert
 				for (int i = 0; i < 8; i++){
-					temp_cp.id = model_vert_.size(); //添加msh的顶点索引为mshVert的大小
-					temp_cp.x = temp_cu.cen.x - sign[i][0]*temp_cu.dx; //左下底角
+					temp_cp.id = model_vert_.size(); //����msh�Ķ�������ΪmshVert�Ĵ�С
+					temp_cp.x = temp_cu.cen.x - sign[i][0]*temp_cu.dx; //���µ׽�
 					temp_cp.y = temp_cu.cen.y - sign[i][1]*temp_cu.dy;
 					temp_cp.z = temp_cu.cen.z - sign[i][2]*temp_cu.dz;
 					temp_id_str = cpoint_id(temp_cp);
 					imsp = map_str_point.find(temp_id_str);
-					//利用map_vert查到当前顶点是否存在,这里需要注意，如果顶点已经存在则只需要将顶点索引置为已存在顶点的索引，不增加顶点计数
+					//����map_vert�鵽��ǰ�����Ƿ����,������Ҫע�⣬��������Ѿ�������ֻ��Ҫ������������Ϊ�Ѵ��ڶ���������������Ӷ������
 					if(imsp!=map_str_point.end()){
 						temp_cu.ids[i] = imsp->second.id;
 					}
-					//若为新的顶点则将其增加到两个映射和一个链表中
+					//��Ϊ�µĶ����������ӵ�����ӳ���һ��������
 					else{
-						temp_cu.ids[i] = temp_cp.id;//新的顶点索引等于顶点集的数量
-						model_vert_.push_back(temp_cp);//将新产生的顶点保存到顶点链表中
-						map_str_point[temp_id_str] = temp_cp;//将新产生的顶点保存到顶点位置映射中
+						temp_cu.ids[i] = temp_cp.id;//�µĶ����������ڶ��㼯������
+						model_vert_.push_back(temp_cp);//���²����Ķ��㱣�浽����������
+						map_str_point[temp_id_str] = temp_cp;//���²����Ķ��㱣�浽����λ��ӳ����
 					}
 				}
 				model_cube_.push_back(temp_cu);
@@ -82,7 +82,7 @@ int GM3D::BuildRegularGrid(char* space_para){
 	else{
 		vert_num_ = model_vert_.size();
 		model_num_ = model_cube_.size();
-		model_block_val_.resize(model_num_,BDL_MAX); //初始化模型块体值为BDL_MAX
+		model_block_val_.resize(model_num_,BDL_MAX); //��ʼ��ģ�Ϳ���ֵΪBDL_MAX
 	}
 	return 0;
 }
\ No newline at end of file
diff --git a/src/forward_delta_t.cpp b/src/forward_delta_t.cpp
index 0ed81eb..8bf9ed2 100644
--- a/src/forward_delta_t.cpp
+++ b/src/forward_delta_t.cpp
@@ -10,18 +10,18 @@ int GM3D::ForwardDeltaT(char* noise_level,char* mag_para){
 	double R222,R122,R212,R112,R221,R121,R211,R111;
 	double G222,G122,G212,G112,G221,G121,G211,G111;
 
-	//初始化正演值和不确定度
+	//��ʼ������ֵ�Ͳ�ȷ����
 	for (int i = 0; i < obs_num_; i++){
 		obs_p_[i].val = obs_p_[i].dev = 0.0;
 	}
 
-	//确定噪声水平
+	//ȷ������ˮƽ
 	double noise_mean, noise_dev, temp_noise;
 	if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
 		noise_mean = noise_dev = 0.0;
 	}
 
-	//确定磁化参数
+	//ȷ���Ż�����
 	if (4 != sscanf(mag_para,"%lf/%lf/%lf/%lf",&I0,&A0,&I,&A)){
 		I0 = I = 90; A0 = A = 0;
 	}
@@ -39,7 +39,7 @@ int GM3D::ForwardDeltaT(char* noise_level,char* mag_para){
 		k6=-sin(I0)*sin(I);
 	}
 
-	//添加高斯噪声值
+	//���Ӹ�˹����ֵ
 	default_random_engine generator;
 	normal_distribution<double> dist(noise_mean, noise_dev);
 
diff --git a/src/forward_delta_tx.cpp b/src/forward_delta_tx.cpp
index f4c32b7..6a188d3 100644
--- a/src/forward_delta_tx.cpp
+++ b/src/forward_delta_tx.cpp
@@ -10,18 +10,18 @@ int GM3D::ForwardDeltaTx(char* noise_level,char* mag_para){
 	double R222,R122,R212,R112,R221,R121,R211,R111;
 	double G222,G122,G212,G112,G221,G121,G211,G111;
 
-	//初始化正演值和不确定度
+	//��ʼ������ֵ�Ͳ�ȷ����
 	for (int i = 0; i < obs_num_; i++){
 		obs_p_[i].val = obs_p_[i].dev = 0.0;
 	}
 
-	//确定噪声水平
+	//ȷ������ˮƽ
 	double noise_mean, noise_dev, temp_noise;
 	if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
 		noise_mean = noise_dev = 0.0;
 	}
 
-	//确定磁化参数
+	//ȷ���Ż�����
 	if (4 != sscanf(mag_para,"%lf/%lf/%lf/%lf",&I0,&A0,&I,&A)){
 		I0 = I = 90; A0 = A = 0;
 	}
@@ -39,7 +39,7 @@ int GM3D::ForwardDeltaTx(char* noise_level,char* mag_para){
 		k6=-sin(I0)*sin(I);
 	}
 
-	//添加高斯噪声值
+	//���Ӹ�˹����ֵ
 	default_random_engine generator;
 	normal_distribution<double> dist(noise_mean, noise_dev);
 
diff --git a/src/forward_delta_ty.cpp b/src/forward_delta_ty.cpp
index db3f5c0..4784d62 100644
--- a/src/forward_delta_ty.cpp
+++ b/src/forward_delta_ty.cpp
@@ -10,18 +10,18 @@ int GM3D::ForwardDeltaTy(char* noise_level,char* mag_para){
 	double R222,R122,R212,R112,R221,R121,R211,R111;
 	double G222,G122,G212,G112,G221,G121,G211,G111;
 
-	//初始化正演值和不确定度
+	//��ʼ������ֵ�Ͳ�ȷ����
 	for (int i = 0; i < obs_num_; i++){
 		obs_p_[i].val = obs_p_[i].dev = 0.0;
 	}
 
-	//确定噪声水平
+	//ȷ������ˮƽ
 	double noise_mean, noise_dev, temp_noise;
 	if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
 		noise_mean = noise_dev = 0.0;
 	}
 
-	//确定磁化参数
+	//ȷ���Ż�����
 	if (4 != sscanf(mag_para,"%lf/%lf/%lf/%lf",&I0,&A0,&I,&A)){
 		I0 = I = 90; A0 = A = 0;
 	}
@@ -39,7 +39,7 @@ int GM3D::ForwardDeltaTy(char* noise_level,char* mag_para){
 		k6=-sin(I0)*sin(I);
 	}
 
-	//添加高斯噪声值
+	//���Ӹ�˹����ֵ
 	default_random_engine generator;
 	normal_distribution<double> dist(noise_mean, noise_dev);
 
diff --git a/src/forward_delta_tz.cpp b/src/forward_delta_tz.cpp
index 5d49a61..7f7ffc3 100644
--- a/src/forward_delta_tz.cpp
+++ b/src/forward_delta_tz.cpp
@@ -11,18 +11,18 @@ int GM3D::ForwardDeltaTz(char* noise_level,char* mag_para){
 	double R222,R122,R212,R112,R221,R121,R211,R111;
 	double G222,G122,G212,G112,G221,G121,G211,G111;
 
-	//初始化正演值和不确定度
+	//��ʼ������ֵ�Ͳ�ȷ����
 	for (int i = 0; i < obs_num_; i++){
 		obs_p_[i].val = obs_p_[i].dev = 0.0;
 	}
 
-	//确定噪声水平
+	//ȷ������ˮƽ
 	double noise_mean, noise_dev, temp_noise;
 	if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
 		noise_mean = noise_dev = 0.0;
 	}
 
-	//确定磁化参数
+	//ȷ���Ż�����
 	if (4 != sscanf(mag_para,"%lf/%lf/%lf/%lf",&I0,&A0,&I,&A)){
 		I0 = I = 90; A0 = A = 0;
 	}
@@ -40,7 +40,7 @@ int GM3D::ForwardDeltaTz(char* noise_level,char* mag_para){
 		k6=-sin(I0)*sin(I);
 	}
 
-	//添加高斯噪声值
+	//���Ӹ�˹����ֵ
 	default_random_engine generator;
 	normal_distribution<double> dist(noise_mean, noise_dev);
 
diff --git a/src/forward_hax.cpp b/src/forward_hax.cpp
index d244bfc..79175a1 100644
--- a/src/forward_hax.cpp
+++ b/src/forward_hax.cpp
@@ -10,18 +10,18 @@ int GM3D::ForwardHax(char* noise_level,char* mag_para){
 	double R222,R122,R212,R112,R221,R121,R211,R111;
 	double G222,G122,G212,G112,G221,G121,G211,G111;
 
-	//初始化正演值和不确定度
+	//��ʼ������ֵ�Ͳ�ȷ����
 	for (int i = 0; i < obs_num_; i++){
 		obs_p_[i].val = obs_p_[i].dev = 0.0;
 	}
 
-	//确定噪声水平
+	//ȷ������ˮƽ
 	double noise_mean, noise_dev, temp_noise;
 	if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
 		noise_mean = noise_dev = 0.0;
 	}
 
-	//确定磁化参数
+	//ȷ���Ż�����
 	if (4 != sscanf(mag_para,"%lf/%lf/%lf/%lf",&I0,&A0,&I,&A)){
 		I = 90; A = 0;
 	}
@@ -34,7 +34,7 @@ int GM3D::ForwardHax(char* noise_level,char* mag_para){
 		Gamma=sin(I);
 	}
 
-	//添加高斯噪声值
+	//���Ӹ�˹����ֵ
 	default_random_engine generator;
 	normal_distribution<double> dist(noise_mean, noise_dev);
 
diff --git a/src/forward_hay.cpp b/src/forward_hay.cpp
index 59db245..ba21719 100644
--- a/src/forward_hay.cpp
+++ b/src/forward_hay.cpp
@@ -10,18 +10,18 @@ int GM3D::ForwardHay(char* noise_level,char* mag_para){
 	double R222,R122,R212,R112,R221,R121,R211,R111;
 	double G222,G122,G212,G112,G221,G121,G211,G111;
 
-	//初始化正演值和不确定度
+	//��ʼ������ֵ�Ͳ�ȷ����
 	for (int i = 0; i < obs_num_; i++){
 		obs_p_[i].val = obs_p_[i].dev = 0.0;
 	}
 
-	//确定噪声水平
+	//ȷ������ˮƽ
 	double noise_mean, noise_dev, temp_noise;
 	if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
 		noise_mean = noise_dev = 0.0;
 	}
 
-	//确定磁化参数
+	//ȷ���Ż�����
 	if (4 != sscanf(mag_para,"%lf/%lf/%lf/%lf",&I0,&A0,&I,&A)){
 		I = 90; A = 0;
 	}
@@ -34,7 +34,7 @@ int GM3D::ForwardHay(char* noise_level,char* mag_para){
 		Gamma=sin(I);
 	}
 
-	//添加高斯噪声值
+	//���Ӹ�˹����ֵ
 	default_random_engine generator;
 	normal_distribution<double> dist(noise_mean, noise_dev);
 
diff --git a/src/forward_vz.cpp b/src/forward_vz.cpp
index 25e85e1..725fe24 100644
--- a/src/forward_vz.cpp
+++ b/src/forward_vz.cpp
@@ -6,18 +6,18 @@ int GM3D::ForwardVz(char* noise_level){
 	double R222,R122,R212,R112,R221,R121,R211,R111;
 	double G222,G122,G212,G112,G221,G121,G211,G111;
 
-	//初始化正演值和不确定度
+	//��ʼ������ֵ�Ͳ�ȷ����
 	for (int i = 0; i < obs_num_; i++){
 		obs_p_[i].val = obs_p_[i].dev = 0.0;
 	}
 
-	//确定噪声水平
+	//ȷ������ˮƽ
 	double noise_mean, noise_dev, temp_noise;
 	if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
 		noise_mean = noise_dev = 0.0;
 	}
 
-	//添加高斯噪声值
+	//���Ӹ�˹����ֵ
 	default_random_engine generator;
 	normal_distribution<double> dist(noise_mean, noise_dev);
 
diff --git a/src/forward_vzx.cpp b/src/forward_vzx.cpp
index 35aaba7..0383633 100644
--- a/src/forward_vzx.cpp
+++ b/src/forward_vzx.cpp
@@ -6,7 +6,7 @@ int GM3D::ForwardVzx(char* noise_level){
 	double R222,R122,R212,R112,R221,R121,R211,R111;
 	double G222,G122,G212,G112,G221,G121,G211,G111;
 
-	//初始化正演值和不确定度
+	//��ʼ������ֵ�Ͳ�ȷ����
 	for (int i = 0; i < obs_num_; i++){
 		obs_p_[i].val = obs_p_[i].dev = 0.0;
 	}
@@ -16,7 +16,7 @@ int GM3D::ForwardVzx(char* noise_level){
 		noise_mean = noise_dev = 0.0;
 	}
 
-	//添加高斯噪声值
+	//���Ӹ�˹����ֵ
 	default_random_engine generator;
 	normal_distribution<double> dist(noise_mean, noise_dev);
 
diff --git a/src/forward_vzy.cpp b/src/forward_vzy.cpp
index d41552a..dd4c6df 100644
--- a/src/forward_vzy.cpp
+++ b/src/forward_vzy.cpp
@@ -6,7 +6,7 @@ int GM3D::ForwardVzy(char* noise_level){
 	double R222,R122,R212,R112,R221,R121,R211,R111;
 	double G222,G122,G212,G112,G221,G121,G211,G111;
 
-	//初始化正演值和不确定度
+	//��ʼ������ֵ�Ͳ�ȷ����
 	for (int i = 0; i < obs_num_; i++){
 		obs_p_[i].val = obs_p_[i].dev = 0.0;
 	}
@@ -16,7 +16,7 @@ int GM3D::ForwardVzy(char* noise_level){
 		noise_mean = noise_dev = 0.0;
 	}
 
-	//添加高斯噪声值
+	//���Ӹ�˹����ֵ
 	default_random_engine generator;
 	normal_distribution<double> dist(noise_mean, noise_dev);
 
diff --git a/src/forward_vzz.cpp b/src/forward_vzz.cpp
index e547e7e..6a24b10 100644
--- a/src/forward_vzz.cpp
+++ b/src/forward_vzz.cpp
@@ -6,7 +6,7 @@ int GM3D::ForwardVzz(char* noise_level){
 	double R222,R122,R212,R112,R221,R121,R211,R111;
 	double G222,G122,G212,G112,G221,G121,G211,G111;
 
-	//初始化正演值和不确定度
+	//��ʼ������ֵ�Ͳ�ȷ����
 	for (int i = 0; i < obs_num_; i++){
 		obs_p_[i].val = obs_p_[i].dev = 0.0;
 	}
@@ -16,7 +16,7 @@ int GM3D::ForwardVzz(char* noise_level){
 		noise_mean = noise_dev = 0.0;
 	}
 
-	//添加高斯噪声值
+	//���Ӹ�˹����ֵ
 	default_random_engine generator;
 	normal_distribution<double> dist(noise_mean, noise_dev);
 
diff --git a/src/forward_za.cpp b/src/forward_za.cpp
index d70ceaf..ca43305 100644
--- a/src/forward_za.cpp
+++ b/src/forward_za.cpp
@@ -10,18 +10,18 @@ int GM3D::ForwardZa(char* noise_level,char* mag_para){
 	double R222,R122,R212,R112,R221,R121,R211,R111;
 	double G222,G122,G212,G112,G221,G121,G211,G111;
 
-	//初始化正演值和不确定度
+	//��ʼ������ֵ�Ͳ�ȷ����
 	for (int i = 0; i < obs_num_; i++){
 		obs_p_[i].val = obs_p_[i].dev = 0.0;
 	}
 
-	//确定噪声水平
+	//ȷ������ˮƽ
 	double noise_mean, noise_dev, temp_noise;
 	if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
 		noise_mean = noise_dev = 0.0;
 	}
 
-	//确定磁化参数
+	//ȷ���Ż�����
 	if (4 != sscanf(mag_para,"%lf/%lf/%lf/%lf",&I0,&A0,&I,&A)){
 		I = 90; A = 0;
 	}
@@ -34,7 +34,7 @@ int GM3D::ForwardZa(char* noise_level,char* mag_para){
 		Gamma=sin(I);
 	}
 
-	//添加高斯噪声值
+	//���Ӹ�˹����ֵ
 	default_random_engine generator;
 	normal_distribution<double> dist(noise_mean, noise_dev);
 
diff --git a/src/gm3d.h b/src/gm3d.h
index db991c1..c14028f 100644
--- a/src/gm3d.h
+++ b/src/gm3d.h
@@ -7,23 +7,24 @@ class GM3D{
 public:
 	GM3D(){}
 	~GM3D(){}
-	int BuildRegularGrid(char*); //初始化反演模型空间
-	int AddModels(char*); //读取模型块体参数文件
-	void AddModels_GUI(); //添加模型块体参数
-	int AddRegularBlock(modelist); //添加普通模型块体
-	int AddTiltedBlock(modelist); //添加倾斜模型块体
-	int AddSphereBlock(modelist); //添加球体椭球体块体
-	int AddInterfaceBlock(modelist); //添加密度界面
-	//模型操作
+	int BuildRegularGrid(char*); //��ʼ������ģ�Ϳռ�
+	int AddModels(char*); //��ȡģ�Ϳ�������ļ�
+	void AddModels_GUI(); //����ģ�Ϳ������
+	int AddRegularBlock(modelist); //������ͨģ�Ϳ���
+	int AddTiltedBlock(modelist); //������бģ�Ϳ���
+	int AddSphereBlock(modelist); //�����������������
+	int AddInterfaceBlock(modelist); //�����ܶȽ���
+	//ģ�Ͳ���
 	int ReadModel(char*,char*);
-	//输出模型
-	int RegisteredOuput(bool); //注册输出的块体模型
-	int OutMshFile(char*,string); //输出模型文件
-	int OutNeighborFile(char*,char*); //输出模型块体或顶点的相邻关系 暂缓
-	//观测数据
+	int ReadModel_mst(char*);
+	//���ģ��
+	int RegisteredOuput(bool); //ע������Ŀ���ģ��
+	int OutMshFile(char*,string); //���ģ���ļ�
+	int OutNeighborFile(char*,char*); //���ģ�Ϳ���򶥵�����ڹ�ϵ �ݻ�
+	//�۲�����
 	int InitObs(char*);
 	int OutObs(char*);
-	//正演函数
+	//���ݺ���
 	int ForwardVz(char*);
 	int ForwardVzx(char*);
 	int ForwardVzy(char*);
@@ -35,16 +36,16 @@ public:
 	int ForwardHax(char*,char*);
 	int ForwardHay(char*,char*);
 	int ForwardZa(char*,char*);
-	//gm3d_gui函数
+	//gm3d_gui����
 	void get_model_list(modelistArray);
 private:
 	int obs_num_, model_num_, vert_num_;
-	//正演数组
+	//��������
 	obspointArray obs_p_;
 	_2dArray input_models_;
 	_1sArray input_model_names_;
 	_1dArray forward_model_;
-	//模型数据
+	//ģ������
 	cubeArray model_cube_;
 	cpointArray model_vert_;
 	_1dArray model_block_val_;
@@ -55,7 +56,7 @@ private:
 	_int2intMap vert_out_map_;
 	_int2intMap ele_data_out_map_;
 
-	_2iArray model_vert_neighbor_; //暂缓
-	_2iArray model_cube_neighbor_; //暂缓
+	_2iArray model_vert_neighbor_; //�ݻ�
+	_2iArray model_cube_neighbor_; //�ݻ�
 };
 #endif
\ No newline at end of file
diff --git a/src/gm3d_gui.cpp b/src/gm3d_gui.cpp
index c827ae6..8e6ad28 100644
--- a/src/gm3d_gui.cpp
+++ b/src/gm3d_gui.cpp
@@ -64,7 +64,7 @@ void cb_mesh_file_btn(Fl_Widget*, void*){
   // Block until user picks something.
   while(chooser.shown()) { Fl::wait(); }
 
-  //从文件的绝对路径中分离文件名并从输出框输出
+  //���ļ��ľ���·���з����ļ���������������
   if (chooser.value() != NULL)
   {
     strcpy(mesh_filename,chooser.value());
@@ -92,7 +92,7 @@ void cb_mod_para_file_btn(Fl_Widget*, void*){
   // Block until user picks something.
   while(chooser.shown()) { Fl::wait(); }
 
-  //下面我们读入文件并添加模型参数到列表显示
+  //�������Ƕ����ļ�������ģ�Ͳ������б���ʾ
   std::ifstream modpara_in;
   modpara_in.open(chooser.value());
   if (!modpara_in)
@@ -149,7 +149,7 @@ void cb_mod_file_out_btn(Fl_Widget*, void*){
   // Block until user picks something.
   while(chooser.shown()) { Fl::wait(); }
 
-  //从文件的绝对路径中分离文件名并从输出框输出
+  //���ļ��ľ���·���з����ļ���������������
   if (chooser.value() != NULL)
   {
     strcpy(out_msh_filename,chooser.value());
@@ -170,17 +170,17 @@ void cb_mod_file_out_btn(Fl_Widget*, void*){
 }
 
 void cb_build_mod_btn(Fl_Widget*, void*){
-  //首先构建一个GM3D实例
+  //���ȹ���һ��GM3Dʵ��
   GM3D gm3d_instance;
-  //设置初始化变量
+  //���ó�ʼ������
   char dimension[1024] = "10/20/990/10/20/990/10/20/490";
   char out_mshname[1024] = "Untitled.msh";
   char elename[1024] = "Untitled";
 
-  //将mesh_para_output拷贝到dimension
+  //��mesh_para_output������dimension
   strcpy(dimension, mesh_filename);
   //std::cout << dimension << std::endl;
-  //将mod_para_brw中的模型参数列表拷贝到一个临时参数列表
+  //��mod_para_brw�е�ģ�Ͳ����б�������һ����ʱ�����б�
   std::stringstream temp_ss;
   modelist temp_list;
   modelistArray brw_model_list;
@@ -196,21 +196,21 @@ void cb_build_mod_btn(Fl_Widget*, void*){
       }
   }
 
-  //将brw_model_list拷贝到GM3D中
+  //��brw_model_list������GM3D��
   gm3d_instance.get_model_list(brw_model_list);
-  //将mod_out_file_input拷贝到out_mshname
+  //��mod_out_file_input������out_mshname
   strcpy(out_mshname, out_msh_filename);
   if (!strcmp(out_mshname,"")){
     fl_message("Output file's name can't be empty !");
     return;
   }
-  //将mod_ele_input_build拷贝到elename
+  //��mod_ele_input_build������elename
   strcpy(elename, mod_ele_input_build->value());
   if (!strcmp(elename,"")){
     fl_message("Element data's name can't be empty !");
     return;
   }
-  //构建模型网络
+  //����ģ������
   if(gm3d_instance.BuildRegularGrid(dimension)){
     fl_message("Mesh Parameters Load Error !");
     return;
@@ -219,12 +219,12 @@ void cb_build_mod_btn(Fl_Widget*, void*){
   gm3d_instance.RegisteredOuput(remove_null);
   if (gm3d_instance.OutMshFile(out_mshname,elename)) return;
 
-  //这里需要一个输出信息窗口
+  //������Ҫһ�������Ϣ����
   fl_message("Model Construction Completed !");
   return;
 }
 
-/***********************************以下是正演计算GUI函数******************************/
+/***********************************���������ݼ���GUI����******************************/
 char in_msh_filename[1024];
 char in_obs_filename[1024];
 char out_res_filename[1024];
@@ -238,13 +238,13 @@ void cb_mod_file_input(Fl_Input*, void*){
 void cb_mod_file_btn(Fl_Widget*, void*){
   // Create the file chooser, and show it
   // directory, filter, chooser type, title
-  Fl_File_Chooser chooser(".", "*.msh", Fl_File_Chooser::SINGLE, "Create a model file.");
+  Fl_File_Chooser chooser(".", "*.msh, *.mst", Fl_File_Chooser::SINGLE, "Create a model file.");
   chooser.show();
 
   // Block until user picks something.
   while(chooser.shown()) { Fl::wait(); }
 
-  //从文件的绝对路径中分离文件名并从输出框输出
+  //���ļ��ľ���·���з����ļ���������������
   if (chooser.value() != NULL)
   {
     strcpy(in_msh_filename,chooser.value());
@@ -279,7 +279,7 @@ void cb_obs_file_btn(Fl_Widget*, void*){
   // Block until user picks something.
   while(chooser.shown()) { Fl::wait(); }
 
-  //从文件的绝对路径中分离文件名并从输出框输出
+  //���ļ��ľ���·���з����ļ���������������
   if (chooser.value() != NULL)
   {
     strcpy(in_obs_filename,chooser.value());
@@ -314,7 +314,7 @@ void cb_res_file_btn(Fl_Widget*, void*){
   // Block until user picks something.
   while(chooser.shown()) { Fl::wait(); }
 
-  //从文件的绝对路径中分离文件名并从输出框输出
+  //���ļ��ľ���·���з����ļ���������������
   if (chooser.value() != NULL)
   {
     strcpy(out_res_filename,chooser.value());
@@ -369,9 +369,9 @@ void cb_mag_data_check(Fl_Check_Button*, void*){
 }
 
 void cb_cal_btn(Fl_Button*, void*){
-  //首先构建一个GM3D实例
+  //���ȹ���һ��GM3Dʵ��
   GM3D gm3d_instance;
-  //设置初始化变量
+  //���ó�ʼ������
   char in_mshname[1024] = "Untitled.msh";
   char in_obspara[1024] = "Untitled.txt";
   char res_outfile[1024] = "Untitled";
@@ -386,11 +386,18 @@ void cb_cal_btn(Fl_Button*, void*){
   strcpy(ele_name,mod_ele_input->value());
   strcpy(noise_para,noise_para_input->value());
   strcpy(mag_para,mag_para_input->value());
-
-  if (gm3d_instance.ReadModel(in_mshname,ele_name))
-  {
-    fl_message("Model Load Error !");
-    return;
+  if(in_mshname[strlen(in_mshname)-1] == 'h'){
+    if (gm3d_instance.ReadModel(in_mshname,ele_name))
+    {
+      fl_message("Model Load Error !");
+      return;
+    }
+  }else if(in_mshname[strlen(in_mshname)-1] == 't'){
+    if (gm3d_instance.ReadModel_mst(in_mshname))
+    {
+      fl_message("Model Load Error !");
+      return;
+    }
   }
 
   if (gm3d_instance.InitObs(in_obspara))
@@ -533,6 +540,7 @@ void cb_cal_btn(Fl_Button*, void*){
       }
     }
   }
+  
 
   fl_message("Forward calculation completed !");
   return;
diff --git a/src/head_func.cpp b/src/head_func.cpp
index c754c77..0d65715 100644
--- a/src/head_func.cpp
+++ b/src/head_func.cpp
@@ -1,6 +1,6 @@
 #include "head_func.h"
-/*************************数据结构函数********************************/
-//cpoint减法
+/*************************���ݽṹ����********************************/
+//cpoint����
 cpoint operator -(cpoint a, cpoint b){
 	cpoint m;
 	m.x=a.x-b.x;
@@ -8,25 +8,25 @@ cpoint operator -(cpoint a, cpoint b){
 	m.z=a.z-b.z;
 	return m;
 }
-//cpoint模长
+//cpointģ��
 double modCpoint(cpoint v){
 	return sqrt(v.x*v.x+v.y*v.y+v.z*v.z);
 }
-/*************************全局函数********************************/
-//正负分离的atan函数 正数返回atan 负数返回atan+pi
+/*************************ȫ�ֺ���********************************/
+//���������atan���� ��������atan ��������atan+pi
 double arctg(double v){
 	double ang;
 	if(v>=0) ang=atan(v);
 	else if(v<0) ang=atan(v)+Pi;
 	return ang;
 }
-//将string转换为stringstream
+//��stringת��Ϊstringstream
 stringstream str2ss(string s){
 	stringstream sstr;
 	sstr.str(""); sstr.clear(); sstr.str(s);
 	return sstr;
 }
-//返回一个cpoint的位置id字符串
+//����һ��cpoint��λ��id�ַ���
 string cpoint_id(cpoint c){
 	string vert_id, mid_id;
 	stringstream sstemp;
@@ -40,7 +40,7 @@ string cpoint_id(cpoint c){
 	vert_id = vert_id + " " + mid_id;
 	return vert_id;
 }
-//测试打开输入文件 如果成功则返回0并输出信息 否则返回1
+//���Դ������ļ� ����ɹ��򷵻�0�������Ϣ ���򷵻�1
 int open_infile(ifstream &infile,char* filename){
 	infile.open(filename);
 	if (!infile){
@@ -49,7 +49,7 @@ int open_infile(ifstream &infile,char* filename){
 	}
 	return 0;
 }
-//测试打开输出文件 如果成功则返回0并输出信息 否则返回1
+//���Դ�����ļ� ����ɹ��򷵻�0�������Ϣ ���򷵻�1
 int open_outfile(ofstream &outfile,char* filename){
 	outfile.open(filename);
 	if (!outfile){
@@ -59,8 +59,8 @@ int open_outfile(ofstream &outfile,char* filename){
 	return 0;
 }
 
-//规则网络插值 长方形内数据插值 距离平方反比
-/*长方体示意图*/
+//���������ֵ �����������ݲ�ֵ ����ƽ������
+/*������ʾ��ͼ*/
 // y
 // |
 // |
@@ -69,10 +69,10 @@ int open_outfile(ofstream &outfile,char* filename){
 // |            |
 // |            |
 // 0------------1--->x
-// 左下角坐标x0 y0
-// 块体尺寸dx dy
-// 插值点坐标x y
-// 四个角点值
+// ���½�����x0 y0
+// ����ߴ�dx dy
+// ��ֵ������x y
+// �ĸ��ǵ�ֵ
 double grid_interpolate(double x0,double y0,double dx,double dy,double x,double y,
 					double d0,double d1,double d2,double d3)
 {
diff --git a/src/head_func.h b/src/head_func.h
index 79c9a94..7e22463 100644
--- a/src/head_func.h
+++ b/src/head_func.h
@@ -8,7 +8,7 @@
 #include "iomanip"
 #include "stdio.h"
 #include "stdlib.h"
-#include "unistd.h"
+//#include "unistd.h"
 #include "vector"
 #include "map"
 #include "algorithm"
@@ -18,42 +18,42 @@
 
 using namespace std;
 
-//数学常量
-#define BDL_MAX 1e+30 ///< 定义变量的最大值
-#define BDL_MIN -1e+30 ///< 定义变量的最小值
-#define PRECISION 16 ///< 定义小数点后需要使用的位数
-#define ZERO 1e-16 ///< 定义零值
+//��ѧ����
+#define BDL_MAX 1e+30 ///< ������������ֵ
+#define BDL_MIN -1e+30 ///< �����������Сֵ
+#define PRECISION 16 ///< ����С�������Ҫʹ�õ�λ��
+#define ZERO 1e-16 ///< ������ֵ
 
-//物理常量
-#define Pi (4.0*atan(1.0)) ///< 圆周率
-#define G0 6.67408e-3 ///< 万有引力常数。注意这里的量级本来应该是e-11，考虑到单位转换，取维度单位为m，密度单位为g/cm^3，乘以G0则重力单位即为mGal
-#define T0 5.0e+4 ///< 地磁场平均强度
-//宏函数
-#define MAX(a,b) (a>b?a:b) ///< 返回a与b的最大值
-#define MIN(a,b) (a<b?a:b) ///< 返回a与b的最小值
-#define SetToBox(a,b,in) (MAX(a,MIN(b,in))) ///< 返回a与b之间的一个值，若in在a与b之间则直接返回，否则返回较近的边界值
-//终端显示控制符
-#define BOLDRED "\033[1m\033[31m" ///< 设置后续字符字体为红色加粗
-#define BOLDGREEN "\033[1m\033[32m" ///< 设置后续字符字体为绿色加粗
-#define BOLDYELLOW "\033[1m\033[33m" ///< 设置后续字符字体为黄色加粗
-#define BOLDBLUE "\033[1m\033[34m" ///< 设置后续字符字体为蓝色加粗
-#define UNDERLINE "\033[1m\033[4m" ///< 设置后续字符为添加下划线
-#define RESET "\033[0m" ///< 重置字符设置
-#define MOVEUP(x) printf("\033[%dA", (x)) ///< 将光标向上挪x行
-#define MOVEDOWN(x) printf("\033[%dB", (x)) ///< 将光标向下娜x行
-#define MOVELEFT(x) printf("\033[%dD", (x)) ///< 将光标向左娜x字符
-#define MOVERIGHT(x) printf("\033[%dC", (x)) ///< 将光标向右娜x字符
-#define MOVETO(y,x) printf("\033[%d;%dH", (y), (x)) ///< 将光标向右娜动y字符,向上挪动x字符
-#define CLEARLINE "\033[K" ///< 清除本行
-#define CLEARALL "\033[2J" ///< 清除终端满屏
-//数据结构
-typedef vector<int> _1iArray; ///< 整形一维向量
-typedef vector<double> _1dArray; ///< 双精度浮点一维向量
-typedef vector<string> _1sArray; ///< 字符串一维向量
-typedef vector<vector<int> > _2iArray; ///< 整形浮点二维向量
-typedef vector<vector<double> > _2dArray; ///< 双精度浮点二维向量
-typedef map<int,int> _int2intMap; ///< 整型到整形的映射
-//模型块体参数
+//��������
+#define Pi (4.0*atan(1.0)) ///< Բ����
+#define G0 6.67408e-3 ///< ��������������ע���������������Ӧ����e-11�����ǵ���λת����ȡά�ȵ�λΪm���ܶȵ�λΪg/cm^3������G0��������λ��ΪmGal
+#define T0 5.0e+4 ///< �شų�ƽ��ǿ��
+//�꺯��
+#define MAX(a,b) (a>b?a:b) ///< ����a��b�����ֵ
+#define MIN(a,b) (a<b?a:b) ///< ����a��b����Сֵ
+#define SetToBox(a,b,in) (MAX(a,MIN(b,in))) ///< ����a��b֮���һ��ֵ����in��a��b֮����ֱ�ӷ��أ����򷵻ؽϽ��ı߽�ֵ
+//�ն���ʾ���Ʒ�
+#define BOLDRED "\033[1m\033[31m" ///< ���ú����ַ�����Ϊ��ɫ�Ӵ�
+#define BOLDGREEN "\033[1m\033[32m" ///< ���ú����ַ�����Ϊ��ɫ�Ӵ�
+#define BOLDYELLOW "\033[1m\033[33m" ///< ���ú����ַ�����Ϊ��ɫ�Ӵ�
+#define BOLDBLUE "\033[1m\033[34m" ///< ���ú����ַ�����Ϊ��ɫ�Ӵ�
+#define UNDERLINE "\033[1m\033[4m" ///< ���ú����ַ�Ϊ�����»���
+#define RESET "\033[0m" ///< �����ַ�����
+#define MOVEUP(x) printf("\033[%dA", (x)) ///< ���������Ųx��
+#define MOVEDOWN(x) printf("\033[%dB", (x)) ///< �����������x��
+#define MOVELEFT(x) printf("\033[%dD", (x)) ///< �����������x�ַ�
+#define MOVERIGHT(x) printf("\033[%dC", (x)) ///< �����������x�ַ�
+#define MOVETO(y,x) printf("\033[%d;%dH", (y), (x)) ///< ����������ȶ�y�ַ�,����Ų��x�ַ�
+#define CLEARLINE "\033[K" ///< �������
+#define CLEARALL "\033[2J" ///< ����ն�����
+//���ݽṹ
+typedef vector<int> _1iArray; ///< ����һά����
+typedef vector<double> _1dArray; ///< ˫���ȸ���һά����
+typedef vector<string> _1sArray; ///< �ַ���һά����
+typedef vector<vector<int> > _2iArray; ///< ����������
+typedef vector<vector<double> > _2dArray; ///< ˫���ȸ����ά����
+typedef map<int,int> _int2intMap; ///< ���͵����ε�ӳ��
+//ģ�Ϳ������
 struct modelist{
 	char mod_type[1024];
 	char val_type[1024];
@@ -61,33 +61,33 @@ struct modelist{
 	double mod_value;
 };
 typedef vector<modelist> modelistArray;
-//直角坐标系点
+//ֱ������ϵ��
 struct cpoint{
 	int id = -1;
 	double x = BDL_MAX; double y = BDL_MAX; double z = BDL_MAX;
 };
 typedef vector<cpoint> cpointArray;
 typedef map<string,cpoint> _str2pointMap;
-//观测点
+//�۲��
 struct obspoint : public cpoint{
 	double val = BDL_MAX; double dev = BDL_MAX;
 };
 typedef vector<obspoint> obspointArray;
-//块体
+//����
 struct cube{
 	cpoint cen;
 	int ids[8] = {-1,-1,-1,-1,-1,-1,-1,-1};
 	double dx = BDL_MAX; double dy = BDL_MAX; double dz = BDL_MAX;
 };
 typedef vector<cube> cubeArray;
-/*************************数据结构函数********************************/
-cpoint operator -(cpoint,cpoint); //矢量减法
-double modCpoint(cpoint); //矢量模
-/*************************全局函数********************************/
-double arctg(double); //正负分离的atan函数 正数返回atan 负数返回atan+pi
-stringstream str2ss(string); //将string转换为stringstream
-string cpoint_id(cpoint); //返回一个cpoint的位置id
-int open_infile(ifstream&,char*); //测试打开输入文件 如果成功则返回0并输出信息 否则返回1
-int open_outfile(ofstream&,char*); //测试打开输出文件 如果成功则返回0并输出信息 否则返回1
-double grid_interpolate(double,double,double,double,double,double,double,double,double,double); //规则网络插值
+/*************************���ݽṹ����********************************/
+cpoint operator -(cpoint,cpoint); //ʸ������
+double modCpoint(cpoint); //ʸ��ģ
+/*************************ȫ�ֺ���********************************/
+double arctg(double); //���������atan���� ��������atan ��������atan+pi
+stringstream str2ss(string); //��stringת��Ϊstringstream
+string cpoint_id(cpoint); //����һ��cpoint��λ��id
+int open_infile(ifstream&,char*); //���Դ������ļ� ����ɹ��򷵻�0�������Ϣ ���򷵻�1
+int open_outfile(ofstream&,char*); //���Դ�����ļ� ����ɹ��򷵻�0�������Ϣ ���򷵻�1
+double grid_interpolate(double,double,double,double,double,double,double,double,double,double); //���������ֵ
 #endif
\ No newline at end of file
diff --git a/src/init_obs.cpp b/src/init_obs.cpp
index 6d3e280..81a2006 100644
--- a/src/init_obs.cpp
+++ b/src/init_obs.cpp
@@ -8,7 +8,7 @@ int GM3D::InitObs(char* obs_para){
 	double xmin,xmax,ymin,ymax;
 	double xs,xe,ys,ye,eleva,dx,dy;
 
-	//按格式解析参数 初始化观测位置 用于正演计算
+	//����ʽ�������� ��ʼ���۲�λ�� �������ݼ���
 	if (7 == sscanf(obs_para,"%lf/%lf/%lf/%lf/%lf/%lf/%lf",&xs,&dx,&xe,&ys,&dy,&ye,&eleva)){
 		xmin = MIN(xs,xe); xmax = MAX(xs,xe);
 		ymin = MIN(ys,ye); ymax = MAX(ys,ye);
@@ -26,7 +26,7 @@ int GM3D::InitObs(char* obs_para){
 			x += dx;
 		}
 	}
-	//解析失败 按文件读入 用于反演使用或者正演计算
+	//����ʧ�� ���ļ����� ���ڷ���ʹ�û������ݼ���
 	else{
 		ifstream infile;
 		if (open_infile(infile,obs_para)) return -1;
@@ -34,7 +34,7 @@ int GM3D::InitObs(char* obs_para){
 		while(getline(infile,temp_str)){
 			if (*(temp_str.begin()) == '#') continue;
 			else{
-				//按每行3个数据解析 初始化为用于正演的观测点
+				//��ÿ��3�����ݽ��� ��ʼ��Ϊ�������ݵĹ۲��
 				if (3 == sscanf(temp_str.c_str(),"%lf %lf %lf",&temp_obs.y,&temp_obs.x,&temp_obs.z)){
 					temp_obs.z *= -1.0;
 					temp_obs.id = obs_p_.size();
diff --git a/src/out_msh_file.cpp b/src/out_msh_file.cpp
index 732682c..6c955ec 100644
--- a/src/out_msh_file.cpp
+++ b/src/out_msh_file.cpp
@@ -6,9 +6,9 @@ int GM3D::OutMshFile(char* filename,string data_name){
 	ofstream outfile;
 	if (open_outfile(outfile,filename)) return -1;
 
-	//好啦 我们这里输出的模型类型应该是块体
+	//���� �������������ģ������Ӧ���ǿ���
 	outfile<<"$MeshFormat"<<endl<<"2.2 0 8"<<endl<<"$EndMeshFormat"<<endl<<"$Nodes"<<endl<< out_vert_ids_.size() <<endl;
-	//输出顶点 最后一个不输出
+	//������� ���һ�������
 	for (int i = 0; i < out_vert_ids_.size(); i++){
 		outfile << i << " " 
 		<< setprecision(16) << model_vert_[out_vert_ids_[i]].x << " " 
diff --git a/src/out_neighbor_file.cpp b/src/out_neighbor_file.cpp
index 982edb8..4049c4d 100644
--- a/src/out_neighbor_file.cpp
+++ b/src/out_neighbor_file.cpp
@@ -3,22 +3,22 @@
 int GM3D::OutNeighborFile(char* v_name,char* b_name){
 	/*
 	if (strcmp(v_name,"NULL") || strcmp(b_name,"NULL")){
-		//整理块体间的相邻关系 先初始化顶点相邻数组
+		//�������������ڹ�ϵ �ȳ�ʼ��������������
 		model_vert_neighbor_.resize(vert_num_);
 		for (int i = 0; i < vert_num_; i++)
 			model_vert_neighbor_[i].resize(8,-1);
 
-		//遍历所有块体整理顶点相邻关系
+		//�������п��������������ڹ�ϵ
 		for (int i = 0; i < model_num_; i++){
 			for (int j = 0; j < 8; j++){
 				model_vert_neighbor_[model_cube_[i].ids[j]][j] = model_cube_[i].cen.id;
 			}
 		}
 
-		//遍历所有顶点相邻关系 所有共点的块体都被认定为相邻块体
+		//�������ж������ڹ�ϵ ���й���Ŀ��嶼���϶�Ϊ���ڿ���
 		model_cube_neighbor_.resize(model_num_);
 
-		//循环顶点相邻列表 都不为-1则相互添加
+		//ѭ�����������б� ����Ϊ-1���໥����
 		for (int i = 0; i < vert_num_; i++){
 			for (int n = 0; n < 8; n++){
 				if (model_vert_neighbor_[i][n] != -1){
@@ -32,14 +32,14 @@ int GM3D::OutNeighborFile(char* v_name,char* b_name){
 			}
 		}
 
-		vector<int>::iterator pos; //整型向量的迭代器
+		vector<int>::iterator pos; //���������ĵ�����
 		for (int i = 0; i < model_num_; i++){
-			sort(model_cube_neighbor_[i].begin(),model_cube_neighbor_[i].end()); //对顶点序列由小到大排序
-			pos = unique(model_cube_neighbor_[i].begin(),model_cube_neighbor_[i].end()); //获取重复序列开始的位置
-			model_cube_neighbor_[i].erase(pos,model_cube_neighbor_[i].end()); //删除重复点
+			sort(model_cube_neighbor_[i].begin(),model_cube_neighbor_[i].end()); //�Զ���������С��������
+			pos = unique(model_cube_neighbor_[i].begin(),model_cube_neighbor_[i].end()); //��ȡ�ظ����п�ʼ��λ��
+			model_cube_neighbor_[i].erase(pos,model_cube_neighbor_[i].end()); //ɾ���ظ���
 		}
 
-		//清理数组
+		//��������
 		for (int i = 0; i < vert_num_; i++){
 			model_vert_neighbor_[i].clear();
 			vector <int>().swap(model_vert_neighbor_[i]);
diff --git a/src/progress_bar.cpp b/src/progress_bar.cpp
index 8218d3b..277957a 100644
--- a/src/progress_bar.cpp
+++ b/src/progress_bar.cpp
@@ -47,7 +47,7 @@ int ProgressBar::GetConsoleWidth(){
         width = csbi.srWindow.Right - csbi.srWindow.Left;
     #else
         struct winsize win;
-        //注意！当我们使用pipe here-doc等通道获取程序参数时无法正确的获取窗口大小 此时我们将使用预定值
+        //ע�⣡������ʹ��pipe here-doc��ͨ����ȡ�������ʱ�޷���ȷ�Ļ�ȡ���ڴ�С ��ʱ���ǽ�ʹ��Ԥ��ֵ
         if (ioctl(0, TIOCGWINSZ, &win) != -1)
             width = win.ws_col;
         else width = 100;
diff --git a/src/progress_bar.h b/src/progress_bar.h
index 5dc7273..3e813c5 100644
--- a/src/progress_bar.h
+++ b/src/progress_bar.h
@@ -1,6 +1,6 @@
 #ifndef _PROGRESS_BAR_
 #define _PROGRESS_BAR_
-#include <sys/ioctl.h>
+//#include </*sys/ioctl.h*/>
 #include <iostream>
 #include <iomanip>
 #include <cstring>
diff --git a/src/read_model.cpp b/src/read_model.cpp
index 3b340c1..df3c8e2 100644
--- a/src/read_model.cpp
+++ b/src/read_model.cpp
@@ -10,14 +10,14 @@ int GM3D::ReadModel(char* filename,char* input_forward_model_name){
 	stringstream temp_ss;
 
 	ifstream mshin;
-	if (open_infile(mshin,filename)) return -1; //检查并打开模型文件
+	if (open_infile(mshin,filename)) return -1; //��鲢��ģ���ļ�
 	while(getline(mshin,temp_str)){
-		//读入模型空间顶点集 msh文件版本为2.2
+		//����ģ�Ϳռ䶥�㼯 msh�ļ��汾Ϊ2.2
 		if (temp_str == "$Nodes"){
 			getline(mshin,temp_str);
 			temp_ss = str2ss(temp_str);
-			temp_ss >> vert_num_; //第一个数为顶点的个数
-			model_vert_.resize(vert_num_); //开辟空间
+			temp_ss >> vert_num_; //��һ����Ϊ����ĸ���
+			model_vert_.resize(vert_num_); //���ٿռ�
 			for (int i = 0; i < vert_num_; i++){
 				getline(mshin,temp_str);
 				temp_ss = str2ss(temp_str);
@@ -25,19 +25,19 @@ int GM3D::ReadModel(char* filename,char* input_forward_model_name){
 				model_vert_[i] = temp_vert;
 			}
 		}
-		//读入模型空间单元体
+		//����ģ�Ϳռ䵥Ԫ��
 		else if (temp_str == "$Elements"){
 			getline(mshin,temp_str);
 			temp_ss = str2ss(temp_str);
-			temp_ss >> model_num_; //第一个数为总元素的个数 包含了所有类型的元素 比如三角形 四边形 块体等
+			temp_ss >> model_num_; //��һ����Ϊ��Ԫ�صĸ��� �������������͵�Ԫ�� ���������� �ı��� �����
 			model_cube_.resize(model_num_);
 			for (int i = 0; i < model_num_; i++){
 				getline(mshin,temp_str);
 				temp_ss = str2ss(temp_str);
 				temp_ss >> temp_cu.cen.id >> ele_type;
-				//只读入块体
+				//ֻ�������
 				if (ele_type == 5){
-					temp_ss >> attri_num; //跳过模型单元的几何组与物理组等信息 以后可能会有用
+					temp_ss >> attri_num; //����ģ�͵�Ԫ�ļ����������������Ϣ �Ժ���ܻ�����
 					for (int a = 0; a < attri_num; a++)
 						temp_ss >> temp_attri;
 					for (int a = 0; a < 8; a++)
@@ -46,44 +46,44 @@ int GM3D::ReadModel(char* filename,char* input_forward_model_name){
 				}
 			}
 		}
-		else continue; //不能识别的单元都被忽略了
+		else continue; //����ʶ��ĵ�Ԫ����������
 	}
 	mshin.close();
 
-	//第二次读入模型文件 初始化模型单元属性
-	if (open_infile(mshin,filename)) return -1; //检查并打开模型文件
+	//�ڶ��ζ���ģ���ļ� ��ʼ��ģ�͵�Ԫ����
+	if (open_infile(mshin,filename)) return -1; //��鲢��ģ���ļ�
 	while(getline(mshin,temp_str)){
-		//读入模型单元属性 注意因为msh文件中$ElementData并未注明所属元素类型
-		//所以可能会将其他元素类型的属性值也读入 但因为其在pyIdMap中并未注册 所以属性值会全为0 在后续使用时我们需要通过名称辨别
+		//����ģ�͵�Ԫ���� ע����Ϊmsh�ļ���$ElementData��δע������Ԫ������
+		//���Կ��ܻὫ����Ԫ�����͵�����ֵҲ���� ����Ϊ����pyIdMap�в�δע�� ��������ֵ��ȫΪ0 �ں���ʹ��ʱ������Ҫͨ�����Ʊ��
 		if (temp_str == "$ElementData"){
-			temp_model.resize(model_num_,0.0); //初始化temp_model 为读入模型单元属性做准备
-			for (int i = 0; i < 2; i++) //先读入元素块的名称 添加到数组
+			temp_model.resize(model_num_,0.0); //��ʼ��temp_model Ϊ����ģ�͵�Ԫ������׼��
+			for (int i = 0; i < 2; i++) //�ȶ���Ԫ�ؿ������ ���ӵ�����
 				getline(mshin,temp_str);
 			input_model_names_.push_back(temp_str);
-			for (int i = 0; i < 6; i++) //跳过元素属性前面的值 最后一次为当前元素块的个数
+			for (int i = 0; i < 6; i++) //����Ԫ������ǰ���ֵ ���һ��Ϊ��ǰԪ�ؿ�ĸ���
 				getline(mshin,temp_str);
 			temp_ss = str2ss(temp_str);
 			temp_ss >> temp_int;
 			for (int i = 0; i < temp_int; i++){
 				getline(mshin,temp_str);
 				temp_ss = str2ss(temp_str);
-				temp_ss >> temp_id >> temp_val; //读入单元体索引与属性值
+				temp_ss >> temp_id >> temp_val; //���뵥Ԫ������������ֵ
 				temp_model[temp_id] = temp_val;
 			}
 			input_models_.push_back(temp_model);
 			temp_model.clear();
 		}
-		else continue; //不能识别的单元都被忽略了
+		else continue; //����ʶ��ĵ�Ԫ����������
 	}
 	mshin.close();
 
-	//清理映射
+	//����ӳ��
 	temp_model.clear();
 	vector <double>().swap(temp_model);
 
-	//初始化数组
+	//��ʼ������
 	forward_model_.resize(model_num_,0.0);
-	//匹配数据名称
+	//ƥ����������
 	char forward_model_name[1024] = "\"";
 	strcat(forward_model_name,input_forward_model_name);
 	strcat(forward_model_name,"\"");
@@ -94,7 +94,7 @@ int GM3D::ReadModel(char* filename,char* input_forward_model_name){
 		}
 	}
 
-	//计算块体的中心位置和尺寸
+	//������������λ�úͳߴ�
 	cpoint corner[8];
 	for (int i = 0; i < model_num_; i++){
 		for (int j = 0; j < 8; j++){
diff --git a/src/read_model_mst.cpp b/src/read_model_mst.cpp
new file mode 100644
index 0000000..ee47e80
--- /dev/null
+++ b/src/read_model_mst.cpp
@@ -0,0 +1,51 @@
+#include "gm3d.h"
+
+int GM3D::ReadModel_mst(char* filename){
+	int temp_int,ele_type,attri_num,temp_attri,temp_id;
+	double temp_val;
+	_1dArray temp_model;
+	cpoint temp_vert;
+	cube temp_cu;
+	string temp_str;
+	stringstream temp_ss;
+
+	ifstream mstin;
+	if (open_infile(mstin,filename)) return -1; //��鲢��ģ���ļ�
+	double ynum, xnum, znum;
+    double ystart, xstart, zstart;
+    double ylength, xlength, zlength;
+    char temp;
+    mstin>>ynum>>xnum>>znum;
+    mstin>>ystart>>xstart>>zstart;
+    mstin>>ynum>>temp>>ylength;
+	mstin>>xnum>>temp>>xlength;
+	mstin>>znum>>temp>>zlength;
+
+    model_num_ = ynum*xnum*znum;
+    model_cube_.resize(model_num_);
+	//��ʼ������
+	forward_model_.resize(model_num_,0.0);
+
+    double f_mdl_tmp;
+
+    int i, j, k;
+    int tmp_num;
+    for(i = 0; i < xnum; i++){
+        for(j = 0; j < ynum; j++){
+            for(k = 0; k < znum; k++){
+                tmp_num = i*ynum*znum + j*znum + k;
+                mstin>>f_mdl_tmp;
+                forward_model_[tmp_num] = f_mdl_tmp;
+                model_cube_[tmp_num].cen.x = xstart + xlength*i + xlength*0.5;
+                model_cube_[tmp_num].cen.y = ystart + ylength*j + ylength*0.5;
+                model_cube_[tmp_num].cen.z = zstart + zlength*k + zlength*0.5;
+                model_cube_[tmp_num].dx = xlength;
+                model_cube_[tmp_num].dy = ylength;
+                model_cube_[tmp_num].dz = zlength;
+            }
+        }
+    }
+    mstin.close();
+
+	return 0;
+}
\ No newline at end of file
diff --git a/src/registered_output.cpp b/src/registered_output.cpp
index 52bc477..e34f83b 100644
--- a/src/registered_output.cpp
+++ b/src/registered_output.cpp
@@ -2,10 +2,10 @@
 
 int GM3D::RegisteredOuput(bool remove_empty_element){
 	int count;
-	//统计输出模型单元块体和顶点 以及输出的块体数据列表
+	//ͳ�����ģ�͵�Ԫ����Ͷ��� �Լ�����Ŀ��������б�
 	if (remove_empty_element){
 		count = 0;
-		//遍历所有块体数据 注册有值的块体
+		//�������п������� ע����ֵ�Ŀ���
 		for (int i = 0; i < model_num_; i++){
 			if (model_block_val_[i] != BDL_MAX){
 				out_ele_data_ids_.push_back(i);
@@ -15,38 +15,38 @@ int GM3D::RegisteredOuput(bool remove_empty_element){
 			}
 		}
 
-		//遍历所有注册的块体 添加顶点
+		//��������ע��Ŀ��� ���Ӷ���
 		for (int i = 0; i < out_ele_ids_.size(); i++){
 			for (int j = 0; j < 8; j++){
 				out_vert_ids_.push_back(model_cube_[out_ele_ids_[i]].ids[j]);
 			}
 		}
 
-		//去除输出顶点中的重复部分
-		vector<int>::iterator pos; //整型向量的迭代器
-		sort(out_vert_ids_.begin(),out_vert_ids_.end()); //对顶点序列由小到大排序
-		pos = unique(out_vert_ids_.begin(),out_vert_ids_.end()); //获取重复序列开始的位置
-		out_vert_ids_.erase(pos,out_vert_ids_.end()); //删除重复点
+		//ȥ����������е��ظ�����
+		vector<int>::iterator pos; //���������ĵ�����
+		sort(out_vert_ids_.begin(),out_vert_ids_.end()); //�Զ���������С��������
+		pos = unique(out_vert_ids_.begin(),out_vert_ids_.end()); //��ȡ�ظ����п�ʼ��λ��
+		out_vert_ids_.erase(pos,out_vert_ids_.end()); //ɾ���ظ���
 
-		//将需要输出的模型顶点序号与它们的排序做一个对应 保证在输出文件中顶点索引号始终是从0开始的连续的序列
+		//����Ҫ�����ģ�Ͷ�����������ǵ�������һ����Ӧ ��֤������ļ��ж���������ʼ���Ǵ�0��ʼ������������
 		for (int i = 0; i < out_vert_ids_.size(); i++){
 			vert_out_map_[out_vert_ids_[i]] = i;
 		}
 	}
 	else{
-		//输出的模型块体为所有
+		//�����ģ�Ϳ���Ϊ����
 		out_ele_ids_.resize(model_num_);
 		for (int i = 0; i < model_num_; i++){
 			out_ele_ids_[i] = i;
 		}
 
-		//输出的模型顶点为所有
+		//�����ģ�Ͷ���Ϊ����
 		out_vert_ids_.resize(vert_num_);
 		for (int i = 0; i < vert_num_; i++){
 			vert_out_map_[i] = out_vert_ids_[i] = i;
 		}
 
-		//注册所有有值的块体数据
+		//ע��������ֵ�Ŀ�������
 		count = 0;
 		for (int i = 0; i < model_num_; i++){
 			if (model_block_val_[i] != BDL_MAX){
@@ -56,6 +56,6 @@ int GM3D::RegisteredOuput(bool remove_empty_element){
 			}
 		}
 	}
-	//输出所有模型单元块体和顶点 只统计输出的块体数据列表
+	//�������ģ�͵�Ԫ����Ͷ��� ֻͳ������Ŀ��������б�
 	return 0;
 }
\ No newline at end of file