source: issm/trunk/test/Par/ISMIPB.par@ 8631

%Ok, start defining model parameters here

disp('      creating thickness');
md.surface=-md.x*tan(0.5*pi/180);
md.bed=md.surface-1000+500*sin(md.x*2*pi/max(md.x));
md.thickness=md.surface-md.bed;
md.firn_layer=0*ones(md.numberofnodes,1);

disp('      creating drag');
md.drag_type=2; %0 none 1 plastic 2 viscous
md.drag_coefficient=200*ones(md.numberofnodes,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 flow law paramter');
md.rheology_B=6.8067*10^7*ones(md.numberofnodes,1);
md.rheology_n=3*ones(md.numberofelements,1);

disp('      boundary conditions for diagnostic model');
%Create node on boundary fist (because we cannot use mesh)
md=SetIceSheetBC(md);