anchor.diff:  
      CALL AllocAry(p%Stif0_GL,6,6,p%ngp*p%elem_total,'Stif0_GL',ErrStat2,ErrMsg2)
         CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  +   p%Stif0_GL(:,:,:) = 0.0D0
      CALL AllocAry(p%Mass0_GL,6,6,p%ngp*p%elem_total,'Mass0_GL',ErrStat2,ErrMsg2)
         CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  +   p%Mass0_GL(:,:,:) = 0.0D0
      DO i=1,p%elem_total
          DO j=1,p%ngp
              temp_id = (i-1)*p%ngp+j
  
positive.code:    temp_L2(1:3,1) = p%uuN0(1:3,1)
     temp_L2(1:3,p%ngp*p%elem_total+2) = p%uuN0(temp_int-5:temp_int-3,p%elem_total)
     DEALLOCATE(temp_GLL)
  
  ! 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 )
     DO i=1,p%ngp
         temp_GL(i) = (temp_GL(i) + 1.0D0)/2.0D0
     ENDDO
  
     DO i=1,p%elem_total
         IF(i .EQ. 1) THEN
             DO j=1,p%ngp
  
negative.code:            ENDDO
         ENDDO
     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 )
     p%beta(:)  = 0.0D0
     p%beta(:)  = InputFileData%InpBl%beta(:)
  
     IF( ErrStat >= AbortErrLev ) RETURN
  