Changeset 11190
- Timestamp:
- 01/24/12 11:43:35 (13 years ago)
- Location:
- issm/trunk-jpl/test
- Files:
-
- 6 added
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test1106.m
r11031 r11190 6 6 7 7 for i=1:length(L_list), 8 L=3*L_list{i}; %in m (3 times the desired lenght for BC problems) 9 nx=60; %number of nodes in x direction 10 ny=60; 11 md=model; 12 md=squaremesh(md,L,L,nx,ny); 8 L=L_list{i}; 9 md=triangle(model,['../Exp/Square_' num2str(L) '.exp'],L/10); %size 3*L 13 10 md=setmask(md,'',''); %ice sheet test 14 11 md=parameterize(md,'../Par/ISMIPC.par'); 15 md=extrude(md,6,1); 12 md.friction.coefficient=sqrt(md.constants.yts.*(1000+1000*sin(md.mesh.x*2*pi/L).*sin(md.mesh.y*2*pi/L))); 13 md=extrude(md,10,1); 14 15 %Add spc on the borders 16 pos=find((md.mesh.x==0 | md.mesh.x==max(md.mesh.x) | md.mesh.y==0 | md.mesh.y==max(md.mesh.y))); 17 md.diagnostic.spcvx(pos)=0; 18 md.diagnostic.spcvy(pos)=0; 19 if(L==5000), 20 md.diagnostic.spcvx(pos)=15.66; 21 md.diagnostic.spcvy(pos)=-0.1967; 22 elseif(L==10000), 23 md.diagnostic.spcvx(pos)=16.04; 24 md.diagnostic.spcvy(pos)=-0.1977; 25 elseif(L==20000), 26 md.diagnostic.spcvx(pos)=16.53; 27 md.diagnostic.spcvy(pos)=-1.27; 28 elseif(L==40000), 29 md.diagnostic.spcvx(pos)=17.23; 30 md.diagnostic.spcvy(pos)=-3.17; 31 elseif(L==80000), 32 md.diagnostic.spcvx(pos)=16.68; 33 md.diagnostic.spcvy(pos)=-2.69; 34 elseif(L==160000), 35 md.diagnostic.spcvx(pos)=16.03; 36 md.diagnostic.spcvy(pos)=-1.27; 37 end 16 38 17 39 md=setflowequation(md,'stokes','all'); 18 40 19 %Create MPCs to have periodic boundary conditions20 21 %posx=find(md.mesh.x==0);22 %posx2=find(md.mesh.x==L);23 24 %posy=find(md.mesh.y==0 & md.mesh.x~=0 & md.mesh.x~=L); %Don't take the same nodes two times25 %posy2=find(md.mesh.y==L & md.mesh.x~=0 & md.mesh.x~=L);26 27 %md.diagnostic.vertex_pairing=[posx,posx2;posy,posy2];28 29 41 %Compute the diagnostic 30 md.verbose=verbose('convergence',true,'solution',true);31 42 md.cluster=generic('name',oshostname(),'np',8); 32 43 md=solve(md,DiagnosticSolutionEnum); … … 49 60 'Vx80km','Vy80km','Vz80km',... 50 61 'Vx160km','Vy160km','Vz160km' 51 } 62 }; 52 63 field_tolerances={... 53 64 1e-12,1e-12,1e-12,...
Note:
See TracChangeset
for help on using the changeset viewer.