Changeset 13971 for issm/trunk-jpl/test/NightlyRun/test1601.m
- Timestamp:
- 11/15/12 15:49:34 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test1601.m
r13670 r13971 3 3 md=parameterize(md,'../Par/SquareShelf.par'); 4 4 md=setflowequation(md,'macayeal','all'); 5 md.diagnostic.spcvx(find(md.mesh.y>0 ))=NaN;6 md.initialization.vx(:)=0 ;7 md.initialization.vy(:)=0 ;8 md.initialization.vel(:)=0 ;5 md.diagnostic.spcvx(find(md.mesh.y>0.))=NaN; 6 md.initialization.vx(:)=0.; 7 md.initialization.vy(:)=0.; 8 md.initialization.vel(:)=0.; 9 9 10 md.cluster=generic('name',oshostname ,'np',2);10 md.cluster=generic('name',oshostname(),'np',2); 11 11 md=solve(md,DiagnosticSolutionEnum()); 12 12 vel0=md.results.DiagnosticSolution.Vel; 13 13 14 theta=30 *pi/180;14 theta=30.*pi/180.; 15 15 x=md.mesh.x; 16 16 y=md.mesh.y; … … 27 27 28 28 %Now, put CS back to normal except on the side where the spc are applied 29 pos=find(x==0 | x==1000000);29 pos=find(x==0. | x==1000000.); 30 30 md.diagnostic.referential(:)=NaN; 31 31 md.diagnostic.referential(pos,1:3)=repmat([cos(theta),sin(theta),0],size(pos,1),1);
Note:
See TracChangeset
for help on using the changeset viewer.