Index: /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp	(revision 26480)
+++ /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp	(revision 26481)
@@ -208,9 +208,9 @@
 	femmodel->parameters->FindParam(&arorder,SmbAutoregressiveOrderEnum);
 	IssmDouble tinit_ar;
-	IssmDouble* beta0                = xNew<IssmDouble>(numbasins);
-	IssmDouble* beta1                = xNew<IssmDouble>(numbasins);
-	IssmDouble* phi                  = xNew<IssmDouble>(numbasins*arorder);
-	IssmDouble* noiseterms           = xNew<IssmDouble>(numbasins);
-	IssmDouble* covmat               = xNew<IssmDouble>(numbasins*numbasins);
+	IssmDouble* beta0      = xNew<IssmDouble>(numbasins);
+	IssmDouble* beta1      = xNew<IssmDouble>(numbasins);
+	IssmDouble* phi        = xNew<IssmDouble>(numbasins*arorder);
+	IssmDouble* noiseterms = xNew<IssmDouble>(numbasins);
+	IssmDouble* covmat     = xNew<IssmDouble>(numbasins*numbasins);
 	femmodel->parameters->FindParam(&tinit_ar,SmbAutoregressionInitialTimeEnum);
 	femmodel->parameters->FindParam(&beta0,&M,SmbBeta0Enum);    _assert_(M==numbasins);
@@ -218,12 +218,17 @@
 	femmodel->parameters->FindParam(&phi,&M,&Nphi,SmbPhiEnum);  _assert_(M==numbasins); _assert_(Nphi==arorder);
 	femmodel->parameters->FindParam(&covmat,&M,&N,SmbCovmatEnum); _assert_(M==numbasins); _assert_(N==numbasins);
-	IssmDouble telapsed_ar = time-tinit_ar; //time elapsed with respect to AR model initial time
+
+	/*time elapsed with respect to AR model initial time*/
+	IssmDouble telapsed_ar = time-tinit_ar; 
+
 	/*Before looping through elements: compute noise term specific to each basin from covmat*/
 	multivariateNormal(&noiseterms,numbasins,0.0,covmat);
+
 	/*Loop over each element to compute SMB at vertices*/
 	for(Object* &object:femmodel->elements->objects){
-		Element* element      = xDynamicCast<Element*>(object); //generate element object
+		Element* element = xDynamicCast<Element*>(object);
 		element->Smbautoregression(isstepforar,arorder,telapsed_ar,beta0,beta1,phi,noiseterms);
 	}
+
 	/*Cleanup*/
 	xDelete<IssmDouble>(beta0);