Changeset 5112

Timestamp:

08/10/10 09:32:52 (15 years ago)

Author:

seroussi

Message:

more ISMIP tests

Location:

issm/trunk/test

Files:

• Archives/Archive1107.mat (added)
• NightlyRun/test1106.m (modified) (1 diff)
• NightlyRun/test1107.m (modified) (1 diff)
• NightlyRun/test1108.m (modified) (1 diff)

issm/trunk/test/NightlyRun/test1106.m

@@ -2 +2 @@
 %Pattyn and Payne 2006
 
-L_list={5000,10000,20000,40000,80000,160000};
+L_list={5000};
+%L_list={5000,10000,20000,40000,80000,160000};
 results={};
 

issm/trunk/test/NightlyRun/test1107.m

@@ -18 +18 @@
 
         %We need one grd on dirichlet: the 4 corners are set to zero
-        md.spcvelocity=zeros(md.numberofgrids,6);
-        pos=find((md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
-        md.spcvelocity(pos,1:3)=1;
+        %md.spcvelocity=zeros(md.numberofgrids,6);
+        %pos=find((md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
+        %md.spcvelocity(pos,1:3)=1;
 
-        %Create MPCs to have periodic boundary conditions
-        posx=find(md.x==0);
-        posx2=find(md.x==max(md.x));
+        %%Create MPCs to have periodic boundary conditions
+        %posx=find(md.x==0);
+        %posx2=find(md.x==max(md.x));
 
-        posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times
-        posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x));
+        %posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times
+        %posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x));
 
-        md.penalties=[posx,posx2;posy,posy2];
+        %md.penalties=[posx,posx2;posy,posy2];
 
         %Compute the diagnostic

issm/trunk/test/NightlyRun/test1108.m

@@ -18 +18 @@
 
         %We need one grd on dirichlet: the 4 corners are set to zero
-        md.spcvelocity=zeros(md.numberofgrids,6);
-        pos=find((md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
-        md.spcvelocity(pos,1:3)=1;
+        %md.spcvelocity=zeros(md.numberofgrids,6);
+        %pos=find((md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
+        %md.spcvelocity(pos,1:3)=1;
 
-        %Create MPCs to have periodic boundary conditions
-        posx=find(md.x==0);
-        posx2=find(md.x==max(md.x));
+        %%Create MPCs to have periodic boundary conditions
+        %posx=find(md.x==0);
+        %posx2=find(md.x==max(md.x));
 
-        posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times
-        posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x));
+        %posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times
+        %posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x));
 
-        md.penalties=[posx,posx2;posy,posy2];
+        %md.penalties=[posx,posx2;posy,posy2];
 
         %Compute the diagnostic