Changeset 5112


Ignore:
Timestamp:
08/10/10 09:32:52 (15 years ago)
Author:
seroussi
Message:

more ISMIP tests

Location:
issm/trunk/test
Files:
1 added
3 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/test/NightlyRun/test1106.m

    r5106 r5112  
    22%Pattyn and Payne 2006
    33
    4 L_list={5000,10000,20000,40000,80000,160000};
     4L_list={5000};
     5%L_list={5000,10000,20000,40000,80000,160000};
    56results={};
    67
  • issm/trunk/test/NightlyRun/test1107.m

    r5098 r5112  
    1818
    1919        %We need one grd on dirichlet: the 4 corners are set to zero
    20         md.spcvelocity=zeros(md.numberofgrids,6);
    21         pos=find((md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
    22         md.spcvelocity(pos,1:3)=1;
     20        %md.spcvelocity=zeros(md.numberofgrids,6);
     21        %pos=find((md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
     22        %md.spcvelocity(pos,1:3)=1;
    2323
    24         %Create MPCs to have periodic boundary conditions
    25         posx=find(md.x==0);
    26         posx2=find(md.x==max(md.x));
     24        %%Create MPCs to have periodic boundary conditions
     25        %posx=find(md.x==0);
     26        %posx2=find(md.x==max(md.x));
    2727
    28         posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times
    29         posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x));
     28        %posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times
     29        %posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x));
    3030
    31         md.penalties=[posx,posx2;posy,posy2];
     31        %md.penalties=[posx,posx2;posy,posy2];
    3232
    3333        %Compute the diagnostic
  • issm/trunk/test/NightlyRun/test1108.m

    r5098 r5112  
    1818
    1919        %We need one grd on dirichlet: the 4 corners are set to zero
    20         md.spcvelocity=zeros(md.numberofgrids,6);
    21         pos=find((md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
    22         md.spcvelocity(pos,1:3)=1;
     20        %md.spcvelocity=zeros(md.numberofgrids,6);
     21        %pos=find((md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
     22        %md.spcvelocity(pos,1:3)=1;
    2323
    24         %Create MPCs to have periodic boundary conditions
    25         posx=find(md.x==0);
    26         posx2=find(md.x==max(md.x));
     24        %%Create MPCs to have periodic boundary conditions
     25        %posx=find(md.x==0);
     26        %posx2=find(md.x==max(md.x));
    2727
    28         posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times
    29         posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x));
     28        %posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times
     29        %posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x));
    3030
    31         md.penalties=[posx,posx2;posy,posy2];
     31        %md.penalties=[posx,posx2;posy,posy2];
    3232
    3333        %Compute the diagnostic
Note: See TracChangeset for help on using the changeset viewer.