Proposal for damping flags

[brendanjmeade]

@aseyboldt Thanks for the idea to reuse existing fields in segment files for damping/regularization. To start with, the potential fields are:

• ss_rate: relevant to all faults
• ss_rate_sig: relevant to all faults
• ss_rate_flag: relevant to all faults
• ds_rate: relevant to non-vertical faults
• ds_rate_sig: relevant to non-vertical faults
• ds_rate_flag: relevant to non-vertical faults
• ts_rate: relevant to vertical faults
• ts_rate_sig: relevant to vertical faults
• ts_rate_flag: relevant to vertical faults

Proposal 1:

If a fault slip component is to be included in the damping set, *s_rate_flag = 2. Why 2? 1 still has a meaning for old-school soft constraints.

Thoughts? Alternatives?

[aseyboldt]

I think that's perfect.

[brendanjmeade]

@aseyboldt There is now a segment file with these damping flags (2) applied in the WNA folder:

https://github.com/brendanjmeade/wna/blob/main/data/segment/graham_nshm23_merge_connected_fixed_segments_3_constraints_1_borderlands_constraint_fix_locking_depths_damped.csv

I'd ask you to implement the indexing and damping of these segments only!

[brendanjmeade]

@aseyboldt This is running successfully. How is the damping weight set? Is this hardcoded or configurable via the <config.json> file? If it's the config file, what's the parameter name? If it's not yet configurable this way, can you/ @maresb / @elc45 implement the config level setting of this?

Speaking of config parameters this might also be a good opportunity to add config level support for the Matern prior (#269) and the annealing (#157). I think @maresb that config support for setting these parameters is all that those to other PRs are weighting for?

[aseyboldt]

It's configurable in the config with these two keys:

https://github.com/brendanjmeade/celeri/blob/main/celeri/config.py#L79
https://github.com/brendanjmeade/celeri/blob/main/celeri/config.py#L108

We should update the docs for the damping flags though, it still says that it refers to all of them...
(@elc45 if you have time? :-) )

[elc45]

Sure, can do!