<p><b>duda</b> 2011-03-01 15:50:04 -0700 (Tue, 01 Mar 2011)</p><p>BRANCH COMMIT<br>
<br>
Bug fix to del4 diffusion for momentum in non-hydrostatic core:<br>
add parentheses so that the entire diffusion term is multiplied<br>
by rho_edge.<br>
<br>
M src/core_nhyd_atmos/module_time_integration.F<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_nhyd_atmos/module_time_integration.F
===================================================================
--- branches/atmos_physics/src/core_nhyd_atmos/module_time_integration.F        2011-03-01 22:46:59 UTC (rev 749)
+++ branches/atmos_physics/src/core_nhyd_atmos/module_time_integration.F        2011-03-01 22:50:04 UTC (rev 750)
@@ -2378,8 +2378,9 @@
! Compute diffusion, computed as </font>
<font color="black">abla divergence - k \times </font>
<font color="red">abla vorticity
! only valid for h_mom_eddy_visc4 == constant
!
- u_diffusion = rho_edge(k,iEdge) * ( delsq_divergence(k,cell2) - delsq_divergence(k,cell1) ) / dcEdge(iEdge) &
- -( delsq_vorticity(k,vertex2) - delsq_vorticity(k,vertex1) ) / dvEdge(iEdge)
+ u_diffusion = rho_edge(k,iEdge) * ( ( delsq_divergence(k,cell2) - delsq_divergence(k,cell1) ) / dcEdge(iEdge) &
+ -( delsq_vorticity(k,vertex2) - delsq_vorticity(k,vertex1) ) / dvEdge(iEdge) &
+ )
tend_u(k,iEdge) = tend_u(k,iEdge) - h_mom_eddy_visc4 * u_diffusion
end do
</font>
</pre>