I used the following code to retrieve data. If some stations in my study area have been equipped with accelerographs or seismographs, it does affect the output. I wrote “output velocity” because I need velocity data; however, I know that some stations are accelerographs. For example, the below station is an accelerograph, so the output has to be velocity or acceleration?
client = Client(‘NCEDC’)
t = UTCDateTime(“2010-01-10T00:27:31”)
st = client.get_waveforms(“BK”, “JCC”, “”, “*H?”, t, t + 3 * 60, attach_response=True)
st.detrend(type=“demean”)
st.detrend(type=“linear”)
st.detrend(type=“polynomial”, order=2)
st.remove_response(output=‘VEL’, water_level=60)
st.filter(type=“bandpass”, freqmin=0.1, freqmax=5)
st.write(“Number17”, format=“TSPAIR”, Units=‘m/s’)
st.plot()
And how do I determine the type of NHERP(Vs30) of stations ?
Thanks alot