%---------------------------------
%     INITIAL.m
%
%     Specify boundary conditions 
%----------------------------------------------
     for I=1:NUMNP
        if NPcode(I) == 2
           NPBC(I) = 1;
           PHI(I)  = 3;
           Q(I)    = 0;
        elseif NPcode(I) == 4
           NPBC(I) = 1;
           PHI(I)  = 9;
           Q(I)    = 0;
        end
     end
