anchor.diff:        totx = x+ddx
         xhat_new = xhat+totx
         call dotprod(totx,d,d2)
  -
       end subroutine update_solution_1DVar
   
   end module mod_1DVar
  
positive.code:   private
  
    public :: init_1DVar
    public :: update_solution_1DVar
  
    real , allocatable , dimension(:,:) :: fmKT
  
negative.code:   interface dotprod
      module procedure dotprod_matrix_matrix
      module procedure dotprod_matrix_array
      module procedure dotprod_array_matrix
      module procedure dotprod_array_array
    end interface
  
  