Skip to content

Commit 8f3f89d

Browse files
ajakefajakef
committed
changed to lowercase panel labels
1 parent 3db8047 commit 8f3f89d

7 files changed

+40
-49
lines changed

‎code/Fig3_array_response.py

Lines changed: 6 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -37,7 +37,7 @@ def circle(r):
3737
plt.axis('equal')
3838
plt.xlabel('Easting (m)')
3939
plt.ylabel('Northing (m)')
40-
plt.title('A. Station Map for XP.PARK', loc='left')
40+
plt.title('a. Station Map for XP.PARK', loc='left')
4141
plt.axis('equal')
4242
plt.legend(['Sub-array (N=3)', 'Full Array (N=22)'], loc = 'upper right', fontsize = 'small')
4343

@@ -54,7 +54,7 @@ def circle(r):
5454

5555
plt.plot(*circle(2*np.pi * 2 / 0.333/100), 'k--')
5656
plt.plot(*circle(2*np.pi * 10 / 0.333/100), 'k--')
57-
plt.title('B. Array Response (N=3)', loc = 'left')
57+
plt.title('b. Array Response (N=3)', loc = 'left')
5858
plt.xlabel('$k_x$ (×$10^2$ rad/km)', labelpad = -1)
5959
plt.ylabel('$k_y$ (×$10^2$ rad/km)', labelpad = -1)
6060

@@ -72,7 +72,7 @@ def circle(r):
7272
cbar.set_label('Semblance')
7373
plt.plot(*circle(2*np.pi * 2 / 0.333/100), 'k--')
7474
plt.plot(*circle(2*np.pi * 10 / 0.333/100), 'k--')
75-
plt.title('C. Array Response (N=22)', loc = 'left')
75+
plt.title('c. Array Response (N=22)', loc = 'left')
7676
plt.xlabel('$k_x$ (×$10^2$ rad/km)', labelpad = -1)
7777
plt.ylabel('$k_y$ (×$10^2$ rad/km)', labelpad = -1)
7878

@@ -89,7 +89,7 @@ def circle(r):
8989
plt.plot(*circle(3), 'k--')
9090
plt.xlabel('$s_x$ (s/km)', labelpad = -1)
9191
plt.ylabel('$s_y$ (s/km)', labelpad = -1)
92-
plt.title('D. Slowness Spectrum (N=3)', loc = 'left')
92+
plt.title('d. Slowness Spectrum (N=3)', loc = 'left')
9393
#%% E: PARK slowness spectrum
9494
plt.subplot(gs[5:7,1])
9595
with open('data/pkl/clean_aftershock_full.pkl', 'rb') as file:
@@ -104,7 +104,7 @@ def circle(r):
104104

105105
plt.xlabel('$s_x$ (s/km)', labelpad = -1)
106106
plt.ylabel("""$s_y$ (s/km)""", labelpad = -1)
107-
plt.title('E. Slowness Spectrum (N=22)', loc = 'left')
107+
plt.title('e. Slowness Spectrum (N=22)', loc = 'left')
108108

109109

110110

@@ -113,6 +113,6 @@ def circle(r):
113113
plt.gcf().set_size_inches(6.5, 9, forward=True) # max fig size for villarrica screen: 12.94x6.65 inch
114114
plt.tight_layout()
115115
#%%
116-
plt.savefig('figures/Fig3_array_response.png')
116+
plt.savefig('figures/Fig3_array_response.png', dpi=300)
117117

118118
####################

‎code/Fig4_clean_aftershock.py

Lines changed: 9 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -70,7 +70,7 @@
7070
# set_xlim leaves some padding; axis('image') changes the aspect ratio. This sets tight axis limits without forcing a bad aspect ratio.
7171
plt.gca().set_xbound(t[0],t[-1])
7272
plt.ylabel('Pa')
73-
plt.title('A. Infrasound Waveform', loc = 'left', fontsize = 11)
73+
plt.title('a. Infrasound Waveform', loc = 'left', fontsize = 11)
7474
plt.xticks(xticks, labels = [])
7575
plt.yticks([-2, 0, 2])
7676

