`geogrid()` fails to find lat/lon values when subsetting Arrays in DAP4

[Mikejmnez]

geogrid(Var, ..., ) works as expected when Var is a Grid array. Var contains a "copy" of the coordinate values lat and lon. Grids are part of the DAP2 model but were dropped out of the DAP4 model.

When Var is an Array geogrid(Var, ...) fails to find the latitude and longitude arrays it needs.

Minimal example

The following errs (not a very helpful error)

https://oceandata.sci.gsfc.nasa.gov/opendap/MODISA/L3SMI/2002/0101/AQUA_MODIS.20020101_20021231.L3m.YR.RRS.Rrs_555.4km.nc.dap.nc4?dap4.function=geogrid(/Rrs_555,35,-121,32,-116)

@jgallagher59701 reached to this conclusion regarding in OPeNDAP's Discourse Thread