Correct an indexing error in the computation of horizontal mixing for w.

Thanks to Matus Martini for finding this bug.
mperego · Sep 4, 2013 · 2686617 · 2686617
1 parent 1bf1146
commit 2686617
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/core_atmosphere/dynamics/mpas_atm_time_integration.F b/src/core_atmosphere/dynamics/mpas_atm_time_integration.F
@@ -2581,7 +2581,7 @@ subroutine atm_compute_dyn_tend(tend, s, diag, grid, rk_step, dt)
                   if (cell1 <= nCellsSolve .or. cell2 <= nCellsSolve) then
    
                      do k=2,nVertLevels
-                        theta_turb_flux = 0.25*(kdiff(k,cell1)+kdiff(k,cell2)+kdiff(k-1,cell1)+kdiff(k,cell2))  &
+                        theta_turb_flux = 0.25*(kdiff(k,cell1)+kdiff(k,cell2)+kdiff(k-1,cell1)+kdiff(k-1,cell2))  &
                                               *(w(k,cell2) - w(k,cell1))/dcEdge(iEdge)
                         theta_turb_flux = theta_turb_flux * meshScalingDel2(iEdge)
                         flux = 0.5*dvEdge (iEdge) * (rho_edge(k,iEdge)+rho_edge(k-1,iEdge)) * theta_turb_flux