Index: /issm/trunk-jpl/src/jl/solve/analyses/stressbalanceanalysis.jl
===================================================================
--- /issm/trunk-jpl/src/jl/solve/analyses/stressbalanceanalysis.jl	(revision 26716)
+++ /issm/trunk-jpl/src/jl/solve/analyses/stressbalanceanalysis.jl	(revision 26717)
@@ -13,9 +13,9 @@
 	for i in 1:md.mesh.numberofvertices
 		if ~isnan(spcvx[i])
-			push!(constraints,Constraint(count,i,1,spcvx[i]))
+			push!(constraints,Constraint(count,i,1,spcvx[i]/md.constants.yts))
 			count+=1
 		end
 		if ~isnan(spcvy[i])
-			push!(constraints,Constraint(count,i,2,spcvy[i]))
+			push!(constraints,Constraint(count,i,2,spcvy[i]/md.constants.yts))
 			count+=1
 		end
Index: /issm/trunk-jpl/src/jl/test.jl
===================================================================
--- /issm/trunk-jpl/src/jl/test.jl	(revision 26717)
+++ /issm/trunk-jpl/src/jl/test.jl	(revision 26717)
@@ -0,0 +1,14 @@
+#!/Applications/Julia-1.6.app/Contents/Resources/julia/bin/julia
+include("issm.jl")
+using .ISSM
+using MAT
+
+#Load model from MATLAB file
+#file = matopen("Pig.mat")
+file = matopen("/Users/mmorligh/Desktop/issmuci/trunk-jpl/examples/Pig2/Pig.mat")
+mat = read(file, "md")
+close(file)
+md=ISSM.model(mat)
+
+#Solve stress balance
+md=ISSM.solve(md,"Stressbalance")