PROGRAM SHIFT_GSRM2p2x Program to create several versions (x = a, b, c, d, e, f,) of a long-term, indefinite-window global shallow seismicity map based on the Global Strain Rate Map, version 2.2 (GSRM2.2) created in 2016 by Corne Kreemer. The methods of scaling strain-rate to seismicity progress from simple to complex in the various versions, with the more complex versions using the set of methods known as Seismic Hazard Inferred From Tectonics (SHIFT), as presented by [Bird & Liu, 2007, Seismol. Res. Lett.] and adapted to global strain-rate maps by Bird et al. [2010, Seismol. Res. Lett.] and by Bird & Kreemer [2015, Bull. Seismol. Soc. Amer.]. By Peter Bird, UCLA, August 2019, for the SHIFT_GSRM2.2x project. Output will be in Peter Bird's .GRD format, in units of epicenters/m^2/s. Now, opening log file SHIFT_GSRM2.2x_log.txt. Creating global grids of 899 rows * 1440 columns = 1294560 points for input strain rate tensors, output epicentral rate densities, & mean cz. Total grid area = 5.10064E+14 square meters, compared to ideal sphere area of 5.10064E+14 square meters. Difference is due to 2 polar caps of 3.88E+08 square meters each. According to the GSRM2.2 data: Plate boundary area = 7.09330E+13 m^2 (13.907%). Intraplate area = 4.39131E+14 m^2 (86.093%). and we can assume the 2 polar caps are intraplate, as well. Enter name of (global, shallow) seismic catalog .EQC file for calibration: GCMT_shallow_m5p767_1977-2013.eqc Here are the first and last earthquakes in the catalog: GlobalCMT 1977.01.01 19:01:40.7 126.150 -2.880 43 6.35 3 46 83 157 7 316 GlobalCMT 2013.12.29 15:37:40.5 126.330 1.410 29 5.78 2 122 17 31 73 217 What is the duration of this catalog, in years? 37.0000000000000 What is the completeness-threshold magnitude of this catalog? 5.76700000000000 Catalog threshold (minimum) scalar moment = 5.018E+17 N m This threshold will be used for calibration of the forecast, regardless of the requested threshold for forecast reporting. PLEASE SELECT DESIRED MODEL VERSION FROM THE FOLLOWING LIST: 1(a): Uniform (calibrated) coupled thickness in all boundaries [2 DOF]. 2(b): Add isotropic smoothing of offshore non-SUB boundaries [2 DOF]. 3(c): Add seaward smoothing of SUB boundaries to outer rises [2 DOF]. 4(d): Add separate calibrations in tectonic zones #0~4 [5 DOF]. 5(e): Add classification & partitioning of strain rate tensors [5 DOF]. 6(f): Add velocity-dependence of SUB and CCB coupling [5 DOF]. CHOOSE INTEGER CODE: 2 Enter desired magnitude threshold for output forecast: 5.76700000000000 Forecast threshold (minimum) scalar moment = 5.018E+17 N m Create a .GRD file-name for the output: SHIFT_GSRM2.2b_m5.767_2014plus.grd File underwater.grd must be available. Please check this. File PB2002_steps.dat must be available. Please check this. PB2002_steps.dat has been converted to a list of 81114 internal points. Finding PB2002 plate boundary nearest to each active GSRM2 grid cell... Search completed. Isotropically smoothing rates adjacent to offshore non-SUB plate boundaries... Isotropic smoothing of offshore non-SUB plate boundaries completed. Earthquakes of magnitude >= 5.767 in plate boundaries: 6239 ( 96.34%) Earthquakes of magnitude >= 5.767 in plate interiors: 237 ( 3.66%) Intraplate seismicity rate (per m^2, per s) = 5.21456E-22 Global plate-boundary seismicity scaling = 1.79694 (relative to SUB). Writing seismicity (EQs/m^2/s) .grd file and computing area-integral... Area integral of seismicity = 5.5463E-06 earthquakes/s = 1.7503E+04 earthquakes/century. Writing effective_cz_meters .grd file... Writing moment_factors_in_SI.dat file... Total seismic moment rate in SI = 7.5009E+14 N m /s. Job completed.