直立六面体正演程序设计实验汇编.docVIP

地球物理重磁勘探实验 实验名称：直立六面体正演程序设计实验 学院名称： 专业名称： 姓名： 学生学号： 指导老师： 一、基本原理 （1）模型示意及有关变量的意义，并分别给出重力异常和磁力异常的计算公式，化极。  （2）方向描述及磁化强度分量计算 （3）计算点输入坐标格式设计 (4)计算点输入坐标格式设计 三，输入/输出数据格式设计 （1）参数 场源系数设计：density，Jmay，ai_source,ax_source,ay_source,x1_source,x2_source,y1_source,y2_source,z1_source,z2_source 场源参数文件：sourcefilename 输入平面规则网文件：Inputfilename 重力异常输出文件：gravityfilename 磁力异常输出文件：magnaticfilename 化极磁力异常：RTPfilename 地磁的倾角：ai，ax，ay 点数线数：mpoint，nline 平面规则网数据：xmin，xmax，ymin，ymax，zmin，zmax （2）类型 曲面规则网：输出用GRD格式 平面规则网：输出也用GRD格式 曲面非规则网：输出采用DAT格式，形式如下： 平面非规则网 I P O 重力异常，磁力异常，化极磁力异常由surfer得到的图形 重力异常： 磁力异常： 化极磁力异常： 附录：源程序代码 !***********************************************************************(主程序) program Forward_Recfangle character*(80) sourcefilename,Inputfilename,gravityfilename,magneticfilename,huajifilename real xmin,xmax,ymin,ymax,zmin,zmax integer number_source,mpoint,nline real,allocatable::density(:),X_source(:,:),Y_source(:,:),Z_source(:,:),Jmag(:),ai_source(:),ax_source(:),ay_source(:),AI_S(:),AX_S(:),AY_S(:),AI_SS(:),AX_SS(:),AY_SS(:),JMAY(:) real,allocatable::Dg(:,:),mg(:,:),HG(:,:) call Read_cmd(Inputfilename,sourcefilename,gravityfilename,magneticfilename,huajifilename) call get_number(sourcefilename,number_source) allocate (density(1:number_source),x_source(1:2,1:number_source),Y_source(1:2,1:number_source),Z_source(1:2,1:number_source),jmag(1:number_source),ai_source(1:number_source),ax_source(1:number_source),ay_source(1:number_source),AI_S(1:number_source),AX_S(1:number_source),AY_S(1:number_source),AI_SS(1:number