Example #1
0
 def distepi(lat_epi, lon_epi, lat, lon):
     """
         Returns the distance, in km, from the epicenter whose latitude, longitude 
         are lat_epi and lon_epi, respectively, and a generic location 
         whose latitude and longitude are lat and lon, respectively.
         All arguments are assumed in degrees. The returned distance is the arc
         of the great circle  passing for the epicenter and the given location
     """
     return deg2km(greatarc_dist(lat_epi, lon_epi, lat, lon))
Example #2
0
 def __call__(self, lat, lon):
     #TO BE IMPLEMENTED
     #getting constants:
     xEP1, xEP2, xEP3, xEP4, xEP5 = self.__constants
     
     dist_epi_sta = deg2km(greatarc_dist(self.lat, self.lon, lat, lon))
     
     M = self.m #magnitude
     depth = self.depth
     
     return xEP1*M + xEP2*log10(depth) + xEP3 - xEP4*0.5*log10((dist_epi_sta/depth)**2 + 1) - xEP5*(sqrt(dist_epi_sta**2 + depth**2)-depth)
Example #3
0
def centralasia_emca_1(M, lat_hyp, lon_hyp, depth_hyp, lat_sta, lon_sta): #h=Depth, R=distanza in R2 tra epicentro e stazione. DISTANZE IN KM!!!!
    """
        Computed intensity prediction equation for central asia.
        
        Parameters:
            M magnitude
            lat_hyp, lon_hyp, depth_hyp: coordinates of the hypocenter (latitude, longitude, depth, respectively). depth must be positive
            lat_sta, lon_sta: coordinates of the station (latitude and longitude, respectively)
            
        Returns:
            a number
    """
    
    dist_epi_sta = deg2km(greatarc_dist(lat_hyp,lon_hyp,lat_sta, lon_sta))
    return _centralasia_emca_1(M, depth_hyp,dist_epi_sta)
Example #4
0
def centralasia_emca_1(
    M, lat_hyp, lon_hyp, depth_hyp, lat_sta, lon_sta
):  #h=Depth, R=distanza in R2 tra epicentro e stazione. DISTANZE IN KM!!!!
    """
        Computed intensity prediction equation for central asia.
        
        Parameters:
            M magnitude
            lat_hyp, lon_hyp, depth_hyp: coordinates of the hypocenter (latitude, longitude, depth, respectively). depth must be positive
            lat_sta, lon_sta: coordinates of the station (latitude and longitude, respectively)
            
        Returns:
            a number
    """

    dist_epi_sta = deg2km(greatarc_dist(lat_hyp, lon_hyp, lat_sta, lon_sta))
    return _centralasia_emca_1(M, depth_hyp, dist_epi_sta)