plotScatter and grouping renders error

[martin-gutwin]

c=Cube()
cube=getCubeData(c,variable=["net_ecosystem_exchange","gross_primary_productivity","terrestrial_ecosystem_respiration"],
longitude=(30.0,32.0),latitude=(50.0,52.0))
p=plotScatter(cube,alongaxis=TimeAxis,xaxis="net_ecosystem_exchange",yaxis="gross_primary_productivity",group="lat", lon=30)

This shows:

MethodError: no method matching grouped_xy(::LinearAlgebra.Transpose{Union{Missing, Float32},Array{Union{Missing, Float32},1}}, ::LinearAlgebra.Transpose{Union{Missing, Float32},Array{Union{Missing, Float32},1}})
Closest candidates are:
  grouped_xy(!Matched::AbstractArray{T,1} where T, ::AbstractArray) at /Net/Groups/BGI/scratch/mgutwin/julia_atacama_depots/packages/StatsPlots/vXig7/src/bar.jl:7
  grouped_xy(::AbstractArray) at /Net/Groups/BGI/scratch/mgutwin/julia_atacama_depots/packages/StatsPlots/vXig7/src/bar.jl:8

I guess the array needs to be transposed and then collected or copied? Example is adapted from the documentation.