@@ -119,7 +119,7 @@ def imageAdj(x): return(np.log(trim(x, ac_max, ac_max*r)))
119119
#plt.axvline(0, color = 'black', lw = 0.5)
120120
plt.axhline(slowness_threshold, color = 'black', lw = 0.5, ls = '--')
121121
plt.ylabel('s/km')
122-
plt.title('B. Power vs. Time, Slowness', loc = 'left')
122+
plt.title('b. Power vs. Time, Slowness', loc = 'left')
123123
w = output['original_sh'] < 4
124124
plt.plot(output['t'][w]-(event - loop_start), output['original_sh'][w], 'k.', markersize=5)
125125
plt.plot(output['t'][w]-(event - loop_start), output['original_sh'][w], 'w.', markersize=2)
@@ -132,7 +132,7 @@ def imageAdj(x): return(np.log(trim(x, ac_max, ac_max*r)))
132132
cleanbf.image(imageAdj(spec_baz), output['t']-(event - loop_start), baz, crosshairs = False)
133133
plt.ylabel('degrees')
134134
#plt.xlabel('Time after earthquake (seconds)')
135-
plt.title(f'C. Power vs Time, Backazimuth (Slowness > {slowness_threshold} s/km)', loc = 'left')
135+
plt.title(f'c. Power vs Time, Backazimuth (Slowness > {slowness_threshold} s/km)', loc = 'left')
136136
plt.plot(output['t']-(event - loop_start), (output['original_az']) % 360 - 180, 'k.', markersize=5)
137137
plt.plot(output['t']-(event - loop_start), (output['original_az']) % 360 - 180, 'w.', markersize=2)
138138
plt.yticks(baz_ticks)
@@ -198,7 +198,7 @@ def imageAdj(x): return(np.log(trim(x, ac_max, ac_max*r)))
198198
cleanbf.image(np.log(trim(spec_s, spec_s[:,w].max(), r*spec_s[:,w].max())), output['t']-(event - loop_start), sh, crosshairs = False)
199199
plt.axhline(slowness_threshold, color = 'black', lw = 0.5, ls = '--')
200200
plt.ylabel('s/km')
201-
plt.title('D. Power vs. Time, Slowness (3 sensors)', loc = 'left')
201+
plt.title('d. Power vs. Time, Slowness (3 sensors)', loc = 'left')
202202
w = output['original_sh'] < 4
203203
plt.plot(output['t'][w]-(event - loop_start), output['original_sh'][w], 'k.', markersize=5)
204204
plt.plot(output['t'][w]-(event - loop_start), output['original_sh'][w], 'w.', markersize=2)
@@ -210,7 +210,7 @@ def imageAdj(x): return(np.log(trim(x, ac_max, ac_max*r)))
210210
im = cleanbf.image(imageAdj(spec_baz), output['t']-(event - loop_start), baz, crosshairs = False) # im used later for colorbar
211211
plt.ylabel('degrees')
212212
#plt.xlabel('Time after earthquake (seconds)')
213-
plt.title(f'E. Power vs. Time, Backazimuth (Slowness > {slowness_threshold} s/km; 3 sensors)', loc = 'left')
213+
plt.title(f'e. Power vs. Time, Backazimuth (Slowness > {slowness_threshold} s/km; 3 sensors)', loc = 'left')
214214
plt.plot(output['t']-(event - loop_start), (output['original_az']) % 360 - 180, 'k.', markersize=5)
215215
plt.plot(output['t']-(event - loop_start), (output['original_az']) % 360 - 180, 'w.', markersize=2)
216216
plt.yticks(baz_ticks)
@@ -224,7 +224,7 @@ def imageAdj(x): return(np.log(trim(x, ac_max, ac_max*r)))
224224
plt.subplot(7,1,6)
225225
plt.axhline(slowness_threshold, color = 'black', lw = 0.5, ls = '--')
226226
plt.ylabel('s/km')
227-
plt.title('F. Slowness, traditional beamforming (3 sensors vs. 20 sensors)', loc = 'left')
227+
plt.title('f. Slowness, traditional beamforming (3 sensors vs. 20 sensors)', loc = 'left')
228228
w = output['original_sh'] < 4
229229
plt.plot(output_full['t'][w]-(event - loop_start), output_full['original_sh'][w], 'b.', markersize=5)
230230
plt.plot(output_3['t'][w]-(event - loop_start), output_3['original_sh'][w], 'r.', markersize=5)
@@ -235,7 +235,7 @@ def imageAdj(x): return(np.log(trim(x, ac_max, ac_max*r)))
235235
plt.subplot(7,1,7)
236236
plt.ylabel('degrees')
237237
plt.xlabel('Time after earthquake (seconds)')
238-
plt.title('G. Backazimuth, traditional beamforming (3 sensors vs. 20 sensors)', loc = 'left')
238+
plt.title('g. Backazimuth, traditional beamforming (3 sensors vs. 20 sensors)', loc = 'left')
239239
plt.plot(output_full['t']-(event - loop_start), (output_full['original_az']) % 360 - 180, 'b.', markersize=5, label = 'N=20')
240240
plt.plot(output_3['t']-(event - loop_start), (output_3['original_az']) % 360 - 180, 'r.', markersize=5, label = 'N=3')
241241
plt.legend(loc='upper right', framealpha=1)
@@ -256,7 +256,7 @@ def imageAdj(x): return(np.log(trim(x, ac_max, ac_max*r)))
256256

