Changeset 9725 for issm/trunk/test/Par/Pig.par

Timestamp:

09/09/11 08:35:15 (14 years ago)

Author:

Mathieu Morlighem

Message:

Some more objects in mesh

File:

• issm/trunk/test/Par/Pig.par (modified) (1 diff)

issm/trunk/test/Par/Pig.par

@@ -11 +11 @@
 md.initialization.vx=md.inversion.vx_obs;
 md.initialization.vy=md.inversion.vy_obs;
-md.initialization.vz=zeros(md.numberofnodes,1);
-md.initialization.pressure=zeros(md.numberofnodes,1);
+md.initialization.vz=zeros(md.mesh.numberofvertices,1);
+md.initialization.pressure=zeros(md.mesh.numberofvertices,1);
 
 %Materials
-md.initialization.temperature=(273-20)*ones(md.numberofnodes,1);
+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.numberofelements,1);
+md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
 md.initialization.temperature=md.initialization.temperature;
 
 %Friction
 pos=find(md.mask.elementonfloatingice);
-md.friction.coefficient=50*ones(md.numberofnodes,1);
+md.friction.coefficient=50*ones(md.mesh.numberofvertices,1);
 md.friction.coefficient(md.elements(pos,:))=0;
-md.friction.p=ones(md.numberofelements,1);
-md.friction.q=ones(md.numberofelements,1);
+md.friction.p=ones(md.mesh.numberofelements,1);
+md.friction.q=ones(md.mesh.numberofelements,1);
 
 %Numerical parameters