add nc_crop

src/NetCDFTools.jl

@@ -35,9 +35,7 @@ include("ncdim_def.jl")
 include("ncatt_put.jl")
 include("tools_rbase.jl")
 
-include("utilize/nc_aggregate.jl")
-include("utilize/nc_combine.jl")
-include("utilize/nc_subset.jl")
+include("utilize/utilize.jl")
 
 include("MFDataset/MFDataset.jl")
 # CMIP

src/utilize/nc_aggregate.jl

@@ -77,4 +77,4 @@ function nc_aggregate_dir(indir;
 end
 
 
-export nc_aggregate, nc_aggregate_dir
+

src/utilize/nc_crop.jl

@@ -0,0 +1,21 @@
+"""
+  nc_crop(f::String, range, fout::String; 
+    compress=false, overwrite=false, kw...)
+
+只有一个变量的情况。
+"""
+function nc_crop(f::String, range::Vector, fout::String; 
+  compress=false, overwrite=false, kw...)
+
+  lon, lat = st_dims(f)
+  ilon = findall(range[1] .<= lon .<= range[2]) |> zip
+  ilat = findall(range[3] .<= lat .<= range[4]) |> zip
+
+  A = nc_read(f, ind=(ilon, ilat, :));
+  var = nc_bands(f)[1]
+  _dims = ncvar_dim(f)
+  _dims["lon"] = _dims["lon"][ilon]
+  _dims["lat"] = _dims["lat"][ilat]
+
+  nc_write(fout, var, A, _dims; compress, overwrite, kw...)
+end

src/utilize/nc_subset.jl

@@ -114,4 +114,3 @@ function nc_subset(d::AbstractDataFrame, range;
 end
 
 
-export nc_subset

src/utilize/utilize.jl

@@ -0,0 +1,8 @@
+include("nc_aggregate.jl")
+include("nc_combine.jl")
+include("nc_subset.jl")
+include("nc_crop.jl")
+
+export nc_aggregate, nc_aggregate_dir
+export nc_combine
+export nc_subset, nc_crop