Changeset 4953 for issm/trunk/test/Par/79North.par
- Timestamp:
- 08/03/10 17:23:56 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/trunk/test/Par/79North.par ΒΆ
r4949 r4953 1 1 %Start defining model parameters here 2 2 3 %Create geometry 4 load('79Northdata.mat'); 5 md.surface=surface; 6 md.bed=bed; 7 md.thickness=thickness; 8 md.vx=vx; 9 md.vy=vy; 3 %Geometry and observation 4 load('./../Data/79North.data','-mat'); 5 md.vx =InterpFromMeshToMesh2d(index,x,y,vx,md.x,md.y,0); 6 md.vy =InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y,0); 7 md.surface =InterpFromMeshToMesh2d(index,x,y,surface,md.x,md.y,0); 8 md.thickness=InterpFromMeshToMesh2d(index,x,y,thickness,md.x,md.y,0); 9 md.bed=md.surface-md.thickness; 10 clear surface thickness vx vy x y index; 10 11 11 %Put melting on ice shelves 12 %Materials 13 md.observed_temperature=(273-20)*ones(md.numberofgrids,1); 14 md.rheology_B=paterson(md.observed_temperature); 15 md.rheology_n=3*ones(md.numberofelements,1); 16 17 %Friction 18 pos=find(md.elementoniceshelf); 19 md.drag_type=2; 20 md.drag_coefficient=50*ones(md.numberofgrids,1); 21 md.drag_coefficient(md.elements(pos,:))=0; 22 md.drag_p=ones(md.numberofelements,1); 23 md.drag_q=ones(md.numberofelements,1); 24 25 %Ice shelf melting 12 26 md.melting_rate=zeros(md.numberofgrids,1); 13 27 pos=find(md.gridoniceshelf); 14 28 md.melting_rate(pos)=10; 15 29 16 %Create drag 17 md.drag_type=2; 18 md.drag_coefficient=50*ones(md.numberofgrids,1); 19 %Take care of iceshelves: no basal drag 20 pos=find(md.elementoniceshelf); 21 md.drag_coefficient(md.elements(pos,:))=0; 22 md.drag_p=ones(md.numberofelements,1); 23 md.drag_q=ones(md.numberofelements,1); 30 %Numerical parameters 24 31 md.viscosity_overshoot=0.3; 32 md.artificial_diffusivity=1; 33 md.verbose=0; 34 md.waitonlock=30; 35 md.np=3; 36 md.eps_res=0.05; 37 md.eps_rel=0.005; 38 md.eps_abs=NaN; 25 39 26 %Create temperature and rheology 27 md.observed_temperature=(273-20)*ones(md.numberofgrids,1); 28 md.rheology_B=paterson(md.observed_temperature); 29 md.rheology_n=3*ones(md.numberofelements,1); 30 31 %Deal with boundary conditions: 40 %Boundary conditions: 32 41 md=SetMarineIceSheetBC(md); 33 42 pos=find(md.gridonboundary); 34 43 md.spcthickness(pos,1)=1; 35 44 md.spcthickness(pos,2)=md.thickness(pos); 36 37 %Parallel options38 md.np=7;39 md.time=50;40 md.waitonlock=90;41 42 %Misc43 md.verbose=0;44 md.eps_abs=NaN;45 md.eps_rel=0.05;
Note:
See TracChangeset
for help on using the changeset viewer.