Index: /issm/trunk-jpl/src/m/solve/loadresultsfrombuffer.js
===================================================================
--- /issm/trunk-jpl/src/m/solve/loadresultsfrombuffer.js	(revision 19816)
+++ /issm/trunk-jpl/src/m/solve/loadresultsfrombuffer.js	(revision 19816)
@@ -0,0 +1,26 @@
+function loadresultsfrombuffer(md,buffer,buffersize){
+//LOADRESULTSFROMBUFFER - load results of solution sequence from memory buffer 
+//
+//   Usage:
+//      loadresultsfrombuffer(md,buffer,buffersize);
+
+	//check number of inputs/outputs
+	if (arguments.length!=3) throw Error('loadresultsfrombuffer error message: wrong usage!');
+
+	if (!md.qmu.isdakota){
+
+		//initialize md.results if not a structure yet
+		if (MapIsEmpty(md.results)) md.results={};
+
+		//load results onto model
+		var array=parseresultsfrombuffer(buffer,buffersize);
+		if (array.length==0) throw Error('No result found in binary buffer. Check for solution crash.');
+
+
+		md.results[array[0]['SolutionType']] = array;
+		
+		//recover solution_type from results
+		md.private.solution=array['SolutionType'];
+	}
+	else throw Error('loadresultsfrombuffer error message: qmu results not supported yet!');
+}
Index: /issm/trunk-jpl/src/m/solve/parseresultsfrombuffer.js
===================================================================
--- /issm/trunk-jpl/src/m/solve/parseresultsfrombuffer.js	(revision 19816)
+++ /issm/trunk-jpl/src/m/solve/parseresultsfrombuffer.js	(revision 19816)
@@ -0,0 +1,116 @@
+function parseresultsfrombuffer(buffer,buffersize){ //{{{
+//PARSERESULTSFROMDISK - ...
+//
+//   Usage:
+//      results=parseresultsfrombuffer(buffer,buffersize);
+
+	//Open file
+	results=[];
+	
+	var fid = new fileptr('mode','r','buffer',buffer,'buffersize',buffersize);
+
+	//Read fields until the end of the file.
+	result  = ReadData(fid);
+
+	if (MapIsEmpty(result))throw Error('no results found in binary memory buffer');
+	
+	var counter = 0;
+	var check_nomoresteps=0;
+	var step    = result['step'];
+
+	while (!MapIsEmpty(result)){
+
+		if (check_nomoresteps){
+			//check that the new result does not add a step, which would be an error: 
+			if (result['step']>=1)throw Error('parsing results for a steady-state core, which incorporates transient results!');
+		}
+
+		//Check step, increase counter if this is a new step
+		if(step!=result['step'] & result['step']>1){
+			counter = counter + 1;
+			step    = result['step'];
+		}
+
+		//Add result
+		if(result['step']==0){
+			//if we have a step = 0, this is a steady state solutoin, don't expect more steps. 
+			index = 0;
+			check_nomoresteps=1;
+		}
+		else if(result['step']==1){
+			index = 0;
+		}
+		else index = counter;
+
+		if(index>results.length-1)results.push({});
+		for(var i=results.length-1;i<index-1;i++)results[i]={};
+		results[index][result['fieldname']]=result['field'];
+		
+		//Get time and step
+		if(result['time']!=-9999){
+			results[index]['step']=result['step'];
+		}
+		if(result['time']!=-9999){
+			results[index]['time']=result['time'];
+		}
+
+		//read next result
+		result  = ReadData(fid);
+	}
+	return results;
+} // }}}
+function ReadData(fid){ //{{{
+//READDATA - ...
+//
+//   Usage:
+//      field=ReadData(fid)
+
+	//read field
+	var length=fid.fread(1,'int');
+	
+	var result={};
+
+	if (length!==-1){
+		fieldname=fid.fread(length,'char');
+		time=fid.fread(1,'double');
+		step=fid.fread(1,'int');
+		type=fid.fread(1,'int');
+		M=fid.fread(1,'int');
+
+		if (type==1) field=fid.fread(M,'double');
+		else if (type==2) field=fid.fread(M,'char');
+		else if (type==3) {
+			N=fid.fread(1,'int');
+			field=fid.fread(N*M,'double');
+		}
+		else throw Error(sprintf("%s%i",'ReadData error message: cannot read data of type ',type));
+
+		//Process units here FIXME: this should not be done here!
+		var yts=365.0*24.0*3600.0;
+		if (fieldname == 'BalancethicknessThickeningRate') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'HydrologyWaterVx') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'HydrologyWaterVy') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'Vx') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'Vy') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'Vz') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'Vel') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'BasalforcingsGroundediceMeltingRate') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'BasalforcingsFloatingiceMeltingRate') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'TotalSmb') for (var i=0;i<field.length;i++)field[i]= field[i]/Math.pow(10,12)*yts; //(GigaTon/year)
+		else if (fieldname == 'SmbMassBalance') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'SmbPrecipitation') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'SmbRunoff') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'SmbCondensation') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'SmbAccumulation') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'SmbMelt') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+		else if (fieldname == 'CalvingCalvingrate') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
+
+		result['fieldname']=fieldname;
+		result['time']=time;
+		if (result['time']!=-9999) result['time']=time/yts;
+		result['step']=step;
+		result['field']=field;
+	}
+	return result;
+
+} // }}}
Index: /issm/trunk-jpl/src/m/solve/solve.js
===================================================================
--- /issm/trunk-jpl/src/m/solve/solve.js	(revision 19815)
+++ /issm/trunk-jpl/src/m/solve/solve.js	(revision 19816)
@@ -183,5 +183,12 @@
 		
 		//Call issm:
