Changeset 5112
- Timestamp:
- 08/10/10 09:32:52 (15 years ago)
- Location:
- issm/trunk/test
- Files:
-
- 1 added
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/NightlyRun/test1106.m
r5106 r5112 2 2 %Pattyn and Payne 2006 3 3 4 L_list={5000,10000,20000,40000,80000,160000}; 4 L_list={5000}; 5 %L_list={5000,10000,20000,40000,80000,160000}; 5 6 results={}; 6 7 -
issm/trunk/test/NightlyRun/test1107.m
r5098 r5112 18 18 19 19 %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; 23 23 24 % Create MPCs to have periodic boundary conditions25 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)); 27 27 28 posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times29 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)); 30 30 31 md.penalties=[posx,posx2;posy,posy2];31 %md.penalties=[posx,posx2;posy,posy2]; 32 32 33 33 %Compute the diagnostic -
issm/trunk/test/NightlyRun/test1108.m
r5098 r5112 18 18 19 19 %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; 23 23 24 % Create MPCs to have periodic boundary conditions25 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)); 27 27 28 posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times29 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)); 30 30 31 md.penalties=[posx,posx2;posy,posy2];31 %md.penalties=[posx,posx2;posy,posy2]; 32 32 33 33 %Compute the diagnostic
Note:
See TracChangeset
for help on using the changeset viewer.