<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>