Index: /issm/trunk-jpl/test/Par/79North.par
===================================================================
--- /issm/trunk-jpl/test/Par/79North.par	(revision 13674)
+++ /issm/trunk-jpl/test/Par/79North.par	(revision 13675)
@@ -9,6 +9,6 @@
 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,15 +23,14 @@
 %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,6 +40,6 @@
 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;
Index: /issm/trunk-jpl/test/Par/79North.py
===================================================================
--- /issm/trunk-jpl/test/Par/79North.py	(revision 13675)
+++ /issm/trunk-jpl/test/Par/79North.py	(revision 13675)
@@ -0,0 +1,69 @@
+import os.path
+import inspect
+import netCDF4
+from numpy import *
+from verbose import *
+from InterpFromMeshToMesh2d import InterpFromMeshToMesh2d
+from paterson import *
+from SetMarineIceSheetBC import *
+
+#Start defining model parameters here
+
+#Geometry and observation
+f = netCDF4.Dataset('../Data/79North.nc','r')
+x         = reshape(f.variables['x'][:],(-1))
+y         = reshape(f.variables['y'][:],(-1))
+vx= f.variables['vx'][:]
+vy= f.variables['vy'][:]
+index     = f.variables['index'][:]
+surface   = f.variables['surface'][:]
+thickness = f.variables['thickness'][:]
+f.close()
+
+[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)
+md.geometry.bed         = md.geometry.surface-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))
+md.initialization.temperature=md.initialization.temperature
+
+#Friction
+pos=numpy.nonzero(md.mask.elementonfloatingice)
+md.friction.coefficient=50.*ones((md.mesh.numberofvertices,1))
+md.friction.coefficient[md.mesh.elements[pos].astype(int)-1]=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=numpy.nonzero(md.mask.elementonfloatingice)
+md.basalforcings.melting_rate=zeros((md.mesh.numberofvertices,1))
+md.basalforcings.melting_rate[md.mesh.elements[pos].astype(int)-1]=10.
+md.surfaceforcings.mass_balance=15*numpy.ones((md.mesh.numberofvertices,1))
+
+#Numerical parameters
+md.diagnostic.viscosity_overshoot=0.3
+md.prognostic.stabilization=1
+md.thermal.stabilization=1
+md.verbose=verbose(0)
+md.settings.waitonlock=30
+md.timestepping.time_step=1.
+md.timestepping.final_time=3.
+md.diagnostic.restol=0.05
+md.diagnostic.reltol=0.005
+md.steadystate.reltol=0.005
+md.diagnostic.abstol=float('NaN')
+
+#Boundary conditions:
+md=SetMarineIceSheetBC(md)
+pos=numpy.nonzero(md.mesh.vertexonboundary)
+md.balancethickness.spcthickness[pos]=md.geometry.thickness[pos]
+md.prognostic.spcthickness[pos]=md.geometry.thickness[pos]
+
+#Change name so that no test have the same name
+if len(inspect.stack()) > 2:
+	md.miscellaneous.name = os.path.basename(inspect.stack()[2][1]).split('.')[0]