c calculate source depth index

	   do iz=1,nz-1
                depth1=rmax-vr2(iz)
                depth2=rmax-vr2(iz+1)
                if(sdepth.gt.depth1.and.sdepth.lt.depth2)then
                        if(sdepth-depth1.le.depth2-sdepth)then
                        isz=iz
                	else
                        isz=iz+1
                	endif
		endif
           enddo
