1. A computer-implemented method for obtaining a seismic risk analysis, the method comprising an act of: causing a computer processor to execute instructions specifically encoded on a memory, such that upon execution, the computer processor performs operations of: a. receiving, by the computer processor, as input a planning period t EBE , where the planning period is the period an investor uses in a financial analysis of an economic-basis earthquake for a particular facility; b. calculating, by the computer processor, a risk analysis parameter based upon the planning period t EBE ; and c. calculating, by the computer processor, an expected annualized loss (EAL) according to the following: 
 EAL= H ×PFL, wherein PFL represents a probable frequent loss, an average loss conditioned on a seismic intensity associated with the economic-basis earthquake, and H represents an economic hazard coefficient; and d. determining, by the computer processor, the EAL in order to provide an expected annualized value of repair cost to the particular facility in a particular location; wherein in the act of calculating a risk analysis parameter, the risk analysis parameter is hazard coefficient H, and the act of calculating hazard coefficient H further comprising acts of: a. calculating G EBE , where G EBE is an average exceedance frequency of s EBE , and where s EBE is the seismic intensity associated with an economic-basis earthquake; b. identifying grid points closest to ô for location ô=({circumflex over (φ)}, {circumflex over (λ)}), where {circumflex over (φ)} denotes a latitude designation and {circumflex over (λ)} denotes a longitude designation, the grid points closest to ô being referred to as the bounding grid points o j ; c. considering the period of interest {circumflex over (T)}, determining a period of at least one hazard curve, referred to as a bounding period T i ; d. letting s r,k denote the k th value (k=1, 2,. .. n s ) of s r in the hazard curve G(s r |o j ,c=c r ,T=T i ,ζ=ζ r ) and adjusting s r,k for each bounding grid point o j and bounding period T i to account for site classification ĉ and for damping ratio {circumflex over (ζ)}, where s is seismic intensity and where adjusting s r,k is done by multiplying s r,k by a site-classification coefficient F c and a damping coefficient F ζ , where by denoting by s k the k th value of s at the same location and period after adjusting for site classification and damping ratio, 
 G ( s k |o j ,c=ĉ,T=T i ,ζ={circumflex over (ζ)})= G ( F c F ζ s r,k |o j ,c=c r ,T=T i ,ζ=ζ r ): k= 1, 2 ,. .. n s ; e. letting G i,j,NZ denote the damping- and site-class-adjusted average exceedance frequency of s NZ given period T i and location o j , and for each T i and o j , interpolating each curve from “Act d” to determine G i,j,NZ ; f. letting G j,NZ denote the damping- and site-class-adjusted value of G NZ at the period of interest {circumflex over (T)} and location o j , and for each j, calculating G j,NZ ; g. letting G NZ denote the damping- and site-class-adjusted value of G NZ at the period of interest {circumflex over (T)} and location of interest ô, and calculating G NZ ; h. calculating risk analysis parameter H (i.e., hazard coefficient) according to the following: H ≡ G NZ ln ⁡ ( G NZ / G EBE ) , whereby through H, a user may calculate an expected annualized loss (EAL) according to the following: 
 EAL= H ×PFL, wherein PFL represents a probable frequent loss, mean loss conditioned on an occurrence of s EBE .