直立六面体正演实验1.docVIP

实验四 直立六面体正演实验 姓名： 学号： 专业： 指导老师： 一、基本原理 输入/输出数据格式设计 输入数据: 场源文件 XYZ输入文件 读入各文件名 sourcefilename.txt inputfilename.txt CMD.txt 输出文件: gravityfilename.grd (重力异常输出) 重要变量说明 density:剩余密度 Jmay:磁化强度 Ai_source:磁化方向倾角 Ax_source:磁化方向与X轴的夹角 Ay_source:磁化方向与Y轴的倾角 X_source:场源中的X坐标 Y_source:场源中的Y坐标 Z_source:场源中的Z坐标 mpoint:点数 nline:条数 三、总体设计 测试结果 重力异常： 磁力异常： 结论 由已有的重力位表达式，借助于联系重力和磁场的泊松公式，直接给出磁位，磁场表达式，从而得出结果。 附录:源程序代码: PROGRAM forward_main PARAMETER(number_source=3) CHARACTER*(80)sourcefilename,inputfilename,gravityfilename REAL Xmin,Xmax,Ymin,Ymax ,Zmin,Zmax INTEGER mpoint,nline REAL,ALLOCATABLE::density(:),Jmay(:),Ai_source(:),Ax_source(:), * Ay_source(:), X_source(:,:),Y_source(:,:),Z_source(:,:) REAL,ALLOCATABLE::DG(:,:) CALL read_cmd(sourcefilename,inputfilename,gravityfilename) CALL get_number_source(sourcefilename,number_source) ALLOCATE(density(1:number_source),Jmay(1:number_source), *Ai_source(1:number_source),Ax_source(1:number_source), *Ay_source(1:number_source),X_source(1:2,1:number_source) Y_source(1:2,1:number_source),Z_source(1:2,1:number_source)) CALL read_source(sourcefilename,number_source,density,Jmay, *Ai_source,Ax_source,Ay_source,X_source,Y_source,Z_source) CALL input_xyz(inputfilename,mpoint,nline,Xmin,Xmax, *Ymin,Ymax,Zmin,Zmax) ALLOCATE(DG(1:mpoint,1:nline)) CALL gravity_sub(number_source,density,X_source,Y_source, * Z_source,mpoint,nline,Xmin,Xmax,Ymin,Ymax,Zmin,Zm