wind.for

Fortran编的波动方程正演模拟,一维的数值算法

      subroutine wcos(w,nw,taper)
	real w(0:nw-1)
	pi=4*atan(1.0)

	do i=0,nw-1
	w(i)=1
	enddo

	m=int( (taper/100)*(nw/2) )
	do i=0,m-1
	w(i)=.5*(1-cos((i+1)*pi/(M+1)))
	enddo

	do i=nw-m-1,nw-1
	w(i)= .5*(1-cos((nw-i)*pi/(M+1)) )
	enddo 
	WRITE(*,*)w
	end


	subroutine wgao(w,nw)
	real w(0:nw-1)
	pi=4*atan(1.0)

	dlt=nw/8
	do i=0,nw-1
	w(i)=exp(-0.5*(i-nw/2)**2/dlt**2)
	enddo
	WRITE(*,*)w
	end