<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(:,:)) > 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(:,:)) > 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>