@@ -14,9 +14,12 @@ function spacing(lon,lat)
1414 dlon = zero (size (lon. val)[1 : 2 ])
1515 dlat = zero (size (lat. val)[1 : 2 ])
1616 @views dlon[2 : end - 1 ,:] = (lon. val[3 : end ,:,1 ] - lon. val[1 : end - 2 ,:,1 ])/ 2
17- @views dlon[1 ,:], dlon[end ,:] = dlon[2 ,:], dlon[end - 1 ,:]
18- @views dlat[:,2 : end - 1 ] = (lat. val[:,3 : end ,1 ] - lat. val[:,1 : end - 2 ,1 ])/ 2
19- @views dlat[:,1 ], dlat[:,end ] = dlat[:,2 ], dlat[:,end - 1 ]
17+ dlon[1 ,:] = dlon[2 ,:]
18+ dlon[end ,:] = dlon[end - 1 ,:]
19+ dlat[:,2 : end - 1 ] = (lat. val[:,3 : end ,1 ] - lat. val[:,1 : end - 2 ,1 ])/ 2
20+ dlat[:,1 ] = dlat[:,2 ]
21+ dlat[:,end ] = dlat[:,end - 1 ]
22+
2023 return dlon, dlat
2124end
2225
@@ -84,7 +87,7 @@ function waterflows(Topo::GeoData, flowdir_fn= WhereTheWaterFlows.d8dir_feature;
8487 nin = zeros (Int8, ni)
8588 c = zeros (Int64, ni)
8689
87- @views area[:,:,1 ], slen[:,:,1 ], dir[:,:,1 ], nout[:,:,1 ], nin[:,:,1 ], sinks, pits, c[:,:,1 ], bnds = waterflows (dem, cellarea, flowdir_fn;
90+ area[:,:,1 ], slen[:,:,1 ], dir[:,:,1 ], nout[:,:,1 ], nin[:,:,1 ], sinks, pits, c[:,:,1 ], bnds = waterflows (dem, cellarea, flowdir_fn;
8891 feedback_fn= feedback_fn, drain_pits= drain_pits, bnd_as_sink= bnd_as_sink)
8992
9093 catchment_large = prune_catchments (c, minsize; val= 0 )
0 commit comments