# VERSION: 1.32	DO NOT REMOVE OR CHANGE THIS LINE
###########################################################################
# ampfactor_file
# File containing the amplitude corrections for velocity and acceleration
# as a function of Q/T/M and amplitude; the path is relative to $shake_home;
# there is no default -- this field is require d if -qtm is specified
###########################################################################
ampfactor_file		: lib/sitecorr/site_corr_cdmg.dat

###########################################################################
# qtm_file
# File of site shear velocity as a function of position; the path is 
# relative to $shake_home; there is no default -- this field is required
# if -qtm is specified
###########################################################################
#qtm_file		: lib/sitecorr/usgs_cgs_geology_60s.bin
# New file from Eli 1/04
#qtm_file                : lib/sitecorr/waor45.bin

qtm_file                : lib/sitecorr/sea_hires.bin

###########################################################################
# stavel_file
# File of site shear velocity for seismic stations; the path is 
# relative to $shake_home; there is no default -- this field is required
# if -qtm is specified (unless you are running the Global ShakeMap system);
# the file format is:
#
#	"Latitude","Longitude","Site","Ptype"
#	lat,lon,station_code,shear_velocity
#	...
#
# Note that the first line of the file is ignored, and could be blank.
###########################################################################
stavel_file             : lib/grind/dig_pnw.txt
#stavel_file		: lib/grind/socalmasterstatlist.txt

###########################################################################
# Culling radii:
#	cthresh - epicenter to station distance within which the computed
#                 epicenter will not be used as a data point
#	pthresh - phantom (richter) station to real station distance within 
#                 which the phantom will not be used as a data point
#
# Defaults:
#	cthresh		 0.0
#	pthresh		 0.0
###########################################################################
#cthresh		: 15.0
#pthresh		: 10.0
# JJL -- One of the old parameters seems to have gone away here.  I'm going
# to put the old one here, but comment it out.
# creject is in the old file, but wasn't initially in this one.
# WTG - 17 Apr 06 - ditto
cthresh         : 15.0
pthresh         : 10.0
#creject         : 25.0



###########################################################################
# The input and output grid intervals specify the sample rate of the 
# interpolated grids manipulated (input_?_grid_interval) or output 
# (output_?_grid_interval) by grind; the default is 1.5 minutes
#
# input_x_grid_interval  : <floating point value><units>
# output_x_grid_interval : <floating point value><units>
# input_y_grid_interval  : <floating point value><units>
# output_y_grid_interval : <floating point value><units>
#
# Where <units> is one of:
#	m => minutes
#	c => seconds
#
# e.g.:
#	input_x_grid_interval	: 30c
#	output_x_grid_interval	: 2.5m
#
# (and both input_ and output_y_grid_interval would default to 1.5m).
###########################################################################
#input_x_grid_interval	: 1m
#output_x_grid_interval	: 1m
#input_y_grid_interval	: 1m
#output_y_grid_interval	: 1m

# Copied from old files by JJL 1-29-2004
# Copied from old files by WTG 17 Aug 06 - comments are JJL's
# Original input
# input_x_grid_interval : 1.5m
# changed for eli at 45 seconds
#input_x_grid_interval   : 45c
input_x_grid_interval   : 7.2c
#Original output changed by JJL because nearneighbor complains about the larger grid not dividing evenly into the map size.
#output_x_grid_interval  : 1.5m
#output_x_grid_interval  : 45c
output_x_grid_interval  : 7.2c
# Original
#input_y_grid_interval  : 1.5m
# changed for eli
#input_y_grid_interval   : 45c
input_y_grid_interval   : 7.2c
#Original output changed by JJL (see above comment)
#output_y_grid_interval  : 1.5m
#output_y_grid_interval  : 45c
output_y_grid_interval  : 7.2c



###########################################################################
# surface_tension provides the value passed to the GMT 'surface' program
# via the -T option; the value must be between zero and one; the default
# is 0.4, but we've changed it to 0.9 here because that seems to pull
# the amplitude at the edges of the map down to the data point better.
###########################################################################
# The one in the old file was 0.4, I'm leaving it here fo rhe moment. JJL 1-29-2004
# WTG 17 Apr 06 - and then must have changed it bacl - it was 0.9 when
#I copied it over
surface_tension		: 0.9

###########################################################################
# mapbuf determines the amount of padding to add to the edges of the
# mapped region to avoid quirkiness at the map edges; the default is
# 0.4 degrees to each edge
###########################################################################
mapbuf	: 0.0
#mapbuf	: 0.1
#mapbuf	: 0.4

###########################################################################
# lonspan specifies the width of the map in degrees of longitude; this 
#   value can be overridden on the command line; the default is 2.5
# latspan specifies the height of the map in degrees of latitude; this
#   value can be overridden on the command line; the default value is
#   'llratio * lonspan'
# llratio specifies the latspan/lonspan ratio to use when setting latspan
#   when latspan is not specified either here or on the command line; the 
#   default is the cosine of the event latitude.  Most regions will want
#   to specify a fixed ratio in order to generate maps of a consistent
#   shape.  Regions that span very large areas of latitude may find
#   it easier to not specify llratio below and let the default system
#   take over.
###########################################################################
# WTG - 17 Apr 06 - not specified in version on Magro

