Changeset 5432 for issm/trunk/test/NightlyRun/test207.m

Timestamp:

08/19/10 17:00:09 (15 years ago)

Author:

seroussi

Message:

some new tests for new Stokes

File:

• issm/trunk/test/NightlyRun/test207.m (modified) (1 diff)

issm/trunk/test/NightlyRun/test207.m

@@ -3 +3 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'macayeal','all','stokes','all');
 md.cluster='none';
+md=setelementstype(md,'macayeal','all');
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
+vx=PatchToVec(md.results.DiagnosticSolution.Vx);
+vy=PatchToVec(md.results.DiagnosticSolution.Vy);
+vz=PatchToVec(md.results.DiagnosticSolution.Vz);
+pos=find(md.spcvelocity(:,1));
+md.spcvelocity(pos,3)=1;
+md=setelementstype(md,'stokes','all');
+md.vx=vx;
+md.vy=vy;
+md.vz=vz;
+md.spcvelocity(:,4)=vx;
+md.spcvelocity(:,5)=vy;
+md.spcvelocity(:,6)=vz;
+pos=find(md.borderstokes);
+md.spcvelocity(pos,1)=1;
+md.spcvelocity(pos,2)=1;
+md.spcvelocity(pos,3)=1;
 md=solve(md,'analysis_type',DiagnosticSolutionEnum);