一维可控源电磁法模型计算程序与变量分析

! Program to calculated Electric Dipole's all Resistivy.
PROGRAM main
character*80 cmdfile
character*80 oldfile
cmdfile='cmd.txt'
open(10,file=cmdfile)
do 
    read (10,*,end=100) oldfile
    call femf(oldfile)
end do
100   end program

! 参数定义

      subroutine femf(oldfile)
CHARACTER OUT1*80, oldfile*80
COMPLEX KK1                      !波数的平方
DIMENSION FF(1000),XS(100)       !ff―频率，XS―点的坐标
  COMMON /PARA/RL,RO,SIN1,COS1,KK1 !
COMMON /CONS/PI,UO,PE            !UO―磁导率，PE―电偶矩
COMMON /LAYER/NL,RD(10),RS(10)   !NL―层数，RD―层厚度，RS―层电阻率
COMMON /HANK/CH0(112),CH1(112)   !汉克尔系数
OPEN(2,FILE='CHYH.DAT',STATUS='OLD')
READ(2,*)(CH0(I),I=1,112)
READ(2,*)(CH1(I),I=1,112)
CLOSE(2)
C  OPEN(23,FILE='CH0112.TXT',STATUS='UNKNOWN')
C  DO I = 1, 112
C  WRITE(23,*)CH0(I)
C  ENDDO

C  READ(2,*)(CH0(I),I=1,112)
C  READ(2,*)(CH1(I),I=1,112)
CLOSE(2)
PI=3.141592653
UO=4.*PI*1.E-7
C 输入数据
! write(*,*)'请输入模型参数文件：'
! read(*,*)oldfile
out1='res_'//trim(oldfile)
OPEN(2,FILE=oldfile,STATUS='OLD')
READ(2,*) AB,NAB,CI              !读入AB长度，NAB分段数，CI电流
PE=AB*CI/2./PI/NAB               !给电偶矩赋值
C 输入频率 
C &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
 READ(2,*) KFRQ 
IF(KFRQ.EQ.1) THEN
 READ(2,*) NF,FMAX,DELTAF
 FF(1)=FMAX
END IF
IF(KFRQ.EQ.2) THEN
 READ(2,*) NF,FMAX,FMIN
 DELTAF=(ALOG10(FMAX)-ALOG10(FMIN))/(NF-1)
 DELTAF=10**DELTAF
 FF(1)=FMIN
END IF
IF(KFRQ.EQ.3) THEN
 OPEN(3,FILE='FREQz.txt',STATUS='OLD')
 READ(3,*) NF
 READ(3,*)(FF(I),I=1,NF)
 CLOSE(3)