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

% If needed force FE to subdivide (increase number FE's) here
% if realtime==280+dt
% m = leaf_subdivide( m, 'morphogen','id_vent',...
%       'min',0.5,'max',1,...
%       'mode','mid','levels','all');
% end
% Cut the mesh along the seams (see above)
% if m.userdata.CutOpen==1
%    m=leaf_dissect(m);
%    m.userdata.CutOpen=2;
%    Relax accumulated stresses slowly i.e. 0.95 to 0.999
%    m = leaf_setproperty( m, 'freezing', 0.999 );
% end