257257

258258
#%% save figure
259-
plt.savefig('figures/Fig4_Aftershock_20_3_paper.png')
259+
plt.savefig('figures/Fig4_Aftershock_20_3_paper.png', dpi = 300)
260260
#%% calculate 6-station clean result for use in supplementary figure
261261

262262
calculate_new_beamform_result_6 = False # set to False after running the calculation once to save time
@@ -283,4 +283,4 @@ def imageAdj(x): return(np.log(trim(x, ac_max, ac_max*r)))
283283
with open(pkl_file6, 'rb') as file:
284284
d = pickle.load(file)
285285
locals().update(d)
286-
output_6 = output
286+
output_6 = output

‎code/Fig5_clean_MonteCristo.py

Lines changed: 8 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -65,7 +65,7 @@ def image_show_outliers(x, n=2):
6565
# set_xlim leaves some padding; axis('image') changes the aspect ratio. This sets tight axis limits without forcing a bad aspect ratio.
6666
plt.gca().set_xbound(-start_offset,3600-start_offset)
6767
plt.ylabel('Pa')
68-
plt.title('A. Infrasound Waveform', loc = 'left', fontsize = 'small')
68+
plt.title('a. Infrasound Waveform', loc = 'left', fontsize = 'small')
6969
plt.xticks(np.arange(0, 3001, 500), labels = [])
7070
plt.yticks(np.arange(-1,2)*0.025)
7171

@@ -105,7 +105,7 @@ def image_show_outliers(x, n=2):
105105
plt.ylabel('s/km')
106106
plt.xticks(np.arange(0, 3001, 500), labels = [])
107107
plt.yticks([0,1,2,3])
108-
plt.title('B. Power vs. Time, Slowness', loc = 'left', fontsize = 'small')
108+
plt.title('b. Power vs. Time, Slowness', loc = 'left', fontsize = 'small')
109109
w = output['original_sh'] < 4
110110

111111
plt.subplot(7,1,3)
@@ -114,15 +114,15 @@ def image_show_outliers(x, n=2):
114114
plt.xticks(np.arange(0, 3001, 500), labels = [])
115115
plt.yticks(baz_ticks)
116116
for i in baz_ticks: plt.axhline(i, color = 'gray', lw = 0.25)
117-
plt.title('C. Power vs. Time, Backazimuth (Slowness > 2 s/km)', loc = 'left', fontsize = 'small')
117+
plt.title('c. Power vs. Time, Backazimuth (Slowness > 2 s/km)', loc = 'left', fontsize = 'small')
118118

119119
plt.subplot(7,1,4)
120120
cleanbf.image(image_show_outliers(spec_baz), output['t']-(event - trace_start), baz, crosshairs = False)
121121
plt.ylabel('degrees')#plt.xlabel('Time after earthquake (seconds)')
122122
plt.xticks(np.arange(0, 3001, 500), labels = [])
123123
plt.yticks(baz_ticks)
124124
for i in baz_ticks: plt.axhline(i, color = 'gray', lw = 0.25)
125-
plt.title('D. Power vs. Time, Backazimuth (Slowness > 2 s/km; above-ambient)', loc = 'left', fontsize = 'small')
125+
plt.title('d. Power vs. Time, Backazimuth (Slowness > 2 s/km; above-ambient)', loc = 'left', fontsize = 'small')
126126

127127
plt.tight_layout()
128128

@@ -171,7 +171,7 @@ def image_show_outliers(x):
171171
plt.ylabel('s/km')
172172
plt.xticks(np.arange(0, 3001, 500), labels = [])
173173
plt.yticks([0,1,2,3])
174-
plt.title('E. Power vs. Time, Slowness (3 sensors)', loc = 'left', fontsize = 'small')
174+
plt.title('e. Power vs. Time, Slowness (3 sensors)', loc = 'left', fontsize = 'small')
175175
w = output['original_sh'] < 4
176176

