source: issm/trunk/examples/ISMIP/runme.m

%which 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
        % Find the indices of the sides of the domain, for x and then for y
        % for x
        % create maxX, list of indices where x is equal to max of x (use >> help find)
        %->

        % create minX, list of indices where x is equal to min of x
        %->

        % for y
        % create maxY, list of indices where y is equal to max of y
        % but not where x is equal to max or min of x
        % (i.e, indices in maxX and minX should be excluded from maxY and minY) 
        %->

        % create minY, list of indices where y is equal to max of y
        % but not where x is equal to max or min of x
        %->

        % set the node that should be paired together, minX with maxX and minY with maxY
        % #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