Changeset 23978 for issm/trunk-jpl/test/NightlyRun/test355.m

Timestamp:

05/31/19 20:08:38 (6 years ago)

Author:

Mathieu Morlighem

Message:

CHG: updated NR test

File:

• issm/trunk-jpl/test/NightlyRun/test355.m (modified) (5 diffs)

issm/trunk-jpl/test/NightlyRun/test355.m

@@ -1 +1 @@
 %Test Name: SquareSheetHydrologyGlaDS
-cluster=generic('name',oshostname(),'np',1);
 
 %create model: 
@@ -11 +10 @@
 md=setmask(md,'',''); %everywhere grounded
 md=setflowequation(md,'SSA','all');
+md.stressbalance.maxiter=2; %Make sure it runs quickly...
 
 %Some constants
 md.constants.g=9.8;
 md.materials.rho_ice=910;
-md.materials.rho_water=1000;
 
 %Geometry
@@ -31 +30 @@
 
 %Materials
-md.materials.rheology_B    = paterson(md.initialization.temperature);
-md.materials.rheology_B(:) = (5e-25)^(-1/3);
-md.materials.rheology_n    = 3.*ones(md.mesh.numberofelements,1);
+md.materials.rheology_B = (5e-25)^(-1/3) * ones(md.mesh.numberofvertices,1);
+md.materials.rheology_n = 3.*ones(md.mesh.numberofelements,1);
 
 %Friction
@@ -39 +37 @@
 md.friction.p=ones(md.mesh.numberofelements,1);
 md.friction.q=ones(md.mesh.numberofelements,1);
-md.friction.coupling=0;
+%md.friction.coupling=0;
 
 %Boundary conditions:
@@ -66 +64 @@
 md.hydrology.spcphi(pos) = md.materials.rho_ice * md.constants.g * md.geometry.thickness(pos);
 
-md.cluster=cluster;
+md.cluster=generic('np',2);
 md=solve(md,'Transient'); %or 'tr'