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