Changeset 4950

Timestamp:

08/03/10 16:53:50 (15 years ago)

Author:

Mathieu Morlighem

Message:

First test

Location:

issm/trunk/test

Files:

• Archives/Archive201.mat (added)
• Data/SquareShelf.data (added)
• NightlyRun/Id2Name.m (added)
• NightlyRun/runme.m (added)
• NightlyRun/test201.m (added)
• NightlyRun/test201_nightly.m (added)
• Par/SquareShelf.par (modified) (2 diffs)
• Validation/ISMIP/TestA/Square.par (modified) (1 diff)

issm/trunk/test/Par/SquareShelf.par

@@ -1 +1 @@
 %Start defining model parameters here
 
-%dynamics
-md.verbose=0;
-md.artificial_diffusivity=1;
-
+%Geometry
 hmin=300;
 hmax=1000;
@@ -13 +10 @@
 md.surface=md.bed+md.thickness;
 
-md.drag_type=2; %0 none 1 plastic 2 viscous
-md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1.
-%Take care of iceshelves: no basal drag
-pos=find(md.elementoniceshelf);
-md.drag_coefficient(md.elements(pos,:))=0;
-md.drag_p=ones(md.numberofelements,1);
-md.drag_q=ones(md.numberofelements,1);
-md.viscosity_overshoot=0.3;
+%Initial velocity 
+load('./../Data/SquareShelf.data','-mat');
+md.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.x,md.y,0);
+md.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y,0);
+clear vx vy x y index;
 
-%Create temperature and rheology
+%Materials
 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
 md.rheology_B=paterson(md.observed_temperature);
 md.rheology_n=3*ones(md.numberofelements,1);
 
-%Deal with boundary conditions:
-md=SetIceShelfBC(md,'Front.exp');
+%Friction
+pos=find(md.elementoniceshelf);
+md.drag_coefficient=20*ones(md.numberofgrids,1);
+md.drag_coefficient(md.elements(pos,:))=0;
+md.drag_p=ones(md.numberofelements,1);
+md.drag_q=ones(md.numberofelements,1);
 
-%Parallel options
+%Numerical parameters
+md.viscosity_overshoot=0.3;
+md.verbose=0;
+md.artificial_diffusivity=1;
 md.np=3;
-md.time=50;
 md.waitonlock=30;
+
+%Boundary conditions:
+md=SetIceShelfBC(md,'./../Exp/SquareFront.exp');

issm/trunk/test/Validation/ISMIP/TestA/Square.par

@@ -24 +24 @@
 md=SetIceSheetBC(md);
 
-md.np=3;
+md.np=10;
+md.cluster=oshostname();
+md.verbose=10;