%wich steps to perform steps are from 1 to 8 
%step 7 is specific to ISMIPA
%step 8 is specific to ISMIPF

steps=[1];

% parameter file to be used, choose between IsmipA.par or IsmipF.par
ParamFile='IsmipA.par'

%Run Steps

% {{{ Mesh Generation #1
if any(steps==1) 

	%initialize md as a new model #help model
	%->

	% generate a squaremesh #help squaremesh
	% Side is 80 km long with 20 points
	%->

	% plot the given mesh #plotdoc
	%->

	% save the given model
	%->

end
% }}}

% {{{ Masks #2
if any(steps==2) 

	% load the preceding step #help loadmodel
	% path is given by the organizer with the name of the given step
	%->

	% set the mask #help setmask
	% all MISMIP nodes are grounded
	%->

	% plot the given mask #md.mask to locate the field
	%->

	% save the given model
	%->

end
% }}}

% {{{ Parameterization #3
if any(steps==3) 

	% load the preceding step #help loadmodel
	% path is given by the organizer with the name of the given step
	%->

	% parametrize the model # help parameterize
	% you will need to fil-up the parameter file (given by the
  % ParamFile variable)
	%->

	% save the given model
	%->

end
% }}}

% {{{ Extrusion #4
if any(steps==4)
	
	% load the preceding step #help loadmodel
	% path is given by the organizer with the name of the given step
	%->

	% vertically extrude the preceding mesh #help extrude
	% only 5 layers exponent 1
	%->

	% plot the 3D geometry #plotdoc
	%->

	% save the given model
	%->

end
% }}}

% {{{ Set the flow computing method #5

if any(steps==5)

	% load the preceding step #help loadmodel
	% path is given by the organizer with the name of the given step
	%->

	% set the approximation for the flow computation #help setflowequation
	% We will be using the Higher Order Model (HO)
	%->

	% save the given model
	%->

end
% }}}

% {{{ Set Boundary Conditions #6
if any(steps==6)

	% load the preceding step #help loadmodel
	% path is given by the organizer with the name of the given step
	%->

	% dirichlet boundary condition are known as SPCs
	% ice frozen to the base, no velocity	#md.stressbalance
	% SPCs are initialized at NaN one value per vertex
	%->

	%->

	%->

	% extract the nodenumbers at the base #md.mesh.vertexonbase
	%->

  % set the sliding to zero on the bed (Vx and Vy)
	%->

	%->

	% periodic boundaries have to be fixed on the sides
	% create tabs with the side of the domain
	% for x
	% create maxX #help find
	%->

	% create minX
	%->

	% for y, max X and minX should be excluded
	% create maxY
	%->

	% create minY
	%->

	% set the node that should be paired together
	% #md.stressbalance.vertex_pairing
	%->

	if (ParamFile=='IsmipF.par')
		% if we are dealing with IsmipF the solution is in masstransport
		md.masstransport.vertex_pairing=md.stressbalance.vertex_pairing;
  end
	% save the given model
	%->

end
% }}}

% {{{ Solving #7
if any(steps==7)
	% load the preceding step #help loadmodel
	% path is given by the organizer with the name of the given step
	%->

	% Set cluster #md.cluster
	% generic parameters #help generic
	% set only the name and number of process
	%->

	% Set which control message you want to see #help verbose
	%->

	% Solve #help solve
	% we are solving a StressBalanc
	%->

	% save the given model
	%->

	% plot the surface velocities #plotdoc
	%->
end
% }}}

% {{{ Solving #8
if any(steps==8)
	% load the preceding step #help loadmodel
	% path is given by the organizer with the name of the given step
	%->

	% Set cluster #md.cluster
	% generic parameters #help generic
	% set only the name and number of process
	%->

	% Set which control message you want to see #help verbose
	%->

	% set the transient model to ignore the thermal model
	% #md.transient 
	%->

	% define the timestepping scheme
	% everything here should be provided in years #md.timestepping
	% give the length of the time_step (4 years)
	%->

	% give final_time (20*4 years time_steps)
	%->

	% Solve #help solve
	% we are solving a TransientSolution
	%->

	% save the given model
	%->

	% plot the surface velocities #plotdoc
	%->

end
% }}}
