source: issm/trunk/test/Par/Pig.par@ 9681

Last change on this file since 9681 was 9681, checked in by Mathieu Morlighem, 14 years ago

Added inversion

  • Property svn:executable set to *
File size: 1.5 KB
RevLine 
[4949]1%Start defining model parameters here
2
[4953]3%Geometry and observation
[8824]4load('../Data/Pig.data','-mat');
[9681]5md.inversion.vx_obs =InterpFromMeshToMesh2d(index,x,y,vx_obs,md.x,md.y);
6md.inversion.vy_obs =InterpFromMeshToMesh2d(index,x,y,vy_obs,md.x,md.y);
[5035]7md.surface =InterpFromMeshToMesh2d(index,x,y,surface,md.x,md.y);
8md.thickness=InterpFromMeshToMesh2d(index,x,y,thickness,md.x,md.y);
[4953]9md.bed=md.surface-md.thickness;
10clear surface thickness vx_obs vy_obs x y index;
[9681]11md.vx=md.inversion.vx_obs;
12md.vy=md.inversion.vy_obs;
[8305]13md.vz=zeros(md.numberofnodes,1);
14md.pressure=zeros(md.numberofnodes,1);
[4949]15
[4953]16%Materials
[9425]17md.temperature=(273-20)*ones(md.numberofnodes,1);
[9636]18md.materials.rheology_B=paterson(md.temperature);
19md.materials.rheology_n=3*ones(md.numberofelements,1);
[9425]20md.temperature=md.temperature;
[4953]21
22%Friction
[9641]23pos=find(md.mask.elementonfloatingice);
[9610]24md.friction.coefficient=50*ones(md.numberofnodes,1);
25md.friction.coefficient(md.elements(pos,:))=0;
26md.friction.p=ones(md.numberofelements,1);
27md.friction.q=ones(md.numberofelements,1);
[4953]28
29%Numerical parameters
[9679]30md.diagnostic.viscosity_overshoot=0.3;
[9642]31md.prognostic.stabilization=1;
[6307]32md.verbose=verbose(0);
[4953]33md.waitonlock=30;
[9628]34md.timestepping.time_step=1;
35md.timestepping.final_time=2;
[9679]36md.diagnostic.restol=0.05;
37md.diagnostic.reltol=1;
[9677]38md.steadystate.reltol=1;
[9679]39md.diagnostic.abstol=NaN;
[4949]40
[4953]41%Boundary conditions:
[4949]42md=SetMarineIceSheetBC(md);
[4999]43
44%Change name so that no test have the same name
45A=dbstack;
[9625]46if (length(A)>2), md.miscellaneous.name=A(3).file(1:end-2); end
Note: See TracBrowser for help on using the repository browser.