Commit 4b637692 authored by tilman.metz's avatar tilman.metz

added overnightbuilt benchmarks

parent 3b76f7a1
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
#!/bin/sh
# Seismogramme als eps erstellen
# input from command line
MODEL=$1
AUTOCLIP=$2
PERCENTAGE=99
#example command line call:
#./seismogramme.sh fullspace
#./seismogramme.sh fullspace clip
#if the command line argument AUTOCLIP='' then NO clip is specifically set
#instead the supswigb option perc=98 is used
#if the command line argument AUTOCLIP='clip' then a there will be a command line
#prompt that asked for constant clip for each individual component
exec 5<&0
#loop over components
while read value
do
case "$AUTOCLIP" in
clip)
#manual clip ( input of constant clip required)
echo "Clipwert: "
read CLIP <&5
supswigb < ../../../overnightbuilt/reference_value/$MODEL/su/$MODEL"_ref_"$value.su key=tracl hbox=9 wbox=9 \
label1="Time / s" label2="Tracenumber" title="Model $MODEL, $value component; clip=$CLIP" clip=$clip > $MODEL"_"$value.eps
;;
*)
#auto clip (no input of constant clip required)
#testrun in order to determine the clip at perc=$PERCENTAGE
#stderr will be redirected to file clip.txt
supswigb < ../../../overnightbuilt/reference_value/$MODEL/su/$MODEL"_ref_"$value.su key=tracl hbox=9 wbox=9 \
label1="Time / s" label2="Tracenumber" title="Model $MODEL, $value component" perc=$PERCENTAGE 1>$MODEL"_"$value.eps 2>clip.txt
#extract clip from file clip.txt
while read line;
do
if [ "$line" != "" ]; then
CLIP=`echo $line | cut -c 32-45`
fi
done < clip.txt
#use CLIP to add this CLIP into the title of the eps figure
supswigb < ../../../overnightbuilt/reference_value/$MODEL/su/$MODEL"_ref_"$value.su key=tracl hbox=9 wbox=9 \
label1="Time / s" label2="Tracenumber" title="Model $MODEL, $value component; clip=$CLIP" clip=$CLIP 1>$MODEL"_"$value.eps
;;
esac
done < components.txt
#!/bin/bash
./seismogramme.sh fullspace
./seismogramme.sh fullspace_el
./seismogramme.sh halfspace
./seismogramme.sh halfspace_el
./seismogramme.sh tunnel_el
./seismogramme_rec_array.sh rec_array
./seismogramme_rec_array.sh rec_array_el
#!/bin/sh
# Seismogramme als eps erstellen
# input from command line
MODEL=$1
AUTOCLIP=$2
PERCENTAGE=99
#example command line call:
#./seismogramme.sh fullspace
#./seismogramme.sh fullspace clip
#if the command line argument AUTOCLIP='' then NO clip is specifically set
#instead the supswigb option perc=98 is used
#if the command line argument AUTOCLIP='clip' then a there will be a command line
#prompt that asked for constant clip for each individual component
exec 5<&0
#loop over components
while read value
do
case "$AUTOCLIP" in
clip)
#manual clip ( input of constant clip required)
echo "Clipwert: "
read CLIP <&5
PROFILE="_gx_28_"
suwind < ../../../overnightbuilt/reference_value/$MODEL/su/$MODEL"_ref_"$value.su \
key=gelev min=2800000 max=2800000 | \
supswigb key=tracl hbox=9 wbox=9 \
label1="Time / s" label2="Receiver position / m" title="Model: $MODEL, Profile at 2800m, Component: $value, Clip: $clip" clip=$clip > $MODEL$PROFILE$value.eps
read clip <&5
PROFILE="_gx_38_"
suwind < ../../../overnightbuilt/reference_value/$MODEL/su/$MODEL"_ref_"$value.su \
key=gelev min=3800000 max=3800000 | \
supswigb key=tracl hbox=9 wbox=9 \
label1="Time / s" label2="Receiver position / m" title="Model: $MODEL, Profile at 3800m, Component: $value, Clip: $clip" clip=$clip > $MODEL$PROFILE$value.eps
;;
*)
#auto clip (no input of constant clip required)
#---- 2D receiver array in 2800 m depth
#testrun in order to determine the clip at perc=$PERCENTAGE
#stderr will be redirected to file clip.txt
PROFILE="_gx_28_"
suwind < ../../../overnightbuilt/reference_value/$MODEL/su/$MODEL"_ref_"$value.su \
key=gelev min=2800000 max=2800000 | \
supswigb key=tracl hbox=9 wbox=9 \
label1="Time / s" label2="Receiver position / m" title="Model $MODEL, Profile at 2800m, Component: $value" perc=$PERCENTAGE 1> $MODEL$PROFILE$value.eps 2>clip.txt
#extract clip from file clip.txt
while read line;
do
if [ "$line" != "" ]; then
CLIP=`echo $line | cut -c 32-45`
fi
done < clip.txt
#use CLIP to add this CLIP into the title of the eps figure
suwind < ../../../overnightbuilt/reference_value/$MODEL/su/$MODEL"_ref_"$value.su \
key=gelev min=2800000 max=2800000 | \
supswigb key=tracl hbox=9 wbox=9 \
label1="Time / s" label2="Receiver position / m" title="Model: $MODEL, 2800m depth, $value component; clip=$CLIP" clip=$CLIP 1> $MODEL$PROFILE$value.eps
#---- 2D receiver array in 3800 m depth
#testrun in order to determine the clip at perc=$PERCENTAGE
#stderr will be redirected to file clip.txt
PROFILE="_gx_38_"
suwind < ../../../overnightbuilt/reference_value/$MODEL/su/$MODEL"_ref_"$value.su \
key=gelev min=3800000 max=3800000 | \
supswigb key=tracl hbox=9 wbox=9 \
label1="Time / s" label2="Receiver position / m" title="Model: $MODEL, Profile at 3800m, Component: $value" perc=$PERCENTAGE 1> $MODEL$PROFILE$value.eps 2>clip.txt
#extract clip from file clip.txt
while read line;
do
if [ "$line" != "" ]; then
CLIP=`echo $line | cut -c 32-45`
fi
done < clip.txt
#use CLIP to add this CLIP into the title of the eps figure
suwind < ../../../overnightbuilt/reference_value/$MODEL/su/$MODEL"_ref_"$value.su \
key=gelev min=3800000 max=3800000 | \
supswigb key=tracl hbox=9 wbox=9 \
label1="Time / s" label2="Receiver position / m" title="Model $MODEL, 3800m depth, $value component; clip=$CLIP" clip=$CLIP 1> $MODEL$PROFILE$value.eps
;;
esac
done < components.txt
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This diff is collapsed.
\subsection{Homogeneous Fullspace}
As you can see in figure \ref{fullspace}, it consists of only one homogeneous material with a density of 2100 kg/m$^3$. The velocities in this model are $v_{P} = 5100$ m/s for the P-wave and $v_{S} = 2800$ m/s for the S-wave. \\
\begin{figure}[!h]
\centering
\includegraphics[width=0.475\textwidth]{fullspace.eps}
\caption{Homogeneous Fullspace}
\label{fullspace}
\end{figure}
The Euclidean distance between each gridpoint is equally spaced in each direction and is 54 meter (DH). The number of gridpoints (NX, NY) amounts in each direction 200 gridpoints. The model dimension thus account for a cube with a edge length of 10800 meters.\\
The source for the seismic waves is an explosive point source with a center source frequency of 5 Hz which is located in the middle of the cube at the coordinates (5400, 5400). To determine the waves, several geophones are used, which are arranged as a chain. The first receiver is located at (1620, 5400) and the last at (9180, 5400). Between each receiver there is an offset of 4 gridpoint, which is equal to 216.0 meters. So the whole chain consists of 35 geophones/receivers with a total length of 7506 meters. Furthermore, no free surface is included in the model which means that the model domain is sourrounded by an absorbing frame.\\
For the computation 4 processor cores are used. 2 cores are computing the expansion of the waves in the x-direction (NPROCX) and 2 in the y-direction (NPROCY). The total time of wave propagation is 4 seconds (TIME) and a single timestep is 0.004 seconds (DT).\\
The following SOFI3D features will be tested during this modeling run:
\begin{itemize}
\item The computational accuracy of your system with a with a equal number of processors for two directions
\item The accuracy of the explosive source
\end{itemize}
In the following, tables of all parameters used for this modeling are listed.\\
The first two tables contain the parameters for the size of the model
\vspace{0.5cm}
\begin{table}[h]
\begin{minipage}{0.4\textwidth}
\begin{center}
\begin{tabular}{c|c}
Direction & Value\\
\hline
NX & 200\\
NY & 200\\
\end{tabular}
\end{center}
\captionsetup{labelformat=empty}
\caption{Number of gridpoints in each direction}
\end{minipage}
\hfill
\begin{minipage}{0.4\textwidth}
\begin{center}
\begin{tabular}{c|c}
Direction & Value \\
\hline
DX & 54 m\\
DY & 54 m\\
\end{tabular}
\end{center}
\captionsetup{labelformat=empty}
\caption{Distance between gridpoints in each direction}
\end{minipage}
\end{table}
\vspace{1cm}
The next two tables contain both the parameters for the model like the velocities and the density and the number of cores used for each direction during the modeling.
\vspace{0.5cm}
\begin{table}[h]
\begin{minipage}{0.45\textwidth}
\begin{center}
\begin{tabular}{c|c}
Parameter & Value \\
\hline
$v_P$ & 5100 m/s\\
$v_S$ & 2800 m/s\\
$\rho$ & 2100 kg/$m^3$\\
\end{tabular}
\end{center}
\captionsetup{labelformat=empty}
\caption{Parameters of the model}
\end{minipage}
\hfill
\begin{minipage}{0.45\textwidth}
\begin{center}
\begin{tabular}{c|c}
Direction & Number of cores\\
\hline
NPROCX & 4\\
NPROCY & 1\\
\end{tabular}
\end{center}
\captionsetup{labelformat=empty}
\caption{Number of cores used}
\end{minipage}
\end{table}
\newpage
In the following figures, seismograms of the \textbf{viscoelastic} calculation of all components (curl, divergenc, pressure, $v_x$ and $v_y$) are shown.
\begin{figure*}[h]
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/fullspace_curl}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/fullspace_div}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/fullspace_p}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/fullspace_vx}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/fullspace_vy}
\end{minipage}
\end{figure*}
\newpage
In the following figures, seismograms of the \textbf{elastic} calculation of all components (curl, divergenc, pressure, $v_x$ and $v_y$) are shown.
\begin{figure*}[h]
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/fullspace_el_curl}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/fullspace_el_div}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/fullspace_el_p}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/fullspace_el_vx}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/fullspace_el_vy}
\end{minipage}
\end{figure*}
This diff is collapsed.
\subsection{Layered Half-Space}
As you can see in figure \ref{halfspace}, the model consists of two homogeneous material arranged as a layered half-space. The density of the materials are 1500 kg/m$^3$ for the upper layer and 2100 $\frac{kg}{m^3}$ for the half-space. The velocities in this model are $v_{P1} = 400 \frac{m}{s}$ for the P-wave and $v_{S2} = 210 \frac{m}{s}$ for the S-wave in the upper layer and $v_{P2} = 5100 $ m/s for the P-wave and $v_{S2} = 2800 $ m/s for the S-wave in the half-space\\
\begin{figure}[htbp]
\centering
\includegraphics[width=0.6\textwidth]{halfspace.eps}
\caption{Picture of the model \textbf{Layerd Halfspace} including the location of the source and receivers}
\label{halfspace}
\end{figure}
The Euclidean distance between each gridpoint is equally spaced in each direction and is 0.5 meter (DH). The number of gridpoints (NX, NY) amounts in each direction 100 gridpoints. The model dimension thus account for a cube with a edge length of 50 meters.\\
The source for the seismic waves is a force in z (in the vertical direction) with a center source frequency of 50 Hz which is located at the coordinates (15.5, 0.5). To determine the waves, several geophones are used, which are arranged as a chain. The first receiver is located at (15, 0.5) and the last at (35, 0.5). Between each receiver there is an offset of 1 gridpoint, which is equal to 0.5 meters. So the whole chain consists of 41 geophones/receivers with a total length of 20 meters. Furthermore, a free surface is included in the model which means that the model domain is sourrounded by an absorbing frame except the upper plane of the cube.\\
For the computation 4 processor cores are used. 2 cores are computing the expansion of the waves in the x-direction (NPROCX) and 2 in the y-direction (NPROCY). The total time of wave propagation is 0.15 seconds (TIME) and a single timestep is $1.4^{-4}$ seconds (DT).\\
In this model we specifie different values for FL and TAU for the viscoelastic calculation. We assume FL = 50 and TAU = 0.05.
The following SOFI2D features will be tested during this modeling run:
\begin{itemize}
\item The computational accuracy of your system with a equal number of processors for two directions
\item The accuracy of the free surface
\end{itemize}
In the following, tables of all parameters used for this modeling are listed.\\
The first two tables contain the parameters for the size of the model.
\vspace{0.5cm}
\begin{table}[h]
\begin{minipage}{0.4\textwidth}
\begin{center}
\begin{tabular}{c|c}
Direction & Value\\
\hline
NX & 100\\
NY & 100\\
\end{tabular}
\end{center}
\captionsetup{labelformat=empty}
\caption{Number of gridpoints in each direction}
\end{minipage}
\hfill
\begin{minipage}{0.4\textwidth}
\begin{center}
\begin{tabular}{c|c}
Direction & Value \\
\hline
DX & 0.5 m\\
DY & 0.5 m\\
\end{tabular}
\end{center}
\captionsetup{labelformat=empty}
\caption{Distance between gridpoints in each direction}
\end{minipage}
\end{table}
\vspace{1cm}
The next two tables contain both the parameters for the model like the velocities and the density and the number and the spreading of cores used for the modeling.
\vspace{0.5cm}
\begin{table}[h]
\begin{minipage}{0.45\textwidth}
\begin{center}
\begin{tabular}{c|c}
Parameter & Value \\
\hline
$v_{P1}$ & 500 m/s\\
$v_{S1}$ & 290 m/s\\
$\rho_{1}$ & 400 kg/$m^3$\\
$v_{P2}$ & 1500 m/s\\
$v_{S2}$ & 860 m/s\\
$\rho_{2}$ & 1200 kg/$m^3$\\
\end{tabular}
\end{center}
\captionsetup{labelformat=empty}
\caption{Parameters of the model}
\end{minipage}
\hfill
\begin{minipage}{0.45\textwidth}
\begin{center}
\begin{tabular}{c|c}
Direction & Number of cores\\
\hline
NPROCX & 2\\
NPROCY & 2\\
\end{tabular}
\end{center}
\captionsetup{labelformat=empty}
\caption{Number of cores used}
\end{minipage}
\end{table}
\newpage
In the following figures, seismograms of the \textbf{viscoelastic} calculation of all components (curl, divergence, pressure, $v_x$ and $v_y$) are shown.
\begin{figure*}[h]
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/halfspace_curl}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/halfspace_div}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/halfspace_p}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/halfspace_vx}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/halfspace_vy}
\end{minipage}
%
\end{figure*}
\newpage
In the following figures, seismograms of the \textbf{elastic} calculation of all components (curl, divergence, pressure, $v_x$ and $v_y$) are shown.
\begin{figure*}[h]
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/halfspace_el_curl}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/halfspace_el_div}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/halfspace_el_p}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/halfspace_el_vx}
\end{minipage}
%
\hfill
%
\begin{minipage}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{eps/halfspace_el_vy}
\end{minipage}
%
\end{figure*}
\ No newline at end of file
This diff is collapsed.
<
\subsection{Receiver Array}
As you can see in figure \ref{rec_array}, the model consists of one homogeneous material. The density of the material is $\rho =$ 2500 kg/m$^3$. The velocities in this model are $v_{P} = 5100$ m/s for the P-wave and $v_{S} = 3100$ m/s for the S-wave.\\
\begin{figure}[htbp]
\centering
\includegraphics[width=7cm]{rec_array.eps}
\caption{bla}
\label{rec_array}
\end{figure}
The Euclidean distance between each gridpoint is equally spaced in each direction and is 40 meter (DH). The number of gridpoints (NX, NY) differs in each direction. NX amounts 300 gridpoints and NY 200 gridpoints. The model dimension thus account for a cuboid with a edge length of 12.000 meters in x and 8.000 meters in y.\\
The source for the seismic waves is a plane wave excitation in a depth of 800 m. The center source frequency is 10 Hz. To determine the waves, a lot of geophones are used, which are arranged in two horizontally layered arrays (REC\_ARRAY). The first receiver array is located in a depth of 2.800 meters (REC\_ARRAY\_DEPTH). The vertical distance between each receiver array is 25 gridpoints or 1.000 meters (REC\_ARRAY\_DIST). According to this, the second array is located in a depth of 3.800 meters.
In each array the distance between the receivers are 7 gridpoint in the x-direction (DRX), hence the Euclidean distance between each geophone is 240 meters in x. Aditionally, a free surface is included in the model which means that the model domain is sourrounded by an absorbing frame except the upper plane of the cuboid.\\
For the computation 4 processor cores are used. 4 cores are computing the expansion of the waves in the x-direction (NPROCX) and 1 in the y-direction (NPROCY). The total time of wave propagation is 4 seconds (TIME) and a single timestep is $3.2\cdot 10^{-3}$ seconds (DT).\\
The following SOFI2D features will be tested during this modeling run:
\begin{itemize}
\item The computational accuracy of your system with a scattered number of processors for each direction
\item The accuracy of the free surface
\item The snapshot output
\end{itemize}
In the following, tables of all parameters used for this modeling are listed.\\
The first two tables contain the parameters for the size of the model
\vspace{0.5cm}
\begin{table}[!h]
\begin{minipage}{0.4\textwidth}
\begin{center}
\begin{tabular}{c|c}
Direction & Value\\
\hline
NX & 300\\
NY & 200\\
\end{tabular}
\end{center}
\captionsetup{labelformat=empty}
\caption{Number of gridpoints in each direction}
\end{minipage}
\hfill
\begin{minipage}{0.4\textwidth}
\begin{center}
\begin{tabular}{c|c}
Direction & Value \\
\hline
DX & 40 m\\
DY & 40 m\\
\end{tabular}