%Start defining model parameters here

disp('      creating geometry');
if md.numberofgrids~=1619,
	error('Pig.par error message: Model can not be parameterize if numgrids != 1619');
end
load('Pigdata.mat');
md.surface=surface;
md.bed=bed;
md.thickness=thickness;
md.vx_obs=vx_obs;
md.vy_obs=vy_obs;

disp('      creating drag');
md.drag_type=2;
md.drag_coefficient=50*ones(md.numberofgrids,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);
md.viscosity_overshoot=0.3;

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:
md=SetMarineIceSheetBC(md);

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

%Misc
md.verbose=0;
md.eps_abs=NaN;
md.eps_rel=0.05;

%control methods
md.control_type='drag_coefficient';
md.cm_min=2;
md.cm_max=120;
md.nsteps=2;
md.fit=0*ones(md.nsteps,1);
md.weights=ones(md.numberofgrids,1);
md.optscal=150*ones(md.nsteps,1);
md.maxiter=10*ones(md.nsteps,1);
md.cm_jump=0.3*ones(md.nsteps,1);