Index: /issm/trunk-jpl/src/c/solutions/hydrology_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutions/hydrology_core.cpp	(revision 14552)
+++ /issm/trunk-jpl/src/c/solutions/hydrology_core.cpp	(revision 14553)
@@ -18,5 +18,5 @@
 
 	/*intermediary*/
-	int        nsteps;
+	int        step,nsteps;
 	int        output_frequency;
 	bool       save_results;
@@ -25,5 +25,5 @@
 	IssmDouble time,dt;
 
-	//first recover parameters common to all solutions
+	/*first recover parameters common to all solutions*/
 	femmodel->parameters->FindParam(&starttime,TimesteppingStartTimeEnum);
 	femmodel->parameters->FindParam(&final_time,TimesteppingFinalTimeEnum);
@@ -31,4 +31,5 @@
 	femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
 	femmodel->parameters->FindParam(&output_frequency,SettingsOutputFrequencyEnum);
+
 	/*first compute slopes: */
 	surfaceslope_core(femmodel);
@@ -41,4 +42,7 @@
 	}
 	else nsteps=reCast<int,IssmDouble>((final_time-starttime)/dt);
+
+	/*initialize: */
+	step=0;
 	time=starttime;
 
@@ -48,5 +52,7 @@
 		if(nsteps)if(VerboseSolution()) _pprintLine_("time step:" << i+1 << "/" << nsteps);
 		time+=dt;
+		step+=1;
 		femmodel->parameters->SetParam(time,TimeEnum);
+		femmodel->parameters->SetParam(step,StepEnum);
 
 		/*Compute hydrology solution: */