代码可用于若干直立六面体重力和磁异常正演,输入观测面(平面或者曲面)和场源文件,输出网格化重力异常,可用