Fixes after JETLS warnings (#1967)

joa-quim · web-flow · commit 72a229cae62b · 2026-03-14T18:19:10.000Z
* Fixes after JETLS warnings

* Comment unintended include

* Try reset the GMT_SERVER

* Correct sever name

* Last try
diff --git a/.github/workflows/ci.yml b/.github/workflows/ci.yml
@@ -28,7 +28,7 @@ jobs:
           #- '1.7' # Leave this line unchanged. '1' will automatically expand to the latest stable 1.x release of Julia.
           - '~1.10.0-0'
           #- '~1.11.0-0'
-          - '~1.12.0-0'
+          #- '~1.12.0-0'
           #- '~1.13.0-0'
           - 'nightly'
         os:
diff --git a/ci/config-gmt-unix.sh b/ci/config-gmt-unix.sh
@@ -11,6 +11,7 @@ set (DCW_ROOT "$ENV{COASTLINEDIR}/dcw")
 
 set (GMT_USE_THREADS TRUE)
 set (GMT_ENABLE_OPENMP TRUE)
+set (GMT_DATA_SERVER https://www.star.nesdis.noaa.gov/data/socd3/lsa/gmtdata/)
 
 # recommended even for release build
 set (CMAKE_C_FLAGS "-Wall -Wdeclaration-after-statement ${CMAKE_C_FLAGS}")
diff --git a/src/GMT.jl b/src/GMT.jl
@@ -407,6 +407,7 @@ using .Laszip
 	grd2cpt(Glix)
 	grd2xyz(Glix)
 	viz(Glix, show=false)
+	viz(I, show=false)
 	histogram(randn(100), T=0.1, center=true, Z=:counts)
 	histogram(I)
 	grdlandmask(R="-10/4/37/45", res=:c, inc=0.1)
diff --git a/src/PrettyTables.jl b/src/PrettyTables.jl
@@ -698,7 +698,7 @@ function _print_table_with_text_back_end(
 
 	# Get size information from the processed table.
 	num_rows, num_columns = _data_size(ptable)
-	num_header_rows, ~ = _header_size(ptable)
+	num_header_rows, _ = _header_size(ptable)
 
 	# == Input Variables Verification and Initial Setup ====================================
 
@@ -2526,7 +2526,7 @@ function _print_custom_text_cell!(display::Display, cell_data::CustomTextCell, c
 	_p!(display, " ", false, 1)
 
 	# Compute the new string given the display size.
-	str, suffix, ~ = _fit_string_in_display(display, cell_processed_str, false, cell_printable_textwidth)
+	str, suffix, _ = _fit_string_in_display(display, cell_processed_str, false, cell_printable_textwidth)
 
 	new_lstr = textwidth(str)
 
diff --git a/src/custom_symb_funs.jl b/src/custom_symb_funs.jl
@@ -17,7 +17,7 @@ function flower_minho(; name::String="flower_minho", width=14, fmt="", format=""
 		figsize=(width,0),             # Fig size. Second arg = 0 means compute the height keeping aspect ratio
 		frame=:none,                   # Do not plot the frame
 		p=opt_p)
-	plot!(x, y, fill=(pattern=TESTSDIR * "tiling2.jpg", dpi=dpi), name = name * _fmt)
+	plot!(x, y, fill=(pattern=tests("tiling2.jpg"), dpi=dpi), name = name * _fmt)
 	(show == 1) && showfig()
 end
 
diff --git a/src/gdal_utils.jl b/src/gdal_utils.jl
@@ -965,7 +965,7 @@ Write a MxNxP `cube` object to disk as a multilayered file.
 function gdalwrite(fname::AbstractString, data, optsP=String[]; opts=String[], kw...)
 	(fname == "") && error("Output file name is missing.")
 	(isempty(optsP) && !isempty(opts)) && (optsP = opts)			# Accept either Positional or KW argument
-	isa(optsP, AbstractString) && (optsP = gdal_opts2vec(optsP))	# Guarantied to return a Vector{String}
+	isa(optsP, AbstractString) && (optsP = Gdal.gdal_opts2vec(optsP))	# Guarantied to return a Vector{String}
 	ds = Gdal.get_gdaldataset(data, optsP)[1]
 	if (Gdal.GDALGetRasterCount(ds.ptr) >= 1)
 		gdaltranslate(ds, optsP; dest=fname, kw...)
diff --git a/src/gmt_main.jl b/src/gmt_main.jl
@@ -614,7 +614,7 @@ function get_PS(API::Ptr{Nothing}, object::Ptr{Nothing})::GMTps
 end
 
 # ---------------------------------------------------------------------------------------------------
-function get_dataset(API::Ptr{Nothing}, object::Ptr{Nothing})::GDtype
+function get_dataset(object::Ptr{Nothing})::GDtype
 # Given a GMT DATASET D, build an array of segment structure and assign values.
 # Each segment will have 6 items:
 # header:	Text string with the segment header (could be empty)
@@ -799,7 +799,7 @@ function GMTJL_Get_Object(API::Ptr{Nothing}, X::GMT_RESOURCE)
 	elseif (X.family == GMT_IS_CUBE)       	# A GMT cube; make it the pos'th output item
 		ptr = get_grid(API, X.object, true)
 	elseif (X.family == GMT_IS_DATASET)		# A GMT table; make it a matrix and the pos'th output item
-		ptr = get_dataset(API, X.object)
+		ptr = get_dataset(X.object)
 	elseif (X.family == GMT_IS_PALETTE)		# A GMT CPT; make it a colormap and the pos'th output item
 		ptr = get_palette(API, X.object)
 	elseif (X.family == GMT_IS_IMAGE)		# A GMT Image; make it the pos'th output item
diff --git a/src/gmtreadwrite.jl b/src/gmtreadwrite.jl
@@ -91,7 +91,7 @@ function gmtread(fname::String; kwargs...)
 end
 function gmtread(fname::String, d::Dict{Symbol, Any})
 
-	endswith(fname, ".xlsx") && return read_xls(_fname; kwargs...)		# Excel extension
+	endswith(fname, ".xlsx") && return read_xls(fname; d...)		# Excel extension
 
 	#fname::String = _fname					# Because args signatures seam to worth shit in body.
 	#d = init_module(false, kwargs...)[1]	# Also checks if the user wants ONLY the HELP mode
diff --git a/src/gmtspatial.jl b/src/gmtspatial.jl
@@ -54,7 +54,7 @@ gmtspatial(arg1, arg2; kw...) = gmtspatial_helper("", arg1, arg2; kw...)
 
 # ---------------------------------------------------------------------------------------------------
 function gmtspatial_helper(cmd0::String, arg1, arg2; kw...)
-	(cmd0 == "" && arg1 === nothing && arg2 === nothing && length(kwargs) == 0) && return gmt("gmtspatial")
+	(cmd0 == "" && arg1 === nothing && arg2 === nothing && length(kw) == 0) && return gmt("gmtspatial")
 	d = init_module(false, kw...)[1]			# Also checks if the user wants ONLY the HELP mode
 	isa(arg1, Matrix) && (arg1 = mat2ds(arg1))
 	invokelatest(_gmtspatial_helper, cmd0, arg1, arg2, d)
diff --git a/src/imshow.jl b/src/imshow.jl
@@ -34,14 +34,17 @@ julia> viz(G, zsize=6, facades=(TESTSDIR*"assets/cenora_base.jpg", TESTSDIR*"bun
 See also: [`grdimage`](@ref), [`grdview`](@ref)
 """
 function imshow(arg1, x::AbstractVector{Float64}=Float64[], y::AbstractVector{Float64}=Float64[]; kw...)
+	_imshow(arg1, x, y, KW(kw))
+end
+function _imshow(@nospecialize(arg1), x::AbstractVector{Float64}, y::AbstractVector{Float64}, d::Dict{Symbol,Any})
 	# Take a 2D array of floats and turn it into a GMTgrid or if input is a string assume it's a file name
 	# In this later case try to figure if it's a grid or an image and act accordingly.
 
-	see = ((val = find_in_kwargs(kw, [:show])[2]) === nothing) ? true : (val != 0)	# No explicit 'show' keyword means show=true
+	see = (!haskey(d, :show)) ? true : (d[:show] != 0)	# No explicit 'show' keyword means show=true
 
-	function isplot3(kw)
+	function _isplot3(d::Dict)
 		_call_plot3 = false
-		opt_p = find_in_kwargs(kw, [:p :view :perspective])[1]
+		opt_p = find_in_dict(d, [:p :view :perspective], false)[1]
 		(isa(opt_p, String) && contains(opt_p, '/')) && (_call_plot3 = true)
 		((isa(opt_p, Tuple) || isa(opt_p, VMr)) && length(opt_p) > 1) && (_call_plot3 = true)
 		return _call_plot3
@@ -60,10 +63,9 @@ function imshow(arg1, x::AbstractVector{Float64}=Float64[], y::AbstractVector{Fl
 			ext = lowercase(ext)
 			(ext == ".jpg" || ext == ".tif" || ext == ".tiff" || ext == ".png" || ext == ".bmp" || ext == ".gif") && (is_image = true)
 			!is_image && (Gg = arg1)
-			#snif_GI_set_CTRLlimits(arg1)			# Set CTRL.limits to be eventually used by J=:guess
 		end
 	elseif (isa(arg1, Array{UInt8}) || isa(arg1, Array{UInt16,3}) || isa(arg1, Array{Bool}) || isa(arg1, BitMatrix))
-		Gi = isa(arg1, BitMatrix) ? mat2img(collect(arg1); kw...) : mat2img(arg1; kw...)
+		Gi = isa(arg1, BitMatrix) ? mat2img(collect(arg1); d...) : mat2img(arg1; d...)
 		call_img = true
 	elseif (isa(arg1, GDtype) || isa(arg1, GMTfv) || isa(arg1, Vector{GMTfv}) || (isa(arg1, AbstractVecOrMat{<:Real}) && size(arg1,2) <= 4) || (isa(arg1, Gdal.AbstractDataset) || isa(arg1, Gdal.AbstractGeometry)))
 		(isa(arg1, Gdal.AbstractDataset) || isa(arg1, Gdal.AbstractGeometry)) && (arg1 = gd2gmt(arg1))
@@ -82,13 +84,15 @@ function imshow(arg1, x::AbstractVector{Float64}=Float64[], y::AbstractVector{Fl
 		call_plot3 = ((isa(arg1, GMTdataset) && arg1.geom == wkbPolygonZM) ||
 		              (isa(arg1, Vector{<:GMTdataset}) && arg1[1].geom == wkbPolygonZM) ||
 					  isFV3D) ? true : false			# Should evolve into a fun that detects the several plot3d cases.
-		!call_plot3 && (call_plot3 = isplot3(kw))
-		return (call_plot3) ? plot3d(arg1; show=see, kw...) : plot(arg1; show=see, kw...)
+		!call_plot3 && (call_plot3 = _isplot3(d))
+		d[:show] = see
+		return (call_plot3) ? plot3d(arg1; d...) : plot(arg1; d...)
 	elseif (isa(arg1, GMTcpt))
-		return (find_in_kwargs(kw, [:D :pos :position])[1] === nothing) ?
-		        psscale(arg1; show=true, D="x0/0+w7+h", kw...) : psscale(arg1; show=true, kw...)
+		return (find_in_dict(d, [:D :pos :position], false)[1] === nothing) ?
+		        psscale(arg1; show=true, D="x0/0+w7+h", d...) : psscale(arg1; show=true, d...)
 	elseif (isdataframe(arg1) || isODE(arg1))
-		return isplot3(kw) ? plot3(arg1; show=see, kw...) :  plot(arg1; show=see, kw...)
+		d[:show] = see
+		return _isplot3(d) ? plot3(arg1; d...) : plot(arg1; d...)
 	else
 		Gg = mat2grid(arg1, x, y, reg=1)					# For displaying, pixel registration is more appropriate
 		call_grd = true
@@ -98,17 +102,18 @@ function imshow(arg1, x::AbstractVector{Float64}=Float64[], y::AbstractVector{Fl
 	call_grd && snif_GI_set_CTRLlimits(Gg)
 
 	if (is_image)
-		grdimage(arg1; show=see, kw...)
+		grdimage(arg1; show=see, d...)
 	else
-		if (isa(Gg, String))		# Guess also if call grdview or grdimage 
-			if (get(kw, :JZ, 0) != 0 || get(kw, :Jz, 0) != 0 || get(kw, :zscale, 0) != 0 || get(kw, :zsize, 0) != 0)
-				(get(kw, :Q, "") == "" && get(kw, :surf, "") == "" && get(kw, :surftype, "") == "") && (kw = (kw..., Q="s"))
-				grdview(arg1; show=see, kw...)				# String when fname is @xxxx
+		if (isa(Gg, String))		# Guess also if call grdview or grdimage
+			d[:show] = see
+			if (get(d, :JZ, 0) != 0 || get(d, :Jz, 0) != 0 || get(d, :zscale, 0) != 0 || get(d, :zsize, 0) != 0)
+				(get(d, :Q, "") == "" && get(d, :surf, "") == "" && get(d, :surftype, "") == "") && (d[:Q] = "s")
+				grdview(arg1; d...)				# String when fname is @xxxx
 			else
-				grdimage(arg1; show=see, kw...)
+				grdimage(arg1; d...)
 			end
 		else
-			imshow(call_img ? Gi : Gg; kw...)				# Call the specialized method
+			_imshow_GI(call_img ? Gi : Gg, d)				# Call the specialized method directly
 		end
 	end
 end
@@ -122,12 +127,14 @@ end
 # - `T, no_interp, tiles`: -T option for grdview
 # - `facades, cubeplot`: Call cubeplot.
 function imshow(arg1::GItype; kw...)
-	# Here the default is to show, but if a 'show' was used let it rule
-	d = KW(kw)
+	_imshow_GI(arg1, KW(kw))
+end
+function _imshow_GI(arg1::GItype, d::Dict{Symbol,Any})
 	see::Bool = (!haskey(d, :show)) ? true : (d[:show] != 0)	# No explicit 'show' keyword means show=true
 
-	if (isa(arg1, GMTimage) && (size(arg1, 3) <= 3 || arg1.layout[4] == 'A'))	# Rest of the work is done in grdiamge
-		return grdimage("", arg1; show=see, kw...)
+	if (isa(arg1, GMTimage) && (size(arg1, 3) <= 3 || arg1.layout[4] == 'A'))	# Rest of the work is done in grdimage
+		d[:show] = see
+		return grdimage_helper(wrapGrids("", arg1), nothing, nothing, true, false, d)
 	end
 
 	if ((cont_opts = find_in_dict(d, [:contour])[1]) !== nothing)
@@ -138,7 +145,7 @@ function imshow(arg1::GItype; kw...)
 	have_tilles::Bool = ((til = find_in_dict(d, [:T :no_interp :tiles])[1]) !== nothing)
 	(!have_tilles && opt_p != "" && !contains(opt_p, '/')) && (flat = true)		# If only 'azimuth' and no 'elev'
 	flat::Bool = (find_in_dict(d, [:flat])[1] !== nothing)	# If true, force the use of grdimage
-	docube = is_in_kwargs(kw, [:facades :cubeplot])
+	docube = is_in_dict(d, [:facades :cubeplot]) !== nothing
 	(flat && docube) && (flat = false)
 
 	if (!docube && (flat || (opt_p == "" && !have_tilles)))
@@ -232,15 +239,16 @@ end
 
 # Simple method to show CPTs. (May grow)
 # Ex: imshow(C, xlabel="Bla", ylabel="Blu"), or imshow(:gray)
-imshow(arg1::Symbol; horizontal::Bool=false, kw...) = imshow(makecpt(arg1); horizontal=horizontal, kw...) 
+imshow(arg1::Symbol; horizontal::Bool=false, kw...) = imshow(makecpt(arg1); horizontal=horizontal, kw...)
 function imshow(arg1::GMTcpt; horizontal::Bool=false, kw...)
-	see::Bool = (!haskey(kw, :show)) ? true : kw[:show]		# No explicit 'show' keyword means show=true
-	horizontal ? psscale(arg1; D="x8c/1c+w12c/0.5c+jTC+h", show=see, kw...) : psscale(arg1; J="X15/0", D="x8c/1c+w12c/0.5c+jBC", show=see, kw...)
+	d = KW(kw)
+	see::Bool = (!haskey(d, :show)) ? true : (d[:show] != 0)
+	horizontal ? psscale(arg1; D="x8c/1c+w12c/0.5c+jTC+h", show=see, d...) : psscale(arg1; J="X15/0", D="x8c/1c+w12c/0.5c+jBC", show=see, d...)
 end
 
 function imshow(arg1::Gdal.AbstractDataset; kw...)
 	(Gdal.GDALGetRasterCount(arg1.ptr) == 0) && return plot(gd2gmt(arg1), show=1)
-	imshow(gd2gmt(arg1); kw...)
+	_imshow_GI(gd2gmt(arg1), KW(kw))
 end
 
 imshow(x::AbstractVector{Float64}, y::AbstractVector{Float64}, f::Function; kw...) = imshow(mat2grid(f, x, y); kw...)
diff --git a/src/laszip/laszip_dll.jl b/src/laszip/laszip_dll.jl
@@ -82,11 +82,11 @@ end
 function laszip_add_vlr(pointer::Ptr{Cvoid}, vlr::Ptr{laszip_vlr})
     ccall((:laszip_add_vlr, laszip), Cint, (Ptr{Cvoid}, Ptr{laszip_vlr}), pointer, vlr)
 end
+=#
 
 function laszip_seek_point(pointer::Ptr{Cvoid}, index::Clonglong)
     ccall((:laszip_seek_point, laszip), Cint, (Ptr{Cvoid}, Clonglong), pointer, index)
 end
-=#
 
 function laszip_open_writer(pointer::Ptr{Cvoid}, fname, compress::Int)
     ccall((:laszip_open_writer, laszip), Cint, (Ptr{Cvoid}, Ptr{UInt8}, Cint), pointer, fname, compress)
diff --git a/src/libgmt.jl b/src/libgmt.jl
@@ -442,8 +442,7 @@ function gmt_centroid_area(API::Ptr{Cvoid}, D, geo::Int=0; ca::Int=3)::Matrix{Fl
 		end
 	else
 		ref1, ref2, nr = pointercols(D)
-		t = gmt_centroid_area(API, ref1, ref2, geo, nr)[inds]
-		mat = bla(mat, t, ca, 1)
+		t = gmt_centroid_area(API, ref1, ref2, geo, nr)
 		if (ca == 3) mat[1], mat[2], mat[3] = t[1], t[2], t[3] elseif (ca == 2) mat[1], mat[2] = t[1], t[2] else mat[1] = t[3] end
 	end
 	return mat
diff --git a/src/maregrams.jl b/src/maregrams.jl
@@ -42,7 +42,7 @@ viz(D, title="Tide Gauge at Lagos (Portugal)")
 """
 function maregrams(; list=false, code="", name="", days=2, starttime::String="", printurl::Bool=false)::GMTdataset{Float64, 2}
 	if (list != 0)			# Either list all stations or search those that contains a "list" string
-		Dl = gmtread(TESTSDIR * "/assets/maregs_online.csv")
+		Dl = gmtread(tests("maregs_online.csv"))
 		list == 1 && return Dl
 		inds = findall(contains.(Dl.text,list))
 		isempty(inds) && (println("'$list' was not found in stations list."); return GMTdataset())
@@ -97,7 +97,7 @@ function maregrams(; list=false, code="", name="", days=2, starttime::String="",
 	D.colnames = ["time", var_in_mareg[1]]
 	D.attrib["Timecol"] = "1"
 	if (have_code)			# When we have the code in our small data file
-		ind = (code != "") ? findfirst(code .== d[:code]) : findfirst(stname .== d[:name])
+		ind = (code != "") ? findfirst(code .== d[:code]) : findfirst(name .== d[:name])
 		D.attrib["Country"] = d[:country][ind]
 		D.attrib["ST_name"] = d[:name][ind]
 		D.attrib["ST_code"] = d[:code][ind]
@@ -119,7 +119,7 @@ function maregrams(x::Real, y::Real; days=2, starttime::String="", printurl::Boo
 end
 
 # ---------------------------------------------------------------------------------------------------
-function read_maregrams(fname=TESTSDIR * "/assets/maregs_online.csv")
+function read_maregrams(fname=tests("maregs_online.csv"))
 	mat = Matrix{Float64}(undef, 378, 2)
 	names = Vector{String}(undef, 378)
 	codes = Vector{String}(undef, 378)
diff --git a/src/marker_name.jl b/src/marker_name.jl
@@ -180,6 +180,5 @@ function seek_custom_symb(marca::AbstractString, with_k::Bool=false)::String
 
 	_siz  = split(marca, '/')[2]			# The custom symbol size
 	_marca = (with_k ? "k" : "")  * fullname * "/" * _siz
-	(GMTver <= v"6.4" && (length(_marca) - length(_siz) -2) > 62) && warn("Due to a GMT <= 6.4 limitation the length of full (name+path) custom symbol name cannot be longer than 62 bytes.")
 	return _marca
 end
diff --git a/src/plot.jl b/src/plot.jl
@@ -1353,7 +1353,7 @@ function radar_helper(cmd0::String, arg1; first::Bool=true, axeslimts=Float64[],
 	(cmd0 != "") && (arg1 = read_data(d, cmd0, "", arg1, " ", false, true)[2])	# Make sure we have the data here
 	if (isa(arg1, GMTdataset))                data::Matrix{Float64} = arg1.data
 	elseif (isa(arg1, Vector{<:GMTdataset}))  data = ds2ds(arg1).data
-	elseif (isa(arg1, Vector{<:Real}))        data = reshape(arg, 1, length(arg1))
+	elseif (isa(arg1, Vector{<:Real}))        data = reshape(arg1, 1, length(arg1))
 	else                                      data = arg1
 	end
 
@@ -2032,7 +2032,7 @@ function piechart(x::VecOrMat; first::Bool=true, kw...)
 	do_explode = false
 	non_exploded = ones(Bool, length(X))		# Default to no explosion
 	if ((val = find_in_dict(d, [:explode])[1]) !== nothing)
-		explode::Vector{Int} = isa(val, Bool) ? find(val) : [(val)...]
+		explode::Vector{Int} = isa(val, Bool) ? findall(val) : [(val)...]
 		maximum(explode) > length(X) && error("The explode option must be a vector of Ints with no values larger than the size of data")
 		data_explode = data[explode, :]			# Get the exploded data
 		non_exploded[explode] .= false
diff --git a/src/psxy.jl b/src/psxy.jl
@@ -45,8 +45,8 @@ function _common_plot_xyz(w::wrapDatasets, caller::String, O::Bool, K::Bool, is3
 
 	parse_paper(d)				# See if user asked to temporarily pass into paper mode coordinates
 
+	opt_B::String, opt_R::String = "", ""
 	if (is_ternary)
-		opt_B::String = ""
 		if (haskey(d, :B))								# Not necessarely the case when ternary!
 			cmd, opt_B = string(cmd, d[:B]), d[:B]		# B option was parsed in plot/ternary
 			delete!(d, :B)
diff --git a/src/statplots.jl b/src/statplots.jl
@@ -889,10 +889,10 @@ function colorize_VecVecVec(@nospecialize(D), N_grp::Int, N_in_each_grp::Vector{
 			D = ds2ds(ds2ds(D)); set_dsBB!(D)	# Crazzy op and wasteful but thse Ds are small
 		end
 		for m = 1:N_grp
-			colorize_candles_violins(D, N_grp, vv[m], !ccolor ? m : 0, custom_colors)	# Assign default colors
+			colorize_candles_violins(D, N_grp, vv[m], m, custom_colors)	# Assign default colors
 		end
 	else
-		colorize_candles_violins(D, N_grp, 1:N_grp, !ccolor ? m : 0, custom_colors)
+		colorize_candles_violins(D, N_grp, 1:N_grp, 0, custom_colors)
 	end
 	return D
 end
diff --git a/src/trend1d.jl b/src/trend1d.jl
@@ -77,7 +77,7 @@ function trend1d(cmd0::String, arg1, d::Dict{Symbol, Any})
 
 	if (opt_F != "")		# Extract column names from opt_F
 		if (opt_F[1] == 'p' || opt_F[1] == 'P' || opt_F[1] == 'c')
-			colnames = ["a$i" for i=0:size(D,2)-1]
+			colnames = ["a$i" for i=0:size(arg1, 2)-1]
 		else
 			cnames = (isa(arg1, GDtype)) ? (isa(arg1, GMTdataset) ? arg1.colnames : arg1[1].colnames) : String[]
 			s = collect(opt_F)
diff --git a/src/utils.jl b/src/utils.jl
@@ -1277,7 +1277,7 @@ end
 Return the width and height of the grid or image from the file `fname`.
 """
 function getsize(fname::String)
-	!isfile(fname) && error("File $GI does not exist")
+	!isfile(fname) && error("File $fname does not exist")
 	ds = Gdal.read(fname)
 	return Gdal.width(ds), Gdal.height(ds)
 end
@@ -1806,3 +1806,4 @@ end
 # ---------------------------------------------------------------------------------------------------
 include("makeDCWs.jl")
 include("getdcw.jl")
+#include("tttAPI.jl")
diff --git a/src/utils_types.jl b/src/utils_types.jl
@@ -1948,7 +1948,7 @@ function stackgrids(names::Vector{String}, v=nothing; zcoord=nothing, zdim_name:
 			v = Dates.datetime2unix.(v);			z_unit = "Seconds since 1970-01-01T00:00:00"
 		elseif (_z_unit == "rata")
 			v = Dates.datetime2rata.(v);			z_unit = "Days since 0000-12-31T00:00:00"
-		elseif (_z_unit == "date" || _zunit == "datetime")			# Crashes Mirone
+		elseif (_z_unit == "date" || _z_unit == "datetime")			# Crashes Mirone
 			v = string.(v);							z_unit = "ISO Date Time"
 		end
 	end
