<p><b>qchen3@fsu.edu</b> 2011-12-22 16:17:11 -0700 (Thu, 22 Dec 2011)</p><p>BRANCH COMMIT<br>
<br>
A distinction is made between the z-level and the isopycnal vertical grids. GM for the isopycnal is implemented. <br>
There is nothing now for the z-level grid -- to come later.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/gmvar/src/core_ocean/mpas_ocn_tendency.F
===================================================================
--- branches/ocean_projects/gmvar/src/core_ocean/mpas_ocn_tendency.F        2011-12-22 17:26:41 UTC (rev 1276)
+++ branches/ocean_projects/gmvar/src/core_ocean/mpas_ocn_tendency.F        2011-12-22 23:17:11 UTC (rev 1277)
@@ -916,15 +916,20 @@
       ! Compute Bolus velocity as bolus == - /kappa (/grad h) / h
       !
       uBolus(:,:) = 0.0
-      call ocn_get_h_kappa(s, grid)
-      if (maxval(h_kappa(:,:)) &gt; epsilon(1d0)) then
-         do iEdge = 1,nEdges
-            cell1 = cellsOnEdge(1,iEdge)
-            cell2 = cellsOnEdge(2,iEdge)
-            do k=1,maxLevelEdgeTop(iEdge)
-               uBolus(k,iEdge) = -h_kappa(k,iEdge) * (h(k,cell2) - h(k,cell1))/dcEdge(iEdge)/h_edge(k,iEdge)
+      if (config_vert_grid_type .EQ. 'z-level') then
+         ! Nothing for now
+
+      else if (config_vert_grid_type .EQ. 'isopycnal') then
+         call ocn_get_h_kappa(s, grid)
+         if (maxval(h_kappa(:,:)) &gt; epsilon(1d0)) then
+            do iEdge = 1,nEdges
+               cell1 = cellsOnEdge(1,iEdge)
+               cell2 = cellsOnEdge(2,iEdge)
+               do k=1,maxLevelEdgeTop(iEdge)
+                  uBolus(k,iEdge) = -h_kappa(k,iEdge) * (h(k,cell2) - h(k,cell1))/dcEdge(iEdge)/h_edge(k,iEdge)
+               end do
             end do
-         end do
+         end if
       end if
 
       !

</font>
</pre>