anchor.diff:        p%du(18 + k) = perturb_b(k)         ! u%InflowOnBlade(:,:,k) = 18 + k
      end do      
      p%du(22) = perturb_t                   ! u%InflowOnTower(:,:) = 22
  -  
  -         
  +   do k=1,p%NumBlades 
  +      p%du(22+k) = perturb                     ! u%UserProp(:,:) = 23,24,25
  +   end do      
         !.....................
         ! get names of linearized inputs
         !.....................
  
positive.code:             index = index + 1
           end do            
        end do
        
     END IF
  
     IF ( PRESENT( y_op ) ) THEN
  
negative.code:       end do !j      
     end do !i
     
     !Module/Mesh/Field: u%UserProp(:,:) = 23,24,25;
     
     do k=1,size(u%UserProp,2) ! p%NumBlades         
        do i=1,size(u%UserProp,1) ! numNodes
              p%Jac_u_indx(index,1) =  22 + k
              p%Jac_u_indx(index,2) =  1 !component index:  this is a scalar, so 1, but is never used
              p%Jac_u_indx(index,3) =  i !Node:   i
              index = index + 1     
        end do !i
     end do !k
        !......................................
        ! default perturbations, p%du:
        !......................................
  