anowacki / sac.jl Goto Github PK

function rotate_to_gcp!(s1::SACtr, s2::SACtr, reverse::Bool=false)
    s2_is_clockwise_of_s1 = angle_difference(s1[:cmpaz], s2[:cmpaz]) > 0
    s2_is_clockwise_of_s1 || ((s1, s2) = (s2, s1))
    !isundefined(s1[:baz]) && !isundefined(s2[:baz]) ||
        throw(ArgumentError("Backazimuth is not defined for both traces"))
    s1[:baz] ≈ s2[:baz] || throw(ArgumentError("Backazimuth not the same for both traces"))
    phi = mod(s1[:baz] + 180 - s1[:cmpaz], 360)
    rotate_through!(s2, s1, phi) # Checks for orthogonality
    reverse && flip_component!(s2)
    s1[:kcmpnm] = sacstring("Radial")
    s2[:kcmpnm] = sacstring((reverse ? "-" : "")*"Trans")
    s1, s2
end

    s2_is_clockwise_of_s1 || ((s1, s2) = (s2, s1))

function tshift!(s::SACtr, tshift::Number, wrap=true)
    @info wrap,tshift
    n = round(Int, tshift/s.delta)
    @info wrap,n,tshift
    if n == 0
        sac_verbose && @info("SAC.tshift!: t ($tshift) is less than delta ($(s.delta)) so no shift applied")
        return
    end
    @info wrap,n,tshift,s.t[1:n]
    s.t = circshift(s.t, n)
    @info wrap,n,tshift,s.t[1:n]
    if !wrap
        n > 0 ? s.t[1:n] = 0. : s.t[end+n+1:end] = 0.
    end
    update_headers!(s)
end

s.t[1:n] .= 0.

s.t[1:n]=0.