Changeset 1803

Timestamp:

08/24/09 14:16:35 (16 years ago)

Author:

Eric.Larour

Message:

Some adaptation of Square.par to new rift formulation

Location:

issm/trunk/test/Verification/IceShelfIceFrontM2dRifts_24

Files:

• Square.par (modified) (1 diff)
• buildme.m (modified) (2 diffs)

issm/trunk/test/Verification/IceShelfIceFrontM2dRifts_24/Square.par

@@ -36 +36 @@
 md.time=50;
 md.waitonlock=1;
-                
-                
-md=solversettoasm(md);
+
+%solver
+md.solverstring=' -mat_type aijmumps -ksp_type preonly -pc_type lu -mat_mumps_icntl_14 80 ';
+
+%rifts
+md.rifts.fill=MelangeEnum();
+md.rifts.fraction=0;
+md.penalty_lock=2;
+md.rifts.fractionincrement=.01;

issm/trunk/test/Verification/IceShelfIceFrontM2dRifts_24/buildme.m

@@ -1 +1 @@
 %mesh resolution
-resolution=50000;
-riftres=1000;
+resolution=10000;
+riftres=500;
 
 md=model;
@@ -11 +11 @@
 md=mesh(md,'DomainOutline.exp','Rifts2.exp',resolution);
 md=meshprocessrifts(md);
+md.numberofelements
 md=geography(md,'all','');
 md=parameterize(md,'Square.par');
 md=setelementstype(md,'macayeal','all');
+
+md.rifts.fill=MelangeEnum();
+md.rifts.fraction=0;
+md.solverstring=' -mat_type aijmumps -ksp_type preonly -pc_type lu -mat_mumps_icntl_14 80 ';
+md.cluster='astrid';
+md.np=8;
+md.penalty_lock=2;
+md.rifts.fractionincrement=.01;
+md=solve(md,'analysis_type','diagnostic');md=tres(md,'diagnostic');
+plotmodel(md,'data','riftrelvel','data','riftvel','data','vel');