Wood Anderson response simulation in obspy.signal.invsim

Hi there,

I wanted to clarify a few things in the estimate_wood_anderson_amplitude_using_response function in obspy.signal.invsim.
In the function, the response is removed with an output of velocity, and then the poles and zeroes of the Wood Anderson instrument are applied to simulate WA response:

response = response.get_evalresp_response_for_frequencies(
[freq], output=“VEL”, start_stage=None, end_stage=None)[0]
response_amplitude = np.absolute(response)
wa_ampl /= response_amplitude
wa_ampl *= paz_2_amplitude_value_of_freq_resp(WOODANDERSON, freq) *

Where WOODANDERSON is defined as:
WOODANDERSON = {‘poles’: [-6.283 + 4.7124j, -6.283 - 4.7124j],
‘zeros’: [0 + 0j], ‘gain’: 1.0, ‘sensitivity’: 2080}

However, these poles and zeroes seem to come from New Manual of Seismological Observatory Practice (2002). I can’t find that exact version, but in the revised version published in 2009 they state on page 570:

Note that the poles and zeros given in Tabs. 11.3 and 11.4 are valid only if the input signal to
the considered seismographs is ground displacement (amplitude Ad).

Where Table 11.3 contains the poles and zeroes for WA.

I’m confused here - perhaps this is just me misunderstanding what they wrote - but doesn’t this suggest that the input to the wood anderson simulation should be displacement and not velocity?
Why do we use the velocity in obspy?
Any help you can give me would be much appreciated.

Edit: I note now that the obspy response has one less zero than the response in the paper. I assume this is allowing us to translate directly from velocity into displacement?