Changeset 8304
- Timestamp:
- 05/16/11 15:49:49 (14 years ago)
- Location:
- issm/trunk/test/NightlyRun
- Files:
-
- 64 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/NightlyRun/test1101.m
r8295 r8304 21 21 22 22 %Create dirichlet on the bed only 23 md.spcvelocity=zeros(md.numberof grids,6);24 pos=find(md. gridonbed);23 md.spcvelocity=zeros(md.numberofnodes,6); 24 pos=find(md.nodeonbed); 25 25 md.spcvelocity(pos,1:2)=1; 26 26 … … 29 29 posx2=find(md.x==max(md.x)); 30 30 31 posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times31 posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same nodes two times 32 32 posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x)); 33 33 … … 42 42 vz=PatchToVec(md.results.DiagnosticSolution.Vz); 43 43 results{i}=md.results.DiagnosticSolution; 44 minvx(i)=min(vx(end-md.numberof grids2d+1:end));45 maxvx(i)=max(vx(end-md.numberof grids2d+1:end));44 minvx(i)=min(vx(end-md.numberofnodes2d+1:end)); 45 maxvx(i)=max(vx(end-md.numberofnodes2d+1:end)); 46 46 47 47 %Now plot vx, vy, vz and vx on a cross section -
issm/trunk/test/NightlyRun/test1102.m
r8295 r8304 28 28 29 29 %Create dirichlet on the bed only 30 %md.spcvelocity=zeros(md.numberof grids,6);31 pos=find(md. gridonbed);30 %md.spcvelocity=zeros(md.numberofnodes,6); 31 pos=find(md.nodeonbed); 32 32 md.spcvelocity(pos,1:3)=1; 33 33 … … 35 35 %posx=find(md.x==0); 36 36 %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); 39 39 40 %posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x) & ~md. gridonbed); %Don't take the same grids two times41 %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); 42 42 43 43 %md.penalties=[posx,posx2;posy,posy2]; … … 55 55 pressure=PatchToVec(md.results.DiagnosticSolution.Pressure); 56 56 results{i}=md.results.DiagnosticSolution; 57 minvx(i)=min(vx(end-md.numberof grids2d+1:end));58 maxvx(i)=max(vx(end-md.numberof grids2d+1:end));57 minvx(i)=min(vx(end-md.numberofnodes2d+1:end)); 58 maxvx(i)=max(vx(end-md.numberofnodes2d+1:end)); 59 59 60 60 %Now plot vx, vy, vz and vx on a cross section -
issm/trunk/test/NightlyRun/test1103.m
r8295 r8304 19 19 20 20 %Create dirichlet on the bed only 21 md.spcvelocity=zeros(md.numberof grids,6);22 pos=find(md. gridonbed);21 md.spcvelocity=zeros(md.numberofnodes,6); 22 pos=find(md.nodeonbed); 23 23 md.spcvelocity(pos,1:2)=1; 24 24 … … 27 27 posx2=find(md.x==max(md.x)); 28 28 29 posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times29 posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same nodes two times 30 30 posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x)); 31 31 … … 40 40 vz=PatchToVec(md.results.DiagnosticSolution.Vz); 41 41 results{i}=md.results.DiagnosticSolution; 42 minvx(i)=min(vx(end-md.numberof grids2d+1:end));43 maxvx(i)=max(vx(end-md.numberof grids2d+1:end));42 minvx(i)=min(vx(end-md.numberofnodes2d+1:end)); 43 maxvx(i)=max(vx(end-md.numberofnodes2d+1:end)); 44 44 45 45 %Now plot vx, vy, vz and vx on a cross section -
issm/trunk/test/NightlyRun/test1104.m
r8295 r8304 18 18 19 19 %Create dirichlet on the bed only 20 md.spcvelocity=zeros(md.numberof grids,6);21 pos=find(md. gridonbed);20 md.spcvelocity=zeros(md.numberofnodes,6); 21 pos=find(md.nodeonbed); 22 22 md.spcvelocity(pos,1:2)=1; 23 23 … … 26 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 times28 posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same nodes two times 29 29 posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x)); 30 30 -
issm/trunk/test/NightlyRun/test1105.m
r8295 r8304 21 21 22 22 %Create MPCs to have periodic boundary conditions 23 md.spcvelocity=zeros(md.numberof grids,6);23 md.spcvelocity=zeros(md.numberofnodes,6); 24 24 25 25 posx=find(md.x==0 & md.y~=0 & md.y~=L); 26 26 posx2=find(md.x==L & md.y~=0 & md.y~=L); 27 27 28 posy=find(md.y==0 & md.x~=0 & md.x~=L); %Don't take the same grids two times28 posy=find(md.y==0 & md.x~=0 & md.x~=L); %Don't take the same nodes two times 29 29 posy2=find(md.y==L & md.x~=0 & md.x~=L); 30 30 … … 32 32 33 33 %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); 35 35 md.spcvelocity(pos,1:2)=1; 36 36 if(L==5000), … … 55 55 56 56 %Spc the bed at zero for vz 57 pos=find(md. gridonbed);57 pos=find(md.nodeonbed); 58 58 md.spcvelocity(pos,3)=1; 59 59 … … 66 66 vz=PatchToVec(md.results.DiagnosticSolution.Vz); 67 67 results{i}=md.results.DiagnosticSolution; 68 minvx(i)=min(vx(end-md.numberof grids2d+1:end));69 maxvx(i)=max(vx(end-md.numberof grids2d+1:end));68 minvx(i)=min(vx(end-md.numberofnodes2d+1:end)); 69 maxvx(i)=max(vx(end-md.numberofnodes2d+1:end)); 70 70 71 71 %Now plot vx, vy, vz and vx on a cross section -
issm/trunk/test/NightlyRun/test1106.m
r8295 r8304 18 18 19 19 %Create MPCs to have periodic boundary conditions 20 %md.spcvelocity=zeros(md.numberof grids,6);20 %md.spcvelocity=zeros(md.numberofnodes,6); 21 21 22 22 %posx=find(md.x==0); 23 23 %posx2=find(md.x==L); 24 24 25 %posy=find(md.y==0 & md.x~=0 & md.x~=L); %Don't take the same grids two times25 %posy=find(md.y==0 & md.x~=0 & md.x~=L); %Don't take the same nodes two times 26 26 %posy2=find(md.y==L & md.x~=0 & md.x~=L); 27 27 -
issm/trunk/test/NightlyRun/test1107.m
r8295 r8304 21 21 22 22 %We need one grd on dirichlet: the 4 corners are set to zero 23 md.spcvelocity=zeros(md.numberof grids,6);23 md.spcvelocity=zeros(md.numberofnodes,6); 24 24 25 25 %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)); 28 28 29 posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times29 posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same nodes two times 30 30 posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x)); 31 31 … … 33 33 34 34 %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); 36 36 md.spcvelocity(:,2)=1; 37 37 md.spcvelocity(pos,1:2)=1; … … 51 51 52 52 %Spc the bed at zero for vz 53 pos=find(md. gridonbed);53 pos=find(md.nodeonbed); 54 54 md.spcvelocity(pos,3)=1; 55 55 … … 62 62 vz=PatchToVec(md.results.DiagnosticSolution.Vz); 63 63 results{i}=md.results.DiagnosticSolution; 64 minvx(i)=min(vx(end-md.numberof grids2d+1:end));65 maxvx(i)=max(vx(end-md.numberof grids2d+1:end));64 minvx(i)=min(vx(end-md.numberofnodes2d+1:end)); 65 maxvx(i)=max(vx(end-md.numberofnodes2d+1:end)); 66 66 67 67 %Now plot vx, vy, vz and vx on a cross section -
issm/trunk/test/NightlyRun/test1108.m
r8295 r8304 18 18 19 19 %We need one grd on dirichlet: the 4 corners are set to zero 20 %md.spcvelocity=zeros(md.numberof grids,6);20 %md.spcvelocity=zeros(md.numberofnodes,6); 21 21 %pos=find((md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y))); 22 22 %md.spcvelocity(pos,1:3)=1; … … 26 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 times28 %posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same nodes two times 29 29 %posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x)); 30 30 -
issm/trunk/test/NightlyRun/test1109.m
r8295 r8304 31 31 32 32 %Create spcs on the bed 33 pos=find(md. gridonbed);34 md.spcvelocity=zeros(md.numberof grids,6);33 pos=find(md.nodeonbed); 34 md.spcvelocity=zeros(md.numberofnodes,6); 35 35 md.spcvelocity(pos,1:3)=1; 36 36 -
issm/trunk/test/NightlyRun/test1110.m
r8295 r8304 23 23 if (i==1 | i==3), 24 24 %Create dirichlet on the bed if no slip 25 pos=find(md. gridonbed);25 pos=find(md.nodeonbed); 26 26 md.spcvelocity(pos,1:3)=1; 27 27 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))); 29 29 md.spcvelocity(pos,1:3)=1; 30 30 md.spcvelocity(pos,4)=100; %because we need a dirichlet somewhere 31 31 32 32 end 33 pos=find(~md. gridonbed);33 pos=find(~md.nodeonbed); 34 34 md.spctemperature(pos,2)=255; 35 35 md.spctemperature(pos,1)=1; -
issm/trunk/test/NightlyRun/test1201.m
r8295 r8304 12 12 13 13 disp(' initial velocity'); 14 md.vx=zeros(md.numberof grids,1);15 md.vy=-400*ones(md.numberof grids,1);16 md.vz=zeros(md.numberof grids,1);14 md.vx=zeros(md.numberofnodes,1); 15 md.vy=-400*ones(md.numberofnodes,1); 16 md.vz=zeros(md.numberofnodes,1); 17 17 18 18 %analytical test -
issm/trunk/test/NightlyRun/test1205.m
r8295 r8304 21 21 md=setelementstype(md,'hutter','all'); 22 22 23 %Spc the grids on the bed24 pos=find(md. gridonbed);23 %Spc the nodes on the bed 24 pos=find(md.nodeonbed); 25 25 md.spcvelocity(pos,1:3)=1; 26 26 … … 31 31 vx=PatchToVec(md.results.DiagnosticSolution.Vx); 32 32 vy=PatchToVec(md.results.DiagnosticSolution.Vy); 33 vel=zeros(md.numberof grids2d,1);33 vel=zeros(md.numberofnodes2d,1); 34 34 35 grid_vel=0;36 for i=1:md.numberof grids2d35 node_vel=0; 36 for i=1:md.numberofnodes2d 37 37 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.numberof grids2d,1).^2)+...40 sqrt(vx(i+(j-1)*md.numberof grids2d,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)); 41 41 end 42 vel(i,1)= grid_vel;43 grid_vel=0;42 vel(i,1)=node_vel; 43 node_vel=0; 44 44 end 45 45 … … 63 63 subplot(2,2,3) 64 64 hold on; 65 plot(sqrt((md.x(1:md.numberof grids2d)).^2+(md.y(1:md.numberofgrids2d)).^2),vel,'r.');65 plot(sqrt((md.x(1:md.numberofnodes2d)).^2+(md.y(1:md.numberofnodes2d)).^2),vel,'r.'); 66 66 plot(sqrt((md.x2d).^2+(md.y2d).^2),vel_obs,'b.'); 67 67 title('Analytical vs calculated velocity','FontSize',14,'FontWeight','bold'); -
issm/trunk/test/NightlyRun/test1206.m
r8295 r8304 21 21 md=setelementstype(md,'pattyn','all'); 22 22 23 %Spc the grids on the bed24 pos=find(md. gridonbed);23 %Spc the nodes on the bed 24 pos=find(md.nodeonbed); 25 25 md.spcvelocity(pos,1:3)=1; 26 26 … … 31 31 vx=PatchToVec(md.results.DiagnosticSolution.Vx); 32 32 vy=PatchToVec(md.results.DiagnosticSolution.Vy); 33 vel=zeros(md.numberof grids2d,1);33 vel=zeros(md.numberofnodes2d,1); 34 34 35 grid_vel=0;36 for i=1:md.numberof grids2d35 node_vel=0; 36 for i=1:md.numberofnodes2d 37 37 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.numberof grids2d,1).^2)+...40 sqrt(vx(i+(j-1)*md.numberof grids2d,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)); 41 41 end 42 vel(i,1)= grid_vel;43 grid_vel=0;42 vel(i,1)=node_vel; 43 node_vel=0; 44 44 end 45 45 … … 62 62 subplot(2,2,3) 63 63 hold on; 64 plot(sqrt((md.x(1:md.numberof grids2d)).^2+(md.y(1:md.numberofgrids2d)).^2),vel,'r.');64 plot(sqrt((md.x(1:md.numberofnodes2d)).^2+(md.y(1:md.numberofnodes2d)).^2),vel,'r.'); 65 65 plot(sqrt((md.x2d).^2+(md.y2d).^2),vel_obs,'b.'); 66 66 title('Analytical vs calculated velocity','FontSize',14,'FontWeight','bold'); -
issm/trunk/test/NightlyRun/test1207.m
r8295 r8304 21 21 md=setelementstype(md,'stokes','all'); 22 22 23 %Spc the grids on the bed24 pos=find(md. gridonbed);23 %Spc the nodes on the bed 24 pos=find(md.nodeonbed); 25 25 md.spcvelocity(pos,1:3)=1; 26 26 … … 31 31 vx=PatchToVec(md.results.DiagnosticSolution.Vx); 32 32 vy=PatchToVec(md.results.DiagnosticSolution.Vy); 33 vel=zeros(md.numberof grids2d,1);33 vel=zeros(md.numberofnodes2d,1); 34 34 35 grid_vel=0;36 for i=1:md.numberof grids2d35 node_vel=0; 36 for i=1:md.numberofnodes2d 37 37 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.numberof grids2d,1).^2)+...40 sqrt(vx(i+(j-1)*md.numberof grids2d,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)); 41 41 end 42 vel(i,1)= grid_vel;43 grid_vel=0;42 vel(i,1)=node_vel; 43 node_vel=0; 44 44 end 45 45 … … 62 62 subplot(2,2,3) 63 63 hold on; 64 plot(sqrt((md.x(1:md.numberof grids2d)).^2+(md.y(1:md.numberofgrids2d)).^2),vel,'r.');64 plot(sqrt((md.x(1:md.numberofnodes2d)).^2+(md.y(1:md.numberofnodes2d)).^2),vel,'r.'); 65 65 plot(sqrt((md.x2d).^2+(md.y2d).^2),vel_obs,'b.'); 66 66 title('Analytical vs calculated velocity','FontSize',14,'FontWeight','bold'); -
issm/trunk/test/NightlyRun/test1208.m
r8295 r8304 12 12 md=setelementstype(md,'hutter','all'); 13 13 14 %Spc the grids on the bed15 pos=find(md. gridonbed);14 %Spc the nodes on the bed 15 pos=find(md.nodeonbed); 16 16 md.spcvelocity(pos,1:3)=1; 17 17 -
issm/trunk/test/NightlyRun/test1301.m
r8295 r8304 13 13 14 14 %Some conditions specific to melting test 15 md.pressure=zeros(md.numberof grids,1);16 md.observed_temperature=273.15*ones(md.numberof grids,1);15 md.pressure=zeros(md.numberofnodes,1); 16 md.observed_temperature=273.15*ones(md.numberofnodes,1); 17 17 md.temperature=md.observed_temperature; 18 18 md.spctemperature(:,2)=md.observed_temperature; … … 28 28 29 29 %plot results 30 comp_melting=zeros(md.numberof grids,1);30 comp_melting=zeros(md.numberofnodes,1); 31 31 comp_melting(md.results.ThermalSolution.MeltingRate.index)=md.results.ThermalSolution.MeltingRate.value; 32 32 relative=abs((comp_melting-melting)./melting)*100; -
issm/trunk/test/NightlyRun/test1302.m
r8295 r8304 14 14 pos1=find(md.elementonbed); md.spctemperature(md.elements(pos1,1:3),1)=1; md.spctemperature(md.elements(pos1,1:3),2)=10; 15 15 pos2=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.numberof grids,1);16 md.vz=0.1*ones(md.numberofnodes,1); 17 17 md.vel=sqrt( md.vx.^2+ md.vy.^2+ md.vz.^2); 18 md.pressure=zeros(md.numberof grids,1);18 md.pressure=zeros(md.numberofnodes,1); 19 19 20 20 %analytical results … … 29 29 30 30 %plot results 31 comp_temp=zeros(md.numberof grids,1);31 comp_temp=zeros(md.numberofnodes,1); 32 32 comp_temp(md.results.ThermalSolution.Temperature.index)=md.results.ThermalSolution.Temperature.value; 33 33 relative=abs((comp_temp-md.observed_temperature)./md.observed_temperature)*100; -
issm/trunk/test/NightlyRun/test1303.m
r8295 r8304 13 13 pos1=find(md.elementonbed); md.spctemperature(md.elements(pos1,1:3),1)=1; md.spctemperature(md.elements(pos1,1:3),2)=10; 14 14 pos2=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.numberof grids,1);15 md.pressure=zeros(md.numberofnodes,1); 16 16 17 17 %analytical results … … 24 24 25 25 %plot results 26 comp_temp=zeros(md.numberof grids,1);26 comp_temp=zeros(md.numberofnodes,1); 27 27 comp_temp(md.results.ThermalSolution.Temperature.index)=md.results.ThermalSolution.Temperature.value; 28 28 relative=abs((comp_temp-md.observed_temperature)./md.observed_temperature)*100; -
issm/trunk/test/NightlyRun/test1304.m
r8295 r8304 13 13 14 14 pos2=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.numberof grids,1);15 md.pressure=zeros(md.numberofnodes,1); 16 16 md.geothermalflux(:)=0.1; %100mW/m^2 17 17 … … 25 25 26 26 %plot results 27 comp_temp=zeros(md.numberof grids,1);27 comp_temp=zeros(md.numberofnodes,1); 28 28 comp_temp(md.results.ThermalSolution.Temperature.index)=md.results.ThermalSolution.Temperature.value; 29 29 relative=abs((comp_temp-md.observed_temperature)./md.observed_temperature)*100; -
issm/trunk/test/NightlyRun/test131.m
r8295 r8304 3 3 md=parameterize(md,'../Par/SquareShelfConstrained.par'); 4 4 %Add boundary conditions on thickness on the border 5 pos=find(md. gridonboundary);5 pos=find(md.nodeonboundary); 6 6 md.spcthickness(pos,1)=1; 7 7 md.spcthickness(pos,2)=md.thickness(pos); -
issm/trunk/test/NightlyRun/test132.m
r8295 r8304 3 3 md=parameterize(md,'../Par/SquareShelfConstrained.par'); 4 4 %Add boundary conditions on thickness on the border 5 pos=find(md. gridonboundary);5 pos=find(md.nodeonboundary); 6 6 md.spcthickness(pos,1)=1; 7 7 md.spcthickness(pos,2)=md.thickness(pos); -
issm/trunk/test/NightlyRun/test133.m
r8295 r8304 4 4 md=extrude(md,5,1); 5 5 %Add boundary conditions on thickness on the border 6 pos=find(md. gridonboundary);6 pos=find(md.nodeonboundary); 7 7 md.spcthickness(pos,1)=1; 8 8 md.spcthickness(pos,2)=md.thickness(pos); -
issm/trunk/test/NightlyRun/test134.m
r8295 r8304 4 4 md=extrude(md,5,1); 5 5 %Add boundary conditions on thickness on the border 6 pos=find(md. gridonboundary);6 pos=find(md.nodeonboundary); 7 7 md.spcthickness(pos,1)=1; 8 8 md.spcthickness(pos,2)=md.thickness(pos); -
issm/trunk/test/NightlyRun/test223.m
r8295 r8304 7 7 md.control_analysis=1; 8 8 md.control_type=RheologyBbarEnum(); 9 md.cm_min=10^6*ones(md.numberof grids,1);10 md.cm_max=2*10^9*ones(md.numberof grids,1);9 md.cm_min=10^6*ones(md.numberofnodes,1); 10 md.cm_max=2*10^9*ones(md.numberofnodes,1); 11 11 md.nsteps=2; 12 12 md.cm_responses=SurfaceAbsVelMisfitEnum*ones(md.nsteps,1); 13 md.weights=ones(md.numberof grids,1);13 md.weights=ones(md.numberofnodes,1); 14 14 md.optscal=10^7*ones(md.nsteps,1); 15 15 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test224.m
r8295 r8304 7 7 md.control_analysis=1; 8 8 md.control_type=RheologyBbarEnum(); 9 md.cm_min=10^6*ones(md.numberof grids,1);10 md.cm_max=2*10^9*ones(md.numberof grids,1);9 md.cm_min=10^6*ones(md.numberofnodes,1); 10 md.cm_max=2*10^9*ones(md.numberofnodes,1); 11 11 md.nsteps=2; 12 12 md.cm_responses=SurfaceAbsVelMisfitEnum*ones(md.nsteps,1); 13 md.weights=ones(md.numberof grids,1);13 md.weights=ones(md.numberofnodes,1); 14 14 md.optscal=10^7*ones(md.nsteps,1); 15 15 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test225.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=RheologyBbarEnum(); 10 md.cm_min=10^6*ones(md.numberof grids,1);11 md.cm_max=2*10^9*ones(md.numberof grids,1);10 md.cm_min=10^6*ones(md.numberofnodes,1); 11 md.cm_max=2*10^9*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=10^7*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test226.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=RheologyBbarEnum(); 10 md.cm_min=10^6*ones(md.numberof grids,1);11 md.cm_max=2*10^9*ones(md.numberof grids,1);10 md.cm_min=10^6*ones(md.numberofnodes,1); 11 md.cm_max=2*10^9*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=10^7*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test227.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=RheologyBbarEnum(); 10 md.cm_min=10^6*ones(md.numberof grids,1);11 md.cm_max=2*10^9*ones(md.numberof grids,1);10 md.cm_min=10^6*ones(md.numberofnodes,1); 11 md.cm_max=2*10^9*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=10^7*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test228.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=RheologyBbarEnum(); 10 md.cm_min=10^6*ones(md.numberof grids,1);11 md.cm_max=2*10^9*ones(md.numberof grids,1);10 md.cm_min=10^6*ones(md.numberofnodes,1); 11 md.cm_max=2*10^9*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=10^7*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test229.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=RheologyBbarEnum(); 10 md.cm_min=10^6*ones(md.numberof grids,1);11 md.cm_max=2*10^9*ones(md.numberof grids,1);10 md.cm_min=10^6*ones(md.numberofnodes,1); 11 md.cm_max=2*10^9*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=10^7*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test230.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=RheologyBbarEnum(); 10 md.cm_min=10^6*ones(md.numberof grids,1);11 md.cm_max=2*10^9*ones(md.numberof grids,1);10 md.cm_min=10^6*ones(md.numberofnodes,1); 11 md.cm_max=2*10^9*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=10^7*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test233.m
r8295 r8304 21 21 22 22 h=1000; 23 md.thickness=h*ones(md.numberof grids,1);23 md.thickness=h*ones(md.numberofnodes,1); 24 24 md.bed=-md.rho_ice/md.rho_water*md.thickness; 25 25 md.surface=md.bed+md.thickness; 26 26 27 27 %Initial velocity and pressure 28 md.vx=zeros(md.numberof grids,1);29 md.vy=zeros(md.numberof grids,1);30 md.vz=zeros(md.numberof grids,1);31 md.pressure=zeros(md.numberof grids,1);28 md.vx=zeros(md.numberofnodes,1); 29 md.vy=zeros(md.numberofnodes,1); 30 md.vz=zeros(md.numberofnodes,1); 31 md.pressure=zeros(md.numberofnodes,1); 32 32 33 33 %Materials 34 md.observed_temperature=(273-20)*ones(md.numberof grids,1);34 md.observed_temperature=(273-20)*ones(md.numberofnodes,1); 35 35 md.rheology_B=paterson(md.observed_temperature); 36 36 md.rheology_n=3*ones(md.numberofelements,1); … … 38 38 39 39 %Boundary conditions: 40 md.spcvelocity=zeros(md.numberof grids,6);40 md.spcvelocity=zeros(md.numberofnodes,6); 41 41 42 42 %constrain flanks to 0 normal velocity … … 51 51 52 52 %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,:);53 nodeonicefront=zeros(md.numberofnodes,1); 54 pos=find(md.y==ymax); nodeonicefront(pos)=1; 55 pos=find(nodeonicefront(md.segments(:,1)) | nodeonicefront(md.segments(:,2))); pressureload=md.segments(pos,:); 56 56 pressureload=[pressureload WaterEnum*md.elementoniceshelf(pressureload(:,end))]; 57 57 md.pressureload=pressureload; -
issm/trunk/test/NightlyRun/test234.m
r8295 r8304 21 21 22 22 h=1000; 23 md.thickness=h*ones(md.numberof grids,1);23 md.thickness=h*ones(md.numberofnodes,1); 24 24 md.bed=-md.rho_ice/md.rho_water*md.thickness; 25 25 md.surface=md.bed+md.thickness; 26 26 27 27 %Initial velocity and pressure 28 md.vx=zeros(md.numberof grids,1);29 md.vy=zeros(md.numberof grids,1);30 md.vz=zeros(md.numberof grids,1);31 md.pressure=zeros(md.numberof grids,1);28 md.vx=zeros(md.numberofnodes,1); 29 md.vy=zeros(md.numberofnodes,1); 30 md.vz=zeros(md.numberofnodes,1); 31 md.pressure=zeros(md.numberofnodes,1); 32 32 33 33 %Materials 34 md.observed_temperature=(273-20)*ones(md.numberof grids,1);34 md.observed_temperature=(273-20)*ones(md.numberofnodes,1); 35 35 md.rheology_B=paterson(md.observed_temperature); 36 36 md.rheology_n=3*ones(md.numberofelements,1); … … 38 38 39 39 %Boundary conditions: 40 md.spcvelocity=zeros(md.numberof grids,6);40 md.spcvelocity=zeros(md.numberofnodes,6); 41 41 42 42 %constrain flanks to 0 normal velocity … … 51 51 52 52 %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,:);53 nodeonicefront=zeros(md.numberofnodes,1); 54 pos=find(md.y==ymax); nodeonicefront(pos)=1; 55 pos=find(nodeonicefront(md.segments(:,1)) | nodeonicefront(md.segments(:,2))); pressureload=md.segments(pos,:); 56 56 pressureload=[pressureload WaterEnum*md.elementoniceshelf(pressureload(:,end))]; 57 57 md.pressureload=pressureload; -
issm/trunk/test/NightlyRun/test235.m
r8295 r8304 18 18 19 19 h=1000; 20 md.thickness=h*ones(md.numberof grids,1);20 md.thickness=h*ones(md.numberofnodes,1); 21 21 md.bed=-md.rho_ice/md.rho_water*md.thickness; 22 22 md.surface=md.bed+md.thickness; 23 23 24 24 %Initial velocity and pressure 25 md.vx=zeros(md.numberof grids,1);26 md.vy=zeros(md.numberof grids,1);27 md.vz=zeros(md.numberof grids,1);28 md.pressure=zeros(md.numberof grids,1);25 md.vx=zeros(md.numberofnodes,1); 26 md.vy=zeros(md.numberofnodes,1); 27 md.vz=zeros(md.numberofnodes,1); 28 md.pressure=zeros(md.numberofnodes,1); 29 29 30 30 %Materials 31 md.observed_temperature=(273-20)*ones(md.numberof grids,1);31 md.observed_temperature=(273-20)*ones(md.numberofnodes,1); 32 32 md.rheology_B=paterson(md.observed_temperature); 33 33 md.rheology_n=3*ones(md.numberofelements,1); … … 35 35 36 36 %Boundary conditions: 37 md.spcvelocity=zeros(md.numberof grids,6);37 md.spcvelocity=zeros(md.numberofnodes,6); 38 38 39 39 %constrain flanks to 0 normal velocity … … 48 48 49 49 %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,:);50 nodeonicefront=zeros(md.numberofnodes,1); 51 pos=find(md.y==ymax); nodeonicefront(pos)=1; 52 pos=find(nodeonicefront(md.segments(:,1)) | nodeonicefront(md.segments(:,2))); pressureload=md.segments(pos,:); 53 53 pressureload=[pressureload WaterEnum*md.elementoniceshelf(pressureload(:,end))]; 54 54 md.pressureload=pressureload; 55 55 56 56 %partitioning 57 md.npart=md.numberof grids;57 md.npart=md.numberofnodes; 58 58 md=partitioner(md,'package','linear','npart',md.npart); 59 59 md.part=md.part-1; -
issm/trunk/test/NightlyRun/test236.m
r8295 r8304 18 18 19 19 h=1000; 20 md.thickness=h*ones(md.numberof grids,1);20 md.thickness=h*ones(md.numberofnodes,1); 21 21 md.bed=-md.rho_ice/md.rho_water*md.thickness; 22 22 md.surface=md.bed+md.thickness; 23 23 24 24 %Initial velocity and pressure 25 md.vx=zeros(md.numberof grids,1);26 md.vy=zeros(md.numberof grids,1);27 md.vz=zeros(md.numberof grids,1);28 md.pressure=zeros(md.numberof grids,1);25 md.vx=zeros(md.numberofnodes,1); 26 md.vy=zeros(md.numberofnodes,1); 27 md.vz=zeros(md.numberofnodes,1); 28 md.pressure=zeros(md.numberofnodes,1); 29 29 30 30 %Materials 31 md.observed_temperature=(273-20)*ones(md.numberof grids,1);31 md.observed_temperature=(273-20)*ones(md.numberofnodes,1); 32 32 md.rheology_B=paterson(md.observed_temperature); 33 33 md.rheology_n=3*ones(md.numberofelements,1); … … 35 35 36 36 %Boundary conditions: 37 md.spcvelocity=zeros(md.numberof grids,6);37 md.spcvelocity=zeros(md.numberofnodes,6); 38 38 39 39 %constrain flanks to 0 normal velocity … … 48 48 49 49 %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,:);50 nodeonicefront=zeros(md.numberofnodes,1); 51 pos=find(md.y==ymax); nodeonicefront(pos)=1; 52 pos=find(nodeonicefront(md.segments(:,1)) | nodeonicefront(md.segments(:,2))); pressureload=md.segments(pos,:); 53 53 pressureload=[pressureload WaterEnum*md.elementoniceshelf(pressureload(:,end))]; 54 54 md.pressureload=pressureload; 55 55 56 56 %partitioning 57 md.npart=md.numberof grids;57 md.npart=md.numberofnodes; 58 58 md=partitioner(md,'package','linear','npart',md.npart); 59 59 md.part=md.part-1; -
issm/trunk/test/NightlyRun/test337.m
r8295 r8304 7 7 md.control_analysis=1; 8 8 md.control_type=DragCoefficientEnum(); 9 md.cm_min=1*ones(md.numberof grids,1);10 md.cm_max=200*ones(md.numberof grids,1);9 md.cm_min=1*ones(md.numberofnodes,1); 10 md.cm_max=200*ones(md.numberofnodes,1); 11 11 md.nsteps=2; 12 12 md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1); 13 md.weights=ones(md.numberof grids,1);13 md.weights=ones(md.numberofnodes,1); 14 14 md.optscal=3*ones(md.nsteps,1); 15 15 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test338.m
r8295 r8304 7 7 md.control_analysis=1; 8 8 md.control_type=DragCoefficientEnum(); 9 md.cm_min=1*ones(md.numberof grids,1);10 md.cm_max=200*ones(md.numberof grids,1);9 md.cm_min=1*ones(md.numberofnodes,1); 10 md.cm_max=200*ones(md.numberofnodes,1); 11 11 md.nsteps=2; 12 12 md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1); 13 md.weights=ones(md.numberof grids,1);13 md.weights=ones(md.numberofnodes,1); 14 14 md.optscal=3*ones(md.nsteps,1); 15 15 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test339.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=DragCoefficientEnum(); 10 md.cm_min=1*ones(md.numberof grids,1);11 md.cm_max=200*ones(md.numberof grids,1);10 md.cm_min=1*ones(md.numberofnodes,1); 11 md.cm_max=200*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=3*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test340.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=DragCoefficientEnum(); 10 md.cm_min=1*ones(md.numberof grids,1);11 md.cm_max=200*ones(md.numberof grids,1);10 md.cm_min=1*ones(md.numberofnodes,1); 11 md.cm_max=200*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=3*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test341.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=DragCoefficientEnum(); 10 md.cm_min=1*ones(md.numberof grids,1);11 md.cm_max=200*ones(md.numberof grids,1);10 md.cm_min=1*ones(md.numberofnodes,1); 11 md.cm_max=200*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceRelVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=3*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test342.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=DragCoefficientEnum(); 10 md.cm_min=1*ones(md.numberof grids,1);11 md.cm_max=200*ones(md.numberof grids,1);10 md.cm_min=1*ones(md.numberofnodes,1); 11 md.cm_max=200*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceRelVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=3*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test343.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=DragCoefficientEnum(); 10 md.cm_min=1*ones(md.numberof grids,1);11 md.cm_max=200*ones(md.numberof grids,1);10 md.cm_min=1*ones(md.numberofnodes,1); 11 md.cm_max=200*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceLogVxVyMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=3*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test344.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=DragCoefficientEnum(); 10 md.cm_min=1*ones(md.numberof grids,1);11 md.cm_max=200*ones(md.numberof grids,1);10 md.cm_min=1*ones(md.numberofnodes,1); 11 md.cm_max=200*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceLogVxVyMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=3*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test423.m
r8295 r8304 6 6 7 7 %partitioning 8 md.npart=md.numberof grids;8 md.npart=md.numberofnodes; 9 9 md=partitioner(md,'package','linear','npart',md.npart); 10 10 md.part=md.part-1; -
issm/trunk/test/NightlyRun/test424.m
r8295 r8304 6 6 7 7 %partitioning 8 md.npart=md.numberof grids;8 md.npart=md.numberofnodes; 9 9 md=partitioner(md,'package','linear','npart',md.npart); 10 10 md.part=md.part-1; -
issm/trunk/test/NightlyRun/test429.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=DragCoefficientEnum(); 10 md.cm_min=1*ones(md.numberof grids,1);11 md.cm_max=200*ones(md.numberof grids,1);10 md.cm_min=1*ones(md.numberofnodes,1); 11 md.cm_max=200*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=3*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test430.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=DragCoefficientEnum(); 10 md.cm_min=1*ones(md.numberof grids,1);11 md.cm_max=200*ones(md.numberof grids,1);10 md.cm_min=1*ones(md.numberofnodes,1); 11 md.cm_max=200*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=3*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test431.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=DragCoefficientEnum(); 10 md.cm_min=1*ones(md.numberof grids,1);11 md.cm_max=200*ones(md.numberof grids,1);10 md.cm_min=1*ones(md.numberofnodes,1); 11 md.cm_max=200*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceRelVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=3*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test432.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=DragCoefficientEnum(); 10 md.cm_min=1*ones(md.numberof grids,1);11 md.cm_max=200*ones(md.numberof grids,1);10 md.cm_min=1*ones(md.numberofnodes,1); 11 md.cm_max=200*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceRelVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=3*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test435.m
r5955 r8304 11 11 md.part=md.part-1; 12 12 13 vector=(1:1:md.numberof grids)';13 vector=(1:1:md.numberofnodes)'; 14 14 vector_on_partition=AreaAverageOntoPartition(md,vector); 15 vector_on_ grids=vector_on_partition(md.part+1);15 vector_on_nodes=vector_on_partition(md.part+1); 16 16 17 field_names ={'vector_on_ grids'};17 field_names ={'vector_on_nodes'}; 18 18 field_tolerances={1e-11}; 19 19 field_values={... 20 vector_on_ grids,...20 vector_on_nodes,... 21 21 }; -
issm/trunk/test/NightlyRun/test446.m
r8295 r8304 4 4 %mesh {{{1 5 5 md=roundmesh(model,radius,resolution); 6 %fix center gridto 0,06 %fix center node to 0,0 7 7 rad=sqrt((md.x).*md.x+(md.y).*md.y); 8 8 pos=find(rad==min(rad)); -
issm/trunk/test/NightlyRun/test521.m
r8295 r8304 13 13 md.control_analysis=1; 14 14 md.control_type=RheologyBbarEnum(); 15 md.cm_min=10^6*ones(md.numberof grids,1);16 md.cm_max=2*10^9*ones(md.numberof grids,1);15 md.cm_min=10^6*ones(md.numberofnodes,1); 16 md.cm_max=2*10^9*ones(md.numberofnodes,1); 17 17 md.nsteps=2; 18 18 md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1); 19 md.weights=ones(md.numberof grids,1);19 md.weights=ones(md.numberofnodes,1); 20 20 md.optscal=10^8*ones(md.nsteps,1); 21 21 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test522.m
r8295 r8304 13 13 md.control_analysis=1; 14 14 md.control_type=RheologyBbarEnum(); 15 md.cm_min=10^6*ones(md.numberof grids,1);16 md.cm_max=2*10^9*ones(md.numberof grids,1);15 md.cm_min=10^6*ones(md.numberofnodes,1); 16 md.cm_max=2*10^9*ones(md.numberofnodes,1); 17 17 md.nsteps=2; 18 18 md.cm_responses=SurfaceAbsVelMisfitEnum()*ones(md.nsteps,1); 19 md.weights=ones(md.numberof grids,1);19 md.weights=ones(md.numberofnodes,1); 20 20 md.optscal=10^8*ones(md.nsteps,1); 21 21 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test523.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=DragCoefficientEnum(); 10 md.cm_min=1*ones(md.numberof grids,1);11 md.cm_max=200*ones(md.numberof grids,1);10 md.cm_min=1*ones(md.numberofnodes,1); 11 md.cm_max=200*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=3*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test524.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=DragCoefficientEnum(); 10 md.cm_min=1*ones(md.numberof grids,1);11 md.cm_max=200*ones(md.numberof grids,1);10 md.cm_min=1*ones(md.numberofnodes,1); 11 md.cm_max=200*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=3*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test525.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=DragCoefficientEnum(); 10 md.cm_min=1*ones(md.numberof grids,1);11 md.cm_max=200*ones(md.numberof grids,1);10 md.cm_min=1*ones(md.numberofnodes,1); 11 md.cm_max=200*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=3*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test526.m
r8295 r8304 8 8 md.control_analysis=1; 9 9 md.control_type=DragCoefficientEnum(); 10 md.cm_min=1*ones(md.numberof grids,1);11 md.cm_max=200*ones(md.numberof grids,1);10 md.cm_min=1*ones(md.numberofnodes,1); 11 md.cm_max=200*ones(md.numberofnodes,1); 12 12 md.nsteps=2; 13 13 md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1); 14 md.weights=ones(md.numberof grids,1);14 md.weights=ones(md.numberofnodes,1); 15 15 md.optscal=3*ones(md.nsteps,1); 16 16 md.maxiter=2*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test527.m
r5098 r8304 27 27 28 28 %refine existing mesh 3 29 hVertices=NaN*ones(md.numberof grids,1);30 hVertices(find(md. gridonboundary))=500;29 hVertices=NaN*ones(md.numberofnodes,1); 30 hVertices(find(md.nodeonboundary))=500; 31 31 md2=bamg(md,'metric',md.dummy,'hmin',1000,'hmax',20000,'gradation',3,'geometricalmetric',1,'anisomax',1,'hVertices',hVertices); 32 32 x5=md2.x; -
issm/trunk/test/NightlyRun/test621.m
r8295 r8304 10 10 md.control_type=DhDtEnum; 11 11 md.thickness_obs=md.thickness; 12 md.cm_min=-50*ones(md.numberof grids,1);13 md.cm_max=50*ones(md.numberof grids,1);12 md.cm_min=-50*ones(md.numberofnodes,1); 13 md.cm_max=50*ones(md.numberofnodes,1); 14 14 md.cm_responses=ThicknessAbsMisfitEnum*ones(md.nsteps,1); 15 md.weights=ones(md.numberof grids,1);15 md.weights=ones(md.numberofnodes,1); 16 16 md.optscal=10/md.yts*ones(md.nsteps,1); 17 17 md.maxiter=4*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test622.m
r8295 r8304 10 10 md.control_type=DhDtEnum; 11 11 md.thickness_obs=md.thickness; 12 md.cm_min=-50*ones(md.numberof grids,1);13 md.cm_max=50*ones(md.numberof grids,1);12 md.cm_min=-50*ones(md.numberofnodes,1); 13 md.cm_max=50*ones(md.numberofnodes,1); 14 14 md.cm_responses=ThicknessAbsMisfitEnum*ones(md.nsteps,1); 15 md.weights=ones(md.numberof grids,1);15 md.weights=ones(md.numberofnodes,1); 16 16 md.optscal=10/md.yts*ones(md.nsteps,1); 17 17 md.maxiter=4*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test623.m
r8295 r8304 11 11 md.control_type=DhDtEnum; 12 12 md.thickness_obs=md.thickness; 13 md.cm_min=-50*ones(md.numberof grids,1);14 md.cm_max=50*ones(md.numberof grids,1);13 md.cm_min=-50*ones(md.numberofnodes,1); 14 md.cm_max=50*ones(md.numberofnodes,1); 15 15 md.cm_responses=ThicknessAbsMisfitEnum*ones(md.nsteps,1); 16 md.weights=ones(md.numberof grids,1);16 md.weights=ones(md.numberofnodes,1); 17 17 md.optscal=10/md.yts*ones(md.nsteps,1); 18 18 md.maxiter=4*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test624.m
r8295 r8304 11 11 md.control_type=DhDtEnum; 12 12 md.thickness_obs=md.thickness; 13 md.cm_min=-50*ones(md.numberof grids,1);14 md.cm_max=50*ones(md.numberof grids,1);13 md.cm_min=-50*ones(md.numberofnodes,1); 14 md.cm_max=50*ones(md.numberofnodes,1); 15 15 md.cm_responses=ThicknessAbsMisfitEnum*ones(md.nsteps,1); 16 md.weights=ones(md.numberof grids,1);16 md.weights=ones(md.numberofnodes,1); 17 17 md.optscal=10/md.yts*ones(md.nsteps,1); 18 18 md.maxiter=4*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test625.m
r8295 r8304 7 7 %Ice sheet only 8 8 md=modelextract(md,md.elementonicesheet); 9 pos=find(md. gridonboundary);9 pos=find(md.nodeonboundary); 10 10 md.spcthickness(pos,1)=1; 11 11 md.spcthickness(:,2)=md.thickness; … … 19 19 md.prognostic_DG=0; 20 20 md.optscal=[10/md.yts*ones(md.nsteps,1) 10/md.yts*ones(md.nsteps,1)]; 21 md.cm_min=[-2000*ones(md.numberof grids,1) -2000*ones(md.numberofgrids,1)];22 md.cm_max=[+2000*ones(md.numberof grids,1) +2000*ones(md.numberofgrids,1)];21 md.cm_min=[-2000*ones(md.numberofnodes,1) -2000*ones(md.numberofnodes,1)]; 22 md.cm_max=[+2000*ones(md.numberofnodes,1) +2000*ones(md.numberofnodes,1)]; 23 23 md.cm_responses=ThicknessAbsMisfitEnum*ones(md.nsteps,1); 24 md.weights=ones(md.numberof grids,1);24 md.weights=ones(md.numberofnodes,1); 25 25 md.maxiter=4*ones(md.nsteps,1); 26 26 md.cm_jump=0.99*ones(md.nsteps,1); -
issm/trunk/test/NightlyRun/test626.m
r8295 r8304 7 7 %Ice sheet only 8 8 md=modelextract(md,md.elementonicesheet); 9 pos=find(md. gridonboundary);9 pos=find(md.nodeonboundary); 10 10 md.spcthickness(pos,1)=1; 11 11 md.spcthickness(:,2)=md.thickness; … … 19 19 md.prognostic_DG=0; 20 20 md.optscal=[10/md.yts*ones(md.nsteps,1) 10/md.yts*ones(md.nsteps,1)]; 21 md.cm_min=[-2000*ones(md.numberof grids,1) -2000*ones(md.numberofgrids,1)];22 md.cm_max=[+2000*ones(md.numberof grids,1) +2000*ones(md.numberofgrids,1)];21 md.cm_min=[-2000*ones(md.numberofnodes,1) -2000*ones(md.numberofnodes,1)]; 22 md.cm_max=[+2000*ones(md.numberofnodes,1) +2000*ones(md.numberofnodes,1)]; 23 23 md.cm_responses=ThicknessAbsMisfitEnum*ones(md.nsteps,1); 24 md.weights=ones(md.numberof grids,1);24 md.weights=ones(md.numberofnodes,1); 25 25 md.maxiter=4*ones(md.nsteps,1); 26 26 md.cm_jump=0.99*ones(md.nsteps,1);
Note:
See TracChangeset
for help on using the changeset viewer.