Index: /issm/trunk/doc/issmdoc.m
===================================================================
--- /issm/trunk/doc/issmdoc.m	(revision 11028)
+++ /issm/trunk/doc/issmdoc.m	(revision 11029)
@@ -8,8 +8,8 @@
 disp(sprintf('%s','	   go to ',ISSM_TIER,'/examples/SquareIceshelf'));
 disp(sprintf('%s','	   md=model;                                %creates a new empty model structure'));
-disp(sprintf('%s','	   md=mesh(md,''DomainOutline.exp'',50000);   %creates a mesh of the domain outline with a resolution of 50000m'));
+disp(sprintf('%s','	   md=triangle(md,''DomainOutline.exp'',50000);   %creates a mesh of the domain outline with a resolution of 50000m'));
 disp(sprintf('%s','	   md=geography(md,''all'','''');               %defines the glacier system as an ice shelf (no island)'));
 disp(sprintf('%s','	   md=parameterize(md,''Square.par'');        %fills all the other fields of the model'));
 disp(sprintf('%s','	   md=setelementstype(md,''macayeal'',''all''); %defines all elements as MacAyeal''s'));
-disp(sprintf('%s','	   md=solve(md,''analysis_type'',''diagnostic'',''package'',''ice'');   %generate the velocity field of the ice shelf using the solver ICE'));
-disp(sprintf('%s','	   plotmodel(md,''data'',md.results.diagnostic.vel);              %displays the velocity (type plotdoc for plotmodel help)'));
+disp(sprintf('%s','	   md=solve(md,DiagnosticSolutionEnum);   %generate the velocity field'));
+disp(sprintf('%s','	   plotmodel(md,''data'',md.results.DiagnosticSolution.Vel);    %displays the velocity (type plotdoc for plotmodel help)'));