Index: /issm/trunk/scripts/historyISSM.m
===================================================================
--- /issm/trunk/scripts/historyISSM.m	(revision 2528)
+++ /issm/trunk/scripts/historyISSM.m	(revision 2529)
@@ -9,10 +9,4 @@
 !vim runme.m
 !vim configuration.m
-%  --  execution --
-runme
-md.analysis_type=DiagnosticAnalysisEnum;
-md.sub_analysis_type=NoneAnalysisEnum;
-md=solve(md,'analysis_type','diagnostic');
-md=loadresultsfromcluster(md);
 %  --  tests  --
 md=runme('analysis_type',{'diagnostic'});
@@ -27,18 +21,22 @@
 md.numberofgrids
 %  --  model set up  --
-md=extrude(md,8,3);
+md=model;
+md=mesh(md,'DomainOutline.exp',5000);
+md=geography(md,'','');
+md=geography(md,'Iceshelves.exp','Islands.exp');
+md=parameterize(md,'Antarctica.par');
+md=extrude(md,8,1);
 md=setelementstype(md,'pattyn','all','stokes','all');
 md=setelementstype(md,'macayeal','all');
+%  --  execution --
+md.analysis_type=DiagnosticAnalysisEnum;
+md.sub_analysis_type=NoneAnalysisEnum;
+md=solve(md,'analysis_type','diagnostic');
 md=solve(md,'analysis_type','diagnostic','batch','yes');
 md=solve(md,'analysis_type','thermal');
 md=solve(md,'analysis_type','thermal','batch','yes');
+md=loadresultsfromcluster(md);
 %  --  plots  --
+plotmodel(md,'data','vel')
 plotmodel(md,'data',md.results.diagnostic.vel)
 print -dpng vel.png
-%  --  setup --
-md=model;
-md=parameterize(md,'Antarctica.par');
-md=geography(md,'','');
-md=geography(md,'Iceshelves.exp','Islands.exp');
-md=extrude(md,8,1);
-md=mesh(md,'DomainOutline.exp',5000);