I need to know the location of a given surface grid point (i,j) in an idealized simulation. The single-domain answer is easy: (i*dx, j*dy) [m] from the (0,0) reference point. However, this can become more complicated in the multi-domain case. Usually WPS filled the XLONG and XLAT values, but this stage is skipped in idealized cases - so XLONG and XLAT are empty. Is there a built-in function / utility that can fill in this role of WPS?