use seissol v1.3/Rusanov fluxes/ vtkhdf/ extend tpv13

Dockerfile

@@ -86,20 +86,15 @@ RUN git clone --recursive https://github.com/TUM-I5/ASAGI.git \
     && CC=mpicc CXX=mpicxx cmake .. -DCMAKE_INSTALL_PREFIX=/home/tools -DSHARED_LIB=off -DSTATIC_LIB=on -DNONUMA=on \
     && make -j$(nproc) && make install
 
-RUN git clone https://github.com/uphoffc/ImpalaJIT.git \
-    && cd ImpalaJIT \
-    && mkdir build && cd build \
-    && cmake .. && make -j $(nproc) install 
-
 RUN git clone https://github.com/SeisSol/easi \
     && cd easi \
     && mkdir build && cd build \
-    && CC=mpicc CXX=mpicxx cmake .. -DEASICUBE=OFF -DLUA=ON -DCMAKE_PREFIX_PATH=/home/tools -DCMAKE_INSTALL_PREFIX=/home/tools -DASAGI=ON -DIMPALAJIT=ON .. \
+    && CC=mpicc CXX=mpicxx cmake .. -DEASICUBE=OFF -DLUA=ON -DCMAKE_PREFIX_PATH=/home/tools -DCMAKE_INSTALL_PREFIX=/home/tools -DASAGI=ON -DIMPALAJIT=OFF .. \
     && make -j$(nproc) && make install
 
 RUN pip install numpy && docker-clean
 
-RUN git clone --recursive --depth 1 --single-branch --branch v1.1.3 https://github.com/SeisSol/SeisSol.git \
+RUN git clone --recursive --depth 1 --single-branch --branch v1.2.0 https://github.com/SeisSol/SeisSol.git \
     && cd SeisSol \
     && mkdir build_hsw && cd build_hsw \
     && export PATH=$PATH:/home/tools/bin \

Dockerfile_jupyterlab

@@ -93,7 +93,7 @@ RUN git clone https://github.com/SeisSol/easi \
     && CC=mpicc CXX=mpicxx cmake .. -DEASICUBE=OFF -DLUA=ON -DCMAKE_PREFIX_PATH=/home/tools -DCMAKE_INSTALL_PREFIX=/home/tools -DASAGI=ON -DIMPALAJIT=OFF .. \
     && make -j$(nproc) && make install
 
-RUN pip install numpy && pip install git+https://github.com/SeisSol/PSpaMM.git@davschneller/compile-fixes && docker-clean
+RUN pip install numpy && pip install git+https://github.com/SeisSol/PSpaMM.git && docker-clean
 
 #RUN git clone --depth 1 --single-branch --branch main https://github.com/libxsmm/libxsmm.git \
 #    && cd libxsmm \
@@ -104,7 +104,7 @@ RUN pip install numpy && pip install git+https://github.com/SeisSol/PSpaMM.git@d
 #       -j$(nproc) generator \
 #    && cp bin/libxsmm_gemm_generator /home/tools/bin
 
-RUN git clone --recursive --depth 1 --single-branch --branch v1.1.3 https://github.com/SeisSol/SeisSol.git \
+RUN git clone --recursive --depth 1 --single-branch --branch v1.2.0 https://github.com/SeisSol/SeisSol.git \
     && cd SeisSol \
     && mkdir build_hsw && cd build_hsw \
     && export PATH=$PATH:/home/tools/bin \

kaikoura/parametersLSW.par

@@ -7,6 +7,8 @@ Tv=0.05
 !Attenuation parameters (ignored if not compiled with attenuation)
 FreqCentral=0.5
 FreqRatio=100
+numflux = 'rusanov' ! The numerical flux. Values: 'godunov', 'rusanov'
+numfluxnearfault = 'rusanov' ! The numerical flux for faces of cells adjacent to the fault (on cells which have at least one DR face), excluding the DR faces themselves. Values: 'godunov', 'rusanov'
 /
 
 &IniCondition

kaikoura/parametersRS.par

@@ -7,6 +7,9 @@ Tv=0.05
 !Attenuation parameters (ignored if not compiled with attenuation)
 FreqCentral=0.5
 FreqRatio=100
