anchor.diff:     p%rot_elem = (p%dof_node/2) * p%nodes_per_elem
   
   
  -   
  +
      !................................
      ! allocate some parameter arrays
      !................................
  
positive.code: 
  
  !FIXME: these need to be computed based on jmj's equations
     
        ! Set values for tip
     idx_node=y%BldMotion%NNodes
  
  
negative.code:    CALL BD_CopyInput(u(1),u_interp,MESH_NEWCOPY,ErrStat2,ErrMsg2)
        call SetErrStat(ErrStat2,ErrMsg2,ErrStat, ErrMsg, RoutineName)
  
     CALL BD_CopyInput(u(1),u_temp,MESH_NEWCOPY,ErrStat2,ErrMsg2) ! this just needs copies of
        call SetErrStat(ErrStat2,ErrMsg2,ErrStat, ErrMsg, RoutineName)
  
        if (ErrStat >= AbortErrLev) then
  