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:    DEALLOCATE(temp_GLL)
     DEALLOCATE(SP_Coef)
  
  ! Compute sectional propertities ( 6 by 6 stiffness and mass matrices)
  ! at Gauss points
     CALL AllocAry(temp_ratio,p%ngp,p%elem_total,'temp_ratio',ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
        if (ErrStat >= AbortErrLev) then
  
negative.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 )
  