Changeset 7831

Timestamp:

04/12/11 07:51:00 (14 years ago)

Author:

seroussi

Message:

ISMIPF some improvements

Location:

issm/trunk/test/NightlyRun

Files:

• test1106.m (modified) (2 diffs)
• test1110.m (modified) (2 diffs)

issm/trunk/test/NightlyRun/test1106.m

@@ -2 +2 @@
 %Pattyn and Payne 2006
 
-L_list={5000};
-%L_list={5000,10000,20000,40000,80000,160000};
+L_list={5000,10000,20000,40000,80000,160000};
 results={};
 
@@ -30 +29 @@
 
         %Compute the diagnostic
+        md.verbose=verbose('convergence',true,'solution',true);
         md=solve(md,'analysis_type',DiagnosticSolutionEnum);
 

issm/trunk/test/NightlyRun/test1110.m

@@ -9 +9 @@
 md=geography(md,'',''); %ice sheet test
 md=parameterize(md,'../Par/ISMIPF.par');
-md=extrude(md,2,1);
+md=extrude(md,3,1);
 
 md=setelementstype(md,'pattyn','all');
@@ -31 +31 @@
 
 md.dt=5;
-md.ndt=2000;
+md.ndt=50;
 
 %Compute the diagnostic
-md.cluster=wilkes;
-md.verbose=10;
-md.rheology_B=md.rheology_B/2;
-%md=solve(md,'analysis_type',DiagnosticSolutionEnum);
-%md=solve(md,'analysis_type',ThermalSolutionEnum);
+md=SetParallel(md,8);
+md.verbose=verbose('convergence',true,'solution',true);
+md.rheology_B=md.rheology_B/2;  %as mu=B/(2 \epsilon_e) and visco=B in matice for linear case
 md=solve(md,'analysis_type',Transient2DSolutionEnum);