source: issm/trunk/test/Validation/PattynStokes/test2_icesheet/Square.par@ 4875

%Ok, start defining model parameters here

disp('      creating thickness');
hmin=300;
hmax=1000;
ymin=min(md.y);
ymax=max(md.y);
md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
md.firn_layer=10*ones(md.numberofgrids,1);
md.bed=-md.rho_ice/md.rho_water*md.thickness+10;
md.surface=md.bed+md.thickness;

disp('      creating drag');
md.drag_type=2; %0 none 1 plastic 2 viscous
md.drag_coefficient=200*ones(md.numberofgrids,1); %q=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);

disp('      creating temperature');
md.observed_temperature=(273-20)*ones(md.numberofgrids,1);

disp('      creating flow law paramter');
md.rheology_B=paterson(md.observed_temperature);
md.rheology_n=3*ones(md.numberofelements,1);

%Deal with boundary conditions:
disp('      boundary conditions for diagnostic model: ');
md=SetMarineIceSheetBC(md,'Front.exp');

%Parallel options
md.np=3;
md.time=50;
md.waitonlock=10;