%
%  qa_ptdambar
%
%  This script runs debug check of 'ptdambar',
%  which plots the damage indicator for bar-like
%  geometry

%
%  Define the geometry and damage index
%
Nodes = [
   1 0 0 0
   2 1 1 0
   3 2 0 0];
      
Bars = [
   1 1 2
   2 2 3
   3 1 3];
     
DIvec = [
   0
   1
   2];
   
inodal = 1;

%
%  Define the plot properties
%

colormap cool

haxes = axes;

ngrad = 24;

naxaz = 330;

naxel = 30;

%
%  Call ptdambar
%

ptdambar(Nodes,Bars,DIvec,haxes,inodal,ngrad,naxaz,naxel);