Changeset 10382
- Timestamp:
- 10/31/11 15:48:23 (13 years ago)
- Location:
- issm/trunk/test
- Files:
-
- 2 added
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/NightlyRun/test1601.m
r10374 r10382 1 1 md=setmesh(model,'../Exp/Square.exp',150000); 2 2 md=setmask(md,'all',''); 3 md.settings.results_on_vertices=1;4 3 md=parameterize(md,'../Par/SquareShelf.par'); 5 4 md=setflowequation(md,'macayeal','all'); … … 9 8 md.initialization.vel(:)=0; 10 9 10 md.settings.results_on_vertices=1; 11 11 md.cluster=generic('name',oshostname,'np',2); 12 12 md=solve(md,DiagnosticSolutionEnum); … … 25 25 26 26 plotmodel(md,'data',vel0,'data',vel1,'data',vel1-vel0,'title','Cartesian CS','title','Rotated CS','title','difference') 27 disp(['Error between Cartesian and rotated CS: ' num2str(max( vel0-vel1))]);27 disp(['Error between Cartesian and rotated CS: ' num2str(max(abs(vel0-vel1))/(max(abs(vel0))+eps)) ]); 28 28 29 29 %Fields and tolerances to track changes
Note:
See TracChangeset
for help on using the changeset viewer.