Changeset 8824


Ignore:
Timestamp:
07/06/11 14:30:08 (14 years ago)
Author:
Eric.Larour
Message:

Split spcvelocity into spcvx, spcvy and spcvz

Location:
issm/trunk/test
Files:
32 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/test/Miscellaneous/Bump/Bump.par

    r8307 r8824  
    2828%Deal with boundary conditions:
    2929md=SetIceSheetBC(md);
    30 md.spcvelocity(:,4)=zeros(md.numberofnodes,1);
    31 md.spcvelocity(:,5)=200*ones(md.numberofnodes,1);
     30md.spcvx=zeros(md.numberofnodes,1);
     31md.spcvy=200*ones(md.numberofnodes,1);
    3232
    3333%Parallel options

  • issm/trunk/test/Miscellaneous/Regularization/Square.par

    r8307 r8824  
    2525md=SetIceSheetBC(md);
    2626pos=find(md.gridonboundary);
    27 md.spcvelocity(pos,1)=1;
    28 md.spcvelocity(pos,2)=1;
    29 md.spcvelocity(pos,3)=1;
    30 md.spcvelocity(pos,4)=0;
    31 md.spcvelocity(pos,5)=200;
    32 md.spcvelocity(pos,6)=0;
     27md.spcvx(pos)=0;
     28md.spcvy(pos)=200;
     29md.spcvz(pos)=0;
    3330
    3431md.eps_res=10^-5;

  • issm/trunk/test/NightlyRun/test1101.m

    r8630 r8824  
    2121
    2222        %Create dirichlet on the bed only
    23         md.spcvelocity=zeros(md.numberofnodes,6);
     23        md.spcvx=NaN*ones(md.numberofnodes,1);
     24        md.spcvy=NaN*ones(md.numberofnodes,1);
     25        md.spcvz=NaN*ones(md.numberofnodes,1);
     26
    2427        pos=find(md.nodeonbed);
    25         md.spcvelocity(pos,1:2)=1;
     28        md.spcvx(pos)=0;
     29        md.spcvy(pos)=0;
    2630
    2731        %Create MPCs to have periodic boundary conditions

  • issm/trunk/test/NightlyRun/test1102.m

    r8630 r8824  
    2828
    2929        %Create dirichlet on the bed only
    30         %md.spcvelocity=zeros(md.numberofnodes,6);
    3130        pos=find(md.nodeonbed);
    32         md.spcvelocity(pos,1:3)=1;
     31        md.spcvx(pos)=0;
     32        md.spcvy(pos)=0;
     33        md.spcvz(pos)=0;
    3334
    3435        %Create MPCs to have periodic boundary conditions

  • issm/trunk/test/NightlyRun/test1103.m

    r8630 r8824  
    1919
    2020        %Create dirichlet on the bed only
    21         md.spcvelocity=zeros(md.numberofnodes,6);
     21        md.spcvx=NaN*ones(md.numberofnodes,1);
     22        md.spcvy=NaN*ones(md.numberofnodes,1);
     23        md.spcvz=NaN*ones(md.numberofnodes,1);
    2224        pos=find(md.nodeonbed);
    23         md.spcvelocity(pos,1:2)=1;
     25        md.spcvx(pos)=0;
     26        md.spcvy(pos)=0;
    2427
    2528        %Create MPCs to have periodic boundary conditions

  • issm/trunk/test/NightlyRun/test1104.m

    r8630 r8824  
    1818
    1919        %Create dirichlet on the bed only
    20         md.spcvelocity=zeros(md.numberofnodes,6);
     20        md.spcvx=NaN*ones(md.numberofnodes,1);
     21        md.spcvy=NaN*ones(md.numberofnodes,1);
     22        md.spcvz=NaN*ones(md.numberofnodes,1);
     23
    2124        pos=find(md.nodeonbed);
    22         md.spcvelocity(pos,1:2)=1;
     25        md.spcvx(pos)=0;
     26        md.spcvy(pos)=0;
    2327
    2428        %Create MPCs to have periodic boundary conditions

  • issm/trunk/test/NightlyRun/test1105.m

    r8630 r8824  
    2121
    2222        %Create MPCs to have periodic boundary conditions
    23         md.spcvelocity=zeros(md.numberofnodes,6);
     23        md.spcvx=NaN*ones(md.numberofnodes,1);
     24        md.spcvy=NaN*ones(md.numberofnodes,1);
     25        md.spcvz=NaN*ones(md.numberofnodes,1);
    2426
    2527        posx=find(md.x==0 & md.y~=0 & md.y~=L);
    … …  
    3335        %Add spc on the corners
    3436        pos=find((md.x==0 | md.x==L) & (md.y==0 | md.y==L) & md.nodeonbed);
    35         md.spcvelocity(pos,1:2)=1;
     37        md.spcvx(pos)=0;
     38        md.spcvy(pos)=0;
    3639        if(L==5000),
    37                 md.spcvelocity(pos,4)=15.66;
    38                 md.spcvelocity(pos,5)=-0.1967;
     40                md.spcvx(pos)=15.66;
     41                md.spcvy(pos)=-0.1967;
    3942        elseif(L==10000),
    40                 md.spcvelocity(pos,4)=16.04;
    41                 md.spcvelocity(pos,5)=-0.1977;
     43                md.spcvx(pos)=16.04;
     44                md.spcvy(pos)=-0.1977;
    4245        elseif(L==20000),
    43                 md.spcvelocity(pos,4)=16.53;
    44                 md.spcvelocity(pos,5)=-1.27;
     46                md.spcvx(pos)=16.53;
     47                md.spcvy(pos)=-1.27;
    4548        elseif(L==40000),
    46                 md.spcvelocity(pos,4)=17.23;
    47                 md.spcvelocity(pos,5)=-3.17;
     49                md.spcvx(pos)=17.23;
     50                md.spcvy(pos)=-3.17;
    4851        elseif(L==80000),
    49                 md.spcvelocity(pos,4)=16.68;
    50                 md.spcvelocity(pos,5)=-2.69;
     52                md.spcvx(pos)=16.68;
     53                md.spcvy(pos)=-2.69;
    5154        elseif(L==160000),
    52                 md.spcvelocity(pos,4)=16.03;
    53                 md.spcvelocity(pos,5)=-1.27;
     55                md.spcvx(pos)=16.03;
     56                md.spcvy(pos)=-1.27;
    5457        end
    5558       
    5659        %Spc the bed at zero for vz
    5760        pos=find(md.nodeonbed);
    58         md.spcvelocity(pos,3)=1;
     61        md.spcvz(pos)=0;
    5962
    6063        %Compute the diagnostic

  • issm/trunk/test/NightlyRun/test1106.m

    r8630 r8824  
    1818
    1919        %Create MPCs to have periodic boundary conditions
    20         %md.spcvelocity=zeros(md.numberofnodes,6);
    2120
    2221        %posx=find(md.x==0);

  • issm/trunk/test/NightlyRun/test1107.m

    r8630 r8824  
    2121
    2222        %We need one grd on dirichlet: the 4 corners are set to zero
    23         md.spcvelocity=zeros(md.numberofnodes,6);
     23        md.spcvx=NaN*ones(md.numberofnodes,1);
     24        md.spcvy=NaN*ones(md.numberofnodes,1);
     25        md.spcvz=NaN*ones(md.numberofnodes,1);
    2426
    2527        %Create MPCs to have periodic boundary conditions
    … …  
    3436        %Add spc on the corners
    3537        pos=find((md.x==0 | md.x==L) & (md.y==0 | md.y==L) & md.nodeonbed);
    36         md.spcvelocity(:,2)=1;
    37         md.spcvelocity(pos,1:2)=1;
     38        md.spcvy(:)=0;
     39        md.spcvx(pos)=0;
    3840        if(L==5000),
    39                 md.spcvelocity(pos,4)=16.0912;
     41                md.spcvx(pos)=16.0912;
    4042        elseif(L==10000),
    41                 md.spcvelocity(pos,4)=16.52;
     43                md.spcvx(pos)=16.52;
    4244        elseif(L==20000),
    43                 md.spcvelocity(pos,4)=17.77;
     45                md.spcvx(pos)=17.77;
    4446        elseif(L==40000),
    45                 md.spcvelocity(pos,4)=19.88;
     47                md.spcvx(pos)=19.88;
    4648        elseif(L==80000),
    47                 md.spcvelocity(pos,4)=18.65;
     49                md.spcvx(pos)=18.65;
    4850        elseif(L==160000),
    49                 md.spcvelocity(pos,4)=16.91;
     51                md.spcvx(pos)=16.91;
    5052        end
    5153       
    5254        %Spc the bed at zero for vz
    5355        pos=find(md.nodeonbed);
    54         md.spcvelocity(pos,3)=1;
     56        md.spcvz(pos,3)=0;
    5557
    5658        %Compute the diagnostic

  • issm/trunk/test/NightlyRun/test1108.m

    r8630 r8824  
    1818
    1919        %We need one grd on dirichlet: the 4 corners are set to zero
    20         %md.spcvelocity=zeros(md.numberofnodes,6);
     20        %md.spcvx=NaN*ones(md.numberofnodes,1);
     21        %md.spcvy=NaN*ones(md.numberofnodes,1);
     22        %md.spcvz=NaN*ones(md.numberofnodes,1);
     23       
    2124        %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;
     25        %md.spcvx(pos)=0;
     26        %md.spcvy(pos)=0;
     27        %md.spcvz(pos)=0;
    2328
    2429        %%Create MPCs to have periodic boundary conditions

  • issm/trunk/test/NightlyRun/test1109.m

    r8665 r8824  
    3232        %Create spcs on the bed
    3333        pos=find(md.nodeonbed);
    34         md.spcvelocity=zeros(md.numberofnodes,6);
    35         md.spcvelocity(pos,1:3)=1;
     34        md.spcvx=NaN*ones(md.numberofnodes,1);
     35        md.spcvy=NaN*ones(md.numberofnodes,1);
     36        md.spcvz=NaN*ones(md.numberofnodes,1);
     37        md.spcvx(pos)=0;
     38        md.spcvy(pos)=0;
     39        md.spcvz(pos)=0;
    3640
    3741        %Remove the spc where there is some sliding (case 3 and 4):
    3842        if i==3 | i==4,
    3943                pos=find(md.y/max(md.y)>=0.44 & md.y/max(md.y)<=0.5);
    40                 md.spcvelocity(pos,1:3)=0;
     44                md.spcvx(pos)=0;
     45                md.spcvy(pos)=0;
     46                md.spcvz(pos)=0;
    4147        end
    4248

  • issm/trunk/test/NightlyRun/test1110.m

    r8589 r8824  
    2020        end
    2121
    22         md.spcvelocity(:,:)=0;
     22        md.spcvx=NaN*ones(md.numberofnodes,1);
     23        md.spcvy=NaN*ones(md.numberofnodes,1);
     24        md.spcvz=NaN*ones(md.numberofnodes,1);
    2325        if (i==1 | i==3),
    2426                %Create dirichlet on the bed if no slip
    2527                pos=find(md.nodeonbed);
    26                 md.spcvelocity(pos,1:3)=1;
     28                md.spcvx(pos)=0;
     29                md.spcvy(pos)=0;
     30                md.spcvz(pos)=0;
    2731        else
    2832                pos=find(md.nodeonbed & (md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
    29                 md.spcvelocity(pos,1:3)=1;
    30                 md.spcvelocity(pos,4)=100; %because we need a dirichlet somewhere
    31 
     33                md.spcvx(pos)=100; %because we need a dirichlet somewhere
     34                md.spcvy(pos)=0;
     35                md.spcvz(pos)=0;
    3236        end
    3337        pos=find(~md.nodeonbed);

  • issm/trunk/test/NightlyRun/test120.m

    r8589 r8824  
    44md=setelementstype(md,'macayeal','all');
    55md.cluster=generic('name',oshostname(),'np',3);
     6
    67md=solve(md,TransientSolutionEnum);
    78

  • issm/trunk/test/NightlyRun/test1201.m

    r8630 r8824  
    1818
    1919        %analytical test
    20         md.spcvelocity(:,1:3)=1;
    21         md.spcvelocity(:,4)=0;
    22         md.spcvelocity(:,5)=-400;
    23         md.spcvelocity(1,1)=0;
     20        md.spcvx=zeros(md.numberofnodes,1);
     21        md.spcvy=zeros(md.numberofnodes,1);
     22        md.spcvz=zeros(md.numberofnodes,1);
     23        md.spcvx(:)=0;
     24        md.spcvy(:)=-400;
     25        md.spcvx(1)=0;
    2426        if type==1,
    2527                md.artificial_diffusivity=1;

  • issm/trunk/test/NightlyRun/test1203.m

    r8630 r8824  
    1111%Impose a non zero velocity on the upper boundary condition (y=max(y))
    1212pos=find(md.y==max(md.y));
    13 md.spcvelocity(pos,5)=400*(((md.x(pos)-100000)/25000).^2-ones(size(pos,1),1)).*heaviside((1+eps)*ones(size(pos,1),1)-((md.x(pos)-100000)/25000).^2);
     13md.spcvy(pos)=400*(((md.x(pos)-100000)/25000).^2-ones(size(pos,1),1)).*heaviside((1+eps)*ones(size(pos,1),1)-((md.x(pos)-100000)/25000).^2);
    1414
    1515%Compute solution for MacAyeal's model

  • issm/trunk/test/NightlyRun/test1204.m

    r8630 r8824  
    1111%Impose a non zero velocity on the upper boundary condition (y=max(y))
    1212pos=find(md.y==max(md.y));
    13 md.spcvelocity(pos,5)=400*(((md.x(pos)-100000)/25000).^2-ones(size(pos,1),1)).*heaviside((1+eps)*ones(size(pos,1),1)-((md.x(pos)-100000)/25000).^2);
     13md.spcvy(pos)=400*(((md.x(pos)-100000)/25000).^2-ones(size(pos,1),1)).*heaviside((1+eps)*ones(size(pos,1),1)-((md.x(pos)-100000)/25000).^2);
    1414
    1515%Compute solution for MacAyeal's model

  • issm/trunk/test/NightlyRun/test1205.m

    r8630 r8824  
    2323%Spc the nodes on the bed
    2424pos=find(md.nodeonbed);
    25 md.spcvelocity(pos,1:3)=1;
     25md.spcvx(pos)=0;
     26md.spcvy(pos)=0;
     27md.spcvz(pos)=0;
    2628
    2729%Now we can solve the problem

  • issm/trunk/test/NightlyRun/test1206.m

    r8630 r8824  
    2323%Spc the nodes on the bed
    2424pos=find(md.nodeonbed);
    25 md.spcvelocity(pos,1:3)=1;
     25md.spcvx(pos)=0;
     26md.spcvy(pos)=0;
     27md.spcvz(pos)=0;
    2628
    2729%Now we can solve the problem

  • issm/trunk/test/NightlyRun/test1207.m

    r8630 r8824  
    2323%Spc the nodes on the bed
    2424pos=find(md.nodeonbed);
    25 md.spcvelocity(pos,1:3)=1;
     25md.spcvx(pos)=0;
     26md.spcvy(pos)=0;
     27md.spcvz(pos)=0;
    2628
    2729%Now we can solve the problem

  • issm/trunk/test/NightlyRun/test1208.m

    r8630 r8824  
    1414%Spc the nodes on the bed
    1515pos=find(md.nodeonbed);
    16 md.spcvelocity(pos,1:3)=1;
     16md.spcvx(pos)=0;
     17md.spcvy(pos)=0;
     18md.spcvz(pos)=0;
    1719
    1820%Adapt the time steps to the resolution

  • issm/trunk/test/NightlyRun/test233.m

    r8304 r8824  
    3838
    3939%Boundary conditions:
    40 md.spcvelocity=zeros(md.numberofnodes,6);
     40md.spcvx=NaN*ones(md.numberofnodes,1);
     41md.spcvy=NaN*ones(md.numberofnodes,1);
     42md.spcvz=NaN*ones(md.numberofnodes,1);
    4143
    4244%constrain flanks to 0 normal velocity
    4345pos=find(md.x==xmin | md.x==xmax);
    44 md.spcvelocity(pos,1)=1;
    45 md.spcvelocity(pos,3)=0;
     46md.spcvx(pos)=0;
     47md.spcvz(pos)=NaN;
    4648
    4749%constrain grounding line to 0 velocity
    4850pos=find(md.y==ymin);
    49 md.spcvelocity(pos,1:2)=1;
    50 md.spcvelocity(pos,3:4)=0;
     51md.spcvx(pos)=0;
     52md.spcvy(pos)=0;
    5153
    5254%icefront

  • issm/trunk/test/NightlyRun/test234.m

    r8589 r8824  
    3838
    3939%Boundary conditions:
    40 md.spcvelocity=zeros(md.numberofnodes,6);
     40md.spcvx=NaN*ones(md.numberofnodes,1);
     41md.spcvy=NaN*ones(md.numberofnodes,1);
     42md.spcvz=NaN*ones(md.numberofnodes,1);
    4143
    4244%constrain flanks to 0 normal velocity
    4345pos=find(md.x==xmin | md.x==xmax);
    44 md.spcvelocity(pos,1)=1;
    45 md.spcvelocity(pos,3)=0;
     46md.spcvx(pos)=0;
     47md.spcvz(pos)=NaN;
    4648
    4749%constrain grounding line to 0 velocity
    4850pos=find(md.y==ymin);
    49 md.spcvelocity(pos,1:2)=1;
    50 md.spcvelocity(pos,3:4)=0;
     51md.spcvx(pos)=0;
     52md.spcvy(pos)=0;
    5153
    5254%icefront

  • issm/trunk/test/NightlyRun/test235.m

    r8304 r8824  
    3535
    3636%Boundary conditions:
    37 md.spcvelocity=zeros(md.numberofnodes,6);
     37md.spcvx=NaN*ones(md.numberofnodes,1);
     38md.spcvy=NaN*ones(md.numberofnodes,1);
     39md.spcvz=NaN*ones(md.numberofnodes,1);
    3840
    3941%constrain flanks to 0 normal velocity
    4042pos=find(md.x==xmin | md.x==xmax);
    41 md.spcvelocity(pos,1)=1;
    42 md.spcvelocity(pos,3)=0;
     43md.spcvx(pos)=0;
     44md.spcvz(pos)=NaN;
    4345
    4446%constrain grounding line to 0 velocity
    4547pos=find(md.y==ymin);
    46 md.spcvelocity(pos,1:2)=1;
    47 md.spcvelocity(pos,3:4)=0;
     48md.spcvx(pos)=0;
     49md.spcvy(pos)=0;
    4850
    4951%icefront

  • issm/trunk/test/NightlyRun/test236.m

    r8589 r8824  
    3535
    3636%Boundary conditions:
    37 md.spcvelocity=zeros(md.numberofnodes,6);
     37md.spcvx=NaN*ones(md.numberofnodes,1);
     38md.spcvy=NaN*ones(md.numberofnodes,1);
     39md.spcvz=NaN*ones(md.numberofnodes,1);
    3840
    3941%constrain flanks to 0 normal velocity
    4042pos=find(md.x==xmin | md.x==xmax);
    41 md.spcvelocity(pos,1)=1;
    42 md.spcvelocity(pos,3)=0;
     43md.spcvx(pos)=0;
     44md.spcvz(pos)=NaN;
    4345
    4446%constrain grounding line to 0 velocity
    4547pos=find(md.y==ymin);
    46 md.spcvelocity(pos,1:2)=1;
    47 md.spcvelocity(pos,3:4)=0;
     48md.spcvx(pos)=0;
     49md.spcvy(pos)=0;
    4850
    4951%icefront

  • issm/trunk/test/NightlyRun/test427.m

    r8295 r8824  
    88
    99%constrain all velocities to 1 m/yr, in the y-direction
    10 md.spcvelocity(:,1:3)=1;
    11 md.spcvelocity(:,4)=0;
    12 md.spcvelocity(:,5)=1;
    13 md.spcvelocity(:,6)=0;
     10md.spcvx(:)=0;
     11md.spcvy(:)=1;
     12md.spcvz(:)=0;
    1413
    1514%Dakota options

  • issm/trunk/test/NightlyRun/test428.m

    r8589 r8824  
    88
    99%constrain all velocities to 1 m/yr, in the y-direction
    10 md.spcvelocity(:,1:3)=1;
    11 md.spcvelocity(:,4)=0;
    12 md.spcvelocity(:,5)=1;
    13 md.spcvelocity(:,6)=0;
     10md.spcvx=zeros(md.numberofnodes,1);
     11md.spcvy=ones(md.numberofnodes,1);
     12md.spcvz=zeros(md.numberofnodes,1);
    1413
    1514%Dakota options

  • issm/trunk/test/NightlyRun/test434.m

    r8589 r8824  
    88
    99%constrain all velocities to 1 m/yr, in the y-direction
    10 md.spcvelocity(:,1:3)=1;
    11 md.spcvelocity(:,4)=0;
    12 md.spcvelocity(:,5)=1;
    13 md.spcvelocity(:,6)=0;
     10md.spcvx(:)=0;
     11md.spcvy(:)=1;
     12md.spcvz(:)=0;
    1413
    1514%Dakota options

  • issm/trunk/test/Par/ISMIPF.par

    r8631 r8824  
    2323%Create node on boundary fist (because we cannot use mesh)
    2424md=SetIceSheetBC(md);
    25 md.spcvelocity(:,4)=100*ones(md.numberofnodes,1);
     25md.spcvx=100*ones(md.numberofnodes,1);
    2626md.vx=zeros(md.numberofnodes,1);
    2727md.vy=zeros(md.numberofnodes,1);

  • issm/trunk/test/Par/Pig.par

    r8305 r8824  
    22
    33%Geometry and observation
    4 load('./../Data/Pig.data','-mat');
     4load('../Data/Pig.data','-mat');
    55md.vx_obs   =InterpFromMeshToMesh2d(index,x,y,vx_obs,md.x,md.y);
    66md.vy_obs   =InterpFromMeshToMesh2d(index,x,y,vy_obs,md.x,md.y);

  • issm/trunk/test/Par/RoundSheetEISMINT.par

    r8631 r8824  
    4848md.x(pos)=0; md.y(pos)=0; %the closest node to the center is changed to be exactly at the center
    4949
    50 md.spcvelocity(pos,1:3)=1;
    51 md.spcvelocity(pos,4:6)=0;
     50md.spcvx(pos)=0;
     51md.spcvy(pos)=0;
     52md.spcvz(pos)=0;
    5253
    5354%parallel options

  • issm/trunk/test/Par/RoundSheetShelf.par

    r8399 r8824  
    7676
    7777%Deal with boundary conditions:
    78 md.spcvelocity=zeros(md.numberofnodes,6);
     78md.spcvx=NaN*ones(md.numberofnodes,1);
     79md.spcvy=NaN*ones(md.numberofnodes,1);
     80md.spcvz=NaN*ones(md.numberofnodes,1);
     81
    7982pos=find(md.x==0 & md.y==0);
    80 md.spcvelocity(pos,1:2)=1;
     83md.spcvx(pos)=0;
     84md.spcvy(pos)=0;
     85
    8186pressureload=md.segments;
    8287pressureload=[pressureload WaterEnum*md.elementoniceshelf(pressureload(:,end))+AirEnum*md.elementonicesheet(pressureload(:,end))];

  • issm/trunk/test/Par/RoundSheetStaticEISMINT.par

    r8631 r8824  
    5454md.x(pos)=0; md.y(pos)=0; %the closest node to the center is changed to be exactly at the center
    5555
    56 md.spcvelocity(pos,1:3)=1;
    57 md.spcvelocity(pos,4:6)=0;
     56md.spcvx(pos)=0;
     57md.spcvy(pos)=0;
     58md.spcvz(pos)=0;
Note: See TracChangeset for help on using the changeset viewer.