source: issm/trunk/test/Validation/ThermalTests/Simpleadvection/Square.par@ 900

Last change on this file since 900 was 900, checked in by Mathieu Morlighem, 16 years ago

simplified Square.par

  • Property svn:executable set to *
File size: 1.5 KB
RevLine 
[16]1%Ok, start defining model parameters here
2
[900]3%parallelization
4md.cluster='none';
[16]5
[900]6disp(' creating thickness');
7h=1000;
8md.thickness=h*ones(md.numberofgrids,1);
9md.firn_layer=10*ones(md.numberofgrids,1);
10md.bed=-1000*ones(md.numberofgrids,1);
11md.surface=md.bed+md.thickness;
[16]12
[900]13disp(' creating velocities');
14md.vx_obs=zeros(md.numberofgrids,1);
15md.vy_obs=zeros(md.numberofgrids,1);
16md.vel_obs=sqrt(md.vx_obs.^2+md.vy_obs.^2);
[16]17
[900]18disp(' creating drag');
19md.drag_type=2; %0 none 1 plastic 2 viscous
20md.drag=200*ones(md.numberofgrids,1); %q=1.
21%Take care of iceshelves: no basal drag
22pos=find(md.elementoniceshelf);
23md.drag(md.elements(pos,:))=0;
24md.p=ones(md.numberofelements,1);
25md.q=ones(md.numberofelements,1);
[16]26
[900]27disp(' creating temperatures');
28md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
[16]29
[900]30disp(' creating flow law paramter');
31md.B=paterson(md.observed_temperature);
32md.n=3*ones(md.numberofelements,1);
33
34disp(' creating accumulation rates');
35md.accumulation=ones(md.numberofgrids,1)/md.yts; %1m/a
36md.melting=0*ones(md.numberofgrids,1)/md.yts; %1m/a
37
38%Deal with boundary conditions:
39
40disp(' boundary conditions for diagnostic model');
41md=SetIceShelfBC(md,'Front.exp');
42
43disp(' boundary conditions for thermal model');
44md.gridondirichlet_thermal=ones(md.numberofgrids,1); %surface temperature
45md.dirichletvalues_thermal=md.observed_temperature;
46md.geothermalflux=zeros(md.numberofgrids,1);
47pos=find(md.elementonicesheet);md.geothermalflux(md.elements(pos,:))=50*10^-3; %50 mW/m^2
Note: See TracBrowser for help on using the repository browser.