<p><b>mperego@fsu.edu</b> 2013-06-07 12:12:16 -0600 (Fri, 07 Jun 2013)</p><p>Land Ice BRANCH: reduced lentghs of lines to comply with 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        2013-06-07 18:09:28 UTC (rev 2791)
+++ branches/land_ice_projects/implement_core/src/core_land_ice/mpas_land_ice_tendency.F        2013-06-07 18:12:16 UTC (rev 2792)
@@ -1130,7 +1130,9 @@
 
             ! Calculate reference height for this cell.
             if (numUpwindNeighbors(iCell) == 0.0_RKIND) then
-               ! It's ok if a cell has no thick ice neighbors now.  We will wait until we see if any of its partially filled neighbors spill into it before deciding what to do about it.  But set the Href to be defined so we don't enter the if-statement below checking for overfilled cells.
+               ! It's ok if a cell has no thick ice neighbors now. 
+               ! We will wait until we see if any of its partially filled neighbors spill into it before deciding what to do about it.
+               !  But set the Href to be defined so we don't enter the if-statement below checking for overfilled cells.
                Href(iCell) = 1.0e36_RKIND
             else
                Href(iCell) = sumUpwindThickness / numUpwindNeighbors(iCell)
@@ -1227,7 +1229,9 @@
                ! check which edges have velocity into this cell   ! + velocity goes from index 1 to 2 in the cellsOnEdge array.  
                ! Note that layerThicknessEdge is no longer valid for these edges because they previously had no flux
 
-               ! these edges should have 0 velocity except for edges adjacent to newly filled cells which will always have all layers fluxing into this cell (what we set in the previous step above), so we don't need to check the entire flux, we can just check a single layer's velocity.
+               ! these edges should have 0 velocity except for edges adjacent to newly filled cells
+               ! which will always have all layers fluxing into this cell (what we set in the previous step above),
+               ! so we don't need to check the entire flux, we can just check a single layer's velocity.
                if (   ( (normalVelocity(1,edgeIndex) &gt; 0.0_RKIND) .and. (cellsOnEdge(2, edgeIndex) == iCell))   &amp;
                  .or. ( (normalVelocity(1,edgeIndex) &lt; 0.0_RKIND) .and. (cellsOnEdge(1, edgeIndex) == iCell))   ) then
 
@@ -1266,7 +1270,9 @@
       end do    ! cell loop
 
 
-      ! One last step.  It is possible for a cell to become dynamically stranded during retreat so that it still has some ice in it but velocity is 0 on all edges feeding it.  In the CFBC subgrid parameterization, this cell is effectively in front of the calving front, and should therefore be calved.
+      ! One last step.  It is possible for a cell to become dynamically stranded during retreat so that it still has some ice in it
+      ! but velocity is 0 on all edges feeding it.  In the CFBC subgrid parameterization, this cell is effectively in front of the calving front,
+      ! and should therefore be calved.
       do iCell = 1, nCells 
          if ( numUpwindNeighbors(iCell) == 0.0_RKIND .and. MASK_IS_THIN_ICE(cellMask(iCell)) .and. MASK_IS_FLOATING(cellMask(iCell)) ) then 
                ! These cells have been dynamically stranded beyond the calving front on retreat, and they don't have any overfilled neighbors to save them - therefore calve them off.

</font>
</pre>