source: issm/trunk/test/Par/RoundSheetEISMINT.par@ 9684

Last change on this file since 9684 was 9684, checked in by seroussi, 14 years ago

added initialization class

  • Property svn:executable set to *
File size: 2.1 KB
Line 
1%Ok, start defining model parameters here
2disp(' creating thickness');
3md.thickness=10*ones(md.numberofnodes,1);
4md.bed=zeros(md.numberofnodes,1);
5md.surface=md.bed+md.thickness;
6
7disp(' creating drag');
8md.friction.coefficient=20*ones(md.numberofnodes,1); %q=1. %no drag is specified in the analytical solution
9md.friction.p=ones(md.numberofelements,1);
10md.friction.q=ones(md.numberofelements,1);
11
12disp(' creating temperatures');
13tmin=238.15; %K
14st=1.67*10^-2/1000; %k/m;
15radius=sqrt((md.x).^2+(md.y).^2);
16md.initialization.temperature=(tmin+st*radius);
17md.basalforcings.geothermalflux=4.2*10^-2*ones(md.numberofnodes,1);
18
19disp(' creating flow law paramter');
20md.materials.rheology_B=6.81*10^(7)*ones(md.numberofnodes,1); %to have the same B as the analytical solution
21md.materials.rheology_n=3*ones(md.numberofelements,1);
22
23disp(' creating surface mass balance');
24smb_max=0.5; %m/yr
25sb=10^-2/1000; %m/yr/m
26rel=450*1000; %m
27md.surfaceforcings.mass_balance=min(smb_max,sb*(rel-radius));
28
29disp(' creating velocities');
30constant=0.3;
31md.inversion.vx_obs=constant/2*md.x.*(md.thickness).^-1;
32md.inversion.vy_obs=constant/2*md.y.*(md.thickness).^-1;
33md.inversion.vel_obs=(sqrt((md.inversion.vx_obs).^2+(md.inversion.vy_obs).^2));
34md.initialization.vx=zeros(md.numberofnodes,1);
35md.initialization.vy=zeros(md.numberofnodes,1);
36md.initialization.vz=zeros(md.numberofnodes,1);
37md.initialization.pressure=zeros(md.numberofnodes,1);
38
39%Deal with boundary conditions:
40disp(' boundary conditions for diagnostic model: ');
41md=SetMarineIceSheetBC(md,'../Exp/RoundFrontEISMINT.exp');
42
43radius=sqrt((md.x).*md.x+(md.y).*md.y);
44pos=find(radius==min(radius));
45md.x(pos)=0; md.y(pos)=0; %the closest node to the center is changed to be exactly at the center
46
47md.diagnostic.spcvx(pos)=0;
48md.diagnostic.spcvy(pos)=0;
49md.diagnostic.spcvz(pos)=0;
50
51%parallel options
52md.timestepping.final_time=50000;
53
54%Constants
55md.materials.rho_ice=910;
56md.materials.thermalconductivity=2.1;
57md.materials.latentheat=3.35*10^5;
58md.materials.beta=8.66*10^-4/(md.materials.rho_ice*md.constants.g); %conversion from K/m to K/Pa
59md.constants.yts=31556926;
Note: See TracBrowser for help on using the repository browser.