Context Navigation

Square.par@ 4883

Last change on this file since 4883 was 4883, checked in by seroussi, 15 years ago
Validation tests ISMIP-HOM
Property svn:executable set to ``*
File size: 910 bytes

Rev	Line
[16]	1	%Ok, start defining model parameters here
	2
[904]	3	disp(' creating thickness');
	4	md.surface=2000-md.xtan(0.1pi/180); %to have z>0
	5	md.bed=md.surface-1000;
	6	md.thickness=md.surface-md.bed;
	7	md.firn_layer=0*ones(md.numberofgrids,1);
[16]	8
[904]	9	disp(' creating drag');
	10	md.drag_type=2; %0 none 1 plastic 2 viscous
[3768]	11	md.drag_coefficient=sqrt(md.yts.(1000+1000sin(md.x2pi/max(md.x)))./(md.g(md.rho_icemd.thickness+md.rho_water*md.bed)));
[904]	12	%Take care of iceshelves: no basal drag
	13	pos=find(md.elementoniceshelf);
[3768]	14	md.drag_coefficient(md.elements(pos,:))=0;
	15	md.drag_p=ones(md.numberofelements,1);
	16	md.drag_q=ones(md.numberofelements,1);
[16]	17
[904]	18	disp(' creating flow law paramter');
[3768]	19	md.rheology_B=6.806710^7ones(md.numberofgrids,1);
	20	md.rheology_n=3*ones(md.numberofelements,1);
[16]	21
[904]	22	disp(' boundary conditions for diagnostic model: ');
	23	%Create grid on boundary fist (because wi can not use mesh)
[4883]	24	md=SetIceSheetBC(md);
	25
	26	md.np=3;

Note: See TracBrowser for help on using the repository browser.