Calculating the grace circle distance from a point and determining maximum value at specific distances

[Carwyn Pelley]

From an old query that came in by a user, I thought I would add an illustrative example here (github gist):
https://gist.github.com/cpelley/6351252

This was not accepted as a suitable example for the iris gallery (not iris enough) so I thought I would post it here since I believe people would benefit from more examples working with the tools at their disposal.

This example demonstrates how we might use an external library (pyproj) to calculate the distance from our points in a cube to a specific location, using our cubes grid mapping. We then associate this information onto our original cube as an auxiliary coordinate, since this relates to the original phenomena.  Following this, a determination of maximum data value as a function of distance is plotted.

I realise I have not chosen example data for this.

[Phil Elson]

So in short it gives you an aux coordinate which represents each data point's distance (spherical distance) from a specified point?

Sounds cool. Thanks for sharing!

--
You received this message because you are subscribed to the Google Groups "Iris" group.
To unsubscribe from this group and stop receiving emails from it, send an email to scitools-iri...@googlegroups.com.
For more options, visit https://groups.google.com/groups/opt_out.

[Carwyn Pelley]

That's right Phil, I then simply put the data in distance bins showing how you can make use of such information.  I have noticed that a number of users define 'distance' by various means which more often than not does not consider what type of sphere the coordinate system describes, so their calculations would likely be radically different to others also making similar calculations.

Thanks for you interest.