%
%  circular_plate_constraint_nodes.m  ver 1.0  September 12, 2012
%
function[con_nodes,num_con_nodes,con_dof,num_con_dof]=...
                    circular_plate_constraint_nodes(nx,ny,nd,idn,nem)
%
%
num_con_nodes=0;
num_con_dof=0;
%
con_nodes(1)=0;
con_dof(1)=0;
%
disp(' ')
disp(' Constrain rotational dof at boundary points?  1=yes 2=no ');
irot=input(' ');
%
if(irot==1)
%
   for i=1:length(nd)
         con_nodes(i)=nd(i);
   end 
%
   ijk=1;
%
   for i=1:length(nd)
         con_dof(ijk)=3*con_nodes(i)-1;
         ijk=ijk+1;
         con_dof(ijk)=3*con_nodes(i);         
         ijk=ijk+1;      
   end
%
   num_con_dof=length(con_dof);
end