qdyn
qdyn copied to clipboard
ydluo
A Quasi-DYNamic earthquake simulator
# What's new in version 3.0.0? This is a major upgrade of the QDYN code base, with (some) compatibility-breaking changes made in both the Fortran code and the Python wrapper....
The restarted simulation differs from the seamless result. I put an example from the QDYN documentation 'Single asperity simulations (3D)' by slightly changing the input. ```python import matplotlib.pyplot as plt...
In `release 3.0.0`, `output_ox_last` is only written at the timestep `N` where `pb%it == pb%itstop` (`N=-1`). It'd be useful to set `N = -1` by default, which signals that a...
Currently QDYN internally generates a mesh, except for 3D MPI simulations. Since the Python wrapper already pre-computes a mesh with default settings, we might as well shift the burden from...
The variable `iot_buf` is intended to be used only by the master processor, but this seems not be the case in the module `output.f90`. This might be the reason behind...
In the single/double asperity example notebooks, `L` is defined as `Length of fault / nucleation length`. This is not correct, and should read `Length of fault / asperity length`. This...
Hi Team, why is number of mesh elements N: # Find next power of two for number of mesh elements N = int(np.power(2, np.ceil(np.log2(resolution * L / Lb)))) ? Where...
See [the documentation](https://ydluo.github.io/qdyn/running_simulations.html#simulation-input-parameters). Under "Discretization and accuracy parameters", the table entry `NPROCS` should be `NPROC`
In [`derivs_all.f90:151`](https://github.com/ydluo/qdyn/blob/master/src/derivs_all.f90#L151) the time derivative of the normal stress needs to be added, i.e.: ``` dmain_var = ( dtau_per + dtau_dt - sigma*dmu_dtheta*dth_dt - dtau_dP*dP_dt ) & / ( sigma*dmu_dv...
# Introduction Regardless of the time frame over which the HDF binary output is scheduled to be implemented, I think it would be useful at some point to explore this...