<p><b>mperego@fsu.edu</b> 2012-06-26 16:37:49 -0600 (Tue, 26 Jun 2012)</p><p>BRANCH COMMIT<br>
land_ice: pass the indexToCellxID array to lifev <br>
</p><hr noshade><pre><font color="gray">Modified: branches/land_ice_projects/implement_core/src/core_land_ice/mpas_land_ice_lifev.F
===================================================================
--- branches/land_ice_projects/implement_core/src/core_land_ice/mpas_land_ice_lifev.F        2012-06-26 19:52:13 UTC (rev 2003)
+++ branches/land_ice_projects/implement_core/src/core_land_ice/mpas_land_ice_lifev.F        2012-06-26 22:37:49 UTC (rev 2004)
@@ -219,6 +219,8 @@
 
       integer, dimension(:,:), pointer :: cellsOnEdge, cellsOnVertex, verticesOnCell, verticesOnEdge
 
+      integer, dimension(:), pointer :: indexToCellID
+
       real (kind=RKIND), dimension(:), pointer :: xCell, yCell, zCell, areaTriangle
 
       ! halo exchange arrays
@@ -249,6 +251,7 @@
       cellsOnVertex =&gt; block % mesh % cellsOnVertex % array
       verticesOnCell =&gt; block % mesh % verticesOnCell % array
       verticesOnEdge =&gt; block % mesh % verticesOnEdge % array
+      indexToCellID =&gt; block % mesh % indexToCellID % array
 
       xCell =&gt; block % mesh % xCell % array
       yCell =&gt; block % mesh % yCell % array
@@ -272,7 +275,7 @@
       !nVertLevels should be equal to nVertLevelsSolve (no split the domain in the vertical direction)
       call mpas_timer_start(&quot;velocity_solver_set_grid_data&quot;)
       call velocity_solver_set_grid_data(nCells, nEdges, nVertices, nVertLevels, nCellsSolve, nEdgesSolve, nVerticesSolve, &amp;
-        cellsOnEdge, cellsOnVertex, verticesOnCell, verticesOnEdge, xCell, yCell, zCell, areaTriangle, &amp;
+        cellsOnEdge, cellsOnVertex, verticesOnCell, verticesOnEdge, indexToCellID, xCell, yCell, zCell, areaTriangle, &amp;
         sendCellsArray, recvCellsArray, sendEdgesArray, recvEdgesArray, sendVerticesArray, recvVerticesArray)
       call mpas_timer_stop(&quot;velocity_solver_set_grid_data&quot;)
 #else
@@ -409,6 +412,7 @@
       case ('FO')
           call mpas_timer_start(&quot;velocity_solver_solve_FO&quot;)
           call velocity_solver_solve_FO(lowerSurface, thickness, beta, tracers(index_temperature,:,:), normalVelocity)
+!          call velocity_solver_estimate_SS_SMB(mesh % sfcMassBal % array)
           call mpas_timer_stop(&quot;velocity_solver_solve_FO&quot;)
           call mpas_timer_start(&quot;velocity_solver export&quot;)
 !          call velocity_solver_init_L1L2(LayerThicknessFractions)          

</font>
</pre>