Changeset 899 for issm/trunk/test/Validation/EISMINT/Transient/Square.par

Timestamp:

06/11/09 12:58:40 (16 years ago)

Author:

Mathieu Morlighem

Message:

fixed runmes.m

File:

• issm/trunk/test/Validation/EISMINT/Transient/Square.par (modified) (2 diffs)

issm/trunk/test/Validation/EISMINT/Transient/Square.par

@@ -85 +85 @@
         disp('      boundary conditions for diagnostic model: ');
         %Build gridonicefront, array of boundary grids belonging to the icefront:
-        gridinsideicefront=ArgusContourToMesh(md.elements,md.x,md.y,expread('Front.exp',1),'node',2);
+        gridinsideicefront=ContourToMesh(md.elements,md.x,md.y,expread('Front.exp',1),'node',2);
         gridonicefront=double(md.gridonboundary & gridinsideicefront);
 
@@ -102 +102 @@
         pos=find(md.y==max(md.y)); %grids on the upper boundary condition
         md.dirichletvalues_diag(pos,2)=400*(((md.x(pos)-100000)/25000).^2-ones(size(pos,1),1)).*heaviside((1+eps)*ones(size(pos,1),1)-((md.x(pos)-100000)/25000).^2);
-        md.segmentonneumann_prog=md.segments(pos,:);
         md.neumannvalues_prog=zeros(size(md.segmentonneumann_prog,1),1);
         md.neumannvalues_prog(:)=NaN; %free radiation