Commit 479998ed authored by Florian Wittkamp's avatar Florian Wittkamp

Added scripts to plot benchmark result

parent ceac2c36
......@@ -33,4 +33,6 @@ ld.Location='northwest';
title('JURECA')
ld.FontSize=16;
ax=gca;
ax.FontSize=16;
\ No newline at end of file
ax.FontSize=16;
grid on
grid minor
\ No newline at end of file
%% Plot a seismogram as imageplot and wiggleplot
% GNU General Public License v3.0
%
% Author: Florian Wittkamp 2016
%
%% Initialitation
clear all;
close all;
%% Parameter
Traces=251;
Samples=1000;
filename='../seismogram_bench_48_p.bin';
%% Read Seismogram
fid=fopen([filename],'r','ieee-le');
Seismogram=fread(fid,[Samples,Traces],'float'); fclose(fid);
%% Wiggleplot
figure
for t=1:Traces;
plot(Seismogram(:,t)./max(Seismogram(:,t))+t,'red')
hold on
end
axis([0 Samples 1 Traces+1])
xlabel('Samples')
ylabel('Traces')
% %% Imageplot
% figure
% imagesc(Seismogram,[-1e-10 1e-10])
% xlabel('Traces')
% ylabel('Samples')
% title(['Imageplot ',filename])
\ No newline at end of file
## plot_seismogram.py
# GNU General Public License v3.0
# Plot seismogram in a wiggle plot
#
# Florian Wittkamp 2016
## Initialisation
from numpy import *
import matplotlib.pyplot as plt
# Seismogram is written out as matrix
# Dimension of Matris is: Traces*Samples
Traces=251
Samples=1000
filename='../seismogram_bench_48_p.bin'
# Read in Seismogram
file = open(filename, 'rb')
Seismogram=fromfile(file,dtype=float32)
Seismogram=reshape(Seismogram,(Traces,Samples))
Seismogram=transpose(Seismogram)
# Trace normalization
for trace in range(Traces):
Seismogram[:,trace]=Seismogram[:,trace]/max(Seismogram[:,trace])
# Plotting
plt.figure(1)
for trace in range(Traces):
plt.plot(Seismogram[:,trace]+trace,color='black')
plt.xlim([0,Samples])
plt.xlabel('Samples')
plt.ylim([-1,Traces])
plt.ylabel('Traces')
plt.savefig('Seismogram.pdf')
plt.show()
# Plotting Seismogram as image plot
# plt.figure(1)
# imgplot=plt.imshow(Seismogram,aspect='auto')
# imgplot.set_clim(-1e-9,1e-9)
# plt.ylabel('Sample')
# plt.xlabel('Tracenumber')
# plt.draw()
......@@ -102,8 +102,8 @@
"Seismograms" : "comment",
"NDT, NDTSHIFT" : "1, 0",
"SEIS_FORMAT" : "1",
"SEIS_FILE" : "../su/seismogram_bench_1536",
"SEIS_FORMAT" : "3",
"SEIS_FILE" : "seismogram_bench_1536",
"Monitoring the simulation" : "comment",
"LOG_FILE" : "../log/bench_1536.log",
......
......@@ -102,8 +102,8 @@
"Seismograms" : "comment",
"NDT, NDTSHIFT" : "1, 0",
"SEIS_FORMAT" : "1",
"SEIS_FILE" : "../su/seismogram_bench_192",
"SEIS_FORMAT" : "3",
"SEIS_FILE" : "seismogram_bench_192",
"Monitoring the simulation" : "comment",
"LOG_FILE" : "../log/bench_192.log",
......
......@@ -102,8 +102,8 @@
"Seismograms" : "comment",
"NDT, NDTSHIFT" : "1, 0",
"SEIS_FORMAT" : "1",
"SEIS_FILE" : "../su/seismogram_bench_3000",
"SEIS_FORMAT" : "3",
"SEIS_FILE" : "seismogram_bench_3000",
"Monitoring the simulation" : "comment",
"LOG_FILE" : "../log/bench_3000.log",
......
......@@ -102,8 +102,8 @@
"Seismograms" : "comment",
"NDT, NDTSHIFT" : "1, 0",
"SEIS_FORMAT" : "1",
"SEIS_FILE" : "../su/seismogram_bench_384",
"SEIS_FORMAT" : "3",
"SEIS_FILE" : "seismogram_bench_384",
"Monitoring the simulation" : "comment",
"LOG_FILE" : "../log/bench_384.log",
......
......@@ -102,8 +102,8 @@
"Seismograms" : "comment",
"NDT, NDTSHIFT" : "1, 0",
"SEIS_FORMAT" : "1",
"SEIS_FILE" : "../su/seismogram_bench_48",
"SEIS_FORMAT" : "3",
"SEIS_FILE" : "seismogram_bench_48",
"Monitoring the simulation" : "comment",
"LOG_FILE" : "../log/bench_48.log",
......
......@@ -102,8 +102,8 @@
"Seismograms" : "comment",
"NDT, NDTSHIFT" : "1, 0",
"SEIS_FORMAT" : "1",
"SEIS_FILE" : "../su/seismogram_bench_768",
"SEIS_FORMAT" : "3",
"SEIS_FILE" : "seismogram_bench_768",
"Monitoring the simulation" : "comment",
"LOG_FILE" : "../log/bench_768.log",
......
......@@ -102,8 +102,8 @@
"Seismograms" : "comment",
"NDT, NDTSHIFT" : "1, 0",
"SEIS_FORMAT" : "1",
"SEIS_FILE" : "../su/seismogram_bench_96",
"SEIS_FORMAT" : "3",
"SEIS_FILE" : "seismogram_bench_96",
"Monitoring the simulation" : "comment",
"LOG_FILE" : "../log/bench_96.log",
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment