anchor.diff:           return
         end if
      
  -! Compute initial position vector uuN0 in blade frame
  +   ! 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 )
  
positive.code:    ENDDO
     DEALLOCATE(temp_GL)
     DEALLOCATE(temp_ratio)
  ! Physical damping flag and 6 damping coefficients
     p%damp_flag  = InputFileData%InpBl%damp_flag
     CALL AllocAry(p%beta,6,'Damping coefficient',ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
negative.code:       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
     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)
  