1dforward.for ! ldforward.for主程序 real ps(30),r(20),p(30),h(9),sd(9) ! ps(30)电测深视电阻率；r(30)极距；p(10)层电阻率；h(9)层厚度，sd(9)深度 character*50 char,name1,name2 integer nalyer,max write(*,*)"请输入模型文件名:" read(*,*)name1 open(2,file=name1) write(*,*)"请输入视电阻率文件名:" read(*,*)name2 open(1,file=name2) read(2,*) char read(2,*)nlayer read(2,*)char do i=1,nlayer read(2,*)p(i) enddo read(2,*)char do i=1,nlayer-1 read(2,*)h(i) enddo !!!!!!!!!!!!!!!!!!!!!!!计算电测深极距共20个，每个数量级5个极距 max=20 do i=1,max r(i)=10.* *(i/5.) enddo !调用一维正演子程序 zyld() call zyld(p,h,nlayer,r ,max,ps) !输出电测深曲线到文件 do i=1,max write(1,*)r(i),ps(i) enddo !!!!!!!!!!!!!!!画模型 sd(1)=0. do i=2,nlayer sd(i)=sd(i-1)+h(i-1) enddo sd(nlayer+1)=10000. open(25,file="model.dat")!画模型的文件 write（25，*)1，p(1) write(25,*)h(1),p(1) do i=2,nlayer write(25,*)sd(i),p(i) write(25,*)sd(i+1),p(i) enddo write(*,*)"模型曲线的文件为model.dat" pause end !!!!!!!!!!!一维正演子程序 subroutine zyld(p,h,nl,r,max,ps) real ps(max),t(nl),ck(20),r(max),m real p(nl),h(nl-1) integer nl,max !max极距个数；nl层数；ps(max)电测深视电阻率；r(max)极距；p(nl)层电阻率； !h(nl-1)层厚度；t(nl)电阻率转换函数；m汉克尔变换积分变量；ck(20)滤波系数； data ck/0.003042,-0.001198,0.01284,0.0235,0.08688,0.2374,0.6194, &. 1.1817,0.4248,-3.4507,2.7044,-1.1324,0.393,-0.1436,0.05812, &. -0.02521,0.01125,-0.004978,0.002072,-0.000318/ t(nl)=p(nl) do i=1,max ps(i)=0. do k=1,20 m=0.11396*10.* *(k/6.)/r(i) do j=nl-1,1,-1 t(j)=p(j)*((t(j+1)+p(j)+(t(j+1)-p(j))*exp(-2.*m*h(j)))& /((t(j+1)+p(j))-(t(j+1)-p(j))*exp((-2.*m*h(j))) enddo ps（i）=ps(i)+t(1)*ck(k) enddo enddo end

Kerry key可控源有限元二维正反演源程序和编译方法，大家可以结合我在csdn的安装教程和使用教程来操作，可以快速入门和使用该可控源程序进行自己所需要的二维可控源或者大地电磁正反演。

大地电磁各种处理的python开源软件，我会在csdn写个详细的安装和使用教程，供大家学习和参考，以便大家能够快速入门和使用，减少大家软件的学习使用成本。

MT一维反演Bostick方法Matlab程序，分为两个Matlaab子程序，一个用于MT一维正演计算，另一个是Bostick反演程序。打开Bostick程序，点击运行，即可自动反演。修改程序中数据，即可完成MT一维反演问题。

大地电磁测深一维正演程序，能够帮助学生初步认识大地电磁，继续研究

程序简单易懂，对于学习模拟退火算法来说简单易懂，方便入门

该程序用于计算水平层状介质大地电磁视电阻率响应。

MT1Dforward 一维正演的FORTRAN程序

一种很巧妙的bostick反演方法. 利用李金铭老师的电法书中的介绍，使用希尔伯特变换做了比较巧妙的处理。

大地电磁MT二维正演模拟软件 含使用教程和例子数据，陈小斌开发。