+
+numflux = 'rusanov' ! The numerical flux. Values: 'godunov', 'rusanov'
+numfluxnearfault = 'rusanov' ! The numerical flux for faces of cells adjacent to the fault (on cells which have at least one DR face), excluding the DR faces themselves. Values: 'godunov', 'rusanov'
 /
 
 &IniCondition

sulawesi/parametersLSW.par

@@ -7,6 +7,8 @@ Tv=0.05
 !Attenuation parameters (ignored if not compiled with attenuation)
 FreqCentral=0.3
 FreqRatio=100
+numflux = 'rusanov' ! The numerical flux. Values: 'godunov', 'rusanov'
+numfluxnearfault = 'rusanov' ! The numerical flux for faces of cells adjacent to the fault (on cells which have at least one DR face), excluding the DR faces themselves. Values: 'godunov', 'rusanov'
 /
 
 &IniCondition

sulawesi/parametersRS.par

@@ -7,6 +7,8 @@ Tv=0.05
 !Attenuation parameters (ignored if not compiled with attenuation)
 FreqCentral=0.3
 FreqRatio=100
+numflux = 'rusanov' ! The numerical flux. Values: 'godunov', 'rusanov'
+numfluxnearfault = 'rusanov' ! The numerical flux for faces of cells adjacent to the fault (on cells which have at least one DR face), excluding the DR faces themselves. Values: 'godunov', 'rusanov'
 /
 
 &IniCondition

tpv13/parameters.par

@@ -4,6 +4,8 @@ MaterialFileName = 'tpv12_13_material.yaml'
 !enable off-fault plasticity (ignored in Plasticity=0)
 Plasticity = 1
 Tv = 0.03
+numflux = 'rusanov' ! The numerical flux. Values: 'godunov', 'rusanov'
+numfluxnearfault = 'rusanov' ! The numerical flux for faces of cells adjacent to the fault (on cells which have at least one DR face), excluding the DR faces themselves. Values: 'godunov', 'rusanov'
 /
 
 &IniCondition

tpv13/tpv13_training.geo

@@ -10,15 +10,15 @@ SetFactory("OpenCASCADE");
 
 // Length of the fault
 l_f = 30e3;
-// Depth of the fault
+// Width of the fault (along-dip)
 w_f = 15e3;
 dip_rad = dip*Pi/180.;
 
 // Explicitly mesh the square nucleation patch, defined here in fault local coordinates 
 // Hypocenter, z_n is here the depth along-dip
 x_n = 0e3;
 z_n = -12e3;
-// size of the nucleaton patch
+// Dimensions of the nucleation patch
 l_n = 3e3;
 w_n = 3e3;
 
@@ -32,7 +32,7 @@ Z0 = -42e3;
 // Create the domain as a box
 domain = newv; Box(domain) = {X0, Y0, Z0, X1-X0, Y1-Y0, -Z0};
 
-// Create the fault as a vertically dipping rectangle, centered in x at the hypocenter
+// Create the fault as a in the x-y plane, centered in x at the hypocenter
 fault = news; Rectangle(fault) = {-l_f/2, -w_f, 0, l_f, w_f};
 
 // Create the nucleation patch as a smaller rectangle
@@ -44,7 +44,7 @@ Rotate{ {1, 0, 0}, {0, 0, 0}, dip_rad } { Surface{fault, nucl}; }
 // Intersect the domain box with the fault rectangle at the free surface
 v() = BooleanFragments{ Volume{domain}; Delete; }{ Surface{fault, nucl}; Delete; };
 
-// Update all coordinates that define important surfaces within the mesh
+// Extract surfaces within the mesh for defining boundary conditions
 eps = 1e-3;
 fault_final[] = Surface In BoundingBox{-l_f/2-eps, -w_f-eps, -w_f-eps, l_f/2+eps, w_f+eps, w_f+eps};
 top[] = Surface In BoundingBox{X0-eps, Y0-eps, -eps, X1+eps, Y1+eps, eps};
@@ -66,4 +66,3 @@ Physical Surface(3) = {fault_final[]};
 Physical Surface(5) = {other[]};
 
 Physical Volume(1) = {domain};
-Mesh.MshFileVersion = 2.2;