zhangyiss/DSurfTomo
1
0
Fork 0
mirror of https://github.com/HongjianFang/DSurfTomo.git synced 2025-05-08 00:01:14 +08:00
Code Issues Actions Packages Projects Releases Wiki Activity

Merge branch 'stable' of https://github.com/HongjianFang/DSurfTomo into stable

Browse Source
This commit is contained in:
Hongjian Fang 2021-02-17 08:43:03 +08:00
commit f35c3d77da
3 changed files with 11 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/CalSurfG.f90
View File

@ -1108,7 +1108,7 @@ subroutine CalSurfG(nx,ny,nz,nparpi,vels,iw,rw,col,dsurf, &
igr=0 igr=0
! print*,kmax ! print*,kmax
call caldespersion(nx,ny,nz,vels,pvRc, & call caldespersion(nx,ny,nz,vels,pvRc, &
iwave,igr,kmax,tRg,depz,minthk) iwave,igr,kmaxRg,tRg,depz,minthk)
igr=1 igr=1
print*,'Rayleigh wave group velocity depth kernel' print*,'Rayleigh wave group velocity depth kernel'
call depthkernel(nx,ny,nz,vels,pvRg,sen_vsRg,sen_vpRg, & call depthkernel(nx,ny,nz,vels,pvRg,sen_vsRg,sen_vpRg, &
@ -1126,7 +1126,7 @@ subroutine CalSurfG(nx,ny,nz,nparpi,vels,iw,rw,col,dsurf, &
iwave=1 iwave=1
igr=0 igr=0
call caldespersion(nx,ny,nz,vels,pvLc, & call caldespersion(nx,ny,nz,vels,pvLc, &
iwave,igr,kmax,tLg,depz,minthk) iwave,igr,kmaxLg,tLg,depz,minthk)
igr=1 igr=1
call depthkernel(nx,ny,nz,vels,pvLg,sen_vsLg,sen_vpLg, & call depthkernel(nx,ny,nz,vels,pvLg,sen_vsLg,sen_vpLg, &
sen_rhoLg,iwave,igr,kmaxLg,tLg,depz,minthk) sen_rhoLg,iwave,igr,kmaxLg,tLg,depz,minthk)

16
src/main.f90
View File

@ -174,9 +174,9 @@ program SurfTomo
if (checkstat > 0) stop 'error allocating RP' if (checkstat > 0) stop 'error allocating RP'
read(10,*)(tRc(i),i=1,kmaxRc) read(10,*)(tRc(i),i=1,kmaxRc)
write(*,*)'Rayleigh wave phase velocity used,periods:(s)' write(*,*)'Rayleigh wave phase velocity used,periods:(s)'
write(*,'(50f7.1)')(tRc(i),i=1,kmaxRc) write(*,'(50f7.2)')(tRc(i),i=1,kmaxRc)
write(66,*)'Rayleigh wave phase velocity used,periods:(s)' write(66,*)'Rayleigh wave phase velocity used,periods:(s)'
write(66,'(50f7.1)')(tRc(i),i=1,kmaxRc) write(66,'(50f7.2)')(tRc(i),i=1,kmaxRc)
endif endif
read(10,*)kmaxRg read(10,*)kmaxRg
if(kmaxRg.gt.0)then if(kmaxRg.gt.0)then
@ -184,9 +184,9 @@ program SurfTomo
if (checkstat > 0) stop 'error allocating RP' if (checkstat > 0) stop 'error allocating RP'
read(10,*)(tRg(i),i=1,kmaxRg) read(10,*)(tRg(i),i=1,kmaxRg)
write(*,*)'Rayleigh wave group velocity used,periods:(s)' write(*,*)'Rayleigh wave group velocity used,periods:(s)'
write(*,'(50f7.1)')(tRg(i),i=1,kmaxRg) write(*,'(50f7.2)')(tRg(i),i=1,kmaxRg)
write(66,*)'Rayleigh wave group velocity used,periods:(s)' write(66,*)'Rayleigh wave group velocity used,periods:(s)'
write(66,'(50f7.1)')(tRg(i),i=1,kmaxRg) write(66,'(50f7.2)')(tRg(i),i=1,kmaxRg)
endif endif
read(10,*)kmaxLc read(10,*)kmaxLc
if(kmaxLc.gt.0)then if(kmaxLc.gt.0)then
@ -194,9 +194,9 @@ program SurfTomo
if (checkstat > 0) stop 'error allocating RP' if (checkstat > 0) stop 'error allocating RP'
read(10,*)(tLc(i),i=1,kmaxLc) read(10,*)(tLc(i),i=1,kmaxLc)
write(*,*)'Love wave phase velocity used,periods:(s)' write(*,*)'Love wave phase velocity used,periods:(s)'
write(*,'(50f7.1)')(tLc(i),i=1,kmaxLc) write(*,'(50f7.2)')(tLc(i),i=1,kmaxLc)
write(66,*)'Love wave phase velocity used,periods:(s)' write(66,*)'Love wave phase velocity used,periods:(s)'
write(66,'(50f7.1)')(tLc(i),i=1,kmaxLc) write(66,'(50f7.2)')(tLc(i),i=1,kmaxLc)
endif endif
read(10,*)kmaxLg read(10,*)kmaxLg
if(kmaxLg.gt.0)then if(kmaxLg.gt.0)then
@ -204,9 +204,9 @@ program SurfTomo
if (checkstat > 0) stop 'error allocating RP' if (checkstat > 0) stop 'error allocating RP'
read(10,*)(tLg(i),i=1,kmaxLg) read(10,*)(tLg(i),i=1,kmaxLg)
write(*,*)'Love wave group velocity used,periods:(s)' write(*,*)'Love wave group velocity used,periods:(s)'
write(*,'(50f7.1)')(tLg(i),i=1,kmaxLg) write(*,'(50f7.2)')(tLg(i),i=1,kmaxLg)
write(66,*)'Love wave group velocity used,periods:(s)' write(66,*)'Love wave group velocity used,periods:(s)'
write(66,'(50f7.1)')(tLg(i),i=1,kmaxLg) write(66,'(50f7.2)')(tLg(i),i=1,kmaxLg)
endif endif
read(10,*)ifsyn read(10,*)ifsyn
read(10,*)noiselevel read(10,*)noiselevel

2
src/surfdisp96.f
View File

@ -308,7 +308,7 @@ c ----- print *, itst,iq,t(k),t1a,t1b,cc0,cc1,gvel
if(iverb(ifunc).eq.0)then if(iverb(ifunc).eq.0)then
iverb(ifunc) = 1 iverb(ifunc) = 1
write(LOT,*)'improper initial value in disper - no zero found' write(LOT,*)'improper initial value in disper - no zero found'
write(*,*)'WARNING:improper initial value in disper - no zero found' !write(*,*)'WARNING:improper initial value in disper - no zero found'
write(LOT,*)'in fundamental mode ' write(LOT,*)'in fundamental mode '
write(LOT,*)'This may be due to low velocity zone ' write(LOT,*)'This may be due to low velocity zone '
write(LOT,*)'causing reverse phase velocity dispersion, ' write(LOT,*)'causing reverse phase velocity dispersion, '