capjoint中的tel3核心代码teleseis3.f90

为了加入更多层的模型 将

teleseis3.f90  /home/capjoint-master/src/tel3/teleseis3.90的地层模型读取部分改为：

program test

      PARAMETER (ND0=2048,NM0=6,LK0=10,NL0=100,PI=3.141593,RAD=.0174533)
      IMPLICIT COMPLEX*8 (Z)
      CHARACTER NAME*40,NAM*4, prefix*80, adel*3, outfile*80
      character*1 wv(9),cm(9), model*80, astr*4,adip*2,arak*4
      character*1 greenV(9)
      character   sta*10
      real*8 gg(500),pp(500),depth,vp00,gs(500),ps(500),vs00
      COMMON /STR0/NL ,VP(NL0),VS(NL0),DEN(NL0),DEP(NL0)
      COMMON /STR1/NL1,VP1(NL0),VS1(NL0),DEN1(NL0),DEP1(NL0)
      COMMON /STR2/NL2,VP2(NL0),VS2(NL0),DEN2(NL0),DEP2(NL0)
 common /sourceRegion/vsrc
      print*, "Please input the model file name:"
      read(*,'(a80)') model
      if(model(1:4) .ne. 'none') then
  print*, model
  open(2,file=model)
         READ(2,'(a40)') name
         print*, name
         READ(2,*) TQP,TQS,NL
         do L = 1, NL
            read(2,*) VP(L),VS(L),DEN(L),DEP(L)
         enddo

         READ(2,*) NL1
         do L = 1, NL1
         read(2,*)VP1(L),VS1(L),DEN1(L),DEP1(L)
         enddo
          READ(2,*)NL2
          do L = 1, NL2
          read(2,*)VP2(L),VS2(L),DEN2(L),DEP2(L)
          enddo
      endif
      print*, "reading down!!"
 end program test

重点部分是标为橘色 的部分，但tel3主要用于远震的格林函数计算，它需要与fk计算的地震波进行比对。

转载于:https://www.cnblogs.com/seisjun/p/7815198.html