anchor.diff:                          end if 
                       end do
                   end do
  +                !$OMP END DO
   
               end if
  +
  +            !$OMP END PARALLEL
   
               ! In the case we use a "not-truely-linear" variant of the linear solver, the energy
               ! calculations do not really make sense
  
positive.code:                         end if
                      end do
                  end do
              else if(.not. domain%is_staggered_grid) then
                  !
                  ! Compute energy statistics using non-staggered grid
                  !
  
                  ! Integrate over the model interior
                  do j = 1 + ecw, (domain%nx(2) - ecw)
                      do i = (1+ecw), (domain%nx(1)-ecw)
                          depth = domain%U(i,j,STG) - domain%U(i,j,ELV)
  
negative.code:                 !
  
                  ! Integrate over the model interior
                  !$OMP DO
                  do j = 1 + ecw, (domain%nx(2) - ecw)
                      do i = (1+ecw), (domain%nx(1)-ecw)
  
  