-		issm(fid, toolkitsstring, solutionstring, md.miscellaneous.name);
+		var outputs = issm(fid, toolkitsstring, solutionstring, md.miscellaneous.name); 
+		
+		//Recover output arguments: 
+		var outputbuffer = outputs[0]; var outputbuffersize = outputs[1];
+			
+		//Load results: 
+		loadresultsfrombuffer(md,outputbuffer,outputbuffersize);
+
 	}
 }
Index: /issm/trunk-jpl/test/NightlyRun/test101.html
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test101.html	(revision 19815)
+++ /issm/trunk-jpl/test/NightlyRun/test101.html	(revision 19816)
@@ -73,4 +73,6 @@
 <script type="text/javascript" src="../../src/m/inversions/marshallcostfunctions.js"></script>
 <script type="text/javascript" src="../../src/m/solve/marshall.js"></script>
+<script type="text/javascript" src="../../src/m/solve/loadresultsfrombuffer.js"></script>
+<script type="text/javascript" src="../../src/m/solve/parseresultsfrombuffer.js"></script>
 <script type="text/javascript" src="../../src/m/solve/WriteData.js"></script>
 <script type="text/javascript" src="../../src/wrappers/TriMesh/TriMesh.js"></script>
@@ -90,18 +92,23 @@
 <body> 
 	
-	<script type="text/javascript" async><!--}}}-->
+<script type="text/javascript" async><!--}}}-->
 
+	console.log('creating model'); 
 	var md = new model();
-	triangle(md,square[0],100000); 
+
+	console.log('meshing');
+	triangle(md,square[0],40000); 
+	
+	console.log('parameterization');
 	setmask(md,'all','');
 	parameterize(md);
 	setflowequation(md,'SSA','all');
-	md.verbose.solution=1;  md.verbose.convergence=1;
-	md=solve(md,StressbalanceSolutionEnum(),'checkconsistency','no');
+	md.verbose.solution=1;  md.verbose.convergence=0;
 
-	/*plotmodel(md,'data',md.initialization.vel,'colorbar','on','colorbarheight',.5,'data',md.initialization.vx,
-	  'data',md.initialization.temperature, 'data',md.basalforcings.floatingice_melting_rate,'data','mesh');*/
-	//var output=issm([5,1,1]);
+	console.log('solve');
+	solve(md,StressbalanceSolutionEnum(),'checkconsistency','no');
 
+	console.log('plotting');
+	plotmodel(md,'data',md.results['StressbalanceSolution'][0]['Vel'],'colorbar#1','on','colorbarheight',.5);
 
 </script> <!--{{{-->