source: issm/trunk/test/Verification/IceShelfIceFrontM2dRifts_24/buildme.m@ 1803

Last change on this file since 1803 was 1803, checked in by Eric.Larour, 16 years ago

Some adaptation of Square.par to new rift formulation

  • Property svn:executable set to *
File size: 731 bytes
RevLine 
[1733]1%mesh resolution
[1803]2resolution=10000;
3riftres=500;
[1733]4
5md=model;
6
7%refine mesh for rift
8expcoarsen('Rifts2.exp','Rifts.exp',riftres);
9
10%md=mesh(md,'DomainOutline.exp',50000);
11md=mesh(md,'DomainOutline.exp','Rifts2.exp',resolution);
12md=meshprocessrifts(md);
[1803]13md.numberofelements
[1733]14md=geography(md,'all','');
15md=parameterize(md,'Square.par');
16md=setelementstype(md,'macayeal','all');
[1803]17
18md.rifts.fill=MelangeEnum();
19md.rifts.fraction=0;
20md.solverstring=' -mat_type aijmumps -ksp_type preonly -pc_type lu -mat_mumps_icntl_14 80 ';
21md.cluster='astrid';
22md.np=8;
23md.penalty_lock=2;
24md.rifts.fractionincrement=.01;
25md=solve(md,'analysis_type','diagnostic');md=tres(md,'diagnostic');
26plotmodel(md,'data','riftrelvel','data','riftvel','data','vel');
Note: See TracBrowser for help on using the repository browser.