docs: def _load_modeling_extent(self):
          """
          # Get extent from GSSHA Grid in LSM coordinates
          # Determine range within LSM Grid
          """
          ####
          # STEP 1: Get extent from GSSHA Grid in LSM coordinates
          ####
          # reproject GSSHA grid and get bounds
          min_x, max_x, min_y, max_y = self.gssha_grid.bounds(as_projection=self.xd.lsm.projection)
  
          # set subset indices
          self._set_subset_indices(min_y,
                                   max_y,
                                   min_x,
                                   max_x)