Index: /issm/trunk-jpl/src/c/cores/stressbalance_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/stressbalance_core.cpp	(revision 23336)
+++ /issm/trunk-jpl/src/c/cores/stressbalance_core.cpp	(revision 23337)
@@ -75,4 +75,8 @@
 	/*Compute vertical velocities*/
 	if (domaintype==Domain3DEnum && (isSIA || isSSA || isL1L2 || isHO)){
+
+		/*We need basal melt rates for vertical velocity*/
+		bmb_core(femmodel);
+
 		analysis = new StressbalanceVerticalAnalysis();
 		analysis->Core(femmodel);