Resolve "Implementation of ZORA"

Description

Add here a short description of what this merge request is about. Please mention if you changed any input option or the defaults in any way. Also say here if this merge request fixes any issue (see the gitlab documentation on how to do this).

News snippet

Add here a one line description to be included in the Changes page. You can delete this section if the changes are trivial and of no consequence for the users.

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 #605 (closed)

src/hamiltonian/hamiltonian_elec.F90

@@ -1348,7 +1348,7 @@ contains
     logical,                  intent(in)    :: accel_copy
 
     FLOAT, allocatable :: zora_scpot(:)
-    integer :: ispin, ip, idim
+    integer :: ispin, ip
 
     PUSH_SUB(hamiltonian_elec_update_pot)
 
@@ -1394,7 +1394,7 @@ contains
         do ip = 1, mesh%np
           this%hm_base%potential(ip, ispin) = M_ZERO
         end do
-        
+
 !       !$omp do simd schedule(static)
 !       continue
 !       do ip = 1, mesh%np
@@ -1433,7 +1433,7 @@ contains
       ! zora_scpot = V
       ! zora_pot   = \frac{c^2}{2c^2 - V}
       do ip = 1, mesh%np
-      ! added extra factor of 2 because of -1/2 in kinetic energy
+        ! added extra factor of 2 because of -1/2 in kinetic energy
         do ispin= 1, this%d%spin_channels
           this%zora_pot(ip,ispin) = M_TWO*P_C**2 / (M_TWO*this%mass*P_C**2 - zora_scpot(ip) )
         end do