Ignore:
Timestamp:
08/19/10 17:00:09 (15 years ago)
Author:
seroussi
Message:

some new tests for new Stokes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/test/NightlyRun/test208.m

    r5098 r5432  
    33md=parameterize(md,'../Par/SquareShelf.par');
    44md=extrude(md,3,2);
    5 md=setelementstype(md,'macayeal','all','stokes','all');
     5md=setelementstype(md,'macayeal','all');
    66md.cluster=oshostname;
     7md=solve(md,'analysis_type',DiagnosticSolutionEnum);
     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;
    724md=solve(md,'analysis_type',DiagnosticSolutionEnum);
    825
Note: See TracChangeset for help on using the changeset viewer.