source: issm/oecreview/Archive/16133-16554/ISSM-16522-16523.diff@ 16556

../trunk-jpl/test/NightlyRun/test1502.m

@@ -4 +4 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
 md=setflowequation(md,'SSA','all');
-md=extrude(md,3,1.);
+%md=extrude(md,3,1.);
 md.cluster=generic('name',oshostname(),'np',1);
 md.transient.isthermal=0;
 
-md.timestepping.time_step=1.;
-md.settings.output_frequency=1;
-md.timestepping.final_time=2000.;
+%md.timestepping.time_step=1.;
+%md.settings.output_frequency=1;
+%md.timestepping.final_time=2000.;
 
 %Solve for thinning rate -> -1 * surface mass balance
 smb= 2.*ones(md.mesh.numberofvertices,1);
-md.surfaceforcings.mass_balance= smb;
-md.basalforcings.melting_rate= smb;
+%md.surfaceforcings.mass_balance= smb;
+%md.basalforcings.melting_rate= smb;
+md.verbose=verbose('all');
 
 md=solve(md,MasstransportSolutionEnum());
 
+error
 for i=1:10
          md=solve(md,MasstransportSolutionEnum());
          md.surfaceforcings.mass_balance= md.surfaceforcings.mass_balance - ((md.results.MasstransportSolution.Thickness)-md.geometry.thickness);