Ignore:
Timestamp:
07/23/13 12:13:24 (12 years ago)
Author:
seroussi
Message:

CHG: no need to redo icefront manually

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TabularUnified issm/trunk-jpl/test/NightlyRun/test218.m

    r15558 r15560  
    11md=squaremesh(model(),1000000,1000000,5,5);
    22md=setmask(md,'all','');
    3 md=parameterize(md,'../Par/SquareShelf2.par');
     3md=parameterize(md,'../Par/SquareShelf.par');
    44md=setflowequation(md,'macayeal','all');
    55md.cluster=generic('name',oshostname(),'np',3);
     
    4848md.diagnostic.spcvy(pos)=0;
    4949
    50 %icefront
    51 nodeonicefront=zeros(md.mesh.numberofvertices,1);
    52 pos=find(md.mesh.y==ymax); nodeonicefront(pos)=1;
    53 pos=find(nodeonicefront(md.mesh.segments(:,1)) | nodeonicefront(md.mesh.segments(:,2))); diagnostic.icefront=md.mesh.segments(pos,:);
    54 diagnostic.icefront=[diagnostic.icefront 1*md.mask.elementonfloatingice(diagnostic.icefront(:,end))];
    55 md.diagnostic.icefront=diagnostic.icefront;
    56 
    5750%partitioning
    5851md.qmu.numberofpartitions=md.mesh.numberofvertices;
Note: See TracChangeset for help on using the changeset viewer.