Skip to content

Commit

Permalink
Author: Dunyu Liu
Browse files Browse the repository at this point in the history
Date:    08/2020-10/08/2020
Contact: [email protected]
Version 1.2.0; Git tag v1.2.0; Parent 1.1.0.
Major changes:
* significant simplification of the code structure.
* files are renamed.
* New files library_bound.f90, library_output.f90, Read_Input_Files.f90, crs.f90, elemcal.f9
*	are created.
* Subroutines bound_load, bound_ft_ku, output_onfault_st, output_offfault_st, output_onfault_transfer,
*	output_timedy, output_globaldat, and output_prof  crs, elemcal, readcurrentcycle,
*	readstations1 and readstations2 are created.
* ./script/ is created to hold MATLAB post processing scripts.
* The license is improved.
* On-sreen printing is improved.
* [Verification] It is applied to SCEC SEAS BP5-QD.
* Detailed changes please refer to the Changelog.md/the google doc.
  • Loading branch information
dunyuliu committed Oct 8, 2020
1 parent 60da844 commit 66958e3
Show file tree
Hide file tree
Showing 47 changed files with 3,217 additions and 5,652 deletions.
68 changes: 68 additions & 0 deletions Changelog.md
Original file line number Diff line number Diff line change
@@ -1,3 +1,71 @@
# Version 1.2.0
*modified: LICENSE
*renamed: batch/AutoEQcycRSF.txt -> batch/ADArun.txt
*deleted: code/Q0/EQQuasi3d.f90
*deleted: code/Q0/Quasi-CRSelem.f90
*deleted: code/Q0/Quasi-c8g2.f90
*deleted: code/Q0/Quasi-faulting.f90
*deleted: code/Q0/Quasi-globalvar.f90
*deleted: code/Q0/Quasi-hrglss.f90
*deleted: code/Q0/Quasi-mesh4num.f90
*deleted: code/Q0/Quasi-meshgen.f90
*deleted: code/Q0/fric.f90
*deleted: code/Q0/makefile
*deleted: code/Qi/EQQuasi3d.f90
*deleted: code/Qi/Quasi-CRSelem.f90
*deleted: code/Qi/Quasi-SHL.f90
*deleted: code/Qi/Quasi-c8g2.f90
*deleted: code/Qi/Quasi-faulting.f90
*deleted: code/Qi/Quasi-globalvar.f90
*deleted: code/Qi/Quasi-hrglss.f90
*deleted: code/Qi/Quasi-mesh4num.f90
*deleted: code/Qi/Quasi-meshgen.f90
*deleted: code/Qi/makefile
*new file: code/Read_Input_Files.f90
*new file: code/c8g2.f90
*new file: code/crs.f90
*new file: code/elemcal.f90
*new file: code/eqquasi.f90
*new file: code/faulting.f90
*renamed: code/Qi/fric.f90 -> code/fric.f90
*new file: code/globalvar.f90
*new file: code/library_bound.f90
*new file: code/library_output.f90
*new file: code/main.f90
*new file: code/makefile
*new file: code/mesh4num.f90
*new file: code/meshgen.f90
*renamed: code/Q0/Quasi-SHL.f90 -> code/shl.f90
*new file: input/FE_Stations.txt
*new file: script/CombineGlobaldat.m
*new file: script/Convert_slip_stress_depth_prof.m
*new file: script/Convert_slip_stress_strike_prof.m
*new file: script/Make_mov_plot_Quasi.m
*new file: script/Plot_Globaldat.m
*new file: script/Plot_Single_onfault_st.m
*new file: script/Plot_rpt.m

#Changes:
New: A license is created.
New: Use one code to simulate all earthquake cycles. For the first cycle, initial conditions
are assigned. For later ones, initial conditions are imported from cplot_EQquasi.txt.
New: Significant simplification of the code structure. No variable is transferred between major subroutines.
New: files are renamed.
New: library_bound.f90 is created. It contains two subroutines bound_load and bound_ft_ku to
handle far_field loading and calculate KU on faults for use in subroutine faulting.
New: library_output.f90 is created. It contains several subroutines output_onfault_st,
output_offfault_st, output_onfault_transfer, output_timedy, output_globaldat, and output_prof.
Future outputs could be implemented here.
New: Read_Input_Files.f90 is copied from EQdyna 3D version 5.2.0 to import station locations and to
read in icstart from currentcycle.txt, which is created by the batch script. Other systematic
parameters should be implemented in a later time.
New: crs.f90 is created to handle the CRS format of stiffness structure construction.
New: elemcal.f90 is created to handle the construction of the stiffness matrix and mass array in
the CRS format.
New: ./script is created for MATLAB scripts for post processing.
Improvement: on-screen information printing is improved.
[Verification] SCEC SEAS benchmark BP-5-QD.

# Version 1.1.0
*modified: REAME.md
*modified: code/Q0/EQQuasi3d.f90
Expand Down
54 changes: 37 additions & 17 deletions LICENSE
Original file line number Diff line number Diff line change
@@ -1,23 +1,43 @@
Copyright (C) 2018-2020 Dunyu Liu, [email protected]
EQquasi Software Licence Agreement

