Changeset 1345


Ignore:
Timestamp:
07/17/09 11:52:41 (16 years ago)
Author:
Mathieu Morlighem
Message:

minor

Location:
issm/trunk/test/Validation/AnisotropicMesh
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/test/Validation/AnisotropicMesh/runme1.m

    r1343 r1345  
    88md=model;
    99md=squaremesh(md,L,L,nx,ny);
    10 md.gridonboundary=zeros(md.numberofgrids,1);
    11 md.gridonboundary([1 ny ny*(nx-1)+1 nx*ny])=1;
    1210
    1311%mesh adaptation loop
    … …  
    1614        md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
    1715        md=YamsCall(md,md.vel_obs,0.001,0.3,1.3,10^-4);
    18         md.gridonboundary=zeros(md.numberofgrids,1);
    19         md.gridonboundary(find((md.x==0 & md.y==0) | (md.x==0 & md.y==L) | (md.x==L & md.y==0) | (md.x==L & md.y==L)));
    2016        plotmodel(md,'data','mesh');pause(0.5);
    2117end
    … …  
    2319        md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
    2420        md=YamsCall(md,md.vel_obs,0.001,0.3,2,10^-4);
    25         md.gridonboundary=zeros(md.numberofgrids,1);
    26         md.gridonboundary(find((md.x==0 & md.y==0) | (md.x==0 & md.y==L) | (md.x==L & md.y==0) | (md.x==L & md.y==L)));
    2721        plotmodel(md,'data','mesh');pause(0.5);
    2822end
    … …  
    3024        md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
    3125        md=YamsCall(md,md.vel_obs,0.001,0.3,2,0.005);
    32         md.gridonboundary=zeros(md.numberofgrids,1);
    33         md.gridonboundary(find((md.x==0 & md.y==0) | (md.x==0 & md.y==L) | (md.x==L & md.y==0) | (md.x==L & md.y==L)));
    3426        plotmodel(md,'data','mesh');pause(0.5);
    3527end

  • issm/trunk/test/Validation/AnisotropicMesh/runme2.m

    r1343 r1345  
    77md=model;
    88md=squaremesh(md,L,L,nx,ny);
    9 md.gridonboundary=zeros(md.numberofgrids,1);
    10 md.gridonboundary(find((md.x==0 & md.y==0) | (md.x==0 & md.y==L) | (md.x==L & md.y==0) | (md.x==L & md.y==L)));
    119
    1210%mesh adaptation loop
    … …  
    1816                +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
    1917        md=YamsCall(md,md.vel_obs,0.005,0.3,1.3,10^-4);
    20         md.gridonboundary=zeros(md.numberofgrids,1);
    21         md.gridonboundary(find((md.x==0 & md.y==0) | (md.x==0 & md.y==L) | (md.x==L & md.y==0) | (md.x==L & md.y==L)));
    2218        plotmodel(md,'data','mesh');pause(0.5);
    2319end
    … …  
    2824                +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
    2925        md=YamsCall(md,md.vel_obs,0.005,0.3,2,10^-4);
    30         md.gridonboundary=zeros(md.numberofgrids,1);
    31         md.gridonboundary(find((md.x==0 & md.y==0) | (md.x==0 & md.y==L) | (md.x==L & md.y==0) | (md.x==L & md.y==L)));
    3226        plotmodel(md,'data','mesh');pause(0.5);
    3327end
    … …  
    3832                +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
    3933        md=YamsCall(md,md.vel_obs,0.005,0.3,2,0.001);
    40         md.gridonboundary=zeros(md.numberofgrids,1);
    41         md.gridonboundary(find((md.x==0 & md.y==0) | (md.x==0 & md.y==L) | (md.x==L & md.y==0) | (md.x==L & md.y==L)));
    4234        plotmodel(md,'data','mesh');pause(0.5);
    4335end
Note: See TracChangeset for help on using the changeset viewer.