anchor.diff:           
               
               ! get central difference:            
  -         call Compute_dY( p, y, y_m, delta, dYdx(:,i) )
  +         call Compute_dY( p, y_p, y_m, delta, dYdx(:,i) )
            
         end do
         
  
positive.code:    !................
     ! names of the columns, InitOut%LinNames_u:
     !................
     
     index = 1
     if (allocated(u%BladePtLoads)) then
        do k=1,p%NumBl
  
negative.code: !----------------------------------------------------------------------------------------------------------------------------------
  !> This subroutine computes the differences of two meshes and packs that value into appropriate locations in the dY array.
  !! Do not change this packing without making sure subroutines elastodyn::ed_init_jacobian and elastodyn::compute_dy are consistant with this routine!
  SUBROUTINE PackMotionMesh_dY(M_p, M_m, dY, indx_first, FieldMask)
     
     TYPE(MeshType)                    , INTENT(IN   ) :: M_p                        !< ED outputs on given mesh at \f$ u + \Delta u \f$ (p=plus)
     TYPE(MeshType)                    , INTENT(IN   ) :: M_m                        !< ED outputs on given mesh at \f$ u - \Delta u \f$ (m=minus)   
     REAL(ReKi)                        , INTENT(INOUT) :: dY(:)                      !< column of dYdu \f$ \frac{\partial Y}{\partial u_i} = \frac{y_p - y_m}{2 \, \Delta u}\f$ 
     INTEGER(IntKi)                    , INTENT(INOUT) :: indx_first                 !< index into dY array; gives location of next array position to fill
     LOGICAL, OPTIONAL                 , INTENT(IN   ) :: FieldMask(FIELDMASK_SIZE)  !< flags to determine if this field is part of the packing
     
        ! local variables:
     INTEGER(IntKi)                :: ErrStat2 ! we're ignoring the errors about small angles
  