Resolve "Enabling free Maxwell propagation in multisystem mode, part 2/4"

Description

Allow Maxwell propagation with external currents.

News snippet

Allow Maxwell propagation with external currents.

Checklist

• I have checked that my code follows the Octopus coding standards
• I have added tests for all the new features added in this request.

Closes #206 (closed)

src/system/xc_kli_inc.F90

@@ -107,10 +107,9 @@ subroutine X(xc_KLI_solve) (namespace, mesh, gr, hm, st, is, oep, first)
   end do
 
   if (oep%has_photons .and. oep%coctranslation_logical) then
-      oep%pt%pol_dipole_array(1:mesh%np, 1) = oep%pt%pol_dipole_array(1:mesh%np,1 ) + &
+    oep%pt%pol_dipole_array(1:mesh%np, 1) = oep%pt%pol_dipole_array(1:mesh%np,1 ) + &
         dmf_dotp(gr%mesh, sum(st%rho(1:mesh%np, :), dim=2), coctranslation(1:mesh%np))/abs(st%qtot)
-      SAFE_DEALLOCATE_A(coctranslation)
-    end if
+    SAFE_DEALLOCATE_A(coctranslation)
   end if
 
   do ip = 1, mesh%np