import NED right away

Author: petrasovaa

content/tutorials/parallelization/SIMWE_parallelization.qmd

@@ -282,7 +282,7 @@ tools.g_extension(extension="r.in.usgs")
 ```
 
 First, we will set the computational region to match the river subwatersheds and then list available datasets
-for that region:
+for that region (the output at the time of creating this tutorial, it may vary later on):
 
 ```{python}
 tools.g_region(vector="river_basins")
@@ -297,41 +297,42 @@ print(
 ```
 
 ```text
-    USGS file(s) to download:
-    -------------------------
-    Total download size: 4.83 GB
-    Tile count: 16
-    USGS SRS: wgs84
-    USGS tile titles:
-    USGS 1/3 Arc Second n36w081 20240510
-    USGS 1/3 Arc Second n36w082 20220504
-    USGS 1/3 Arc Second n36w082 20220512
-    USGS 1/3 Arc Second n36w082 20240510
-    USGS 1/3 Arc Second n37w080 20210305
-    USGS 1/3 Arc Second n37w081 20210305
-    USGS 1/3 Arc Second n37w081 20240510
-    USGS 1/3 Arc Second n37w082 20210305
-    USGS 1/3 Arc Second n37w082 20220512
-    USGS 1/3 Arc Second n37w082 20240510
-    USGS 13 arc-second n36w080 1 x 1 degree
-    USGS 13 arc-second n36w081 1 x 1 degree
-    USGS 13 arc-second n36w082 1 x 1 degree
-    USGS 13 arc-second n37w080 1 x 1 degree
-    USGS 13 arc-second n37w081 1 x 1 degree
-    USGS 13 arc-second n37w082 1 x 1 degree
-    -------------------------
+USGS file(s) to download:
+-------------------------
+Total download size:	5.32 GB
+Tile count:	17
+USGS SRS:	wgs84
+USGS tile titles:
+USGS 1/3 Arc Second n36w080 20240611
+USGS 1/3 Arc Second n36w081 20240611
+USGS 1/3 Arc Second n36w082 20220504
+USGS 1/3 Arc Second n36w082 20220512
+USGS 1/3 Arc Second n36w082 20240611
+USGS 1/3 Arc Second n37w080 20210305
+USGS 1/3 Arc Second n37w081 20210305
+USGS 1/3 Arc Second n37w081 20240611
+USGS 1/3 Arc Second n37w082 20210305
+USGS 1/3 Arc Second n37w082 20220512
+USGS 1/3 Arc Second n37w082 20240611
+USGS 13 arc-second n36w080 1 x 1 degree
+USGS 13 arc-second n36w081 1 x 1 degree
+USGS 13 arc-second n36w082 1 x 1 degree
+USGS 13 arc-second n37w080 1 x 1 degree
+USGS 13 arc-second n37w081 1 x 1 degree
+USGS 13 arc-second n37w082 1 x 1 degree
+-------------------------
 ```
 
 We will select the 2024 data and use the `title_filter` option to filter them.
-We will first only download the datasets and use them later.
+We will download and import the datasets.
 
 ```{python}
 tools.r_in_usgs(
     product="ned",
     ned_dataset="ned13sec",
-    flags="d",
     output_directory=os.getcwd(),
-    title_filter="20240510",
+    title_filter="20240611",
+    output_name="ned",
 )
 ```
 
@@ -362,31 +363,22 @@ Set the computational region with [g.region](https://grass.osgeo.org/grass-devel
 tools.g_region(vector=f"basin_{huc12}")
 ```
 
-Import the NED data for that subwatershed. r.in.usgs will use already downloaded data.
+Display the NED data for that region:
 
 ```{python}
-tools.r_in_usgs(
-    product="ned",
-    ned_dataset="ned13sec",
-    cache_size_tolerance=int(1e10),
-    output_directory=os.getcwd(),
-    title_filter="20240510",
-    output_name=f"ned_{huc12}",
-)
-
-basin_map = gj.Map()
-basin_map.d_rast(map=f"ned_{huc12}")
+basin_map = gj.Map(current_region=True)
+basin_map.d_rast(map="ned")
 basin_map.d_vect(map=f"basin_{huc12}", fill_color="none")
 basin_map.show()
 ```
 
 ![](SIMWE_images/basin_030401010402.webp)
 
-Set the region to match the elevation raster.
+Set the region to match the extent of the basin and align resolution with the elevation raster.
 Convert the selected vector to a raster using [v.to.rast](https://grass.osgeo.org/grass-devel/manuals/v.to.rast.html).
 
 ```{python}
-tools.g_region(raster=f"ned_{huc12}")
+tools.g_region(raster="ned", vector=f"basin_{huc12}")
 tools.v_to_rast(input=f"basin_{huc12}", output=f"basin_{huc12}", use="val")
 ```
 
@@ -397,10 +389,10 @@ MaskManager is new in GRASS version 8.5. Alternatively, you can use [r.mask](htt
 mask = gs.MaskManager(mask_name=f"basin_{huc12}")
 mask.activate()
 
-basin_map = gj.Map()
-basin_map.d_rast(map=f"ned_{huc12}")
+basin_map = gj.Map(current_region=True)
+basin_map.d_rast(map="ned")
 basin_map.d_legend(
-    raster=f"ned_{huc12}",
+    raster="ned",
     flags="t",
     at=[10, 15, 40, 95],
     title="Elevation [m]",
@@ -476,7 +468,7 @@ tools.r_recode(input=f"nlcd_{huc12}", output=f"runoff_{huc12}", rules="runoff.tx
 Then, compute slope components (dx and dy) from the elevation raster to support hydrologic modeling.
 
 ```{python}
-tools.r_slope_aspect(elevation=f"ned_{huc12}", dx=f"dx_{huc12}", dy=f"dy_{huc12}")
+tools.r_slope_aspect(elevation="ned", dx=f"dx_{huc12}", dy=f"dy_{huc12}")
 ```
 
 Run hydrologic and sediment simulations for the selected subwatershed.
@@ -487,7 +479,7 @@ Running the simulation may take a while.
 ```{python}
 simulation_time = 180
 tools.r_sim_water(
-    elevation=f"ned_{huc12}",
+    elevation="ned",
     dx=f"dx_{huc12}",
     dy=f"dy_{huc12}",
     depth=f"depth_{huc12}",
@@ -519,7 +511,7 @@ tools.r_mapcalc(expression=f"transport_capacity_{huc12} = 0.001")
 tools.r_mapcalc(expression=f"detachment_coef_{huc12} = 0.001")
 tools.r_mapcalc(expression=f"shear_stress_{huc12} = 0.01")
 tools.r_sim_sediment(
-    elevation=f"ned_{huc12}",
+    elevation="ned",
     water_depth=f"depth_{huc12}",
     dx=f"dx_{huc12}",
     dy=f"dy_{huc12}",
@@ -645,17 +637,8 @@ def compute(huc12):
         where=f"huc12 == '{huc12}'",
         flags="t",
     )
-    # Set the computational region to match non-null area in the new raster.
-    with gs.RegionManager(vector=f"basin_{huc12}"):
-        tools.r_in_usgs(
-            product="ned",
-            ned_dataset="ned13sec",
-            cache_size_tolerance=int(1e10),
-            output_directory=os.getcwd(),
-            title_filter="20240510",
-            output_name=f"ned_{huc12}",
-        )
-    with gs.RegionManager(raster=f"ned_{huc12}"):
+    # Set the computational region to match the basin
+    with gs.RegionManager(vector=f"basin_{huc12}", raster="ned"):
         tools.v_to_rast(
             input=f"basin_{huc12}",
             output=f"basin_{huc12}",
@@ -680,12 +663,12 @@ def compute(huc12):
                 rules="runoff.txt",
             )
             tools.r_slope_aspect(
-                elevation=f"ned_{huc12}",
+                elevation="ned",
                 dx=f"dx_{huc12}",
                 dy=f"dy_{huc12}",
             )
             tools.r_sim_water(
-                elevation=f"ned_{huc12}",
+                elevation="ned",
                 dx=f"dx_{huc12}",
                 dy=f"dy_{huc12}",
                 depth=f"depth_{huc12}",
@@ -698,7 +681,7 @@ def compute(huc12):
             tools.r_mapcalc(expression=f"shear_stress_{huc12} = 0.01")
             region = tools.g_region(flags="p", format="json")
             tools.r_sim_sediment(
-                elevation=f"ned_{huc12}",
+                elevation="ned",
                 water_depth=f"depth_{huc12}",
                 dx=f"dx_{huc12}",
                 dy=f"dy_{huc12}",