!Control file for program exsim_dmb
! Revision of program involving a change in the control file on this date:
   05/28/09
!Title
  Runs for comparing EXSIM and SMSIM 
!Write acc, psa, husid files for each site?
 Y
!MW, Stress, flag (0=fmax; 1=kappa), fmax or kappa
  7.0 140.0  1  0.005        	
!lat and lon of upper edge of fault
  0.0 0.0  
!strike,dip, depth of fault
  0.0 50.0 2.0            	
!fault type (S=strikeslip; R=reverse; N=normal; U=undifferentiated) 
! (Only used if Wells and Coppersmith is used to obtain FL and FW).
  R                               
!fault length and width, dl, dw, stress_ref
!Note: Force program to use Wells and Coppersmith (WC) for FL and/or FW if either entry = 0.0
! If Wells and Coppersmith are used to obtain FL and/or FW, the WC values are modified
! to account for the scaling implied by differences in the stress specified above and a
! stress that is assumed to be valid for the generic WC relations; this stress is stress_ref.
! The value of 70 bars is an educated guess, but it is not based on a quantitative analysis.
! The WC values of FL and/or FW are multiplied by the factor (stress_ref/stress)^(1/3).
  0  0  1.5 1.5 70.0 ! used for WC scaling (need an entry as a placeholder even if not used)
!y (=vrup/beta)
  0.8              
!hypo location in along fault and down dip distance from the fault 
!reference point (an upper corner)(-1.0, -1.0 for a random location); 
!number of iterations over hypocenter (need an entry, but only used if 
!either of the first two values are -1.0, indicating a random location)
 -1.0 -1.0 10                           
!Enter type of risetime (1=original, 2=1/f0)
 2
!tpadl, tpadt, delta t
 20.0 20.0 0.002               	
!beta , rho
  3.7 2.8                 	
!gsprd: r_ref, nsegs, (rlow(i), a_s, b_s, m_s(i))  (Usually set r_ref = 1.0 km)
    1.0
    3
      1.0 -1.3 0.0 6.5
     70.0 +0.2 0.0 6.5
    140.0 -0.5 0.0 6.5
!q: fr1, Qr1, s1, ft1, ft2, fr2, qr2, s2, c_q
   1.0 1000.0 0.0 1.4242 1.4242 1.0 893.0 0.32 3.7   
!trilinear duration and properties (rmin,rd1,rd2,durmin,b1,b2,b3)
  10. 70. 130. 0.0 0.16 -0.03 0.04 	
!Type of window: 1 for Saragoni-Hart taper windows, 0 for tapered boxcar window
  0                           	 
!low-cut filter corner, nslope (0 ==> no filter)
 0.05 8
! %damping of response spectra
 5.0
!# of f and Min and Max F for response spectra
  100 0.1   100.               	
!no. of frequencies for summary output (10 max):
 4 
!frequency (-1.0, 99.0 for pgv, pga):
 -1.0 99.0 0.5 5.0
!Output file names stem:
  M7.0_dl_1.5_dw_1.5_140b_10_20_pulse50_orig_rt_az045 
!Name of crustal amplification file:
  ab06_hr_crustal_amps.txt
!No. of frequencies in the crustal amplification file
  5                              
!Name of site amplification file:
  siteAmplification.txt
!No. of frequencies in the site amplification file (0 to skip)
  0                               
!DynamicFlag (0=no), PulsingPercent
  1   50.0                  	
!iflagscalefactor (1=vel^2; 2=acc^2; 3=asymptotic acc^2 (dmb))
  2                               
!iflagfas_avg (1=arithmetic; 2=geometric, 3=rms: USE 3!)
  3                               
!iflagpsa_avg (1=arithmetic; 2=geometric: USE 2!, 3=rms)
  2                               
!deterministic flag,gama,nu,t0, impulse peak
  0   1.0  90.0  4.0  10.		
!iseed, # of trials 
  309 20                    	
!Number of Sites, coord flag (1=lat,long; 2=R,Az; 3=N,E)
  54 2                      	
!Coordinates of each site
 1              45
 1.122018454	45
 1.258925412	45
 1.412537545	45
 1.584893192	45
 1.77827941	45
 1.995262315	45
 2.238721139	45
 2.511886432	45
 2.818382931	45
 3.16227766	45
 3.548133892	45
 3.981071706	45
 4.466835922	45
 5.011872336	45
 5.623413252	45
 6.309573445	45
 7.079457844	45
 7.943282347	45
 8.912509381	45
 10             45
 11.22018454	45
 12.58925412	45
 14.12537545	45
 15.84893192	45
 17.7827941	45
 19.95262315	45
 22.38721139	45
 25.11886432	45
 28.18382931	45
 31.6227766	45
 35.48133892	45
 39.81071706	45
 44.66835922	45
 50.11872336	45
 56.23413252	45
 63.09573445	45
 70.79457844	45
 79.43282347	45
 89.12509381	45
 100            45
 112.2018454	45
 125.8925412	45
 141.2537545	45
 158.4893192	45
 177.827941	45
 199.5262315	45
 223.8721139	45
 251.1886432	45
 281.8382931	45
 316.227766	45
 354.8133892	45
 398.1071706	45
 446.6835922	45
!If "Y" below and strike = 0.0:
!  if isitecoordflag = 2, move origin of the radial line to the midpoint of
!                         the top edge of the fault
!  if isitecoordflag = 3 and siteLocation(1) = 0, redefine siteLocation(1) = 0 to be
!                         the midpoint of the top edge of the fault (so that the
!                         sites will be along a line normal to the midpoint
!  if isitecoordflag = 3 and siteLocation(2) = 0, redefine siteLocation(1) = 0 to be the
!                         far end of the fault, so that the sites are along a line
!                         along the strike of the fault 
 Y
!islipweight=-1 -> unity slip for all subfaults,=0 -> specify slips below, =1 -> random weights
  -1                       	
!Matrix of slip weights
  1.0 1.0 1.0 1.0 1.0
  1.1 1.3 1.4 1.1 1.0
  1.1 1.4 1.5 1.2 1.2
  1.2 1.5 1.6 2.2 1.3