source: issm/trunk/test/Par/RoundSheetStaticEISMINT.par@ 9641

Last change on this file since 9641 was 9641, checked in by Mathieu Morlighem, 14 years ago

Added mask object

  • Property svn:executable set to *
File size: 1.9 KB
Line 
1disp(' creating thickness');
2hmin=0.01;
3hmax=2756.7;
4radius=(sqrt((md.x).^2+(md.y).^2));
5radiusmax=max(radius);
6md.thickness=hmin*ones(size(md.x,1),1)+hmax*(4*((1/2)^(4/3)*ones(size(md.x,1),1)-((radius)./(2*radiusmax)).^(4/3))).^(3/8);
7md.bed=0*md.thickness;
8md.surface=md.bed+md.thickness;
9
10disp(' creating drag');
11md.friction.coefficient=20*ones(md.numberofnodes,1); %q=1. %no drag is specified in the analytical solution
12%Take care of iceshelves: no basal drag
13pos=find(md.mask.elementonfloatingice);
14md.friction.coefficient(md.elements(pos,:))=0;
15md.friction.p=ones(md.numberofelements,1);
16md.friction.q=ones(md.numberofelements,1);
17
18disp(' creating temperatures');
19tmin=238.15; %K
20st=1.67*10^-2/1000; %k/m;
21md.temperature=(tmin+st*radius);
22md.basalforcings.geothermalflux=4.2*10^-2*ones(md.numberofnodes,1);
23
24disp(' creating flow law paramter');
25md.materials.rheology_B=6.81*10^(7)*ones(md.numberofnodes,1); %to have the same B as the analytical solution
26md.materials.rheology_n=3*ones(md.numberofelements,1);
27
28disp(' creating surface mass balance');
29smb_max=0.5; %m/yr
30sb=10^-2/1000; %m/yr/m
31rel=450*1000; %m
32md.surfaceforcings.mass_balance=min(smb_max,sb*(rel-radius));
33
34disp(' creating velocities');
35constant=0.3;
36md.vx_obs=constant/2*md.x.*(md.thickness).^-1;
37md.vy_obs=constant/2*md.y.*(md.thickness).^-1;
38md.vel_obs=(sqrt((md.vx_obs).^2+(md.vy_obs).^2));
39md.vx=zeros(md.numberofnodes,1);
40md.vy=zeros(md.numberofnodes,1);
41md.vz=zeros(md.numberofnodes,1);
42md.pressure=zeros(md.numberofnodes,1);
43
44%Deal with boundary conditions:
45disp(' boundary conditions for diagnostic model: ');
46md=SetMarineIceSheetBC(md,'../Exp/RoundFrontEISMINT.exp');
47
48radius=sqrt((md.x).*md.x+(md.y).*md.y);
49pos=find(radius==min(radius));
50md.x(pos)=0; md.y(pos)=0; %the closest node to the center is changed to be exactly at the center
51
52md.spcvx(pos)=0;
53md.spcvy(pos)=0;
54md.spcvz(pos)=0;
Note: See TracBrowser for help on using the repository browser.