177177
plt.subplot(7,1,6)
@@ -180,7 +180,7 @@ def image_show_outliers(x):
180180
plt.xticks(np.arange(0, 3001, 500), labels = [])
181181
plt.yticks(baz_ticks)
182182
for i in baz_ticks: plt.axhline(i, color = 'gray', lw = 0.25)
183-
plt.title('F. Power vs. Time, Backazimuth (Slowness > 2 s/km)', loc = 'left', fontsize = 'small')
183+
plt.title('f. Power vs. Time, Backazimuth (Slowness > 2 s/km)', loc = 'left', fontsize = 'small')
184184

185185
plt.subplot(7,1,7)
186186
im=cleanbf.image(image_show_outliers(spec_baz), output['t']-(event - loop_start), baz, crosshairs = False)
@@ -190,7 +190,7 @@ def image_show_outliers(x):
190190
plt.xticks(np.arange(0, 3001, 500))
191191
plt.yticks(baz_ticks)
192192
for i in baz_ticks: plt.axhline(i, color = 'gray', lw = 0.25)
193-
plt.title('G. Power vs. Time, Backazimuth (Slowness > 2 s/km; above-ambient)', loc = 'left', fontsize = 'small')
193+
plt.title('g. Power vs. Time, Backazimuth (Slowness > 2 s/km; above-ambient)', loc = 'left', fontsize = 'small')
194194

195195

196196
#%%
@@ -204,5 +204,5 @@ def image_show_outliers(x):
204204

205205

206206
#%%
207-
plt.savefig('figures/Fig5_montecristo_17_3_paper.png')
207+
plt.savefig('figures/Fig5_montecristo_17_3_paper.png', dpi = 300)
208208

‎code/Fig6_backprojection_map.py

Lines changed: 3 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -81,11 +81,10 @@
8181
ax.set_aspect(1/np.cos(lat_station * np.pi/180))
8282
fig.colorbar(a, ax = ax, location = 'top', ticks = 10**np.arange(3, 6), shrink = 0.35, label = 'Infrasound Production')
8383
fig.colorbar(b, ax = ax, location = 'right', ticks = np.arange(1500, 3001, 500), label = 'Elevation (m)')
84-
fig.savefig('figures/Fig6a_aftershock_backprojection_map.png')
84+
fig.savefig('figures/Fig6a_aftershock_backprojection_map.png', dpi = 300)
8585

8686

8787
#%% download seismic data for the Monte Cristo EQ backprojection plot
88-
import obspy, shutil
8988
from obspy.clients.fdsn.mass_downloader import GlobalDomain, Restrictions, MassDownloader
9089

9190
## change project_dir to match where you store these files
@@ -191,7 +190,7 @@
191190
plt.ylabel('Epicentral Distance (degrees)')
192191
plt.tight_layout()
193192

194-
plt.savefig('figures/Fig6c_MonteCristo_seismograms.png')
193+
plt.savefig('figures/Fig6c_MonteCristo_seismograms.png', dpi = 300)
195194
#%% set up Monte Cristo backprojection
196195
slowness_threshold = 2
197196
with open('data/pkl/clean_MonteCristo_full.pkl', 'rb') as file:
@@ -318,4 +317,4 @@
318317
plt.text(station_info.Longitude[i], station_info.Latitude[i]+station_info.rms[i]*0.6e4, station_info.Station[i], size = 'small')
319318

320319

321-
fig.savefig('figures/Fig6b_MonteCristo_backprojection_map.png')
320+
fig.savefig('figures/Fig6b_MonteCristo_backprojection_map.png', dpi = 300)

‎code/FigS1_clean_MonteCristo_N6.py

Lines changed: 6 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -114,21 +114,13 @@ def circle(r):
114114
ac_max = np.quantile(spec_baz, 0.999)
115115

116116

