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
|
Rev | Line | |
---|
[1733] | 1 | %mesh resolution
|
---|
[1803] | 2 | resolution=10000;
|
---|
| 3 | riftres=500;
|
---|
[1733] | 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);
|
---|
[1803] | 13 | md.numberofelements
|
---|
[1733] | 14 | md=geography(md,'all','');
|
---|
| 15 | md=parameterize(md,'Square.par');
|
---|
| 16 | md=setelementstype(md,'macayeal','all');
|
---|
[1803] | 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.