<p><b>ringler@lanl.gov</b> 2012-04-18 13:44:08 -0600 (Wed, 18 Apr 2012)</p><p><br>
temporarily estimate dvEdge as dcEdge/sqrt(3) as it is for perfect hexagons<br>
</p><hr noshade><pre><font color="gray">Modified: branches/omp_blocks/io/src/core_ocean/mpas_ocn_vel_hmix_del4.F
===================================================================
--- branches/omp_blocks/io/src/core_ocean/mpas_ocn_vel_hmix_del4.F        2012-04-17 20:58:17 UTC (rev 1791)
+++ branches/omp_blocks/io/src/core_ocean/mpas_ocn_vel_hmix_del4.F        2012-04-18 19:44:08 UTC (rev 1792)
@@ -175,7 +175,7 @@
          do k=1,maxLevelEdgeTop(iEdge)
             ! Compute </font>
<font color="black">abla^2 u = </font>
<font color="black">abla divergence + k \times </font>
<font color="red">abla vorticity
             delsq_u =          ( divergence(k,cell2)  - divergence(k,cell1) ) * invDcEdge  &amp;
-                -viscVortCoef *( vorticity(k,vertex2) - vorticity(k,vertex1)) * invDvEdge 
+                -viscVortCoef *( vorticity(k,vertex2) - vorticity(k,vertex1)) * invDcEdge * sqrt(3.0)   ! TDR
 
             ! vorticity using </font>
<font color="gray">abla^2 u
             r_tmp = dcEdge(iEdge) * delsq_u
@@ -203,7 +203,7 @@
 
          do k=1,maxLevelEdgeTop(iEdge)
             u_diffusion = (delsq_divergence(k,cell2) - delsq_divergence(k,cell1)) * invDcEdge  &amp;
-                -viscVortCoef * (delsq_vorticity(k,vertex2) - delsq_vorticity(k,vertex1) ) * invDvEdge
+                -viscVortCoef * (delsq_vorticity(k,vertex2) - delsq_vorticity(k,vertex1) ) * invDcEdge * sqrt(3.0) ! TDR
 
             tend(k,iEdge) = tend(k,iEdge) - edgeMask(k, iEdge) * u_diffusion * r_tmp
          end do

</font>
</pre>