Index: /issm/trunk-jpl/src/c/classes/Elements/Element.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 27583)
+++ /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 27584)
@@ -4237,8 +4237,4 @@
 	const int NLOOP = 10; // default internal iteration for SEMIC (Ruckamp et al. 2018).
 
-	IssmDouble* s=xNew<IssmDouble>(NUM_VERTICES);
-	IssmDouble* s0gcm=xNew<IssmDouble>(NUM_VERTICES);
-	IssmDouble* st=xNew<IssmDouble>(NUM_VERTICES);
-
 	// daily forcing inputs
 	IssmDouble* dailyrainfall   =xNew<IssmDouble>(NUM_VERTICES);
@@ -4252,4 +4248,10 @@
 	IssmDouble* dailytemperature=xNew<IssmDouble>(NUM_VERTICES);
 
+	// inputs: geometry
+	IssmDouble* s=xNew<IssmDouble>(NUM_VERTICES);
+	IssmDouble* s0gcm=xNew<IssmDouble>(NUM_VERTICES);
+	IssmDouble* st=xNew<IssmDouble>(NUM_VERTICES);
+
+	// inputs
 	IssmDouble* tsurf_in        =xNew<IssmDouble>(NUM_VERTICES); 
 	IssmDouble* mask_in         =xNew<IssmDouble>(NUM_VERTICES); 
@@ -4260,5 +4262,5 @@
 	IssmDouble* hsnow_in        =xNew<IssmDouble>(NUM_VERTICES); 
 
-	// daily outputs
+	// outputs
 	IssmDouble* tsurf_out  =xNew<IssmDouble>(NUM_VERTICES); memset(tsurf_out, 0., NUM_VERTICES*sizeof(IssmDouble));
 	IssmDouble* smb_out    =xNew<IssmDouble>(NUM_VERTICES); memset(smb_out, 0., NUM_VERTICES*sizeof(IssmDouble));
@@ -4429,7 +4431,11 @@
 				&tsurf_out[iv], &smb_out[iv], &smbi_out[iv], &smbs_out[iv], &saccu_out[iv], &smelt_out[iv],&albedo_out[iv], &albedo_snow_out[iv], &hsnow_out[iv], &hice_out[iv]);
 
-		/* unit conversion: water -> ice */
+		/* 
+		 unit conversion: water -> ice
+		 w.e. : water equivalenet.
+		 */
 		smbi_out[iv] = smbi_out[iv]*rho_water/rho_ice;
-		smbs_out[iv] = smbs_out[iv]*yts; // m/sec -> m/yr
+		smbs_out[iv] = smbs_out[iv]*yts; // w.e. m/sec -> m/yr
+		sacc_out[iv] = accu_out[iv]*yts; // w.e. m/sec -> m/yr
 	}
 
@@ -4468,5 +4474,5 @@
 	}
 
-	/*clean-up*/
+	/*clean-up {{{*/
 	delete gauss;
 	xDelete<IssmDouble>(dailysnowfall);
@@ -4481,4 +4487,5 @@
 	xDelete<IssmDouble>(dailytemperature);
 
+	/*for outputs*/
 	xDelete<IssmDouble>(tsurf_out);
 	xDelete<IssmDouble>(smb_out);
@@ -4489,5 +4496,8 @@
 	xDelete<IssmDouble>(albedo_out);
 	xDelete<IssmDouble>(albedo_snow_out);
-
+	xDelete<IssmDouble>(hsnow_out)
+	xDelete<IssmDouble>(hice_out)
+
+	/*for inputs*/
 	xDelete<IssmDouble>(hsnow_in);
 	xDelete<IssmDouble>(hice_in);
@@ -4497,7 +4507,10 @@
 	xDelete<IssmDouble>(albedo_snow_in);
 	xDelete<IssmDouble>(tsurf_in);
+
+	/* for inputs:geometry */
 	xDelete<IssmDouble>(s);
 	xDelete<IssmDouble>(st);
 	xDelete<IssmDouble>(s0gcm);
+	/*}}}*/
 }
 /*}}}*/