add log10(largest_area)

JuliaGeodynamics · Jul 5, 2024 · 23d9d4a · 23d9d4a
1 parent 6e1b1bb
commit 23d9d4a
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/src/WaterFlow.jl b/src/WaterFlow.jl
@@ -98,7 +98,8 @@ function waterflows(Topo::GeoData, flowdir_fn= WhereTheWaterFlows.d8dir_feature;
     largest_area[.!largest_catchment] .= NaN
 
     log10_area = log10.(area)
-
-    Topo_water =  addfield(Topo,(;area, slen, dir, nout, nin, c, cellarea_m2, catchment_large, log10_area, largest_catchment, largest_area))                        
+    log10_largest_area = log10.(largest_area)
+
+    Topo_water =  addfield(Topo,(;area, slen, dir, nout, nin, c, cellarea_m2, catchment_large, log10_area, largest_catchment, largest_area, log10_largest_area))                        
     return Topo_water, sinks, pits, bnds 
 end