WCPS - add non-mentioned axes to domain intervals for scale() implicitly

[Bang Pham Huu]

For example, if a coverage is 3D (time, lat and long axes), and one only specifies scale values for lat and long axes.

for $c in (test_irr_cube_2)

return encode(
   scale(
         $c[ansi:"CRS:1"(0:3)] ,  
         {E:"CRS:1"(0:20),  N:"CRS:1"(0:30) }
        ), "json" )

Then, the grid domains of the ansi (time axis) should be added implicitly to the domain intervals for the scale().

[Bang Pham Huu]

NOTE: currently, there is a feature to scale X and Y axes implicitly: ​https://doc.rasdaman.org/05_geo-services-guide.html#auto-ratio-for-scaling-x-or-y-axis-in-wcps

This ticket should not break this behavior for X and Y axes.