117-
def trim(x, a, b): x[x>a]=a; x[x<b]=b; return x
118-
def image_trim_log(x): return(np.log(trim(x, ac_max, ac_max*r)))
119-
def image_show_outliers(x):
120-
mean = np.einsum('ij->j', x)/x.shape[0]
121-
std = np.sqrt(np.einsum('ij->j', (x - mean)**2)/x.shape[0])
122-
x[x < (mean +2*std)] = 0
123-
return np.log(trim(x, ac_max, r*ac_max))
124-
125117
plt.subplot(5,1,1)
126118
cleanbf.image(image_trim_log(spec_s), output['t']-(event - loop_start), sh, crosshairs = False)
127119
plt.axhline(slowness_threshold, color = 'black', lw = 0.5, ls = '--')
128120
plt.ylabel('s/km')
129121
plt.xticks(np.arange(0, 3001, 500), labels = [])
130122
plt.yticks([0,1,2,3])
131-
plt.title('A. N=6, Slowness (6 sensors)', loc = 'left', fontsize = 11)
123+
plt.title('a. N=6, Slowness (6 sensors)', loc = 'left', fontsize = 11)
132124
w = output['original_sh'] < 4
133125

134126
plt.subplot(5,1,2)
@@ -137,7 +129,7 @@ def image_show_outliers(x):
137129
plt.xticks(np.arange(0, 3001, 500), labels = [])
138130
plt.yticks(baz_ticks)
139131
for i in baz_ticks: plt.axhline(i, color = 'gray', lw = 0.25)
140-
plt.title('B. N=6, Backazimuth (Slowness > 2 s/km)', loc = 'left')
132+
plt.title('b. N=6, Backazimuth (Slowness > 2 s/km)', loc = 'left')
141133

142134
plt.subplot(5,1,3)
143135
im = cleanbf.image(image_show_outliers(spec_baz), output['t']-(event - loop_start), baz, crosshairs = False)
@@ -146,7 +138,7 @@ def image_show_outliers(x):
146138
plt.xticks(np.arange(0, 3001, 500), labels = [])
147139
plt.yticks(baz_ticks)
148140
for i in baz_ticks: plt.axhline(i, color = 'gray', lw = 0.25)
149-
plt.title('C. N=6, Backazimuth (Slowness > 2 s/km; above-ambient)', loc = 'left', fontsize = 11)
141+
plt.title('c. N=6, Backazimuth (Slowness > 2 s/km; above-ambient)', loc = 'left', fontsize = 11)
150142
#%% N = 3
151143
with open(pkl_file3, 'rb') as file:
152144
d = pickle.load(file)
@@ -169,7 +161,7 @@ def image_show_outliers(x):
169161
plt.xticks(np.arange(0, 3001, 500), labels = [])
170162
plt.yticks(baz_ticks)
171163
for i in baz_ticks: plt.axhline(i, color = 'gray', lw = 0.25)
172-
plt.title('D. N=3, Backazimuth (Slowness > 2 s/km; above-ambient)', loc = 'left', fontsize = 11)
164+
plt.title('d. N=3, Backazimuth (Slowness > 2 s/km; above-ambient)', loc = 'left', fontsize = 11)
173165

174166
#%% full array
175167
with open(pkl_file, 'rb') as file:
@@ -192,7 +184,7 @@ def image_show_outliers(x):
192184
plt.xticks(np.arange(0, 3001, 500))
193185
plt.yticks(baz_ticks)
194186
for i in baz_ticks: plt.axhline(i, color = 'gray', lw = 0.25)
195-
plt.title('E. N=17, Backazimuth (Slowness > 2 s/km; above-ambient)', loc = 'left', fontsize = 11)
187+
plt.title('e. N=17, Backazimuth (Slowness > 2 s/km; above-ambient)', loc = 'left', fontsize = 11)
196188
plt.xlabel('Time after earthquake (seconds)')
197189
plt.ylabel('degrees')
198190

@@ -205,4 +197,4 @@ def image_show_outliers(x):
205197
cbar_ax = fig.add_axes([0.93, 0.25, 0.02, 0.5])
206198
fig.colorbar(im,cax=cbar_ax, label='Infrasound Power', ticks = [])
207199
#%%
208-
plt.savefig('figures/FigS1_montecristo_17_6_paper.png')
200+
plt.savefig('figures/FigS1_montecristo_17_6_paper.png', dpi = 300)

‎code/FigS2array_response_N6.py

Lines changed: 8 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -49,7 +49,7 @@ def circle(r):
4949
plt.axis('equal')
5050
plt.xlabel('Easting (m)')
5151
plt.ylabel('Northing (m)')
52-
plt.title('A. Station Map for XP.PARK', loc='left', fontsize='medium')
52+
plt.title('a. Station Map for XP.PARK', loc='left', fontsize='medium')
5353
plt.legend(loc = 'upper right', fontsize = 'small')
5454
#%% B: 3-element array response (wavenumber)
5555
plt.subplot(gs[3:5,0])
@@ -64,7 +64,7 @@ def circle(r):
6464

