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