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:         ! Optionally report integrated energy
          if(report_energy_statistics) then
  
              ecw = domain%exterior_cells_width
              energy_potential = 0.0_force_double
              energy_kinetic = 0.0_force_double
  
              if(domain%is_staggered_grid .and. domain%linear_solver_is_truely_linear) then
                  !
  
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)
  
  