6565
plt.plot(*circle(2*np.pi * 2 / 0.333/100), 'k--')
6666
plt.plot(*circle(2*np.pi * 10 / 0.333/100), 'k--')
67-
plt.title('B. Array Response (N=3)', loc = 'left', fontsize='medium')
67+
plt.title('b. Array Response (N=3)', loc = 'left', fontsize='medium')
6868
plt.xlabel('$k_x$ (×$10^2$ rad/km)', labelpad = -1)
6969
plt.ylabel('$k_y$ (×$10^2$ rad/km)', labelpad = -1)
7070

@@ -81,7 +81,7 @@ def circle(r):
8181
plt.plot(*circle(3), 'k--')
8282
plt.xlabel('$s_x$ (s/km)', labelpad = -1)
8383
plt.ylabel('$s_y$ (s/km)', labelpad = -1)
84-
plt.title('C. Slowness Spectrum (N=3)', loc = 'left', fontsize = 'medium')
84+
plt.title('c. Slowness Spectrum (N=3)', loc = 'left', fontsize = 'medium')
8585
#%% C: 6-element array response (wavenumber)
8686
plt.subplot(gs[3:5,1])
8787
xy = np.append(np.array(coords_sub_6.loc[:,['x', 'y']]),np.zeros([6,1]), axis=1)
@@ -95,7 +95,7 @@ def circle(r):
9595

9696
plt.plot(*circle(2*np.pi * 2 / 0.333/100), 'k--')
9797
plt.plot(*circle(2*np.pi * 10 / 0.333/100), 'k--')
98-
plt.title('D. N=6', loc = 'left', fontsize='medium')
98+
plt.title('d. N=6', loc = 'left', fontsize='medium')
9999
plt.xlabel('$k_x$ (×$10^2$ rad/km)', labelpad = -1)
100100
plt.ylabel('$k_y$ (×$10^2$ rad/km)', labelpad = -1)
101101

@@ -112,7 +112,7 @@ def circle(r):
112112
plt.plot(*circle(3), 'k--')
113113
plt.xlabel('$s_x$ (s/km)', labelpad = -1)
114114
plt.ylabel('$s_y$ (s/km)', labelpad = -1)
115-
plt.title('E. N=6', loc = 'left', fontsize = 'medium')
115+
plt.title('e. N=6', loc = 'left', fontsize = 'medium')
116116
#%% D: PARK response (wavenumber)
117117
plt.subplot(gs[3:5,2])
118118
xy = np.append(np.array(coords_full.loc[:,['x', 'y']]), np.zeros([22,1]), axis=1)
@@ -126,7 +126,7 @@ def circle(r):
126126
#cbar.set_label('Semblance')
127127
plt.plot(*circle(2*np.pi * 2 / 0.333/100), 'k--')
128128
plt.plot(*circle(2*np.pi * 10 / 0.333/100), 'k--')
129-
plt.title('F. Full Array', loc = 'left', fontsize='medium')
129+
plt.title('f. Full Array', loc = 'left', fontsize='medium')
130130
plt.xlabel('$k_x$ (×$10^2$ rad/km)', labelpad = -1)
131131
plt.ylabel('$k_y$ (×$10^2$ rad/km)', labelpad = -1)
132132
###############################
@@ -143,7 +143,7 @@ def circle(r):
143143
plt.plot(*circle(3), 'k--')
144144
plt.xlabel('$s_x$ (s/km)', labelpad = -1)
145145
plt.ylabel('$s_y$ (s/km)', labelpad = -1)
146-
plt.title('G. Full Array', loc = 'left', fontsize = 'medium')
146+
plt.title('g. Full Array', loc = 'left', fontsize = 'medium')
147147

148148

149149
#%%
@@ -156,4 +156,4 @@ def circle(r):
156156
plt.gcf().set_size_inches(9, 9, forward=True) # max fig size for villarrica screen: 12.94x6.65 inch
157157
plt.tight_layout()
158158
#%%
159-
plt.savefig('figures/FigS2_Array_response_3_6_22.png')
159+
plt.savefig('figures/FigS2_Array_response_3_6_22.png', dpi = 300)

‎figures/Fig3_array_response.png

369 KB
Loading

0 commit comments

Comments
 (0)