source: issm/trunk/test/NightlyRun/test408.m@ 5432

Last change on this file since 5432 was 5432, checked in by seroussi, 15 years ago

some new tests for new Stokes

File size: 1.2 KB
RevLine 
[4996]1md=mesh(model,'../Exp/Square.exp',180000);
2md=geography(md,'../Exp/SquareShelf.exp','');
3md=parameterize(md,'../Par/SquareSheetShelf.par');
4md=extrude(md,3,1);
[5432]5md=setelementstype(md,'pattyn','all');
6md.cluster=oshostname;
[4996]7md=solve(md,'analysis_type',DiagnosticSolutionEnum);
[5432]8vx=PatchToVec(md.results.DiagnosticSolution.Vx);
9vy=PatchToVec(md.results.DiagnosticSolution.Vy);
10vz=PatchToVec(md.results.DiagnosticSolution.Vz);
11pos=find(md.spcvelocity(:,1));
12md.spcvelocity(pos,3)=1;
13md=setelementstype(md,'stokes','all');
14md.vx=vx;
15md.vy=vy;
16md.vz=vz;
17md.spcvelocity(:,4)=vx;
18md.spcvelocity(:,5)=vy;
19md.spcvelocity(:,6)=vz;
20pos=find(md.borderstokes);
21md.spcvelocity(pos,1)=1;
22md.spcvelocity(pos,2)=1;
23md.spcvelocity(pos,3)=1;
24md=solve(md,'analysis_type',DiagnosticSolutionEnum);
[5098]25
26%Fields and tolerances to track changes
27field_names ={'Vx','Vy','Vz','Vel','Pressure'};
[5101]28field_tolerances={1e-08,1e-08,1e-07,1e-08,1e-08};
[5098]29field_values={...
30 PatchToVec(md.results.DiagnosticSolution.Vx),...
31 PatchToVec(md.results.DiagnosticSolution.Vy),...
32 PatchToVec(md.results.DiagnosticSolution.Vz),...
33 PatchToVec(md.results.DiagnosticSolution.Vel),...
34 PatchToVec(md.results.DiagnosticSolution.Pressure),...
35 };
Note: See TracBrowser for help on using the repository browser.