search

ydluo / qdyn

A Quasi-DYNamic earthquake simulator
59 stars 28 forks source link

Output velocity time-series for each individual fault #89

Closed crpiceda closed 7 months ago

crpiceda commented 8 months ago

In the latest version of QDYN release 3.0.0, the output_vmax function inside output.f90 provides vmax values cor-responding to the entire mesh. While this is suitable for a single fault scenario, it becomes less practical when dealing with multiple faults, particularly in cases where events overlap across different faults. As a solution, this function was extended to output vmax timeseries for each individual fault. These data points are now stored in two additional columns per fault within the file output_fault: the index of element with peak slip rate (ivmax_fault) and the corresponding value (vmax_fault). This modification allows for more accurate event timing calculations on a per-fault basis. Furthermore, I've updated the wrapper to support the new output_fault format, and additionally, introduced two post-processing functions to extract the catalog and nucleation points.

crpiceda commented 7 months ago

the changes appear on the last commit ff00b63

crpiceda commented 7 months ago

I addressed the comments on the last commit 19a81c8