publicImage/gprMax
1
0
派生 0
镜像自地址 https://gitee.com/sunhf/gprMax.git 已同步 2025-08-07 04:56:51 +08:00
代码 工单 软件包 项目 版本发布 百科 动态

Added commented out settings for y axes limits.

浏览代码
这个提交包含在：
Craig Warren
2016-02-25 17:57:38 +00:00
父节点 a18ea7fe47
当前提交 11bf373b11
共有 1 个文件被更改，包括 8 次插入 和 2 次删除
下载 Patch 文件 下载 Diff 文件 展开所有文件 折叠所有文件

10
tools/plot_Ascan.py
查看文件

@@ -129,7 +129,7 @@ for rx in range(1, nrx + 1):
fig, ax = plt.subplots(subplot_kw=dict(xlabel='Time [s]', ylabel=outputtext + ' field strength [V/m]'), num='rx' + str(rx), figsize=(20, 10), facecolor='w', edgecolor='w') fig, ax = plt.subplots(subplot_kw=dict(xlabel='Time [s]', ylabel=outputtext + ' field strength [V/m]'), num='rx' + str(rx), figsize=(20, 10), facecolor='w', edgecolor='w')
line = ax.plot(time, outputdata,'r', lw=2, label=outputtext) line = ax.plot(time, outputdata,'r', lw=2, label=outputtext)
ax.set_xlim([0, np.amax(time)]) ax.set_xlim([0, np.amax(time)])
#ax.set_ylim([-0.25, 0.20]) #ax.set_ylim([-15, 20])
ax.grid() ax.grid()
if 'H' in output: if 'H' in output:
@@ -139,7 +139,7 @@ for rx in range(1, nrx + 1):
plt.setp(line, color='b') plt.setp(line, color='b')
plt.setp(ax, ylabel=outputtext + ', current [A]') plt.setp(ax, ylabel=outputtext + ', current [A]')
# If multiple fields required, create all nine subplots and populate only the specified ones # If multiple outputs required, create all nine subplots and populate only the specified ones
else: else:
fig, ax = plt.subplots(subplot_kw=dict(xlabel='Time [s]'), num='rx' + str(rx), figsize=(20, 10), facecolor='w', edgecolor='w') fig, ax = plt.subplots(subplot_kw=dict(xlabel='Time [s]'), num='rx' + str(rx), figsize=(20, 10), facecolor='w', edgecolor='w')
gs = gridspec.GridSpec(3, 3, hspace=0.3, wspace=0.3) gs = gridspec.GridSpec(3, 3, hspace=0.3, wspace=0.3)
@@ -164,26 +164,32 @@ for rx in range(1, nrx + 1):
ax = plt.subplot(gs[0, 0]) ax = plt.subplot(gs[0, 0])
ax.plot(time, outputdata,'r', lw=2, label=outputtext) ax.plot(time, outputdata,'r', lw=2, label=outputtext)
ax.set_ylabel(outputtext + ', field strength [V/m]') ax.set_ylabel(outputtext + ', field strength [V/m]')
#ax.set_ylim([-15, 20])
elif output == 'Ey': elif output == 'Ey':
ax = plt.subplot(gs[1, 0]) ax = plt.subplot(gs[1, 0])
ax.plot(time, outputdata,'r', lw=2, label=outputtext) ax.plot(time, outputdata,'r', lw=2, label=outputtext)
ax.set_ylabel(outputtext + ', field strength [V/m]') ax.set_ylabel(outputtext + ', field strength [V/m]')
#ax.set_ylim([-15, 20])
elif output == 'Ez': elif output == 'Ez':
ax = plt.subplot(gs[2, 0]) ax = plt.subplot(gs[2, 0])
ax.plot(time, outputdata,'r', lw=2, label=outputtext) ax.plot(time, outputdata,'r', lw=2, label=outputtext)
ax.set_ylabel(outputtext + ', field strength [V/m]') ax.set_ylabel(outputtext + ', field strength [V/m]')
#ax.set_ylim([-15, 20])
elif output == 'Hx': elif output == 'Hx':
ax = plt.subplot(gs[0, 1]) ax = plt.subplot(gs[0, 1])
ax.plot(time, outputdata,'g', lw=2, label=outputtext) ax.plot(time, outputdata,'g', lw=2, label=outputtext)
ax.set_ylabel(outputtext + ', field strength [A/m]') ax.set_ylabel(outputtext + ', field strength [A/m]')
#ax.set_ylim([-0.03, 0.03])
elif output == 'Hy': elif output == 'Hy':
ax = plt.subplot(gs[1, 1]) ax = plt.subplot(gs[1, 1])
ax.plot(time, outputdata,'g', lw=2, label=outputtext) ax.plot(time, outputdata,'g', lw=2, label=outputtext)
ax.set_ylabel(outputtext + ', field strength [A/m]') ax.set_ylabel(outputtext + ', field strength [A/m]')
#ax.set_ylim([-0.03, 0.03])
elif output == 'Hz': elif output == 'Hz':
ax = plt.subplot(gs[2, 1]) ax = plt.subplot(gs[2, 1])
ax.plot(time, outputdata,'g', lw=2, label=outputtext) ax.plot(time, outputdata,'g', lw=2, label=outputtext)
ax.set_ylabel(outputtext + ', field strength [A/m]') ax.set_ylabel(outputtext + ', field strength [A/m]')
#ax.set_ylim([-0.03, 0.03])
elif output == 'Ix': elif output == 'Ix':
ax = plt.subplot(gs[0, 2]) ax = plt.subplot(gs[0, 2])
ax.plot(time, outputdata,'b', lw=2, label=outputtext) ax.plot(time, outputdata,'b', lw=2, label=outputtext)