Skip to content

Instantly share code, notes, and snippets.

@dsentinel

dsentinel/bug_sac_pz.ipynb

Last active July 24, 2018 17:03
Show Gist options
  • Save dsentinel/f8aeec00f42afc9e32e552ae2a8ebbef to your computer and use it in GitHub Desktop.
Save dsentinel/f8aeec00f42afc9e32e552ae2a8ebbef to your computer and use it in GitHub Desktop.
Download ZIP
Display the source blob
Display the rendered blob
Raw
bug_sac_pz.ipynb
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"from obspy.clients.fdsn import Client\n",
"import io"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"client = Client('IRIS')"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"inv = client.get_stations(network='TA',\n",
" station='034A',\n",
" level='response',\n",
" )"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"scrolled": true
},
"outputs": [],
"source": [
"f = io.StringIO()"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"TA.034A.ACE has no paz. Skipping.\n",
"TA.034A.LCE has no paz. Skipping.\n",
"TA.034A.LCQ has no paz. Skipping.\n",
"TA.034A.LOG has no paz. Skipping.\n",
"TA.034A.OCF has no paz. Skipping.\n",
"TA.034A.VCO has no paz. Skipping.\n",
"TA.034A.VEA has no paz. Skipping.\n",
"TA.034A.VEC has no paz. Skipping.\n",
"TA.034A.VEP has no paz. Skipping.\n",
"TA.034A.VKI has no paz. Skipping.\n",
"TA.034A.VM0 has no paz. Skipping.\n",
"TA.034A.VM1 has no paz. Skipping.\n",
"TA.034A.VM2 has no paz. Skipping.\n",
"TA.034A.VM3 has no paz. Skipping.\n",
"TA.034A.VM4 has no paz. Skipping.\n",
"TA.034A.VM5 has no paz. Skipping.\n",
"TA.034A.VM6 has no paz. Skipping.\n",
"TA.034A.VPB has no paz. Skipping.\n",
"TA.034A.LCE has no paz. Skipping.\n",
"TA.034A.LCO has no paz. Skipping.\n",
"TA.034A.LDM has no paz. Skipping.\n",
"TA.034A.LEP has no paz. Skipping.\n",
"TA.034A.LIM has no paz. Skipping.\n",
"TA.034A.LKM has no paz. Skipping.\n",
"TA.034A.QEP has unrecognized input units in response: COUNTS. Skipping\n"
]
}
],
"source": [
"inv.write(f, format='SACPZ')"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/plain": [
"0"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"f.seek(0)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"scrolled": true
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"* **************************************************\n",
"* NETWORK : TA\n",
"* STATION : 034A\n",
"* LOCATION : \n",
"* CHANNEL : BHE\n",
"* CREATED : 2018-07-24T16:59:54.815787Z\n",
"* START : 2010-01-08T00:00:00.000000Z\n",
"* END : 2011-11-17T17:05:00.000000Z\n",
"* DESCRIPTION : Hebronville, TX, USA\n",
"* LATITUDE : 27.064699\n",
"* LONGITUDE : -98.683296\n",
"* ELEVATION : 155.0\n",
"* DEPTH : 0.0\n",
"* DIP (SEED) : 0.0\n",
"* AZIMUTH : 89.1\n",
"* SAMPLE RATE : 40.0\n",
"* INPUT UNIT : M\n",
"* OUTPUT UNIT : COUNTS\n",
"* INSTTYPE : None\n",
"* INSTGAIN : 1504.2 (M/S)\n",
"* SENSITIVITY : 630907000.0 (M/S)\n",
"* A0 : 571400000.0\n",
"* **************************************************\n",
"ZEROS 3\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
"POLES 5\n",
" -3.701000e-02 +3.701000e-02\n",
" -3.701000e-02 -3.701000e-02\n",
" -1.131000e+03 +0.000000e+00\n",
" -1.005000e+03 +0.000000e+00\n",
" -5.027000e+02 +0.000000e+00\n",
"CONSTANT 3.605003e+17\n",
"\n",
"\n",
"* **************************************************\n",
"* NETWORK : TA\n",
"* STATION : 034A\n",
"* LOCATION : \n",
"* CHANNEL : BHN\n",
"* CREATED : 2018-07-24T16:59:54.815787Z\n",
"* START : 2010-01-08T00:00:00.000000Z\n",
"* END : 2011-11-17T17:05:00.000000Z\n",
"* DESCRIPTION : Hebronville, TX, USA\n",
"* LATITUDE : 27.064699\n",
"* LONGITUDE : -98.683296\n",
"* ELEVATION : 155.0\n",
"* DEPTH : 0.0\n",
"* DIP (SEED) : 0.0\n",
"* AZIMUTH : 359.1\n",
"* SAMPLE RATE : 40.0\n",
"* INPUT UNIT : M\n",
"* OUTPUT UNIT : COUNTS\n",
"* INSTTYPE : None\n",
"* INSTGAIN : 1504.2 (M/S)\n",
"* SENSITIVITY : 630907000.0 (M/S)\n",
"* A0 : 571400000.0\n",
"* **************************************************\n",
"ZEROS 3\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
"POLES 5\n",
" -3.701000e-02 +3.701000e-02\n",
" -3.701000e-02 -3.701000e-02\n",
" -1.131000e+03 +0.000000e+00\n",
" -1.005000e+03 +0.000000e+00\n",
" -5.027000e+02 +0.000000e+00\n",
"CONSTANT 3.605003e+17\n",
"\n",
"\n",
"* **************************************************\n",
"* NETWORK : TA\n",
"* STATION : 034A\n",
"* LOCATION : \n",
"* CHANNEL : BHZ\n",
"* CREATED : 2018-07-24T16:59:54.815787Z\n",
"* START : 2010-01-08T00:00:00.000000Z\n",
"* END : 2011-11-17T17:05:00.000000Z\n",
"* DESCRIPTION : Hebronville, TX, USA\n",
"* LATITUDE : 27.064699\n",
"* LONGITUDE : -98.683296\n",
"* ELEVATION : 155.0\n",
"* DEPTH : 0.0\n",
"* DIP (SEED) : -90.0\n",
"* AZIMUTH : 0.0\n",
"* SAMPLE RATE : 40.0\n",
"* INPUT UNIT : M\n",
"* OUTPUT UNIT : COUNTS\n",
"* INSTTYPE : None\n",
"* INSTGAIN : 1504.2 (M/S)\n",
"* SENSITIVITY : 630907000.0 (M/S)\n",
"* A0 : 571400000.0\n",
"* **************************************************\n",
"ZEROS 3\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
"POLES 5\n",
" -3.701000e-02 +3.701000e-02\n",
" -3.701000e-02 -3.701000e-02\n",
" -1.131000e+03 +0.000000e+00\n",
" -1.005000e+03 +0.000000e+00\n",
" -5.027000e+02 +0.000000e+00\n",
"CONSTANT 3.605003e+17\n",
"\n",
"\n",
"* **************************************************\n",
"* NETWORK : TA\n",
"* STATION : 034A\n",
"* LOCATION : \n",
"* CHANNEL : LHE\n",
"* CREATED : 2018-07-24T16:59:54.815787Z\n",
"* START : 2010-01-08T00:00:00.000000Z\n",
"* END : 2011-11-17T17:05:00.000000Z\n",
"* DESCRIPTION : Hebronville, TX, USA\n",
"* LATITUDE : 27.064699\n",
"* LONGITUDE : -98.683296\n",
"* ELEVATION : 155.0\n",
"* DEPTH : 0.0\n",
"* DIP (SEED) : 0.0\n",
"* AZIMUTH : 89.1\n",
"* SAMPLE RATE : 1.0\n",
"* INPUT UNIT : M\n",
"* OUTPUT UNIT : COUNTS\n",
"* INSTTYPE : None\n",
"* INSTGAIN : 1499.76 (M/S)\n",
"* SENSITIVITY : 629044000.0 (M/S)\n",
"* A0 : 573092000.0\n",
"* **************************************************\n",
"ZEROS 3\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
"POLES 5\n",
" -3.701000e-02 +3.701000e-02\n",
" -3.701000e-02 -3.701000e-02\n",
" -1.131000e+03 +0.000000e+00\n",
" -1.005000e+03 +0.000000e+00\n",
" -5.027000e+02 +0.000000e+00\n",
"CONSTANT 3.605001e+17\n",
"\n",
"\n",
"* **************************************************\n",
"* NETWORK : TA\n",
"* STATION : 034A\n",
"* LOCATION : \n",
"* CHANNEL : LHN\n",
"* CREATED : 2018-07-24T16:59:54.815787Z\n",
"* START : 2010-01-08T00:00:00.000000Z\n",
"* END : 2011-11-17T17:05:00.000000Z\n",
"* DESCRIPTION : Hebronville, TX, USA\n",
"* LATITUDE : 27.064699\n",
"* LONGITUDE : -98.683296\n",
"* ELEVATION : 155.0\n",
"* DEPTH : 0.0\n",
"* DIP (SEED) : 0.0\n",
"* AZIMUTH : 359.1\n",
"* SAMPLE RATE : 1.0\n",
"* INPUT UNIT : M\n",
"* OUTPUT UNIT : COUNTS\n",
"* INSTTYPE : None\n",
"* INSTGAIN : 1499.76 (M/S)\n",
"* SENSITIVITY : 629044000.0 (M/S)\n",
"* A0 : 573092000.0\n",
"* **************************************************\n",
"ZEROS 3\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
"POLES 5\n",
" -3.701000e-02 +3.701000e-02\n",
" -3.701000e-02 -3.701000e-02\n",
" -1.131000e+03 +0.000000e+00\n",
" -1.005000e+03 +0.000000e+00\n",
" -5.027000e+02 +0.000000e+00\n",
"CONSTANT 3.605001e+17\n",
"\n",
"\n",
"* **************************************************\n",
"* NETWORK : TA\n",
"* STATION : 034A\n",
"* LOCATION : \n",
"* CHANNEL : LHZ\n",
"* CREATED : 2018-07-24T16:59:54.815787Z\n",
"* START : 2010-01-08T00:00:00.000000Z\n",
"* END : 2011-11-17T17:05:00.000000Z\n",
"* DESCRIPTION : Hebronville, TX, USA\n",
"* LATITUDE : 27.064699\n",
"* LONGITUDE : -98.683296\n",
"* ELEVATION : 155.0\n",
"* DEPTH : 0.0\n",
"* DIP (SEED) : -90.0\n",
"* AZIMUTH : 0.0\n",
"* SAMPLE RATE : 1.0\n",
"* INPUT UNIT : M\n",
"* OUTPUT UNIT : COUNTS\n",
"* INSTTYPE : None\n",
"* INSTGAIN : 1499.76 (M/S)\n",
"* SENSITIVITY : 629044000.0 (M/S)\n",
"* A0 : 573092000.0\n",
"* **************************************************\n",
"ZEROS 3\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
"POLES 5\n",
" -3.701000e-02 +3.701000e-02\n",
" -3.701000e-02 -3.701000e-02\n",
" -1.131000e+03 +0.000000e+00\n",
" -1.005000e+03 +0.000000e+00\n",
" -5.027000e+02 +0.000000e+00\n",
"CONSTANT 3.605001e+17\n",
"\n",
"\n",
"* **************************************************\n",
"* NETWORK : TA\n",
"* STATION : 034A\n",
"* LOCATION : \n",
"* CHANNEL : UHE\n",
"* CREATED : 2018-07-24T16:59:54.815787Z\n",
"* START : 2010-01-08T00:00:00.000000Z\n",
"* END : 2011-11-17T17:05:00.000000Z\n",
"* DESCRIPTION : Hebronville, TX, USA\n",
"* LATITUDE : 27.064699\n",
"* LONGITUDE : -98.683296\n",
"* ELEVATION : 155.0\n",
"* DEPTH : 0.0\n",
"* DIP (SEED) : 0.0\n",
"* AZIMUTH : 89.1\n",
"* SAMPLE RATE : 0.01\n",
"* INPUT UNIT : M\n",
"* OUTPUT UNIT : COUNTS\n",
"* INSTTYPE : None\n",
"* INSTGAIN : 86.5644 (M/S)\n",
"* SENSITIVITY : 36307700.0 (M/S)\n",
"* A0 : 9929020000.0\n",
"* **************************************************\n",
"ZEROS 3\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
"POLES 5\n",
" -3.701000e-02 +3.701000e-02\n",
" -3.701000e-02 -3.701000e-02\n",
" -1.131000e+03 +0.000000e+00\n",
" -1.005000e+03 +0.000000e+00\n",
" -5.027000e+02 +0.000000e+00\n",
"CONSTANT 3.604999e+17\n",
"\n",
"\n",
"* **************************************************\n",
"* NETWORK : TA\n",
"* STATION : 034A\n",
"* LOCATION : \n",
"* CHANNEL : UHN\n",
"* CREATED : 2018-07-24T16:59:54.815787Z\n",
"* START : 2010-01-08T00:00:00.000000Z\n",
"* END : 2011-11-17T17:05:00.000000Z\n",
"* DESCRIPTION : Hebronville, TX, USA\n",
"* LATITUDE : 27.064699\n",
"* LONGITUDE : -98.683296\n",
"* ELEVATION : 155.0\n",
"* DEPTH : 0.0\n",
"* DIP (SEED) : 0.0\n",
"* AZIMUTH : 359.1\n",
"* SAMPLE RATE : 0.01\n",
"* INPUT UNIT : M\n",
"* OUTPUT UNIT : COUNTS\n",
"* INSTTYPE : None\n",
"* INSTGAIN : 86.5644 (M/S)\n",
"* SENSITIVITY : 36307700.0 (M/S)\n",
"* A0 : 9929020000.0\n",
"* **************************************************\n",
"ZEROS 3\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
"POLES 5\n",
" -3.701000e-02 +3.701000e-02\n",
" -3.701000e-02 -3.701000e-02\n",
" -1.131000e+03 +0.000000e+00\n",
" -1.005000e+03 +0.000000e+00\n",
" -5.027000e+02 +0.000000e+00\n",
"CONSTANT 3.604999e+17\n",
"\n",
"\n",
"* **************************************************\n",
"* NETWORK : TA\n",
"* STATION : 034A\n",
"* LOCATION : \n",
"* CHANNEL : UHZ\n",
"* CREATED : 2018-07-24T16:59:54.815787Z\n",
"* START : 2010-01-08T00:00:00.000000Z\n",
"* END : 2011-11-17T17:05:00.000000Z\n",
"* DESCRIPTION : Hebronville, TX, USA\n",
"* LATITUDE : 27.064699\n",
"* LONGITUDE : -98.683296\n",
"* ELEVATION : 155.0\n",
"* DEPTH : 0.0\n",
"* DIP (SEED) : -90.0\n",
"* AZIMUTH : 0.0\n",
"* SAMPLE RATE : 0.01\n",
"* INPUT UNIT : M\n",
"* OUTPUT UNIT : COUNTS\n",
"* INSTTYPE : None\n",
"* INSTGAIN : 86.5644 (M/S)\n",
"* SENSITIVITY : 36307700.0 (M/S)\n",
"* A0 : 9929020000.0\n",
"* **************************************************\n",
"ZEROS 3\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
"POLES 5\n",
" -3.701000e-02 +3.701000e-02\n",
" -3.701000e-02 -3.701000e-02\n",
" -1.131000e+03 +0.000000e+00\n",
" -1.005000e+03 +0.000000e+00\n",
" -5.027000e+02 +0.000000e+00\n",
"CONSTANT 3.604999e+17\n",
"\n",
"\n",
"* **************************************************\n",
"* NETWORK : TA\n",
"* STATION : 034A\n",
"* LOCATION : \n",
"* CHANNEL : VHE\n",
"* CREATED : 2018-07-24T16:59:54.815787Z\n",
"* START : 2010-01-08T00:00:00.000000Z\n",
"* END : 2011-11-17T17:05:00.000000Z\n",
"* DESCRIPTION : Hebronville, TX, USA\n",
"* LATITUDE : 27.064699\n",
"* LONGITUDE : -98.683296\n",
"* ELEVATION : 155.0\n",
"* DEPTH : 0.0\n",
"* DIP (SEED) : 0.0\n",
"* AZIMUTH : 89.1\n",
"* SAMPLE RATE : 0.1\n",
"* INPUT UNIT : M\n",
"* OUTPUT UNIT : COUNTS\n",
"* INSTTYPE : None\n",
"* INSTGAIN : 1235.82 (M/S)\n",
"* SENSITIVITY : 518339000.0 (M/S)\n",
"* A0 : 695491000.0\n",
"* **************************************************\n",
"ZEROS 3\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
"POLES 5\n",
" -3.701000e-02 +3.701000e-02\n",
" -3.701000e-02 -3.701000e-02\n",
" -1.131000e+03 +0.000000e+00\n",
" -1.005000e+03 +0.000000e+00\n",
" -5.027000e+02 +0.000000e+00\n",
"CONSTANT 3.605001e+17\n",
"\n",
"\n",
"* **************************************************\n",
"* NETWORK : TA\n",
"* STATION : 034A\n",
"* LOCATION : \n",
"* CHANNEL : VHN\n",
"* CREATED : 2018-07-24T16:59:54.815787Z\n",
"* START : 2010-01-08T00:00:00.000000Z\n",
"* END : 2011-11-17T17:05:00.000000Z\n",
"* DESCRIPTION : Hebronville, TX, USA\n",
"* LATITUDE : 27.064699\n",
"* LONGITUDE : -98.683296\n",
"* ELEVATION : 155.0\n",
"* DEPTH : 0.0\n",
"* DIP (SEED) : 0.0\n",
"* AZIMUTH : 359.1\n",
"* SAMPLE RATE : 0.1\n",
"* INPUT UNIT : M\n",
"* OUTPUT UNIT : COUNTS\n",
"* INSTTYPE : None\n",
"* INSTGAIN : 1235.82 (M/S)\n",
"* SENSITIVITY : 518339000.0 (M/S)\n",
"* A0 : 695491000.0\n",
"* **************************************************\n",
"ZEROS 3\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
"POLES 5\n",
" -3.701000e-02 +3.701000e-02\n",
" -3.701000e-02 -3.701000e-02\n",
" -1.131000e+03 +0.000000e+00\n",
" -1.005000e+03 +0.000000e+00\n",
" -5.027000e+02 +0.000000e+00\n",
"CONSTANT 3.605001e+17\n",
"\n",
"\n",
"* **************************************************\n",
"* NETWORK : TA\n",
"* STATION : 034A\n",
"* LOCATION : \n",
"* CHANNEL : VHZ\n",
"* CREATED : 2018-07-24T16:59:54.815787Z\n",
"* START : 2010-01-08T00:00:00.000000Z\n",
"* END : 2011-11-17T17:05:00.000000Z\n",
"* DESCRIPTION : Hebronville, TX, USA\n",
"* LATITUDE : 27.064699\n",
"* LONGITUDE : -98.683296\n",
"* ELEVATION : 155.0\n",
"* DEPTH : 0.0\n",
"* DIP (SEED) : -90.0\n",
"* AZIMUTH : 0.0\n",
"* SAMPLE RATE : 0.1\n",
"* INPUT UNIT : M\n",
"* OUTPUT UNIT : COUNTS\n",
"* INSTTYPE : None\n",
"* INSTGAIN : 1235.82 (M/S)\n",
"* SENSITIVITY : 518339000.0 (M/S)\n",
"* A0 : 695491000.0\n",
"* **************************************************\n",
"ZEROS 3\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
" +0.000000e+00 +0.000000e+00\n",
"POLES 5\n",
" -3.701000e-02 +3.701000e-02\n",
" -3.701000e-02 -3.701000e-02\n",
" -1.131000e+03 +0.000000e+00\n",
" -1.005000e+03 +0.000000e+00\n",
" -5.027000e+02 +0.000000e+00\n",
"CONSTANT 3.605001e+17\n",
"\n",
"\n",
"\n",
"\n"
]
}
],
"source": [
"print(f.read())"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.6.5"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment