RSQSim Restart and Processing Instructions - Revision history

Kmilner: /* Restarting an RSQSim simulation in order to extend a catalog from where it ended */

2019-07-09T01:29:21Z

Restarting an RSQSim simulation in order to extend a catalog from where it ended

2019-07-09T01:28:30Z

Restarting an RSQSim simulation in order to extend a catalog from where it ended

2019-02-23T05:27:51Z

Filtering a catalog by event time and removing events in the simulation runup

2019-02-23T05:27:16Z

Filtering a catalog by which fault(s) ruptured

2019-02-23T05:26:45Z

Plotting the fault model in plotFault3d

2019-02-23T05:24:46Z

2019-02-22T22:14:17Z

Instructions for restarting an RSQSim simulation in order to extend a catalog from where it ended

2019-02-22T22:12:44Z

Instructions for filtering a catalog by magnitude and saving a new filtered catalog

2019-02-22T22:05:34Z

Instructions for filtering a catalog by magnitude and saving a new filtered catalog

2019-02-22T22:04:18Z

Instructions for filtering a catalog by magnitude and saving a new filtered catalog

@@ Line 3: / Line 3: @@
 '''
-Note: You'll need some of the R scripts from RSQSimPostProcess
+Note: You'll need some of the R scripts from RSQSimPostProcess. All commands below are R commands.
 .  Create a new directory for the simulation you are setting up, and copy over the input file, pbs script, fault file (e.g. *.flt, or zfault_Deepen.in), *.KZero, and *.neighbors files. For this example, we'll call the previous catalog "cat1" and the extension "cat2"

@@ Line 5: / Line 5: @@
 Note: You'll need some of the R scripts from RSQSimPostProcess
-.  Create a new directory for the simulation you are setting up, and copy over the cat1.in, cat1.pbs, UCERF3.flt, UCERF3.KZero, and UCERF3.neighbors files.
+.  Create a new directory for the simulation you are setting up, and copy over the input file, pbs script, fault file (e.g. *.flt, or zfault_Deepen.in), *.KZero, and *.neighbors files. For this example, we'll call the previous catalog "cat1" and the extension "cat2"
-. Rename cat1.in and cat1.pbs to match the new catalog name.
+. Rename cat1.in and cat1.pbs to match the new catalog name (cat2.in and cat2.pbs).
 . Update the catalog name and input file name in the new cat2.pbs, and the outNameInfix in the new cat2.in files.
@@ Line 13: / Line 13: @@
 . Read in catalog 1:
 <pre>
-eqs = readEqs("home/user/name/catalog/eqs.catalog1.out")
+eqs = readEqs("/path/to/eqs.cat1.out")
 </pre>
 . Create the new input files for shear stress, normal stress, and theta from the end of catalog 1 and fill those in for the variables initTauFname, initSigmaFname, and initThetaFname in your cat2.in file:
@@ Line 48: / Line 48: @@
 . Update maxT in your cat2.in file to the length you want the extended catalog to be (in seconds).
 =='''Combining extended/restarted catalogs into one long catalog''' ==

@@ Line 125: / Line 125: @@
 </pre>
 . Choose the minimun and maximum catalog times you’re interested in:
 Note: The simulation runup time depends on the slip rates of your faults (you want each one to rupture at least once). If you are running a simulation with the UCERF3 California fault model and you’re interested in the San Andreas Fault, a reasonable runup time is 5,000 years.
 <pre>
@@ Line 141: / Line 142: @@
 inRange = which(eqs$t0yr > tmin & eqs$t0yr <tmax & eqs$M >= Mmin & eqs$M <=Mmax)
 </pre>
 =='''Filtering a catalog by which fault(s) ruptured'''==

@@ Line 150: / Line 150: @@
 </pre>
 . Make a list of the fault segments you are interested in:
 Note: If you are using the UCERF3 fault model, you can pass eqs$fault$faultName to the function unique() to get a list of the options.
 <pre>
@@ Line 164: / Line 165: @@
 use = euse[which(euse %in% inRange)]
 </pre>
 == '''Plotting the fault model in plotFault3d''' ==

@@ Line 171: / Line 171: @@
 . Initialize a variable for the plot colors:
 Note: This will create list of NA (which will show up as clear when you plot it) that is the length of the number of fault patches in the model.
 <pre>

@@ Line 1: / Line 1: @@
 '''
-== Instructions for restarting an RSQSim simulation in order to extend a catalog from where it ended ==
+== Restarting an RSQSim simulation in order to extend a catalog from where it ended ==
 '''
@@ Line 49: / Line 49: @@
 . Update maxT in your cat2.in file to the length you want the extended catalog to be (in seconds).
-=='''Instructions for combining extended/restarted catalogs into one long catalog''' ==
 Note: You'll need some of the R scripts from RSQSimPostProcess
@@ Line 78: / Line 79: @@
 load("combinedCat.RData")
 </pre>
 . If you also saved and need to combine your transitions file:
 <pre>
@@ Line 88: / Line 88: @@
-== '''Instructions for filtering a catalog by magnitude and saving a new filtered catalog''' ==
+== '''Filtering a catalog by magnitude and saving a new filtered catalog''' ==
 . Read the full catalog in (or load the eqs.RData file):
@@ Line 111: / Line 110: @@
 subName = "new_catalog_M7"
 writeEqsAndLists(eqsNew, outFnameInfix = paste(subName, sep = ""))
 </pre>

@@ Line 39: / Line 39: @@
 </pre>
 . Update pin from the list of patches that got locked in cat1 and write out the new pin file for cat2:
 <pre>
 system("grep Eliminating ../*e | awk '{print $4}' > ../cat1.locked")
 pin2 = scan(file="../cat1.locked")

@@ Line 92: / Line 92: @@
 . Read the full catalog in (or load the eqs.RData file):
 <pre>
 eqs = readEqs("eqs.cat.out")
 </pre>
-.Pick a minimum magnitude M for your new catalog (you'll get M=muse & greater):
+. Pick a minimum magnitude M for your new catalog (you'll get M=muse & greater):
 <pre>
 muse = 7

@@ Line 103: / Line 103: @@
 euse = which(eqs$M>muse)
 </pre>
-. Subset/filter the catalog by euse:
+. Subset/filter the catalog by euse (set renumberEvents=FALSE if you want to keep the original event ID's):
 <pre>
 eqsNew = subsetEqs(eqs,euse, renumberEvents=FALSE)