anchor.diff:        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
         IF( ErrStat >= AbortErrLev ) RETURN
   
  -!Read inputs from Driver/Glue code
  -!1 Global position vector
  -!2 Global rotation tensor
  -!3 Gravity vector
  +   !Read inputs from Driver/Glue code
  +   !1 Global position vector
  +   !2 Global rotation tensor
  +   !3 Gravity vector
      CALL AllocAry(p%GlbPos,3,'Global position vector',ErrStat2,ErrMsg2)
         CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
      CALL AllocAry(p%GlbRot,3,3,'Global rotation tensor',ErrStat2,ErrMsg2)
  
positive.code:           CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
         temp_id = temp_id2
     ENDDO
  ! Compute blade/member/segment lengths and the ratios between member/segment and blade lengths
     CALL AllocAry(p%member_length,InputFileData%member_total,2,'member length array',ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
  
negative.code:          return
        end if
     
     ! Compute initial position vector uuN0 in blade frame
     temp_int = p%node_elem * p%dof_node
     CALL AllocAry(p%uuN0,temp_int,p%elem_total,'uuN0 (initial position) array',ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  