<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 => block % mesh % cellsOnVertex % array
verticesOnCell => block % mesh % verticesOnCell % array
verticesOnEdge => block % mesh % verticesOnEdge % array
+ indexToCellID => block % mesh % indexToCellID % array
xCell => block % mesh % xCell % array
yCell => 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("velocity_solver_set_grid_data")
call velocity_solver_set_grid_data(nCells, nEdges, nVertices, nVertLevels, nCellsSolve, nEdgesSolve, nVerticesSolve, &
-        cellsOnEdge, cellsOnVertex, verticesOnCell, verticesOnEdge, xCell, yCell, zCell, areaTriangle, &
+        cellsOnEdge, cellsOnVertex, verticesOnCell, verticesOnEdge, indexToCellID, xCell, yCell, zCell, areaTriangle, &
        sendCellsArray, recvCellsArray, sendEdgesArray, recvEdgesArray, sendVerticesArray, recvVerticesArray)
call mpas_timer_stop("velocity_solver_set_grid_data")
#else
@@ -409,6 +412,7 @@
case ('FO')
call mpas_timer_start("velocity_solver_solve_FO")
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("velocity_solver_solve_FO")
call mpas_timer_start("velocity_solver export")
!         call velocity_solver_init_L1L2(LayerThicknessFractions)
</font>
</pre>