Changeset 5410


Ignore:
Timestamp:
08/19/10 10:41:50 (15 years ago)
Author:
Mathieu Morlighem
Message:

Better runme

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r5182 r5410  
    11%test the anisotropic mesh adaptation
     2printingflag=true;
    23
    34%create square mesh
    45L=1; %in m
    5 nx=70; %numberof nodes in x direction
    6 ny=70;
     6nx=30; %numberof nodes in x direction
     7ny=30;
    78md=model;
    89
    910%mesh adaptation loop YAMS
    1011md=squaremesh(md,L,L,nx,ny);
    11 plotmodel(md,'data','mesh');pause(0.5);
    12 for i=1:1
    13         u=4*md.x-2; v=4*md.y-2;
    14         md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
    15                 +tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
    16                 +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
    17         md=YamsCall(md,md.vel_obs,0.005,0.3,1.3,10^-4);
    18         plotmodel(md,'data','mesh');pause(0.5);
     12u=4*md.x-2; v=4*md.y-2;
     13md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
     14        +tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
     15        +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
     16plotmodel(md,'data','vel_obs','data','vel_obs','nlines',1,'ncols',2,'title','','figposition',[500 500 1000 500],'axis#all','equal','xlim#all',[0 1],'ylim#all',[0 1],'edgecolor#1','w'); pause(0.5);
     17if printingflag,
     18        set(gcf,'Color','w')
     19        printmodel('mesh2_yams1','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
     20        system(['mv mesh2_yams1.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
    1921end
    20 for i=1:1
    21         u=4*md.x-2; v=4*md.y-2;
    22         md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
    23                 +tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
    24                 +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
    25         md=YamsCall(md,md.vel_obs,0.005,0.3,2,10^-4);
    26         plotmodel(md,'data','mesh');pause(0.5);
     22
     23md=YamsCall(md,md.vel_obs,0.005,0.3,2.3,10^-2);
     24u=4*md.x-2; v=4*md.y-2;
     25md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
     26        +tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
     27        +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
     28plotmodel(md,'data','vel_obs','data','vel_obs','nlines',1,'ncols',2,'title','','figposition',[500 500 1000 500],'axis#all','equal','xlim#all',[0 1],'ylim#all',[0 1],'edgecolor#1','w'); pause(0.5);
     29if printingflag,
     30        set(gcf,'Color','w')
     31        printmodel('mesh2_yams2','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
     32        system(['mv mesh2_yams2.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
    2733end
    28 for i=1:3,
    29         u=4*md.x-2; v=4*md.y-2;
    30         md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
    31                 +tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
    32                 +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
    33         md=YamsCall(md,md.vel_obs,0.005,0.3,2,0.001);
    34         plotmodel(md,'data','mesh');pause(0.5);
     34
     35md=YamsCall(md,md.vel_obs,0.005,0.3,3,0.005);
     36u=4*md.x-2; v=4*md.y-2;
     37md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
     38        +tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
     39        +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
     40plotmodel(md,'data','vel_obs','data','vel_obs','nlines',1,'ncols',2,'title','','figposition',[500 500 1000 500],'axis#all','equal','xlim#all',[0 1],'ylim#all',[0 1],'edgecolor#1','w'); pause(0.5);
     41if printingflag,
     42        set(gcf,'Color','w')
     43        printmodel('mesh2_yams3','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
     44        system(['mv mesh2_yams3.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
    3545end
    36 md.vel_obs=tanh(30*((4*md.x-2).^2+(4*md.y-2).^2-0.25))+...
    37         tanh(30*((4*md.x-2.75).^2+(4*md.y-1.25).^2-2.75))+...
    38         tanh(30*((4*md.x-2.75).^2+(4*md.y-2).^2-1.25))+...
    39         tanh(30*((4*md.x-1.25).^2+(4*md.y-2).^2-2.75))+...
    40         tanh(30*((4*md.x-1.25).^2+(4*md.y-2).^2-1.25));
    4146
    4247%mesh adaptation loop BAMG
    4348md=squaremesh(md,L,L,nx,ny);
    44 plotmodel(md,'data','mesh');pause(0.5);
    45 for i=1:1
    46         u=4*md.x-2; v=4*md.y-2;
    47         md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
    48                 +tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
    49                 +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
    50         md.bamg=NaN;
    51         md=bamg(md,'field',md.vel_obs,'hmin',0.005,'hmax',0.3,'gradation',1.3,'err',10^-4);
    52         plotmodel(md,'data','mesh');pause(0.5);
     49u=4*md.x-2; v=4*md.y-2;
     50md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
     51        +tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
     52        +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
     53plotmodel(md,'data','vel_obs','data','vel_obs','nlines',1,'ncols',2,'title','','figposition',[500 500 1000 500],'axis#all','equal','xlim#all',[0 1],'ylim#all',[0 1],'edgecolor#1','w'); pause(0.5);
     54if printingflag,
     55        set(gcf,'Color','w')
     56        printmodel('mesh2_bamg1','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
     57        system(['mv mesh2_bamg1.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
    5358end
    54 for i=1:1
    55         u=4*md.x-2; v=4*md.y-2;
    56         md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
    57                 +tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
    58                 +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
    59         md.bamg=NaN;
    60         md=bamg(md,'field',md.vel_obs,'hmin',0.005,'hmax',0.3,'gradation',2,'err',10^-4);
    61         plotmodel(md,'data','mesh');pause(0.5);
     59
     60md.bamg=NaN;
     61md=bamg(md,'field',md.vel_obs,'hmin',0.005,'hmax',0.3,'gradation',2.3,'err',10^-2);
     62u=4*md.x-2; v=4*md.y-2;
     63md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
     64        +tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
     65        +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
     66plotmodel(md,'data','vel_obs','data','vel_obs','nlines',1,'ncols',2,'title','','figposition',[500 500 1000 500],'axis#all','equal','xlim#all',[0 1],'ylim#all',[0 1],'edgecolor#1','w'); pause(0.5);
     67if printingflag,
     68        set(gcf,'Color','w')
     69        printmodel('mesh2_bamg2','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
     70        system(['mv mesh2_bamg2.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
    6271end
    63 for i=1:3,
    64         u=4*md.x-2; v=4*md.y-2;
    65         md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
    66                 +tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
    67                 +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
    68         md.bamg=NaN;
    69         md=bamg(md,'field',md.vel_obs,'hmin',0.001,'hmax',0.3,'gradation',2,'err',10^-3);
    70         plotmodel(md,'data','mesh');pause(0.5);
     72
     73md.bamg=NaN;
     74md=bamg(md,'field',md.vel_obs,'hmin',0.005,'hmax',0.3,'gradation',3,'err',0.005);
     75u=4*md.x-2; v=4*md.y-2;
     76md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
     77        +tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
     78        +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
     79plotmodel(md,'data','vel_obs','data','vel_obs','nlines',1,'ncols',2,'title','','figposition',[500 500 1000 500],'axis#all','equal','xlim#all',[0 1],'ylim#all',[0 1],'edgecolor#1','w'); pause(0.5);
     80if printingflag,
     81        set(gcf,'Color','w')
     82        printmodel('mesh2_bamg3','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
     83        system(['mv mesh2_bamg3.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
    7184end
    72 md.vel_obs=tanh(30*((4*md.x-2).^2+(4*md.y-2).^2-0.25))+...
    73         tanh(30*((4*md.x-2.75).^2+(4*md.y-1.25).^2-2.75))+...
    74         tanh(30*((4*md.x-2.75).^2+(4*md.y-2).^2-1.25))+...
    75         tanh(30*((4*md.x-1.25).^2+(4*md.y-2).^2-2.75))+...
    76         tanh(30*((4*md.x-1.25).^2+(4*md.y-2).^2-1.25));
     85
     86md.bamg=NaN;
     87md=bamg(md,'field',md.vel_obs,'hmin',0.005,'hmax',0.3,'gradation',1.5,'err',0.003,'anisomax',1);
     88u=4*md.x-2; v=4*md.y-2;
     89md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
     90        +tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
     91        +tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
     92plotmodel(md,'data','vel_obs','data','vel_obs','nlines',1,'ncols',2,'title','','figposition',[500 500 1000 500],'axis#all','equal','xlim#all',[0 1],'ylim#all',[0 1],'edgecolor#1','w'); pause(0.5);
     93if printingflag,
     94        set(gcf,'Color','w')
     95        printmodel('mesh2_bamgiso','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
     96        system(['mv mesh2_bamgiso.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
     97end
Note: See TracChangeset for help on using the changeset viewer.