Changeset 17119


Ignore:
Timestamp:
01/15/14 12:00:26 (11 years ago)
Author:
seroussi
Message:

BUG: fixing level set in NR

Location:
issm/trunk-jpl/test/NightlyRun
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/test/NightlyRun/test217.m

    r15948 r17119  
    5454nodeonicefront=zeros(md.mesh.numberofvertices,1);
    5555pos=find(md.mesh.y==ymax); nodeonicefront(pos)=1;
    56 md.mask.ice_levelset=1-nodeonicefront;
     56md.mask.ice_levelset=-1+nodeonicefront;
    5757
    5858md=solve(md,StressbalanceSolutionEnum());

  • issm/trunk-jpl/test/NightlyRun/test217.py

    r15948 r17119  
    6666pos=numpy.nonzero(md.mesh.y==ymax)
    6767nodeonicefront[pos]=1
    68 md.mask.ice_levelset=1-nodeonicefront
     68md.mask.ice_levelset=-1+nodeonicefront
    6969
    7070md=solve(md,StressbalanceSolutionEnum())

  • issm/trunk-jpl/test/NightlyRun/test701.m

    r16938 r17119  
    1111
    1212%mask
    13 md.mask.ice_levelset  = ones(md.mesh.numberofvertices,1);
     13md.mask.ice_levelset  = - ones(md.mesh.numberofvertices,1);
    1414md.mask.ice_levelset(find(vertexflags(md.mesh,2))) = 0;
    1515md.mask.groundedice_levelset = double(md.mesh.x<0)-.5;

  • issm/trunk-jpl/test/NightlyRun/test702.m

    r16938 r17119  
    1313md.geometry.bed       = interp1(x,b,md.mesh.x);
    1414md.geometry.thickness = md.geometry.surface-md.geometry.bed;
    15 md.mask.ice_levelset = ones(md.mesh.numberofvertices,1);
     15md.mask.ice_levelset = - ones(md.mesh.numberofvertices,1);
    1616md.mask.ice_levelset(find(vertexflags(md.mesh,2))) = 0;
    1717md.mask.groundedice_levelset = double(md.mesh.x<0)-.5;
Note: See TracChangeset for help on using the changeset viewer.