docs: def get_median_area(self, mag, rake):
          """
          Calculates median fault area from magnitude.
          """
          if rake is None:
              # Return average of strike-slip and dip-slip curves
              return power(10.0, (mag - 4.185))
          elif (-45 <= rake <= 45) or (rake >= 135) or (rake <= -135):
              # strike-slip
              return power(10.0, (mag - 4.18))
          else:
              # Dip-slip (thrust or normal), and undefined rake
              return power(10.0, (mag - 4.19))