4: add option to write out standard CMTSOLUTION files

Author: jbrussell

4-writeCMTSOLUTION.py

@@ -0,0 +1,57 @@
+# %% markdown
+# # writeCMTSOLUTION
+# 
+# Write out standard CMTSOLUTION files for each event
+# 
+# ##### JBR - 10/11/2021
+# %% codecell
+# %load_ext autoreload
+# %autoreload
+from setup_parameters import *
+import matplotlib.pyplot as plt
+import obspy
+from obspy.clients.fdsn import Client
+from obspy import UTCDateTime
+from obspy.core import AttribDict
+from obspy.io.sac import SACTrace
+from obspy import read
+import numpy as np
+import matplotlib.pylab as plt
+# %matplotlib inline
+from datetime import datetime
+import os
+import calendar
+import urllib
+
+# %% codecell
+# Load Client
+client = Client(webservice)
+print(client)
+
+if not os.path.exists(CMTSOLUTION_path):
+    os.makedirs(CMTSOLUTION_path)
+# %% codecell
+t1 = UTCDateTime(tstart)
+t2 = UTCDateTime(tend)
+
+# Load events from GCMT catalogue using IRIS SPUD
+url_query = 'http://ds.iris.edu/spudservice/momenttensor/ids?' \
+           +'evtstartdate='+t1.strftime('%Y-%m-%dT%H:%M:%S') \
+           +'&evtenddate='+t2.strftime('%Y-%m-%dT%H:%M:%S') \
+           +'&evtminmag='+str(minmagnitude)
+evids = urllib.request.urlopen(url_query)
+events_str = '&'.join([line.decode("utf-8").replace("\n", "") for line in evids])+'&'
+url_ndk = 'http://ds.iris.edu/spudservice/momenttensor/bundleids/ndk?'+events_str
+cat_evts = obspy.read_events(url_ndk)
+
+
+# %% codecell
+for iev, ev in enumerate(cat_evts) :
+    
+    evname = ev.origins[0].time.strftime('%Y%m%d%H%M');
+    outpath = CMTSOLUTION_path+'/'+evname+'_CMTSOLUTION'
+    print('Working on '+outpath)
+    
+    ev.write(outpath,format='CMTSOLUTION')
+    
+    

setup_parameters.py

@@ -41,4 +41,7 @@
 
 
 ########## 3-CMT2idagrn ##########
-CMT2idagrn_path = search_dir+'CMT2idagrn/' # OUTPUT PATH
+CMT2idagrn_path = search_dir+'CMT2idagrn/' # OUTPUT PATH
+
+########## 4-writeCMTSOLUTION ##########
+CMTSOLUTION_path = search_dir+'CMTSOLUTIONS/' # OUTPUT PATH