<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)  &amp;
-                                                 -( 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)  &amp;
+                                                    -( delsq_vorticity(k,vertex2) - delsq_vorticity(k,vertex1) ) / dvEdge(iEdge) &amp;
+                                                  )
 
                tend_u(k,iEdge) = tend_u(k,iEdge) - h_mom_eddy_visc4 * u_diffusion
             end do

</font>
</pre>