%---------------------------------------------------
%                  INITIAL.m
%
%  The following arrays must be initialized by user
%  if they differ from default values (all zeros):
%
%      NPBC, XBC, YBC, TX, TY, THETA, MAT
%---------------------------------------------------
      for I=1:NUMNP
        if NPcode(I) == 1;
           NPBC(I) = 1;
           XBC(I)  = 0.0;
           YBC(I)  = 0.0;
           TX(I)   = 0.0;
           TY(I)   = 0.0;
           THETA(I)= 0.0;
        elseif NPcode(I) == 3;
           NPBC(I) = 3;
           XBC(I)  = 0.0;
           YBC(I)  = 0.0;
           TX(I)   = 0.0;
           TY(I)   = 0.0;
           THETA(I)= 0.3;
        elseif NPcode(I) == 7;
           angle   = atan(10./6.);
           NPBC(I) = 0;
           XBC(I)  = 0.0;
           YBC(I)  = 0.0;
           TX(I)   = 300/sqrt(136)*cos(angle);
           TY(I)   = 300/sqrt(136)*sin(angle);
           THETA(I)= -angle;
        end
      end

