<p><b>xylar@lanl.gov</b> 2010-04-22 16:56:13 -0600 (Thu, 22 Apr 2010)</p><p>BRANCH COMMIT<br>
<br>
Committed without a comment by mistake last time. The changes then were:<br>
* incorporating the changes to make g95-serial work, including replacing an MPI_abort with dmpar_abort<br>
* adding 4 geometric fields to the registry that are useful for interpolation<br>
* adding significant changes to the interpolation code that uses radial basis functions to compute coefficients for both 2D and 3D vectors and scalars<br>
<br>
A few added comments on the ibInterp_initialize subroutine in this commit<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/IBinterp/mpas/src/operators/module_IB_interpolation.F
===================================================================
--- branches/ocean_projects/IBinterp/mpas/src/operators/module_IB_interpolation.F        2010-04-22 22:29:47 UTC (rev 204)
+++ branches/ocean_projects/IBinterp/mpas/src/operators/module_IB_interpolation.F        2010-04-22 22:56:13 UTC (rev 205)
@@ -31,11 +31,18 @@
subroutine ibInterp_initialize(grid)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
- ! Purpose:
+ ! Purpose: compute geometric fields that will be potentially useful for calling
+ ! the interpolation routines
!
- ! Input:
+ ! Input: the grid
!
- ! Output:
+ ! Output:
+ ! radialUnitVectors - the unit vector in the radial direction at each cell center
+ ! edgeNormalVectors - the unit vector at the center of each edge tangent to the sphere
+ ! edgeLocations - the location of each edge center in cartesian coordinates
+ ! cellTangentPlane - 2 orthoginal unit vectors in the tangent plane of each cell
+ ! The first unit vector is chosen to point toward the center of the first
+ ! edge on the cell.
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
</font>
</pre>