This output of mgdiff demonstrates the difference in the code involved by using ezscint instead of fscint. The ezscint code is on the left, fscint code is on the rigth.

In this example, we take a sample geopotential field at 500 mb that we interpolate from the current version of the GEM model to a Gaussian Grid. We then interpolate the data back on the GEM grid, and compute the differences using the 'statfld' subroutine.

For the sake of simplicity, the dimensions are defined in advance, instead of being dynamically allocated. Real working code would have to be much lengthier, taking care of changing dimensions and error codes.

Click on the image to view it in full size.

Go to the list of functions, or the general index