<p><b>mperego@fsu.edu</b> 2012-06-26 17:01:46 -0600 (Tue, 26 Jun 2012)</p><p>BRANCH COMMIT<br>
land_ice: break long lines into small pieces to please pgi compiler<br>
</p><hr noshade><pre><font color="gray">Modified: branches/land_ice_projects/implement_core/src/core_land_ice/mpas_land_ice_tendency.F
===================================================================
--- branches/land_ice_projects/implement_core/src/core_land_ice/mpas_land_ice_tendency.F        2012-06-26 22:37:49 UTC (rev 2004)
+++ branches/land_ice_projects/implement_core/src/core_land_ice/mpas_land_ice_tendency.F        2012-06-26 23:01:46 UTC (rev 2005)
@@ -466,7 +466,9 @@
             do k=1, nVertLevels
                ! Calculate h on edges using first order
                VelSign = sign(1.0, normalVelocity(k, iEdge))
-               layerThicknessEdge(k,iEdge) = max(VelSign * layerThickness(k, cell1), VelSign * (-1.0) * layerThickness(k, cell2))  ! + velocity goes from index 1 to 2 in the cellsOnEdge array.  !  Doug does the calculation as: h_edge = max(VelSign, 0.0) * h1 - min(VelSign, 0.0) * h2
+               layerThicknessEdge(k,iEdge) = max(VelSign * layerThickness(k, cell1), VelSign * (-1.0) * layerThickness(k, cell2)) 
+                ! + velocity goes from index 1 to 2 in the cellsOnEdge array.  
+                !  Doug does the calculation as: h_edge = max(VelSign, 0.0) * h1 - min(VelSign, 0.0) * h2
             end do
             ! thickness_edge is not currently in registry and not currenly needed.  If it is, uncomment the next line
             !h_edge = max(thickness(cell1), thickness(cell2))
@@ -688,7 +690,9 @@
       do iEdge=1,nEdges
          cell1 = cellsOnEdge(1,iEdge)
          cell2 = cellsOnEdge(2,iEdge)
-         if ( (layerThickness(1,cell1) .gt. 0.0) .or. (layerThickness(1,cell2) .gt. 0.0) ) then  ! Don't calculate for non-ice cells - would result in divide by 0 (just need to check one layer)  \todo: this should use the edgeMask and use the dynamic thickness limit
+         if ( (layerThickness(1,cell1) .gt. 0.0) .or. (layerThickness(1,cell2) .gt. 0.0) ) then  
+            ! Don't calculate for non-ice cells - would result in divide by 0 (just need to check one layer)
+            ! \todo: this should use the edgeMask and use the dynamic thickness limit
 
             VelSignSum = 0.0
             do k=1, nVertLevels

</font>
</pre>