Metric Size Field

[breuera]

The metric mesh refinement based on Python-expressions is broken. This is reproducible in one of the test cases (see below). Possible solution: The size field should be sampled in meters (not deg), then projected from UTM to Lon-lat, and finally the data set should be queried.

https://github.com/shingleproject/Shingle/blob/75c61afcd29a388b9621adced8f9216a6ee9aeee/test/UK_NorthSea_region_opendap_utm_metric/UK_NorthSea_region_opendap_utm_metric.brml#L102

• Observed output in NorthSea_metric.pos samples in lon-lat (ignoring the defined projection):

-13.9833 63.95 0
0.0332776867289 -0.0332714923016 1
599 539 1
27789.6
[...]

• Expected output in NorthSea_metric.pos is a 3D rectangular structured mesh in utm coordinates, which specifies the size field.