source: issm/trunk-jpl/test/Par/SquareNoDyn.par@ 24863

%Start defining model parameters here

%Geometry
md.geometry.thickness=1000.0*ones(md.mesh.numberofvertices,1);
md.geometry.base=zeros(md.mesh.numberofvertices,1);
md.geometry.surface=md.geometry.base+md.geometry.thickness;

%Materials
md.initialization.temperature=(273.-20.)*ones(md.mesh.numberofvertices,1);
md.materials.rheology_B=paterson(md.initialization.temperature);
md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);

%Friction
md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
md.friction.coefficient(find(md.mask.ocean_levelset<0.))=0.;
md.friction.p=ones(md.mesh.numberofelements,1);
md.friction.q=ones(md.mesh.numberofelements,1);

%Some necessary fields to fool checkonsistency
md.initialization.vx=zeros(md.mesh.numberofvertices,1);
md.initialization.vy=zeros(md.mesh.numberofvertices,1);
md.initialization.vz=zeros(md.mesh.numberofvertices,1);
md.initialization.pressure=zeros(md.mesh.numberofvertices,1);

md.stressbalance.spcvx=zeros(md.mesh.numberofvertices,1);
md.stressbalance.spcvy=zeros(md.mesh.numberofvertices,1);
md.stressbalance.spcvz=zeros(md.mesh.numberofvertices,1);

md.stressbalance.referential=NaN(md.mesh.numberofvertices,6);
md.stressbalance.loadingforce=zeros(md.mesh.numberofvertices,3);

md.smb.mass_balance=zeros(md.mesh.numberofvertices,1);

md.basalforcings.groundedice_melting_rate=zeros(md.mesh.numberofvertices,1);
md.basalforcings.floatingice_melting_rate=zeros(md.mesh.numberofvertices,1);


%Numerical parameters
md.verbose=verbose(0);
md.settings.waitonlock=30;
md.groundingline.migration='None';

md.transient=deactivateall(md.transient);

%Change name so that no test have the same name
A=dbstack;
if (length(A)>2), md.miscellaneous.name=A(3).file(1:end-2); end