#llratio : 0.666667

###########################################################################
# bad_station specifies stations to flag as bad under certain circumstances
# the format of the statement is:
#
#	bad_station : code mag start_date-[end_date [mag start_date-[...]]]
#
# Where 'code' is the station code, 'mag' is the event magnitude cutoff 
# below which the station is considered bad, 'start_date' is the event 
# date to begin applying the cutoff, and 'end_date' is the event date 
# at which the cutoff no longer applies; dates are given in the yyyymmdd
# format, and are UTC (i.e. GMT) dates; a missing end date implys dates 
# inclusive of the current date; multiple 'mag start_date-end_date' groups 
# are allowed, e.g.:
#
#	bad_station : BC3 3.8 19990101-19990407 2.2 19990407-
#
# In the above example, the station 'BC3' will be flagged as bad for
# events smaller than 3.8 from January 1, 1999 to April 7, 1999 and
# for events smaller than 2.2 from April 7, 1999 to the present.  The
# station will not be flagged for events before January 1, 1999
###########################################################################
#bad_station             : RIN 9.9 20011001
#bad_station		: DAN 9.9 20010105-20010302

#bad_station		: MWC 4.0 19980101-
#bad_station		: BC3 3.8 19980101-
#bad_station		: PLS 4.1 19980101-
#bad_station		: BAR 3.8 19980101-
#bad_station		: ELC 3.8 19980101-
#bad_station		: MPM 9.9 19980101-

bad_station             : WMC 3.0 19980101-

###########################################################################
# regression specifies which regression relation to use for particular
# magnitude ranges. Regressions are defined by the name of the module that
# encapsulates its rules.
# 
# The format of the statement is:
#
#	regression : regression-module min-mag max-mag [depth-min depth-mag]
#
# Where 'regression-module' is the name of the module for the regression
# relation, 'min-mag' and 'max-mag' are the minimum and maximum magnitudes
# over which the relation should be used.
#
# For example:
#
#	regression : BJF97 0.0 9.9 0 30
#
# In the above example, the module Regression::BJF97 will be used to 
# calculate the regression values for events over the magnitude range of 
# 0.0 <= mag < 9.9 and for depths < 30 km.
#
# Multiple regression: lines can be used, in which case the first line
# that meets the current magnitude and depth parameters.
#
# If the event has no depth parameter then the regression: depth limits are 
# ignored.
#
###########################################################################
regression : HazusPGV 4.3 9.9 0 39
regression : Youngs97 0.0 9.9 39 100
regression : Small 0.0 4.3 0 39
#regression : HazusPGV 5.3 9.9 
#regression : Small 0.0 5.3

###########################################################################
# outlier_deviation_level specifies how many standard deviations define
# an outlier.
# 
# The format of the statement is:
#
#	outlier_deviation_level : level
#
# Where 'level' is the number of standard deviations defining outliers.
# This number can be any positive float. The default is 3.
#
# outlier_max_mag specifies a magnitude above which the automatic flagging
# of outliers will no longer take place (this cutoff will have no effect
# on manually flagged stations placed in the flagged_stations.txt file in
# the input directory (which always supersedes any automatic flagging);
# the purpose of this parameter is to prevent valid data from being 
# flagged because a long fault rupture might cause stations far from
# the epicenter but close to the rupture to show very high amplitudes;
# the default is 7.0; The format of the statement is:
#
#	outlier_max_mag : magnitude
#
# Examples:
#
#	outlier_deviation_level : 2
#	outlier_max_mag         : 8.2
#
###########################################################################
outlier_deviation_level : 3
outlier_max_mag         : 7.0

###########################################################################
# pga_accepted_halfwidth and pgv_accepted_halfwidth specify the minimum 
# regions above and below the pga and pgv regression curves in which data 
# values must be accepted (not flagged as outliers). This halfwidth 
# overrides the outlier bounds based on the standard deviation of the 
# regression curve, which may be very narrow, particularly at large distances.
# 
# The format of the statement is:
#
#	pga_accepted_halfwidth : halfwidth
#	pgv_accepted_halfwidth : halfwidth
#
# Where 'halfwidth' is a positive float. The default for both parameters
# is 0.01.
###########################################################################
# WTG - 17 Apr 06 set to values from Magro
#pga_accepted_halfwidth : 0.01
#pgv_accepted_halfwidth : 0.01
pga_accepted_halfwidth : 0.3
pgv_accepted_halfwidth : 0.3