If you copy, download, install or use EQquasi if you agree to be bound
by the terms and conditions of the following agreement.
/* Copyright (C) 2018-2020, Earthquake Modeling Lab @ Texas A&M University.
* All rights reserved.
* Authors: Benchun Duan ([email protected]), Dunyu Liu ([email protected])/

EQquasi Software Licence Agreement
If you copy, download, install or use EQquasi, you agree to be bound by the
terms and conditions of the following agreement:

'Software' shall mean the EQquasi in the binary excutable
format, together with a user guide, exemplary source code files included
in the package. The author grants to the licensee a non-exclusive, non-
transferable license to download, install and use the Software free of
charge for non-commercial, academic research, including industry-sponsored
research, and education purposes, except for the part of MUMPS, a parallel
iterative solver, which is distributed under CeCILL-C license (please see
https://cecill.info/licences/Licence_CeCILL-C_V1-en.html) abd the part of
AZTEC, a parallel iterative solver, which is distributed under the Aztec
License Agreesment that can be referred to at http://www.cs.sandia.gov/CRF/Aztec_download.html.
'Software' shall mean EQquasi in the binary excutable format, together with
a user guide, exemplary source code files included in the package. Permission
to use, copy, modify, and distribute this software and its documentation for
educational, research, and not-for-profit purposes, without fee and without a
signed licensing agreement, is hereby granted, except for that the MUMPS 5.1.2, the
parallel direct solver used by EQquasi, is released under CeCill-C license, which
can be referred to by the link https://cecill.info/licences/Licence_CeCILL-C_V1-en.html.)
provided that the above copyright notice, this paragraph and the following six
paragraphs appear in all copies, modifications, and distributions, Please contact
Dr. Benchun Duan by [email protected] for commercial licensing opportunities.

Disclaimer and Limitation of Warranty and Liability. The licensee acknowledge
that the Software is a research tool, provided with potential faults and bugs,
wihtout any maintenance, debugging, improvement and other support.
IN NO EVENT SHALL THE AUTHORS BE LIABLE TO ANY PARTY FOR DIRECT, INDIRECT, SPECIAL,
INCIDENTAL, OR CONSEQUENTIAL DAMAGES, INCLUDING LOST PROFITS, ARISING OUT OF
THE USE OF THIS SOFTWARE AND ITS DOCUMENTATION, EVEN IF REGENTS HAS BEEN ADVISED
OF THE POSSIBILITY OF SUCH DAMAGE.

THE AUTHORS SPECIFICALLY DISCLAIMS ANY WARRANTIES, INCLUDING, BUT NOT LIMITED TO,
THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
THE SOFTWARE AND ACCOMPANYING DOCUMENTATION, IF ANY, PROVIDED HEREUNDER IS
PROVIDED "AS IS". REGENTS HAS NO OBLIGATION TO PROVIDE MAINTENANCE, SUPPORT,
UPDATES, ENHANCEMENTS, OR MODIFICATIONS.

You can acknowlege (using references [1], [2], and [3]) the contribution of EQquasi in
any scientific publication dependent on the use of the package.

[1] D. Liu, B. Duan and B. Luo, EQsimu: a 3-D finite element dynamic earthquake simulator
for multicycle dynamics of geometrically complex faults governed by rate- and state-
dependent friction." Geophysical Journal International 220(1): 598-609 (2020).
[2] P. R. Amestoy, I. S. Duff, J. Koster and J.-Y. L'Excellent, A fully asynchronous
multifrontal solver using distributed dynamic scheduling, SIAM Journal on Matrix Analysis
and Applications, Vol 23, No 1, pp 15-41 (2001).
[3] P. R. Amestoy, A. Buttari, J.-Y. L'Excellent and T. Mary, Performance and scalability
of the block low-rank multifrontal factorization on multicore architectures, ACM
Transactions on Mathematical Software, Vol 45, Issue 1, pp 2:1-2:26 (2019)
18 changes: 18 additions & 0 deletions REAME.md
Original file line number Diff line number Diff line change
Expand Up @@ -9,6 +9,24 @@ post-seismic phases of an earthquake cycle. It adopts parallel solvers MUMPS and
handle the heavy computing loads. Rate- and state- friction law with various forms are imple-
mented.

### Author: Dunyu Liu
### Date: 08/2020-10/08/2020
### Contact: [email protected]
## Version 1.2.0; Git tag v1.2.0; Parent 1.1.0.
# Major changes:
* significant simplification of the code structure.
* files are renamed.
* New files library_bound.f90, library_output.f90, Read_Input_Files.f90, crs.f90, elemcal.f9
* are created.
* Subroutines bound_load, bound_ft_ku, output_onfault_st, output_offfault_st, output_onfault_transfer,
* output_timedy, output_globaldat, and output_prof crs, elemcal, readcurrentcycle,
* readstations1 and readstations2 are created.
* ./script/ is created to hold MATLAB post processing scripts.
* The license is improved.
* On-sreen printing is improved.
* [Verification] It is applied to SCEC SEAS BP5-QD.
* Detailed changes please refer to the Changelog.md/the google doc.

### Author: Dunyu Liu
### Date: 11/25/2019-01/06/2020
### Contact: [email protected]
Expand Down
6 changes: 3 additions & 3 deletions batch/AutoEQcycRSF.txt → batch/ADArun.txt
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
#BSUB -J M8AutoSharp-NormCap-NewSwitch
#BSUB -J jobname
#BSUB -L /bin/bash
#BSUB -W 20:00
#BSUB -n 100 # 100: total number of cores/jobslots to allocate for the job
Expand All @@ -13,9 +13,9 @@
##BSUB -B -N # send emails on job begin (-B) and end (-N)
#
export OMP_NUM_THREADS=1

#load intel toolchain
ml intel/2018a
#

for i in `seq 16 17`;
do
echo "Cycle #:" $i
Expand Down
88 changes: 0 additions & 88 deletions code/Q0/EQQuasi3d.f90

This file was deleted.

Loading

0 comments on commit 66958e3

Please sign in to comment.