har = [.2 .3 .4 .5 .6 .7 .8 .9 1. 2. 3. 4. 5.];

colorch = [0.99 0.99 0.95];

sourcedef  = 4;         %mogi(1), penny(2), spheroid(3), dislocation(4)

strikedef  = '90';     %default fault strike

dipdef     = '89.99'; 	%default fault dip

depthdef   = '1e-3'; 	%default fault depth

widthdef   = '1.5e4'; 	%default fault width

lengthdef  = '7e4'; 	%default fault length

undef    = '0'; 	%default fault opening

usdef    = '1'; 	%default strike-slip displacement

uddef    = '0'; 	%default dip-slip displacement

geoddef    = 1;         % default geodectic measurement

rheodef    = 1;         %default rheology

xmaxdef    = '100';      %default xgrid and ygrid size

xmindef    = '-100';

ymaxdef    = '100';

ymindef    = '-100';

shiftxdef     = '0';      % "x" distance from the source epicenter to the origin 

shiftydef     = '0';      % "y" distance from the source epicenter to the origin 

xspacedef  = '1';

yspacedef  = '1';

hdef    = '0.5';            %defines dimensionless crack depth (see array har(hdef))

mudef    = '30';         %shear modulus, GPa

nudef    = '0.27';       %Poisson's ratio

pressuredef  = '1e6';      %Excess magma pressure
 
adef  = '3e3';      	 %Major axis of a Yang-Davis spheroid, m
 
bdef  = '1e3';      	 %Minor axis of a Yang-Davis spheroid, m
 
thetadef  = '89.5';      	 %Plunge of a major axis of a Yang-Davis spheroid, deg

 u1=[];
 u1= char(evalc(['dir SIC_par.mat']));
 if u1(3) == 'e'
  load SIC_par.mat
 end