###########################################################################
# estimate_x_grid_interval and estimate_y_grid_interval specify the size
# of grid elements in the grids of estimated amplitude values. Values are
# in units of degrees.
#
# The estimate grid is normally much coarser than the input or output
# grids. It is used to fill large gaps in the data distribution and is
# applied prior to interpolation to the "input" grid defined above.
# 
# The format of the statements is:
#
#	estimate_x_grid_interval : interval
#	estimate_y_grid_interval : interval
#
# Where 'interval' is a positive float. The default for both parameters
# is 0.3 degrees.
###########################################################################
# WTG - 17 Apr 06 set to values from Magro
#estimate_x_grid_interval : 0.1
#estimate_y_grid_interval : 0.1
estimate_x_grid_interval : 0.1
estimate_y_grid_interval : 0.1

###########################################################################
# bias parameters: because the attenuation relations are averages of events
# of varying mechanisms, they do no necessarily accurately reflect the
# mechanics of any specific event, nor are we guaranteed that the initial
# earthquake magnitude is accurate.  To overcome these deficiencies, we
# compute a factor by which the predicted ground motions are multiplied
# to bring them in line with the event data.  This factor is computed by
# minimizing the difference between the data values at the seismic stations
# and the estimated values at those locations.  The minimization is in 
# either a least squares sense or an absolute deviation sense.  Because
# there is naturally a lot of scatter in seismic data, the absolute
# deviation (i.e. L1 norm) seems like a better bet and, in fact, has
# proved to be so in practice.
#
# bias_norm: acceptable values are 'l2' (for least squares) or 'l1' (for
#   absolute deviation); the default is 'l1'.
# bias_max_range: distance (in kilometers) beyond which seismic stations 
#   will be excluded from the bias calculation; this helps to insure that 
#   the bias is computed using the (hopefully) more accurate near-source 
#   data; the default is 120km.
# bias_min_stations: sets the minimum number of seismic stations (within the 
#   search radius that are required to compute the bias; fewer than this 
#   number will result in the bias being set to 1.0, and a warning message 
#   being issued; the default is 6.
# bias_max_mag: earthquake magnitude above which the bias calculation is
#   not performed; this insures (for areas like southern California) that
#   a long fault rupture will not confuse the bias calculation; the default
#   is 7.0.
# bias_max_bias: the maximum value that the bias is allowed to take; this
#   is the maximum factor by which all estimates (including, possibly, the 
#   centroid) are multiplied; the default is 4.0; this parameter also
#   sets the minimum bias, which is (1.0 / bias_max_bias).
# bias_log_bias: determines whether the bias is calculated with linear
#   or log amplitude.  Valid values are 'true' and 'false'.  Defaults 
#   to false in order to maintain consistency with earlier versions.
#
# Examples:
#
#	bias_norm         : l1
#	bias_max_range    : 120
#	bias_min_stations : 6
#	bias_max_mag      : 7.0
#	bias_max_bias     : 4.0
#	bias_log_amp	  : true
#
###########################################################################
# WTG - 17 Apr 06 set to values from Magro
bias_norm         : l1
#bias_max_range    : 120
bias_max_range    : 150
#bias_min_stations : 6
bias_min_stations : 4
#bias_max_mag      : 7.0
bias_max_mag      : 6.5
bias_max_bias     : 4.0
bias_log_amp	  : false

###########################################################################
# For scenarios: 
# master list of stations for forward modeling (if not specified, 
#	uses stavel_file)
# filename to dump forward estimates
# forward modeling estimate grid intervals and number of grids
###########################################################################
#fwstatlist		: lib/grind/socalmasterstatlist.txt
fwstatlist              : lib/grind/dig_pnw.txt
fwdata_file		: forward.xml

###########################################################################
# These are the defaults for the estimate grid for scenarios
# (much finer than the usual richter coarse grid)
# fw_x_grid_interval : 0.075
# fw_y_grid_interval : 0.075
###########################################################################
fw_x_grid_interval : 0.075
fw_y_grid_interval : 0.075

###########################################################################
# rounding_fix: correction factor ADDED to lat/lon map boundary to fix
# roundoff errors if the map boundary is past the edge of the QTM grid.
# Playing with the rounding_fix parameter in grind.conf may fix this. 
# The order is N S E W
# Default is 0 0 0 0 (no correction)
#
# For Utah:
# rounding_fix: 0.000 0.0001 -0.0001 0.000
#
# (from Kris Pankow)
###########################################################################
#rounding_fix: 0.000 0.0001 -0.0001 0.000


###########################################################################
# strictbound: fixes the map in place and size given by the argument.
# There is no default and this parameter is not required.
# Usage:
#
#	strictbound : W/S/E/N
#
# For example:
#
#	strictbound : -118.5/33.5/-117.5/34.25
#
# Would center a 1.0 degree by 0.75 degree map approximately over Los
# Angeles.
#Seattle Shakemap
#47.34 - 47.92
#   121.9 122.55
#strictbound : -122.55/47.34/-121.9/47.92
###########################################################################
strictbound : -122.55 47.34 -121.9 47.92

