重磁异常正演fortran代码

代码说明：

对给定的密度体和磁性体计算起伏观测面以及平面上的重力异常、磁力异常以及化极磁力异常, 场源参数保存在“source.dat”中。第一列为剩余密度(g/cm3)；第二列为磁化强度(10-6CGSM)；第三列为磁化方向倾角(DEG)；第四列为磁化方向与x轴的夹角(DEG)；第五列为磁化方向与y轴的夹角(DEG)；第六列～第七列为x坐标的起点和终点(km)；第八列～第九列为y坐标的起点和终点(km)；第十列～第十一列为z坐标的起点和终点(km，向下为正)。