
% In this section you may modify the mesh in any way whatsoever.

    m=leaf_plotoptions(m,'layeroffset', 0.2);    
    % set colour of polariser gradient arrows
    m=leaf_plotoptions(m,'highgradcolor',[0,0,0],'lowgradcolor',[1,0,0]);
    m=leaf_plotoptions(m,'decorscale',1.5);
    m=leaf_plotoptions(m,'arrowthickness',1.5);
    m=leaf_plotoptions(m,'sidegrad','AB'); %polariser gradient will be plotted on both sides
if (Steps(m)==0) && m.globalDynamicProps.doinit % First iteration
    
    % Set up names for variant models.  Useful for running multiple models on a cluster.
    
    m.userdata.ranges.modelname.range {1} = 'Figure 9 A-F. Phase I + II of div domes'; %div model with two morphogenetic switch at day 12 when boosting the model4 growth rates in kper and kpar to produce cell files and then reducing then from day 15 to make shape
    m.userdata.ranges.modelname.range {2} = 'Figure 9J. No orthogonal conflict'; %div model with two morphogenetic switch at day 12 when boosting the model4 growth rates in kper and kpar to produce cell files and then reducing then from day 15 to make shape
    m.userdata.ranges.modelname.range {3} = 'Figure 9I. No directional conflict and specified anisotropy'; %div model with two morphogenetic switch at day 12 when boosting the model4 growth rates in kper and kpar to produce cell files and then reducing then from day 15 to make shape
    m.userdata.ranges.modelname.range {4} = 'Figure 9G. No surface conflict'; %div model with two morphogenetic switch at day 12 when boosting the model4 growth rates in kper and kpar to produce cell files and then reducing then from day 15 to make shape
    m.userdata.ranges.modelname.range {5} = 'Figure 9H. No areal conflict'; %div model with two morphogenetic switch at day 12 when boosting the model4 growth rates in kper and kpar to produce cell files and then reducing then from day 15 to make shape
    m.userdata.ranges.modelname.range {6} = 'WT_criss-cross'; %div model with two morphogenetic switch at day 12 when boosting the model4 growth rates in kper and kpar to produce cell files and then reducing then from day 15 to make shape
    m.userdata.ranges.modelname.range {7} = 'Figure 10A-F.   Phase I + II of wild-type wedge'; %WT model with two morphogenetic switch at day 12 when boosting the model4 growth rates in kper and kpar to produce cell files and then reducing then from day 15 to make shape
    m.userdata.ranges.modelname.range {8} = 'WT_criss-cross + sinus_A+B'; %WT model based on model6 from palate&lip modelling
    m.userdata.ranges.modelname.range {9} = 'Figure 10I. No directional conflict'; %WT model with two morphogenetic switch at day 12 when boosting the model4 growth rates in kper and kpar to produce cell files and then reducing then from day 15 to make shape
    m.userdata.ranges.modelname.range {10} = 'Figure 10J. No orthogonal directional conflict'; %WT model with two morphogenetic switch at day 12 when boosting the model4 growth rates in kper and kpar to produce cell files and then reducing then from day 15 to make shape
    m.userdata.ranges.modelname.range {11} = 'WT_criss-cross + sinus_B_no surface conflict'; %WT model with two morphogenetic switch at day 12 when boosting the model4 growth rates in kper and kpar to produce cell files and then reducing then from day 15 to make shape
    m.userdata.ranges.modelname.range {12} = 'Figure 10H. No areal conflict'; %WT model with two morphogenetic switch at day 12 when boosting the model4 growth rates in kper and kpar to produce cell files and then reducing then from day 15 to make shape
    m.userdata.ranges.modelname.range {13} = 'WT_timing all at day 10'; %WT model with two morphogenetic switch at day 12 when boosting the model4 growth rates in kper and kpar to produce cell files and then reducing then from day 15 to make shape
    m.userdata.ranges.modelname.index =11;                       % CLUSTER
    
    
    %m = leaf_setzeroz( m ); % keeps the noise in the canvas to zero
    m = leaf_bowlz(m, 0.01 ); %induce a slightly curve in the medio-lateral axis of the canvas
    m = leaf_setthickness(m, 0.06); %thickness of canvas, relates to the 60 um of the real data
    
    %     %ID_EDGE
    %    m = leaf_mgen_edge( m, 'ID_EDGE', 1);
    % %     id_edge_p = m.morphogens(:,id_edge_i);
    %
    %    m = leaf_fix_vertex( m, 'ID_EDGE', 'dfs', 'z' );
    %  m = leaf_fix_vertex( m, 'vertex', find(m.nodes(:,2)>=min(m.nodes(:,2)+0.01)), 'dfs', 'y' );
    %     m = leaf_fix_vertex( m, 'vertex', find(m.nodes(:,2)>=max(m.nodes(:,2)-0.01)), 'dfs', 'z' );
    %
    %m = leaf_fix_vertex( m, 'vertex', find(m.nodes(:,1)<=min(m.nodes(:,1)+0.01)), 'dfs', 'xyz' );
    %m = leaf_fix_vertex( m, 'vertex', find(m.nodes(:,1)>=max(m.nodes(:,1)-0.01)), 'dfs', 'xyz' );
    

    
    m.globalProps.timestep=10; %each step relates to 10 hours in real time
    
    
end
modelname = m.userdata.ranges.modelname.range{m.userdata.ranges.modelname.index};  % CLUSTER
disp(sprintf('model2',mfilename, modelname));
switch modelname
    case {'MODEL1','MODEL2'}
        
    otherwise
        % If you reach here, you probably forgot a case.
end