Changeset 11190


Ignore:
Timestamp:
01/24/12 11:43:35 (13 years ago)
Author:
seroussi
Message:

fixing nightly runs

Location:
issm/trunk-jpl/test
Files:
6 added
4 edited

Legend:

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

    r11031 r11190  
    66
    77for 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
    1310        md=setmask(md,'',''); %ice sheet test
    1411        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
    1638
    1739        md=setflowequation(md,'stokes','all');
    1840
    19         %Create MPCs to have periodic boundary conditions
    20 
    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 times
    25         %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 
    2941        %Compute the diagnostic
    30         md.verbose=verbose('convergence',true,'solution',true);
    3142        md.cluster=generic('name',oshostname(),'np',8);
    3243        md=solve(md,DiagnosticSolutionEnum);
     
    4960        'Vx80km','Vy80km','Vz80km',...
    5061        'Vx160km','Vy160km','Vz160km'
    51 }
     62};
    5263field_tolerances={...
    5364        1e-12,1e-12,1e-12,...
Note: See TracChangeset for help on using the changeset viewer.