Index: /issm/trunk-jpl/src/c/analyses/SamplingAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/SamplingAnalysis.cpp	(revision 26009)
+++ /issm/trunk-jpl/src/c/analyses/SamplingAnalysis.cpp	(revision 26010)
@@ -76,5 +76,5 @@
   iomodel->FetchDataToInput(inputs,elements,"md.sampling.kappa",SamplingKappaEnum);
 	iomodel->FetchDataToInput(inputs,elements,"md.sampling.beta",SamplingBetaEnum,0.);
-	//iomodel->FetchDataToInput(inputs,elements,"md.initialization.sample",SampleEnum,0.);
+	iomodel->FetchDataToInput(inputs,elements,"md.initialization.sample",SampleEnum,0.);
 
 }/*}}}*/
Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 26009)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 26010)
@@ -853,5 +853,5 @@
 		case TransientSolutionEnum:{
 			/*We have multiple analyses here, process one by one*/
-			bool isSIA,isFS,isthermal,isenthalpy,ismasstransport,isgroundingline,isstressbalance,ismovingfront,ishydrology,isdamage,issmb,isslc,isesa,isgia;
+			bool isSIA,isFS,isthermal,isenthalpy,ismasstransport,isgroundingline,isstressbalance,ismovingfront,ishydrology,isdamage,issmb,isslc,isesa,isgia,issampling;
 			iomodel->FindConstant(&isthermal,"md.transient.isthermal");
 			iomodel->FindConstant(&ismovingfront,"md.transient.ismovingfront");
@@ -865,6 +865,7 @@
 			iomodel->FindConstant(&isesa,"md.transient.isesa");
 			iomodel->FindConstant(&isgia,"md.transient.isgia");
-         int* analyses_iter     = NULL;
-         int  num_analyses_iter = 0;
+			iomodel->FindConstant(&issampling,"md.transient.issampling");
+      int* analyses_iter     = NULL;
+      int  num_analyses_iter = 0;
 			if(isstressbalance){
 				SolutionAnalysesList(&analyses_iter,&num_analyses_iter,iomodel,StressbalanceSolutionEnum);
@@ -911,4 +912,7 @@
 			if(isgia){
 				analyses_temp[numanalyses++]=GiaAnalysisEnum;
+			}
+			if(issampling){
+				analyses_temp[numanalyses++]=SamplingAnalysisEnum;
 			}
 
Index: /issm/trunk-jpl/src/c/cores/transient_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/transient_core.cpp	(revision 26009)
+++ /issm/trunk-jpl/src/c/cores/transient_core.cpp	(revision 26010)
@@ -138,5 +138,5 @@
 
 	/*parameters: */
-	bool isstressbalance,ismasstransport,issmb,isthermal,isgroundingline,isgia,isesa;
+	bool isstressbalance,ismasstransport,issmb,isthermal,isgroundingline,isgia,isesa,issampling;;
 	bool isslc,ismovingfront,isdamageevolution,ishydrology,isoceancoupling,save_results;
 	int  step,sb_coupling_frequency;
@@ -160,4 +160,5 @@
 	femmodel->parameters->FindParam(&isdamageevolution,TransientIsdamageevolutionEnum);
 	femmodel->parameters->FindParam(&ishydrology,TransientIshydrologyEnum);
+	femmodel->parameters->FindParam(&issampling,TransientIssamplingEnum);
 	femmodel->parameters->FindParam(&numoutputs,TransientNumRequestedOutputsEnum);
 
@@ -239,4 +240,10 @@
 		if(VerboseSolution()) _printf0_("   computing sea level change\n");
 		sealevelchange_core(femmodel);
+	}
+
+	/*Sampling: */
+	if(issampling){
+		if(VerboseSolution()) _printf0_("   computing Gaussian random field\n");
+		sampling_core(femmodel);
 	}
 
Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp	(revision 26009)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp	(revision 26010)
@@ -86,5 +86,6 @@
 		parameters->AddObject(iomodel->CopyConstantObject("md.transient.isoceancoupling",TransientIsoceancouplingEnum));
 		parameters->AddObject(iomodel->CopyConstantObject("md.transient.amr_frequency",TransientAmrFrequencyEnum));
-
+		parameters->AddObject(iomodel->CopyConstantObject("md.transient.issampling",TransientIssamplingEnum));
+		
 		/*For stress balance only*/
 		parameters->AddObject(iomodel->CopyConstantObject("md.flowequation.isFS",FlowequationIsFSEnum));