docs: def lambert_window(self, radius, lat0, long0):
          ''' Square Lambert Azimuthal equal area projection of
          a window centered at (lat0, long0) with a given radius (km).
  
          Args:
              radius(float): Radius of the window (km).
              lat0(float): Latitude at the center (degree).
              long0(float): Longitude at the center (degree).
  
          Returns:
              A tuple ``(longll, longtr, latll, lattr)` with ``longll``
              the longitude of the lower left corner, ``longtr`` the
              longitude of the top right corner, ``latll`` the latitude
              of the lower left corner and ``lattr`` the latitude of the
              top right corner.
  
          Note:
              All return coordinates are in degree
  
          '''
  
          radius = radius * 360.0 / (np.pi * 2 * 1734.4)
          radius = radius * np.pi / 180.0
          lat0 = lat0 * np.pi / 180.0
          long0 = long0 * np.pi / 180.0
  
          bot = self._kp_func(lat0 - radius, long0, lat0, long0)
          bot = bot * (np.cos(lat0) * np.sin(lat0 - radius) -
                       np.sin(lat0) * np.cos(lat0 - radius))
          x = bot
          y = bot
          rho = np.sqrt(x**2 + y**2)
          c = 2.0 * np.arcsin(rho / float(2.0))
          latll = np.arcsin(np.cos(c) * np.sin(lat0) + y * np.sin(c)
                            * np.cos(lat0) / rho) * float(180.0) / np.pi
          lon = long0 + np.arctan2(x * np.sin(c), rho * np.cos(lat0)
                                   * np.cos(c) - y * np.sin(lat0) * np.sin(c))
          longll = lon * 180.0 / np.pi
  
          x = -bot
          y = -bot
          rho = np.sqrt(x**2 + y**2)
          c = 2.0 * np.arcsin(rho / 2.0)
          lattr = np.arcsin(np.cos(c) * np.sin(lat0) + y * np.sin(c)
                            * np.cos(lat0) / rho) * float(180.0) / np.pi
          lon = long0 + np.arctan2(x * np.sin(c), rho * np.cos(lat0)
                                   * np.cos(c) - y * np.sin(lat0) * np.sin(c))
          longtr = lon * 180.0 / np.pi
  
          return longll, longtr, latll, lattr