Index: /issm/trunk-jpl/src/m/classes/diagnostic.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/diagnostic.py	(revision 14529)
+++ /issm/trunk-jpl/src/m/classes/diagnostic.py	(revision 14530)
@@ -138,5 +138,5 @@
 		md = checkfield(md,'diagnostic.maxiter','size',[1],'>=',1)
 		md = checkfield(md,'diagnostic.referential','size',[md.mesh.numberofvertices,6])
-		md = checkfield(md,'diagnostic.loadingforce','size',[md.mesh.loadingforce,3])
+		md = checkfield(md,'diagnostic.loadingforce','size',[md.mesh.numberofvertices,3])
 		if not md.diagnostic.requested_outputs:
 			md = checkfield(md,'diagnostic.requested_outputs','size',[float('NaN'),1])
Index: /issm/trunk-jpl/src/m/classes/model/model.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/model/model.py	(revision 14529)
+++ /issm/trunk-jpl/src/m/classes/model/model.py	(revision 14530)
@@ -634,4 +634,5 @@
 		md.balancethickness.spcthickness=project3d(md,'vector',md.balancethickness.spcthickness,'type','node')
 		md.diagnostic.referential=project3d(md,'vector',md.diagnostic.referential,'type','node')
+		md.diagnostic.loadingforce=project3d(md,'vector',md.diagnostic.loadingforce,'type','node')
 
 		#in 3d, pressureload: [node1 node2 node3 node4 element]