%---------------------------------------------------
%                     INITIAL.m
%
%   Must specify initial condtions for PHI
%   Must specify:
%       MNI      Maximum number of increments
%       DTIME    Time step
%       INCsv    First increment to save PHI values
%       INTsv    Subsequent intervals to save PHI values
%       IRZ      Axi-symmetric designation
%---------------------------------------------------

%---------------------------------------------------
%    For Test Problem
%---------------------------------------------------
     for I=1:NUMNP
        PHI(I)  = 3*XORD(I) -5*YORD(I);
     end
     PHI(4)=0.0;    % Perturb node 4 from steady state

%-----------------------------------------------------
%    Calculate time step.
%
%    For Code Check: 
%
%    R=3, mu=2, dx = 5, HS=0
%
%    dT = 0.5*(mu)*(dx^2)/R
%    dT = 0.5*(2)*(5^2)/3   = 8.34
%-----------------------------------------------------

     DTIME = 8.00;             % Stable but some oscillation
     DTIME = 11.0;             % Unstable 
     DTIME = 10.0;             % Stable but severe oscillation
     DTIME = 4.00;             % Stable 

     DTIME = 0.02;             % Use for smooth analysis

     MNI   = round(200/DTIME)  % Maximum number of increments
               
     IRZ   = 0;                % x-y coordinates

     INCsv = 100;               % First increment to save PHI values
     INTsv = 100;               % Subsequent intervals to save PHI values


