@@ -14,9 +14,12 @@ function spacing(lon,lat)
14
14
dlon = zero (size (lon. val)[1 : 2 ])
15
15
dlat = zero (size (lat. val)[1 : 2 ])
16
16
@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
+
20
23
return dlon, dlat
21
24
end
22
25
@@ -84,7 +87,7 @@ function waterflows(Topo::GeoData, flowdir_fn= WhereTheWaterFlows.d8dir_feature;
84
87
nin = zeros (Int8, ni)
85
88
c = zeros (Int64, ni)
86
89
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;
88
91
feedback_fn= feedback_fn, drain_pits= drain_pits, bnd_as_sink= bnd_as_sink)
89
92
90
93
catchment_large = prune_catchments (c, minsize; val= 0 )
0 commit comments