Changeset 4953 for issm/trunk/test/Par/79North.par

Timestamp:

08/03/10 17:23:56 (15 years ago)

Author:

Mathieu Morlighem

Message:

updated par files

File:

• issm/trunk/test/Par/79North.par (modified) (1 diff)

issm/trunk/test/Par/79North.par

@@ -1 +1 @@
 %Start defining model parameters here
 
-%Create geometry
-load('79Northdata.mat');
-md.surface=surface;
-md.bed=bed;
-md.thickness=thickness;
-md.vx=vx;
-md.vy=vy;
+%Geometry and observation
+load('./../Data/79North.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);
+md.surface  =InterpFromMeshToMesh2d(index,x,y,surface,md.x,md.y,0);
+md.thickness=InterpFromMeshToMesh2d(index,x,y,thickness,md.x,md.y,0);
+md.bed=md.surface-md.thickness;
+clear surface thickness vx vy x y index;
 
-%Put melting on ice shelves
+%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);
+
+%Friction
+pos=find(md.elementoniceshelf);
+md.drag_type=2;
+md.drag_coefficient=50*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);
+
+%Ice shelf melting
 md.melting_rate=zeros(md.numberofgrids,1);
 pos=find(md.gridoniceshelf);
 md.melting_rate(pos)=10;
 
-%Create drag
-md.drag_type=2;
-md.drag_coefficient=50*ones(md.numberofgrids,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);
+%Numerical parameters
 md.viscosity_overshoot=0.3;
+md.artificial_diffusivity=1;
+md.verbose=0;
+md.waitonlock=30;
+md.np=3;
+md.eps_res=0.05;
+md.eps_rel=0.005;
+md.eps_abs=NaN;
 
-%Create temperature and rheology
-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:
+%Boundary conditions:
 md=SetMarineIceSheetBC(md);
 pos=find(md.gridonboundary);
 md.spcthickness(pos,1)=1;
 md.spcthickness(pos,2)=md.thickness(pos);
-
-%Parallel options
-md.np=7;
-md.time=50;
-md.waitonlock=90;
-
-%Misc
-md.verbose=0;
-md.eps_abs=NaN;
-md.eps_rel=0.05;