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/grid/derivatives_inc.F90

@@ -674,7 +674,7 @@ subroutine X(derivatives_partial)(der, ff, op_ff, dir, ghost_update, set_bc)
   POP_SUB(X(derivatives_partial))
 end subroutine X(derivatives_partial)
 
-subroutine X(derivatives_batch_grad)(der, ffb, opffb, ghost_update, set_bc, to_cartesian, metric)
+subroutine X(derivatives_batch_grad)(der, ffb, opffb, ghost_update, set_bc, to_cartesian, metric, factor)
   type(derivatives_t), intent(in)    :: der
   class(batch_t),      intent(inout) :: ffb
   class(batch_t),      intent(inout) :: opffb(:)
@@ -682,6 +682,7 @@ subroutine X(derivatives_batch_grad)(der, ffb, opffb, ghost_update, set_bc, to_c
   logical, optional,   intent(in)    :: set_bc
   logical, optional,   intent(in)    :: to_cartesian
   FLOAT,   optional,   intent(in)    :: metric(:,:)
+  FLOAT,   optional,   intent(in)    :: factor
 
   integer :: idir
   logical :: set_bc_, ghost_update_
@@ -697,7 +698,7 @@ subroutine X(derivatives_batch_grad)(der, ffb, opffb, ghost_update, set_bc, to_c
 
   do idir = 1, der%dim
     call X(derivatives_batch_perform)(der%grad(idir), der, ffb, opffb(idir), &
-      ghost_update=ghost_update_, set_bc=set_bc_)
+      ghost_update=ghost_update_, set_bc=set_bc_, factor=factor)
     set_bc_       = .false. ! there is no need to update again
     ghost_update_ = .false. ! the boundary or ghost points
   end do