There are 3 sets of IDL procedures, one for each method and each has two variants a) to process GHCN V3C and b) to process Hadcrut4.5 using Crutem4 station files, both using Hadsst3 Ocean temperature data. You can download each one as follows
<ol>
 	<li>2-D triangulation
<ol>
 	<li><a href="http://clivebest.com/data/t-study/V3-2D-tri-average.pro">V3-2D-tri-average.pro</a></li>
 	<li><a href="http://clivebest.com/data/t-study/Had4-2D-tri-average.pro">Had4-2D-tri-average.pro</a></li>
</ol>
</li>
 	<li>3-D triangulation
<ol>
 	<li><a href="http://clivebest.com/data/t-study/V3-tri-sphere-average.pro">V3-tri-sphere-average.pro</a></li>
 	<li><a href="http://clivebest.com/data/t-study/H4-tri-sphere-average.pro">H4-tri-sphere-average.pro</a></li>
</ol>
</li>
 	<li>Icosahedron Binning
<ol>
 	<li><a href="http://clivebest.com/data/t-study/V3-Isoc.pro">V3-Isoc.pro</a></li>
 	<li><a href="http://clivebest.com/data/t-study/H4-Isoc.pro">H4-Isoc.pro</a></li>
</ol>
</li>
</ol>
The Icosahedral binning uses <a href="http://clivebest.com/data/t-study/Zicos.zip">the following routines</a> which can be placed at the end of the above procedures or in your IDL library. In addition you may need this routine <a href="http://clivebest.com/data/t-study/nint.pro">nint.pro</a> which just calculates the nearest integer.
<p> It should not be too difficult to adapt them to other programming languages.