Commit bb05c704 authored by Alex Breuer's avatar Alex Breuer

Added HSP2a input.

parent b36c7767
HSP2a
=====
This directory contains input and output for the three-dimensional HSP2a benchmark.
The benchmark tests the incorporation of viscoelastic attenuation.
The velocity model is homogeneous.
The source is a point source.
Further details are available from http://www.sismowine.org/.
* Contact: Alex Breuer
Config for the HSP2a Benchmark
==============================
This directory contains the config for the HSP2a benchmark.
* Contact: Alex Breuer
<?xml version="1.0" encoding="UTF-8"?>
<!--
@file This file is part of EDGE.
@author Alexander Breuer (anbreuer AT ucsd.edu)
@section LICENSE
Copyright (c) 2019, Alexander Breuer
Copyright (c) 2016-2017, Regents of the University of California
All rights reserved.
Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
3. Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission.
THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
@section DESCRIPTION
Config for the HSP2a benchmark.
-->
<edge>
<cfr>
<mesh>
<options>
<read>;</read>
</options>
<files>
<in>meshes/tet4_300.msh</in>
</files>
<boundary>
<outflow>105</outflow>
</boundary>
</mesh>
<velocity_model>
<domain>
<half_space>
<origin>
<x> 0</x>
<y> 0</y>
<z>-50000</z>
</origin>
<normal>
<x>0</x>
<y>0</y>
<z>1</z>
</normal>
</half_space>
<rho>2700</rho>
<lambda>32.4E9</lambda>
<mu>32.4E9</mu>
<qp>60</qp>
<qs>30</qs>
</domain>
</velocity_model>
<setups>
<attenuation>
<central_frequency>2.5</central_frequency>
<frequency_ratio>100</frequency_ratio>
</attenuation>
<point_sources>
<file>sources/hsp2a_uns_src.h5</file>
</point_sources>
<end_time>5.0</end_time>
</setups>
<output>
<receivers>
<path_to_dir>hsp2a_uns_300</path_to_dir>
<freq>0.005</freq>
<receiver>
<name>hsp2a_1</name>
<coords>
<x> 0</x>
<y> 693</y>
<z> 0</z>
</coords>
</receiver>
<receiver>
<name>hsp2a_2</name>
<coords>
<x> 0</x>
<y> 5543</y>
<z> 0</z>
</coords>
</receiver>
<receiver>
<name>hsp2a_3</name>
<coords>
<x> 0</x>
<y>10392</y>
<z> 0</z>
</coords>
</receiver>
<receiver>
<name>hsp2a_4</name>
<coords>
<x> 490</x>
<y> 490</y>
<z> 0</z>
</coords>
</receiver>
<receiver>
<name>hsp2a_5</name>
<coords>
<x> 3919</x>
<y> 3919</y>
<z> 0</z>
</coords>
</receiver>
<receiver>
<name>hsp2a_6</name>
<coords>
<x> 7348</x>
<y> 7348</y>
<z> 0</z>
</coords>
</receiver>
<receiver>
<name>hsp2a_7</name>
<coords>
<x> 400</x>
<y> 400</y>
<z> 400</z>
</coords>
</receiver>
<receiver>
<name>hsp2a_8</name>
<coords>
<x> 3200</x>
<y> 3200</y>
<z> 3200</z>
</coords>
</receiver>
<receiver>
<name>hsp2a_9</name>
<coords>
<x> 6000</x>
<y> 6000</y>
<z> 6000</z>
</coords>
</receiver>
<receiver>
<name>hsp2a_10</name>
<coords>
<x> 555</x>
<y> 370</y>
<z> 185</z>
</coords>
</receiver>
<receiver>
<name>hsp2a_11</name>
<coords>
<x> 4443</x>
<y> 2962</y>
<z> 1481</z>
</coords>
</receiver>
<receiver>
<name>hsp2a_12</name>
<coords>
<x> 8331</x>
<y> 5554</y>
<z> 2777</z>
</coords>
</receiver>
</receivers>
</output>
</cfr>
</edge>
.. image:: http://i.creativecommons.org/p/zero/1.0/88x31.png
:target: http://creativecommons.org/publicdomain/zero/1.0/
To the extent possible under law, `Alexander Nikolas Breuer <https://dial3343.org/>`_ has waived all copyright and related or neighboring rights to the following work:
* `tet4_300.tar.xz`
This work is published from: Germany.
Meshes for the HSP2a Benchmark
==============================
This directory contains the meshes for the HSP2a benchmark.
* Contact: Alex Breuer
Used Commands
=============
.. code-block:: bash
gmsh -format msh2 -3 tet4_300.geo 2>&1 | tee tet4_300.msh.log
/**
* @file This file is part of EDGE.
*
* @author Alexander Breuer (anbreuer AT ucsd.edu)
*
* @section LICENSE
* Copyright (c) 2019, Alexander Breuer
* Copyright (c) 2017, Regents of the University of California
* All rights reserved.
*
* Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
*
* 1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
*
* 2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
*
* 3. Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*
* @section DESCRIPTION
* Tet4 setup of HSP2a benchmark with a char length of 300 in the roi.
**/
// characteristic length at the source
cl1 = 300;
// remainder of the domain
cl2 = cl1*5;
Include "tet4_base.geo";
This source diff could not be displayed because it is stored in LFS. You can view the blob instead.
/**
* @file This file is part of EDGE.
*
* @author Alexander Breuer (anbreuer AT ucsd.edu)
*
* @section LICENSE
* Copyright (c) 2019, Alexander Breuer
* Copyright (c) 2016-2017, Regents of the University of California
* All rights reserved.
*
* Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
*
* 1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
*
* 2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
*
* 3. Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*
* @section DESCRIPTION
* Tet4 setup of HSP2a benchmark.
*
* The minimal circumscribing sphere of the receivers and source has:
* Radius: 5938.8
* Center: 2621.5 5196 1182.5
*
* Source is located at (0,0,0).
* Mesh is refined at (2622 5196 1183) with radius 8908 (1.5x the minimum) to characteristic length cl1.
* Remainder of sphere with radius 11878 (twice the minimum) gradually coarsens to cl2.
**/
// use OpenCascade Kernel for construction of the sphere
SetFactory("OpenCASCADE");
// mesh optimizations
Mesh.Optimize=1;
Sphere(1) = { 2622, 5196, 1183, 26000 };
// increase mesh resolution in sphere containing source and receivers
Field[1] = Ball;
Field[1].VIn = cl1;
Field[1].VOut = cl2;
Field[1].Radius = 8908;
Field[1].XCenter = 2622;
Field[1].YCenter = 5196;
Field[1].ZCenter = 1183;
// graduate coarsening around center of refinement
Point(101) = {2622,5196,1183};
Field[2] = Attractor;
Field[2].NodesList = {101};
Field[3] = Threshold;
Field[3].IField = 2;
Field[3].LcMin = cl1;
Field[3].LcMax = cl2;
Field[3].DistMin = 8908;
Field[3].DistMax = 11878;
Field[4] = Min;
Field[4].FieldsList = {1,3};
Background Field = 4;
// set outflow boundary conditions
Physical Surface(105) = {1};
// set volume
Physical Volume(1) = {1};
.. image:: http://i.creativecommons.org/p/zero/1.0/88x31.png
:target: http://creativecommons.org/publicdomain/zero/1.0/
To the extent possible under law, `Alexander Nikolas Breuer <https://dial3343.org/>`_ has waived all copyright and related or neighboring rights to the following work:
* `hsp2a_uns_src.h5`
This work is published from: Germany.
Source for the HSP2a Benchmark
==============================
This directory contains the point source for the HSP2a benchmark.
It was converted from the respective HSP1a source to HDF5.
EDGE-version for the script: ``64d041e185b67dd7e58e7855eba55457cfe9cc7e``
* Contact: Alex Breuer
Used Commands
-------------
.. code-block:: bash
wget https://opt.dial3343.org/3343/edge_opt/raw/master/bench/seismic/wp/hsp1a/hsp1a_uns_src.nc
nrfToHdf5.py -l 32.4E9 -m 32.4E9 -i hsp1a_uns_src.nc -o hsp2a_uns_src.h5 2>&1 | tee hsp2a_uns_src.h5.log
This source diff could not be displayed because it is stored in LFS. You can view the blob instead.
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