Changeset 13675

Timestamp:

10/15/12 14:06:03 (12 years ago)

Author:

Mathieu Morlighem

Message:

NEW: added 79North.py

Location:

issm/trunk-jpl/test/Par

Files:

• 79North.par (modified) (3 diffs)
• 79North.py (added)

issm/trunk-jpl/test/Par/79North.par

@@ -9 +9 @@
 surface=transpose(ncread('../Data/79North.nc','surface'));
 thickness=transpose(ncread('../Data/79North.nc','thickness'));
-md.initialization.vx       =InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y);
-md.initialization.vy       =InterpFromMeshToMesh2d(index,x,y,vy,md.mesh.x,md.mesh.y);
+md.initialization.vx =InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y);
+md.initialization.vy =InterpFromMeshToMesh2d(index,x,y,vy,md.mesh.x,md.mesh.y);
 md.geometry.surface  =InterpFromMeshToMesh2d(index,x,y,surface,md.mesh.x,md.mesh.y);
 md.geometry.thickness=InterpFromMeshToMesh2d(index,x,y,thickness,md.mesh.x,md.mesh.y);
@@ -23 +23 @@
 %Friction
 pos=find(md.mask.elementonfloatingice);
-md.friction.coefficient=50*ones(md.mesh.numberofvertices,1);
-md.friction.coefficient(md.mesh.elements(pos,:))=0;
+md.friction.coefficient=50.*ones(md.mesh.numberofvertices,1);
+md.friction.coefficient(md.mesh.elements(pos,:))=0.;
 md.friction.p=ones(md.mesh.numberofelements,1);
 md.friction.q=ones(md.mesh.numberofelements,1);
 
 %Ice shelf melting and surface mass balance 
+pos=find(md.mask.elementonfloatingice);
 md.basalforcings.melting_rate=zeros(md.mesh.numberofvertices,1);
-pos=zeros(md.mesh.numberofvertices,1);
-pos(md.mesh.elements(find(md.mask.elementonfloatingice),:))=1;
-md.basalforcings.melting_rate(find(pos))=10;
-md.surfaceforcings.mass_balance=15*ones(md.mesh.numberofvertices,1);
+md.basalforcings.melting_rate(md.mesh.elements(pos,:))=10;
+md.surfaceforcings.mass_balance=15.*ones(md.mesh.numberofvertices,1);
 
 %Numerical parameters
@@ -41 +40 @@
 md.verbose=verbose(0);
 md.settings.waitonlock=30;
-md.timestepping.time_step=1;
-md.timestepping.final_time=3;
+md.timestepping.time_step=1.;
+md.timestepping.final_time=3.;
 md.diagnostic.restol=0.05;
 md.diagnostic.reltol=0.005;