<p><b>sprice@lanl.gov</b> 2012-01-10 16:37:56 -0700 (Tue, 10 Jan 2012)</p><p>Branch Commit (land ice): More corrections for subbing thck_layer for h.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/land_ice/mpas/src/core_land_ice/mpas_land_ice_time_integration.F
===================================================================
--- branches/land_ice/mpas/src/core_land_ice/mpas_land_ice_time_integration.F        2012-01-10 23:30:39 UTC (rev 1349)
+++ branches/land_ice/mpas/src/core_land_ice/mpas_land_ice_time_integration.F        2012-01-10 23:37:56 UTC (rev 1350)
@@ -271,7 +271,7 @@
       real (kind=RKIND), dimension(:), pointer :: h_s, fVertex, fEdge, dvEdge, dcEdge, areaCell, areaTriangle, &amp;
                                                   meshScalingDel2, meshScalingDel4
       real (kind=RKIND), dimension(:,:), pointer :: vh, weightsOnEdge, kiteAreasOnVertex, &amp;
-                                                    h_edge, thck_layer, unorm, tend_h, tend_u, &amp;
+                                                    h_edge, thck_layer, unorm, tend_thck_layer, &amp;
 !!!                                                    h_edge, h, u, v, tend_h, tend_u, &amp;
 !!!                                                    circulation, vorticity, ke, pv_edge, divergence,  &amp; 
                                                     h_vertex
@@ -338,12 +338,12 @@
       !
       ! Compute height tendency for each cell
       !
-      tend_h(:,:) = 0.0
+      tend_thck_layer(:,:) = 0.0
       do iEdge=1,nEdges
          cell1 = cellsOnEdge(1,iEdge)
          cell2 = cellsOnEdge(2,iEdge)
          do k=1,nVertLevels
-            flux = u(k,iEdge) * dvEdge(iEdge) * h_edge(k,iEdge)
+            flux = unorm(k,iEdge) * dvEdge(iEdge) * h_edge(k,iEdge)
             tend_thck_layer(k,cell1) = tend_thck_layer(k,cell1) - flux
             tend_thck_layer(k,cell2) = tend_thck_layer(k,cell2) + flux
          end do

</font>
</pre>