%
%  qa_ptdamshl
%
%  This script runs debug check of 'ptdamshl',
%  which plots the damage indicator for shell-like
%  geometry

%
%  Define the geometry and damage index
%
Nodes = [
   1 0 0 0
   2 1 0 0
   3 1 1 0
   4 0 1 0
   5 2 0 0
   6 2 1 0
];
      
Shells = [
   1 1 2 3 4
   2 2 5 6 3
    ];
    
%inodal = 1;

%DIvec = [
%   0
%   25
%   8
%   4
%   7
%   6];

inodal = 0;

DIvec(1,:,:) = [
   0  12
   25 16
   ];

DIvec(2,:,:) = [
   4 16
   8 30
   ];
%
%  Define the plot properties
%

colormap cool

haxes = axes;

ngrad = 2;

naxaz = 330;

naxel = 30;

%
%  Call ptdamshl
%

ptdamshl(Nodes,Shells,DIvec,haxes,inodal,ngrad,naxaz,naxel);