Changeset 8304


Ignore:
Timestamp:
05/16/11 15:49:49 (14 years ago)
Author:
seroussi
Message:

remove grid in nightly runs

Location:
issm/trunk/test/NightlyRun
Files:
64 edited

Legend:

Unmodified
Added
Removed

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

    r8295 r8304  
    2121
    2222        %Create dirichlet on the bed only
    23         md.spcvelocity=zeros(md.numberofgrids,6);
    24         pos=find(md.gridonbed);
     23        md.spcvelocity=zeros(md.numberofnodes,6);
     24        pos=find(md.nodeonbed);
    2525        md.spcvelocity(pos,1:2)=1;
    2626
    … …  
    2929        posx2=find(md.x==max(md.x));
    3030
    31         posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times
     31        posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same nodes two times
    3232        posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x));
    3333
    … …  
    4242        vz=PatchToVec(md.results.DiagnosticSolution.Vz);
    4343        results{i}=md.results.DiagnosticSolution;
    44         minvx(i)=min(vx(end-md.numberofgrids2d+1:end));
    45         maxvx(i)=max(vx(end-md.numberofgrids2d+1:end));
     44        minvx(i)=min(vx(end-md.numberofnodes2d+1:end));
     45        maxvx(i)=max(vx(end-md.numberofnodes2d+1:end));
    4646
    4747        %Now plot vx, vy, vz and vx on a cross section

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

    r8295 r8304  
    2828
    2929        %Create dirichlet on the bed only
    30         %md.spcvelocity=zeros(md.numberofgrids,6);
    31         pos=find(md.gridonbed);
     30        %md.spcvelocity=zeros(md.numberofnodes,6);
     31        pos=find(md.nodeonbed);
    3232        md.spcvelocity(pos,1:3)=1;
    3333
    … …  
    3535        %posx=find(md.x==0);
    3636        %posx2=find(md.x==max(md.x));
    37         %posx=find(md.x==0 & md.y~=0 & md.y~=max(md.y) & ~md.gridonbed);
    38         %posx2=find(md.x==max(md.x) &  md.y~=0 & md.y~=max(md.y) & ~md.gridonbed);
     37        %posx=find(md.x==0 & md.y~=0 & md.y~=max(md.y) & ~md.nodeonbed);
     38        %posx2=find(md.x==max(md.x) &  md.y~=0 & md.y~=max(md.y) & ~md.nodeonbed);
    3939
    40         %posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x) & ~md.gridonbed); %Don't take the same grids two times
    41         %posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x) & ~md.gridonbed);
     40        %posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x) & ~md.nodeonbed); %Don't take the same nodes two times
     41        %posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x) & ~md.nodeonbed);
    4242
    4343        %md.penalties=[posx,posx2;posy,posy2];
    … …  
    5555        pressure=PatchToVec(md.results.DiagnosticSolution.Pressure);
    5656        results{i}=md.results.DiagnosticSolution;
    57         minvx(i)=min(vx(end-md.numberofgrids2d+1:end));
    58         maxvx(i)=max(vx(end-md.numberofgrids2d+1:end));
     57        minvx(i)=min(vx(end-md.numberofnodes2d+1:end));
     58        maxvx(i)=max(vx(end-md.numberofnodes2d+1:end));
    5959
    6060        %Now plot vx, vy, vz and vx on a cross section

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

    r8295 r8304  
    1919
    2020        %Create dirichlet on the bed only
    21         md.spcvelocity=zeros(md.numberofgrids,6);
    22         pos=find(md.gridonbed);
     21        md.spcvelocity=zeros(md.numberofnodes,6);
     22        pos=find(md.nodeonbed);
    2323        md.spcvelocity(pos,1:2)=1;
    2424
    … …  
    2727        posx2=find(md.x==max(md.x));
    2828
    29         posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times
     29        posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same nodes two times
    3030        posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x));
    3131
    … …  
    4040        vz=PatchToVec(md.results.DiagnosticSolution.Vz);
    4141        results{i}=md.results.DiagnosticSolution;
    42         minvx(i)=min(vx(end-md.numberofgrids2d+1:end));
    43         maxvx(i)=max(vx(end-md.numberofgrids2d+1:end));
     42        minvx(i)=min(vx(end-md.numberofnodes2d+1:end));
     43        maxvx(i)=max(vx(end-md.numberofnodes2d+1:end));
    4444
    4545        %Now plot vx, vy, vz and vx on a cross section

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

    r8295 r8304  
    1818
    1919        %Create dirichlet on the bed only
    20         md.spcvelocity=zeros(md.numberofgrids,6);
    21         pos=find(md.gridonbed);
     20        md.spcvelocity=zeros(md.numberofnodes,6);
     21        pos=find(md.nodeonbed);
    2222        md.spcvelocity(pos,1:2)=1;
    2323
    … …  
    2626        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
     28        posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same nodes two times
    2929        posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x));
    3030

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

    r8295 r8304  
    2121
    2222        %Create MPCs to have periodic boundary conditions
    23         md.spcvelocity=zeros(md.numberofgrids,6);
     23        md.spcvelocity=zeros(md.numberofnodes,6);
    2424
    2525        posx=find(md.x==0 & md.y~=0 & md.y~=L);
    2626        posx2=find(md.x==L & md.y~=0 & md.y~=L);
    2727
    28         posy=find(md.y==0 & md.x~=0 & md.x~=L); %Don't take the same grids two times
     28        posy=find(md.y==0 & md.x~=0 & md.x~=L); %Don't take the same nodes two times
    2929        posy2=find(md.y==L & md.x~=0 & md.x~=L);
    3030
    … …  
    3232
    3333        %Add spc on the corners
    34         pos=find((md.x==0 | md.x==L) & (md.y==0 | md.y==L) & md.gridonbed);
     34        pos=find((md.x==0 | md.x==L) & (md.y==0 | md.y==L) & md.nodeonbed);
    3535        md.spcvelocity(pos,1:2)=1;
    3636        if(L==5000),
    … …  
    5555       
    5656        %Spc the bed at zero for vz
    57         pos=find(md.gridonbed);
     57        pos=find(md.nodeonbed);
    5858        md.spcvelocity(pos,3)=1;
    5959
    … …  
    6666        vz=PatchToVec(md.results.DiagnosticSolution.Vz);
    6767        results{i}=md.results.DiagnosticSolution;
    68         minvx(i)=min(vx(end-md.numberofgrids2d+1:end));
    69         maxvx(i)=max(vx(end-md.numberofgrids2d+1:end));
     68        minvx(i)=min(vx(end-md.numberofnodes2d+1:end));
     69        maxvx(i)=max(vx(end-md.numberofnodes2d+1:end));
    7070
    7171        %Now plot vx, vy, vz and vx on a cross section

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

    r8295 r8304  
    1818
    1919        %Create MPCs to have periodic boundary conditions
    20         %md.spcvelocity=zeros(md.numberofgrids,6);
     20        %md.spcvelocity=zeros(md.numberofnodes,6);
    2121
    2222        %posx=find(md.x==0);
    2323        %posx2=find(md.x==L);
    2424
    25         %posy=find(md.y==0 & md.x~=0 & md.x~=L); %Don't take the same grids two times
     25        %posy=find(md.y==0 & md.x~=0 & md.x~=L); %Don't take the same nodes two times
    2626        %posy2=find(md.y==L & md.x~=0 & md.x~=L);
    2727

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

    r8295 r8304  
    2121
    2222        %We need one grd on dirichlet: the 4 corners are set to zero
    23         md.spcvelocity=zeros(md.numberofgrids,6);
     23        md.spcvelocity=zeros(md.numberofnodes,6);
    2424
    2525        %Create MPCs to have periodic boundary conditions
    26         posx=find(md.x==0 & ~(md.y==0 & md.gridonbed) & ~(md.y==L & md.gridonbed));
    27         posx2=find(md.x==max(md.x) & ~(md.y==0 & md.gridonbed) & ~(md.y==L & md.gridonbed));
     26        posx=find(md.x==0 & ~(md.y==0 & md.nodeonbed) & ~(md.y==L & md.nodeonbed));
     27        posx2=find(md.x==max(md.x) & ~(md.y==0 & md.nodeonbed) & ~(md.y==L & md.nodeonbed));
    2828
    29         posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times
     29        posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same nodes two times
    3030        posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x));
    3131
    … …  
    3333
    3434        %Add spc on the corners
    35         pos=find((md.x==0 | md.x==L) & (md.y==0 | md.y==L) & md.gridonbed);
     35        pos=find((md.x==0 | md.x==L) & (md.y==0 | md.y==L) & md.nodeonbed);
    3636        md.spcvelocity(:,2)=1;
    3737        md.spcvelocity(pos,1:2)=1;
    … …  
    5151       
    5252        %Spc the bed at zero for vz
    53         pos=find(md.gridonbed);
     53        pos=find(md.nodeonbed);
    5454        md.spcvelocity(pos,3)=1;
    5555
    … …  
    6262        vz=PatchToVec(md.results.DiagnosticSolution.Vz);
    6363        results{i}=md.results.DiagnosticSolution;
    64         minvx(i)=min(vx(end-md.numberofgrids2d+1:end));
    65         maxvx(i)=max(vx(end-md.numberofgrids2d+1:end));
     64        minvx(i)=min(vx(end-md.numberofnodes2d+1:end));
     65        maxvx(i)=max(vx(end-md.numberofnodes2d+1:end));
    6666
    6767        %Now plot vx, vy, vz and vx on a cross section

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

    r8295 r8304  
    1818
    1919        %We need one grd on dirichlet: the 4 corners are set to zero
    20         %md.spcvelocity=zeros(md.numberofgrids,6);
     20        %md.spcvelocity=zeros(md.numberofnodes,6);
    2121        %pos=find((md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
    2222        %md.spcvelocity(pos,1:3)=1;
    … …  
    2626        %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
     28        %posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same nodes two times
    2929        %posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x));
    3030

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

    r8295 r8304  
    3131
    3232        %Create spcs on the bed
    33         pos=find(md.gridonbed);
    34         md.spcvelocity=zeros(md.numberofgrids,6);
     33        pos=find(md.nodeonbed);
     34        md.spcvelocity=zeros(md.numberofnodes,6);
    3535        md.spcvelocity(pos,1:3)=1;
    3636

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

    r8295 r8304  
    2323        if (i==1 | i==3),
    2424                %Create dirichlet on the bed if no slip
    25                 pos=find(md.gridonbed);
     25                pos=find(md.nodeonbed);
    2626                md.spcvelocity(pos,1:3)=1;
    2727        else
    28                 pos=find(md.gridonbed & (md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
     28                pos=find(md.nodeonbed & (md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
    2929                md.spcvelocity(pos,1:3)=1;
    3030                md.spcvelocity(pos,4)=100; %because we need a dirichlet somewhere
    3131
    3232        end
    33         pos=find(~md.gridonbed);
     33        pos=find(~md.nodeonbed);
    3434        md.spctemperature(pos,2)=255;
    3535        md.spctemperature(pos,1)=1;

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

    r8295 r8304  
    1212
    1313        disp('      initial velocity');
    14         md.vx=zeros(md.numberofgrids,1);
    15         md.vy=-400*ones(md.numberofgrids,1);
    16         md.vz=zeros(md.numberofgrids,1);
     14        md.vx=zeros(md.numberofnodes,1);
     15        md.vy=-400*ones(md.numberofnodes,1);
     16        md.vz=zeros(md.numberofnodes,1);
    1717
    1818        %analytical test

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

    r8295 r8304  
    2121md=setelementstype(md,'hutter','all');
    2222
    23 %Spc the grids on the bed
    24 pos=find(md.gridonbed);
     23%Spc the nodes on the bed
     24pos=find(md.nodeonbed);
    2525md.spcvelocity(pos,1:3)=1;
    2626
    … …  
    3131vx=PatchToVec(md.results.DiagnosticSolution.Vx);
    3232vy=PatchToVec(md.results.DiagnosticSolution.Vy);
    33 vel=zeros(md.numberofgrids2d,1);
     33vel=zeros(md.numberofnodes2d,1);
    3434
    35 grid_vel=0;
    36 for i=1:md.numberofgrids2d
     35node_vel=0;
     36for i=1:md.numberofnodes2d
    3737        for j=1:(md.numlayers-1)
    38                 grid_vel=grid_vel+1/(2*(md.numlayers-1))*(sqrt(vx(i+j*md.numberofgrids2d,1).^2+...
    39                         vy(i+j*md.numberofgrids2d,1).^2)+...
    40                         sqrt(vx(i+(j-1)*md.numberofgrids2d,1).^2+vy(i+(j-1)*md.numberofgrids2d,1).^2));
     38                node_vel=node_vel+1/(2*(md.numlayers-1))*(sqrt(vx(i+j*md.numberofnodes2d,1).^2+...
     39                        vy(i+j*md.numberofnodes2d,1).^2)+...
     40                        sqrt(vx(i+(j-1)*md.numberofnodes2d,1).^2+vy(i+(j-1)*md.numberofnodes2d,1).^2));
    4141        end
    42         vel(i,1)=grid_vel;
    43         grid_vel=0;
     42        vel(i,1)=node_vel;
     43        node_vel=0;
    4444end
    4545
    … …  
    6363subplot(2,2,3)
    6464hold on;
    65 plot(sqrt((md.x(1:md.numberofgrids2d)).^2+(md.y(1:md.numberofgrids2d)).^2),vel,'r.');
     65plot(sqrt((md.x(1:md.numberofnodes2d)).^2+(md.y(1:md.numberofnodes2d)).^2),vel,'r.');
    6666plot(sqrt((md.x2d).^2+(md.y2d).^2),vel_obs,'b.');
    6767title('Analytical vs calculated velocity','FontSize',14,'FontWeight','bold');

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

    r8295 r8304  
    2121md=setelementstype(md,'pattyn','all');
    2222
    23 %Spc the grids on the bed
    24 pos=find(md.gridonbed);
     23%Spc the nodes on the bed
     24pos=find(md.nodeonbed);
    2525md.spcvelocity(pos,1:3)=1;
    2626
    … …  
    3131vx=PatchToVec(md.results.DiagnosticSolution.Vx);
    3232vy=PatchToVec(md.results.DiagnosticSolution.Vy);
    33 vel=zeros(md.numberofgrids2d,1);
     33vel=zeros(md.numberofnodes2d,1);
    3434
    35 grid_vel=0;
    36 for i=1:md.numberofgrids2d
     35node_vel=0;
     36for i=1:md.numberofnodes2d
    3737        for j=1:(md.numlayers-1)
    38                 grid_vel=grid_vel+1/(2*(md.numlayers-1))*(sqrt(vx(i+j*md.numberofgrids2d,1).^2+...
    39                         vy(i+j*md.numberofgrids2d,1).^2)+...
    40                         sqrt(vx(i+(j-1)*md.numberofgrids2d,1).^2+vy(i+(j-1)*md.numberofgrids2d,1).^2));
     38                node_vel=node_vel+1/(2*(md.numlayers-1))*(sqrt(vx(i+j*md.numberofnodes2d,1).^2+...
     39                        vy(i+j*md.numberofnodes2d,1).^2)+...
     40                        sqrt(vx(i+(j-1)*md.numberofnodes2d,1).^2+vy(i+(j-1)*md.numberofnodes2d,1).^2));
    4141        end
    42         vel(i,1)=grid_vel;
    43         grid_vel=0;
     42        vel(i,1)=node_vel;
     43        node_vel=0;
    4444end
    4545
    … …  
    6262subplot(2,2,3)
    6363hold on;
    64 plot(sqrt((md.x(1:md.numberofgrids2d)).^2+(md.y(1:md.numberofgrids2d)).^2),vel,'r.');
     64plot(sqrt((md.x(1:md.numberofnodes2d)).^2+(md.y(1:md.numberofnodes2d)).^2),vel,'r.');
    6565plot(sqrt((md.x2d).^2+(md.y2d).^2),vel_obs,'b.');
    6666title('Analytical vs calculated velocity','FontSize',14,'FontWeight','bold');

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

    r8295 r8304  
    2121md=setelementstype(md,'stokes','all');
    2222
    23 %Spc the grids on the bed
    24 pos=find(md.gridonbed);
     23%Spc the nodes on the bed
     24pos=find(md.nodeonbed);
    2525md.spcvelocity(pos,1:3)=1;
    2626
    … …  
    3131vx=PatchToVec(md.results.DiagnosticSolution.Vx);
    3232vy=PatchToVec(md.results.DiagnosticSolution.Vy);
    33 vel=zeros(md.numberofgrids2d,1);
     33vel=zeros(md.numberofnodes2d,1);
    3434
    35 grid_vel=0;
    36 for i=1:md.numberofgrids2d
     35node_vel=0;
     36for i=1:md.numberofnodes2d
    3737        for j=1:(md.numlayers-1)
    38                 grid_vel=grid_vel+1/(2*(md.numlayers-1))*(sqrt(vx(i+j*md.numberofgrids2d,1).^2+...
    39                         vy(i+j*md.numberofgrids2d,1).^2)+...
    40                         sqrt(vx(i+(j-1)*md.numberofgrids2d,1).^2+vy(i+(j-1)*md.numberofgrids2d,1).^2));
     38                node_vel=node_vel+1/(2*(md.numlayers-1))*(sqrt(vx(i+j*md.numberofnodes2d,1).^2+...
     39                        vy(i+j*md.numberofnodes2d,1).^2)+...
     40                        sqrt(vx(i+(j-1)*md.numberofnodes2d,1).^2+vy(i+(j-1)*md.numberofnodes2d,1).^2));
    4141        end
    42         vel(i,1)=grid_vel;
    43         grid_vel=0;
     42        vel(i,1)=node_vel;
     43        node_vel=0;
    4444end
    4545
    … …  
    6262subplot(2,2,3)
    6363hold on;
    64 plot(sqrt((md.x(1:md.numberofgrids2d)).^2+(md.y(1:md.numberofgrids2d)).^2),vel,'r.');
     64plot(sqrt((md.x(1:md.numberofnodes2d)).^2+(md.y(1:md.numberofnodes2d)).^2),vel,'r.');
    6565plot(sqrt((md.x2d).^2+(md.y2d).^2),vel_obs,'b.');
    6666title('Analytical vs calculated velocity','FontSize',14,'FontWeight','bold');

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

    r8295 r8304  
    1212md=setelementstype(md,'hutter','all');
    1313
    14 %Spc the grids on the bed
    15 pos=find(md.gridonbed);
     14%Spc the nodes on the bed
     15pos=find(md.nodeonbed);
    1616md.spcvelocity(pos,1:3)=1;
    1717

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

    r8295 r8304  
    1313
    1414%Some conditions specific to melting test
    15 md.pressure=zeros(md.numberofgrids,1);
    16 md.observed_temperature=273.15*ones(md.numberofgrids,1);
     15md.pressure=zeros(md.numberofnodes,1);
     16md.observed_temperature=273.15*ones(md.numberofnodes,1);
    1717md.temperature=md.observed_temperature;
    1818md.spctemperature(:,2)=md.observed_temperature;
    … …  
    2828
    2929%plot results
    30 comp_melting=zeros(md.numberofgrids,1);
     30comp_melting=zeros(md.numberofnodes,1);
    3131comp_melting(md.results.ThermalSolution.MeltingRate.index)=md.results.ThermalSolution.MeltingRate.value;
    3232relative=abs((comp_melting-melting)./melting)*100;

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

    r8295 r8304  
    1414pos1=find(md.elementonbed);     md.spctemperature(md.elements(pos1,1:3),1)=1; md.spctemperature(md.elements(pos1,1:3),2)=10;
    1515pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6),1)=1; md.spctemperature(md.elements(pos2,4:6),2)=0;
    16 md.vz=0.1*ones(md.numberofgrids,1);
     16md.vz=0.1*ones(md.numberofnodes,1);
    1717md.vel=sqrt( md.vx.^2+ md.vy.^2+ md.vz.^2);
    18 md.pressure=zeros(md.numberofgrids,1);
     18md.pressure=zeros(md.numberofnodes,1);
    1919
    2020%analytical results
    … …  
    2929
    3030%plot results
    31 comp_temp=zeros(md.numberofgrids,1);
     31comp_temp=zeros(md.numberofnodes,1);
    3232comp_temp(md.results.ThermalSolution.Temperature.index)=md.results.ThermalSolution.Temperature.value;
    3333relative=abs((comp_temp-md.observed_temperature)./md.observed_temperature)*100;

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

    r8295 r8304  
    1313pos1=find(md.elementonbed);     md.spctemperature(md.elements(pos1,1:3),1)=1; md.spctemperature(md.elements(pos1,1:3),2)=10;
    1414pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6),1)=1; md.spctemperature(md.elements(pos2,4:6),2)=0;
    15 md.pressure=zeros(md.numberofgrids,1);
     15md.pressure=zeros(md.numberofnodes,1);
    1616
    1717%analytical results
    … …  
    2424
    2525%plot results
    26 comp_temp=zeros(md.numberofgrids,1);
     26comp_temp=zeros(md.numberofnodes,1);
    2727comp_temp(md.results.ThermalSolution.Temperature.index)=md.results.ThermalSolution.Temperature.value;
    2828relative=abs((comp_temp-md.observed_temperature)./md.observed_temperature)*100;

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

    r8295 r8304  
    1313
    1414pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6),1)=1; md.spctemperature(md.elements(pos2,4:6),2)=0;
    15 md.pressure=zeros(md.numberofgrids,1);
     15md.pressure=zeros(md.numberofnodes,1);
    1616md.geothermalflux(:)=0.1; %100mW/m^2
    1717
    … …  
    2525
    2626%plot results
    27 comp_temp=zeros(md.numberofgrids,1);
     27comp_temp=zeros(md.numberofnodes,1);
    2828comp_temp(md.results.ThermalSolution.Temperature.index)=md.results.ThermalSolution.Temperature.value;
    2929relative=abs((comp_temp-md.observed_temperature)./md.observed_temperature)*100;

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

    r8295 r8304  
    33md=parameterize(md,'../Par/SquareShelfConstrained.par');
    44%Add boundary conditions on thickness on the border
    5 pos=find(md.gridonboundary);
     5pos=find(md.nodeonboundary);
    66md.spcthickness(pos,1)=1;
    77md.spcthickness(pos,2)=md.thickness(pos);

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

    r8295 r8304  
    33md=parameterize(md,'../Par/SquareShelfConstrained.par');
    44%Add boundary conditions on thickness on the border
    5 pos=find(md.gridonboundary);
     5pos=find(md.nodeonboundary);
    66md.spcthickness(pos,1)=1;
    77md.spcthickness(pos,2)=md.thickness(pos);

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

    r8295 r8304  
    44md=extrude(md,5,1);
    55%Add boundary conditions on thickness on the border
    6 pos=find(md.gridonboundary);
     6pos=find(md.nodeonboundary);
    77md.spcthickness(pos,1)=1;
    88md.spcthickness(pos,2)=md.thickness(pos);

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

    r8295 r8304  
    44md=extrude(md,5,1);
    55%Add boundary conditions on thickness on the border
    6 pos=find(md.gridonboundary);
     6pos=find(md.nodeonboundary);
    77md.spcthickness(pos,1)=1;
    88md.spcthickness(pos,2)=md.thickness(pos);

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

    r8295 r8304  
    77md.control_analysis=1;
    88md.control_type=RheologyBbarEnum();
    9 md.cm_min=10^6*ones(md.numberofgrids,1);
    10 md.cm_max=2*10^9*ones(md.numberofgrids,1);
     9md.cm_min=10^6*ones(md.numberofnodes,1);
     10md.cm_max=2*10^9*ones(md.numberofnodes,1);
    1111md.nsteps=2;
    1212md.cm_responses=SurfaceAbsVelMisfitEnum*ones(md.nsteps,1);
    13 md.weights=ones(md.numberofgrids,1);
     13md.weights=ones(md.numberofnodes,1);
    1414md.optscal=10^7*ones(md.nsteps,1);
    1515md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    77md.control_analysis=1;
    88md.control_type=RheologyBbarEnum();
    9 md.cm_min=10^6*ones(md.numberofgrids,1);
    10 md.cm_max=2*10^9*ones(md.numberofgrids,1);
     9md.cm_min=10^6*ones(md.numberofnodes,1);
     10md.cm_max=2*10^9*ones(md.numberofnodes,1);
    1111md.nsteps=2;
    1212md.cm_responses=SurfaceAbsVelMisfitEnum*ones(md.nsteps,1);
    13 md.weights=ones(md.numberofgrids,1);
     13md.weights=ones(md.numberofnodes,1);
    1414md.optscal=10^7*ones(md.nsteps,1);
    1515md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=RheologyBbarEnum();
    10 md.cm_min=10^6*ones(md.numberofgrids,1);
    11 md.cm_max=2*10^9*ones(md.numberofgrids,1);
     10md.cm_min=10^6*ones(md.numberofnodes,1);
     11md.cm_max=2*10^9*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=10^7*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=RheologyBbarEnum();
    10 md.cm_min=10^6*ones(md.numberofgrids,1);
    11 md.cm_max=2*10^9*ones(md.numberofgrids,1);
     10md.cm_min=10^6*ones(md.numberofnodes,1);
     11md.cm_max=2*10^9*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=10^7*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=RheologyBbarEnum();
    10 md.cm_min=10^6*ones(md.numberofgrids,1);
    11 md.cm_max=2*10^9*ones(md.numberofgrids,1);
     10md.cm_min=10^6*ones(md.numberofnodes,1);
     11md.cm_max=2*10^9*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=10^7*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=RheologyBbarEnum();
    10 md.cm_min=10^6*ones(md.numberofgrids,1);
    11 md.cm_max=2*10^9*ones(md.numberofgrids,1);
     10md.cm_min=10^6*ones(md.numberofnodes,1);
     11md.cm_max=2*10^9*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=10^7*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=RheologyBbarEnum();
    10 md.cm_min=10^6*ones(md.numberofgrids,1);
    11 md.cm_max=2*10^9*ones(md.numberofgrids,1);
     10md.cm_min=10^6*ones(md.numberofnodes,1);
     11md.cm_max=2*10^9*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=10^7*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=RheologyBbarEnum();
    10 md.cm_min=10^6*ones(md.numberofgrids,1);
    11 md.cm_max=2*10^9*ones(md.numberofgrids,1);
     10md.cm_min=10^6*ones(md.numberofnodes,1);
     11md.cm_max=2*10^9*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=10^7*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    2121
    2222h=1000;
    23 md.thickness=h*ones(md.numberofgrids,1);
     23md.thickness=h*ones(md.numberofnodes,1);
    2424md.bed=-md.rho_ice/md.rho_water*md.thickness;
    2525md.surface=md.bed+md.thickness;
    2626
    2727%Initial velocity and pressure
    28 md.vx=zeros(md.numberofgrids,1);
    29 md.vy=zeros(md.numberofgrids,1);
    30 md.vz=zeros(md.numberofgrids,1);
    31 md.pressure=zeros(md.numberofgrids,1);
     28md.vx=zeros(md.numberofnodes,1);
     29md.vy=zeros(md.numberofnodes,1);
     30md.vz=zeros(md.numberofnodes,1);
     31md.pressure=zeros(md.numberofnodes,1);
    3232
    3333%Materials
    34 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
     34md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    3535md.rheology_B=paterson(md.observed_temperature);
    3636md.rheology_n=3*ones(md.numberofelements,1);
    … …  
    3838
    3939%Boundary conditions:
    40 md.spcvelocity=zeros(md.numberofgrids,6);
     40md.spcvelocity=zeros(md.numberofnodes,6);
    4141
    4242%constrain flanks to 0 normal velocity
    … …  
    5151
    5252%icefront
    53 gridonicefront=zeros(md.numberofgrids,1);
    54 pos=find(md.y==ymax); gridonicefront(pos)=1;
    55 pos=find(gridonicefront(md.segments(:,1)) | gridonicefront(md.segments(:,2))); pressureload=md.segments(pos,:);
     53nodeonicefront=zeros(md.numberofnodes,1);
     54pos=find(md.y==ymax); nodeonicefront(pos)=1;
     55pos=find(nodeonicefront(md.segments(:,1)) | nodeonicefront(md.segments(:,2))); pressureload=md.segments(pos,:);
    5656pressureload=[pressureload WaterEnum*md.elementoniceshelf(pressureload(:,end))];
    5757md.pressureload=pressureload;

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

    r8295 r8304  
    2121
    2222h=1000;
    23 md.thickness=h*ones(md.numberofgrids,1);
     23md.thickness=h*ones(md.numberofnodes,1);
    2424md.bed=-md.rho_ice/md.rho_water*md.thickness;
    2525md.surface=md.bed+md.thickness;
    2626
    2727%Initial velocity and pressure
    28 md.vx=zeros(md.numberofgrids,1);
    29 md.vy=zeros(md.numberofgrids,1);
    30 md.vz=zeros(md.numberofgrids,1);
    31 md.pressure=zeros(md.numberofgrids,1);
     28md.vx=zeros(md.numberofnodes,1);
     29md.vy=zeros(md.numberofnodes,1);
     30md.vz=zeros(md.numberofnodes,1);
     31md.pressure=zeros(md.numberofnodes,1);
    3232
    3333%Materials
    34 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
     34md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    3535md.rheology_B=paterson(md.observed_temperature);
    3636md.rheology_n=3*ones(md.numberofelements,1);
    … …  
    3838
    3939%Boundary conditions:
    40 md.spcvelocity=zeros(md.numberofgrids,6);
     40md.spcvelocity=zeros(md.numberofnodes,6);
    4141
    4242%constrain flanks to 0 normal velocity
    … …  
    5151
    5252%icefront
    53 gridonicefront=zeros(md.numberofgrids,1);
    54 pos=find(md.y==ymax); gridonicefront(pos)=1;
    55 pos=find(gridonicefront(md.segments(:,1)) | gridonicefront(md.segments(:,2))); pressureload=md.segments(pos,:);
     53nodeonicefront=zeros(md.numberofnodes,1);
     54pos=find(md.y==ymax); nodeonicefront(pos)=1;
     55pos=find(nodeonicefront(md.segments(:,1)) | nodeonicefront(md.segments(:,2))); pressureload=md.segments(pos,:);
    5656pressureload=[pressureload WaterEnum*md.elementoniceshelf(pressureload(:,end))];
    5757md.pressureload=pressureload;

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

    r8295 r8304  
    1818
    1919h=1000;
    20 md.thickness=h*ones(md.numberofgrids,1);
     20md.thickness=h*ones(md.numberofnodes,1);
    2121md.bed=-md.rho_ice/md.rho_water*md.thickness;
    2222md.surface=md.bed+md.thickness;
    2323
    2424%Initial velocity and pressure
    25 md.vx=zeros(md.numberofgrids,1);
    26 md.vy=zeros(md.numberofgrids,1);
    27 md.vz=zeros(md.numberofgrids,1);
    28 md.pressure=zeros(md.numberofgrids,1);
     25md.vx=zeros(md.numberofnodes,1);
     26md.vy=zeros(md.numberofnodes,1);
     27md.vz=zeros(md.numberofnodes,1);
     28md.pressure=zeros(md.numberofnodes,1);
    2929
    3030%Materials
    31 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
     31md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    3232md.rheology_B=paterson(md.observed_temperature);
    3333md.rheology_n=3*ones(md.numberofelements,1);
    … …  
    3535
    3636%Boundary conditions:
    37 md.spcvelocity=zeros(md.numberofgrids,6);
     37md.spcvelocity=zeros(md.numberofnodes,6);
    3838
    3939%constrain flanks to 0 normal velocity
    … …  
    4848
    4949%icefront
    50 gridonicefront=zeros(md.numberofgrids,1);
    51 pos=find(md.y==ymax); gridonicefront(pos)=1;
    52 pos=find(gridonicefront(md.segments(:,1)) | gridonicefront(md.segments(:,2))); pressureload=md.segments(pos,:);
     50nodeonicefront=zeros(md.numberofnodes,1);
     51pos=find(md.y==ymax); nodeonicefront(pos)=1;
     52pos=find(nodeonicefront(md.segments(:,1)) | nodeonicefront(md.segments(:,2))); pressureload=md.segments(pos,:);
    5353pressureload=[pressureload WaterEnum*md.elementoniceshelf(pressureload(:,end))];
    5454md.pressureload=pressureload;
    5555
    5656%partitioning
    57 md.npart=md.numberofgrids;
     57md.npart=md.numberofnodes;
    5858md=partitioner(md,'package','linear','npart',md.npart);
    5959md.part=md.part-1;

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

    r8295 r8304  
    1818
    1919h=1000;
    20 md.thickness=h*ones(md.numberofgrids,1);
     20md.thickness=h*ones(md.numberofnodes,1);
    2121md.bed=-md.rho_ice/md.rho_water*md.thickness;
    2222md.surface=md.bed+md.thickness;
    2323
    2424%Initial velocity and pressure
    25 md.vx=zeros(md.numberofgrids,1);
    26 md.vy=zeros(md.numberofgrids,1);
    27 md.vz=zeros(md.numberofgrids,1);
    28 md.pressure=zeros(md.numberofgrids,1);
     25md.vx=zeros(md.numberofnodes,1);
     26md.vy=zeros(md.numberofnodes,1);
     27md.vz=zeros(md.numberofnodes,1);
     28md.pressure=zeros(md.numberofnodes,1);
    2929
    3030%Materials
    31 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
     31md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
    3232md.rheology_B=paterson(md.observed_temperature);
    3333md.rheology_n=3*ones(md.numberofelements,1);
    … …  
    3535
    3636%Boundary conditions:
    37 md.spcvelocity=zeros(md.numberofgrids,6);
     37md.spcvelocity=zeros(md.numberofnodes,6);
    3838
    3939%constrain flanks to 0 normal velocity
    … …  
    4848
    4949%icefront
    50 gridonicefront=zeros(md.numberofgrids,1);
    51 pos=find(md.y==ymax); gridonicefront(pos)=1;
    52 pos=find(gridonicefront(md.segments(:,1)) | gridonicefront(md.segments(:,2))); pressureload=md.segments(pos,:);
     50nodeonicefront=zeros(md.numberofnodes,1);
     51pos=find(md.y==ymax); nodeonicefront(pos)=1;
     52pos=find(nodeonicefront(md.segments(:,1)) | nodeonicefront(md.segments(:,2))); pressureload=md.segments(pos,:);
    5353pressureload=[pressureload WaterEnum*md.elementoniceshelf(pressureload(:,end))];
    5454md.pressureload=pressureload;
    5555
    5656%partitioning
    57 md.npart=md.numberofgrids;
     57md.npart=md.numberofnodes;
    5858md=partitioner(md,'package','linear','npart',md.npart);
    5959md.part=md.part-1;

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

    r8295 r8304  
    77md.control_analysis=1;
    88md.control_type=DragCoefficientEnum();
    9 md.cm_min=1*ones(md.numberofgrids,1);
    10 md.cm_max=200*ones(md.numberofgrids,1);
     9md.cm_min=1*ones(md.numberofnodes,1);
     10md.cm_max=200*ones(md.numberofnodes,1);
    1111md.nsteps=2;
    1212md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1);
    13 md.weights=ones(md.numberofgrids,1);
     13md.weights=ones(md.numberofnodes,1);
    1414md.optscal=3*ones(md.nsteps,1);
    1515md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    77md.control_analysis=1;
    88md.control_type=DragCoefficientEnum();
    9 md.cm_min=1*ones(md.numberofgrids,1);
    10 md.cm_max=200*ones(md.numberofgrids,1);
     9md.cm_min=1*ones(md.numberofnodes,1);
     10md.cm_max=200*ones(md.numberofnodes,1);
    1111md.nsteps=2;
    1212md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1);
    13 md.weights=ones(md.numberofgrids,1);
     13md.weights=ones(md.numberofnodes,1);
    1414md.optscal=3*ones(md.nsteps,1);
    1515md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=DragCoefficientEnum();
    10 md.cm_min=1*ones(md.numberofgrids,1);
    11 md.cm_max=200*ones(md.numberofgrids,1);
     10md.cm_min=1*ones(md.numberofnodes,1);
     11md.cm_max=200*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=3*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=DragCoefficientEnum();
    10 md.cm_min=1*ones(md.numberofgrids,1);
    11 md.cm_max=200*ones(md.numberofgrids,1);
     10md.cm_min=1*ones(md.numberofnodes,1);
     11md.cm_max=200*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=3*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=DragCoefficientEnum();
    10 md.cm_min=1*ones(md.numberofgrids,1);
    11 md.cm_max=200*ones(md.numberofgrids,1);
     10md.cm_min=1*ones(md.numberofnodes,1);
     11md.cm_max=200*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceRelVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=3*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=DragCoefficientEnum();
    10 md.cm_min=1*ones(md.numberofgrids,1);
    11 md.cm_max=200*ones(md.numberofgrids,1);
     10md.cm_min=1*ones(md.numberofnodes,1);
     11md.cm_max=200*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceRelVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=3*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=DragCoefficientEnum();
    10 md.cm_min=1*ones(md.numberofgrids,1);
    11 md.cm_max=200*ones(md.numberofgrids,1);
     10md.cm_min=1*ones(md.numberofnodes,1);
     11md.cm_max=200*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceLogVxVyMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=3*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=DragCoefficientEnum();
    10 md.cm_min=1*ones(md.numberofgrids,1);
    11 md.cm_max=200*ones(md.numberofgrids,1);
     10md.cm_min=1*ones(md.numberofnodes,1);
     11md.cm_max=200*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceLogVxVyMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=3*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    66
    77%partitioning
    8 md.npart=md.numberofgrids;
     8md.npart=md.numberofnodes;
    99md=partitioner(md,'package','linear','npart',md.npart);
    1010md.part=md.part-1;

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

    r8295 r8304  
    66
    77%partitioning
    8 md.npart=md.numberofgrids;
     8md.npart=md.numberofnodes;
    99md=partitioner(md,'package','linear','npart',md.npart);
    1010md.part=md.part-1;

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=DragCoefficientEnum();
    10 md.cm_min=1*ones(md.numberofgrids,1);
    11 md.cm_max=200*ones(md.numberofgrids,1);
     10md.cm_min=1*ones(md.numberofnodes,1);
     11md.cm_max=200*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=3*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=DragCoefficientEnum();
    10 md.cm_min=1*ones(md.numberofgrids,1);
    11 md.cm_max=200*ones(md.numberofgrids,1);
     10md.cm_min=1*ones(md.numberofnodes,1);
     11md.cm_max=200*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=3*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=DragCoefficientEnum();
    10 md.cm_min=1*ones(md.numberofgrids,1);
    11 md.cm_max=200*ones(md.numberofgrids,1);
     10md.cm_min=1*ones(md.numberofnodes,1);
     11md.cm_max=200*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceRelVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=3*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=DragCoefficientEnum();
    10 md.cm_min=1*ones(md.numberofgrids,1);
    11 md.cm_max=200*ones(md.numberofgrids,1);
     10md.cm_min=1*ones(md.numberofnodes,1);
     11md.cm_max=200*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceRelVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=3*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r5955 r8304  
    1111md.part=md.part-1;
    1212
    13 vector=(1:1:md.numberofgrids)';
     13vector=(1:1:md.numberofnodes)';
    1414vector_on_partition=AreaAverageOntoPartition(md,vector);
    15 vector_on_grids=vector_on_partition(md.part+1);
     15vector_on_nodes=vector_on_partition(md.part+1);
    1616
    17 field_names     ={'vector_on_grids'};
     17field_names     ={'vector_on_nodes'};
    1818field_tolerances={1e-11};
    1919field_values={...
    20          vector_on_grids,...
     20         vector_on_nodes,...
    2121        };

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

    r8295 r8304  
    44%mesh {{{1
    55md=roundmesh(model,radius,resolution);
    6 %fix center grid to 0,0
     6%fix center node to 0,0
    77rad=sqrt((md.x).*md.x+(md.y).*md.y);
    88pos=find(rad==min(rad));

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

    r8295 r8304  
    1313md.control_analysis=1;
    1414md.control_type=RheologyBbarEnum();
    15 md.cm_min=10^6*ones(md.numberofgrids,1);
    16 md.cm_max=2*10^9*ones(md.numberofgrids,1);
     15md.cm_min=10^6*ones(md.numberofnodes,1);
     16md.cm_max=2*10^9*ones(md.numberofnodes,1);
    1717md.nsteps=2;
    1818md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1);
    19 md.weights=ones(md.numberofgrids,1);
     19md.weights=ones(md.numberofnodes,1);
    2020md.optscal=10^8*ones(md.nsteps,1);
    2121md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    1313md.control_analysis=1;
    1414md.control_type=RheologyBbarEnum();
    15 md.cm_min=10^6*ones(md.numberofgrids,1);
    16 md.cm_max=2*10^9*ones(md.numberofgrids,1);
     15md.cm_min=10^6*ones(md.numberofnodes,1);
     16md.cm_max=2*10^9*ones(md.numberofnodes,1);
    1717md.nsteps=2;
    1818md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1);
    19 md.weights=ones(md.numberofgrids,1);
     19md.weights=ones(md.numberofnodes,1);
    2020md.optscal=10^8*ones(md.nsteps,1);
    2121md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=DragCoefficientEnum();
    10 md.cm_min=1*ones(md.numberofgrids,1);
    11 md.cm_max=200*ones(md.numberofgrids,1);
     10md.cm_min=1*ones(md.numberofnodes,1);
     11md.cm_max=200*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=3*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=DragCoefficientEnum();
    10 md.cm_min=1*ones(md.numberofgrids,1);
    11 md.cm_max=200*ones(md.numberofgrids,1);
     10md.cm_min=1*ones(md.numberofnodes,1);
     11md.cm_max=200*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=3*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=DragCoefficientEnum();
    10 md.cm_min=1*ones(md.numberofgrids,1);
    11 md.cm_max=200*ones(md.numberofgrids,1);
     10md.cm_min=1*ones(md.numberofnodes,1);
     11md.cm_max=200*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=3*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r8295 r8304  
    88md.control_analysis=1;
    99md.control_type=DragCoefficientEnum();
    10 md.cm_min=1*ones(md.numberofgrids,1);
    11 md.cm_max=200*ones(md.numberofgrids,1);
     10md.cm_min=1*ones(md.numberofnodes,1);
     11md.cm_max=200*ones(md.numberofnodes,1);
    1212md.nsteps=2;
    1313md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1);
    14 md.weights=ones(md.numberofgrids,1);
     14md.weights=ones(md.numberofnodes,1);
    1515md.optscal=3*ones(md.nsteps,1);
    1616md.maxiter=2*ones(md.nsteps,1);

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

    r5098 r8304  
    2727
    2828%refine existing mesh 3
    29 hVertices=NaN*ones(md.numberofgrids,1);
    30 hVertices(find(md.gridonboundary))=500;
     29hVertices=NaN*ones(md.numberofnodes,1);
     30hVertices(find(md.nodeonboundary))=500;
    3131md2=bamg(md,'metric',md.dummy,'hmin',1000,'hmax',20000,'gradation',3,'geometricalmetric',1,'anisomax',1,'hVertices',hVertices);
    3232x5=md2.x;

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

    r8295 r8304  
    1010md.control_type=DhDtEnum;
    1111md.thickness_obs=md.thickness;
    12 md.cm_min=-50*ones(md.numberofgrids,1);
    13 md.cm_max=50*ones(md.numberofgrids,1);
     12md.cm_min=-50*ones(md.numberofnodes,1);
     13md.cm_max=50*ones(md.numberofnodes,1);
    1414md.cm_responses=ThicknessAbsMisfitEnum*ones(md.nsteps,1);
    15 md.weights=ones(md.numberofgrids,1);
     15md.weights=ones(md.numberofnodes,1);
    1616md.optscal=10/md.yts*ones(md.nsteps,1);
    1717md.maxiter=4*ones(md.nsteps,1);

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

    r8295 r8304  
    1010md.control_type=DhDtEnum;
    1111md.thickness_obs=md.thickness;
    12 md.cm_min=-50*ones(md.numberofgrids,1);
    13 md.cm_max=50*ones(md.numberofgrids,1);
     12md.cm_min=-50*ones(md.numberofnodes,1);
     13md.cm_max=50*ones(md.numberofnodes,1);
    1414md.cm_responses=ThicknessAbsMisfitEnum*ones(md.nsteps,1);
    15 md.weights=ones(md.numberofgrids,1);
     15md.weights=ones(md.numberofnodes,1);
    1616md.optscal=10/md.yts*ones(md.nsteps,1);
    1717md.maxiter=4*ones(md.nsteps,1);

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

    r8295 r8304  
    1111md.control_type=DhDtEnum;
    1212md.thickness_obs=md.thickness;
    13 md.cm_min=-50*ones(md.numberofgrids,1);
    14 md.cm_max=50*ones(md.numberofgrids,1);
     13md.cm_min=-50*ones(md.numberofnodes,1);
     14md.cm_max=50*ones(md.numberofnodes,1);
    1515md.cm_responses=ThicknessAbsMisfitEnum*ones(md.nsteps,1);
    16 md.weights=ones(md.numberofgrids,1);
     16md.weights=ones(md.numberofnodes,1);
    1717md.optscal=10/md.yts*ones(md.nsteps,1);
    1818md.maxiter=4*ones(md.nsteps,1);

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

    r8295 r8304  
    1111md.control_type=DhDtEnum;
    1212md.thickness_obs=md.thickness;
    13 md.cm_min=-50*ones(md.numberofgrids,1);
    14 md.cm_max=50*ones(md.numberofgrids,1);
     13md.cm_min=-50*ones(md.numberofnodes,1);
     14md.cm_max=50*ones(md.numberofnodes,1);
    1515md.cm_responses=ThicknessAbsMisfitEnum*ones(md.nsteps,1);
    16 md.weights=ones(md.numberofgrids,1);
     16md.weights=ones(md.numberofnodes,1);
    1717md.optscal=10/md.yts*ones(md.nsteps,1);
    1818md.maxiter=4*ones(md.nsteps,1);

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

    r8295 r8304  
    77%Ice sheet only
    88md=modelextract(md,md.elementonicesheet);
    9 pos=find(md.gridonboundary);
     9pos=find(md.nodeonboundary);
    1010md.spcthickness(pos,1)=1;
    1111md.spcthickness(:,2)=md.thickness;
    … …  
    1919md.prognostic_DG=0;
    2020md.optscal=[10/md.yts*ones(md.nsteps,1) 10/md.yts*ones(md.nsteps,1)];
    21 md.cm_min=[-2000*ones(md.numberofgrids,1) -2000*ones(md.numberofgrids,1)];
    22 md.cm_max=[+2000*ones(md.numberofgrids,1) +2000*ones(md.numberofgrids,1)];
     21md.cm_min=[-2000*ones(md.numberofnodes,1) -2000*ones(md.numberofnodes,1)];
     22md.cm_max=[+2000*ones(md.numberofnodes,1) +2000*ones(md.numberofnodes,1)];
    2323md.cm_responses=ThicknessAbsMisfitEnum*ones(md.nsteps,1);
    24 md.weights=ones(md.numberofgrids,1);
     24md.weights=ones(md.numberofnodes,1);
    2525md.maxiter=4*ones(md.nsteps,1);
    2626md.cm_jump=0.99*ones(md.nsteps,1);

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

    r8295 r8304  
    77%Ice sheet only
    88md=modelextract(md,md.elementonicesheet);
    9 pos=find(md.gridonboundary);
     9pos=find(md.nodeonboundary);
    1010md.spcthickness(pos,1)=1;
    1111md.spcthickness(:,2)=md.thickness;
    … …  
    1919md.prognostic_DG=0;
    2020md.optscal=[10/md.yts*ones(md.nsteps,1) 10/md.yts*ones(md.nsteps,1)];
    21 md.cm_min=[-2000*ones(md.numberofgrids,1) -2000*ones(md.numberofgrids,1)];
    22 md.cm_max=[+2000*ones(md.numberofgrids,1) +2000*ones(md.numberofgrids,1)];
     21md.cm_min=[-2000*ones(md.numberofnodes,1) -2000*ones(md.numberofnodes,1)];
     22md.cm_max=[+2000*ones(md.numberofnodes,1) +2000*ones(md.numberofnodes,1)];
    2323md.cm_responses=ThicknessAbsMisfitEnum*ones(md.nsteps,1);
    24 md.weights=ones(md.numberofgrids,1);
     24md.weights=ones(md.numberofnodes,1);
    2525md.maxiter=4*ones(md.nsteps,1);
    2626md.cm_jump=0.99*ones(md.nsteps,1);
Note: See TracChangeset for help on using the changeset viewer.