Index: /issm/branches/trunk-larour-SLPS2022/src/c/analyses/SealevelchangeAnalysis.cpp
===================================================================
--- /issm/branches/trunk-larour-SLPS2022/src/c/analyses/SealevelchangeAnalysis.cpp	(revision 27074)
+++ /issm/branches/trunk-larour-SLPS2022/src/c/analyses/SealevelchangeAnalysis.cpp	(revision 27075)
@@ -42,5 +42,4 @@
 	/*external solidearthsolution: solid-Earth model*/
 	iomodel->FetchData(&isexternal,"md.solidearth.isexternal");
-
 	if(isexternal){
 
@@ -632,15 +631,12 @@
 	/*Resolve Mmes using the modelid, if necessary: meaning if we are running a transient model and that UQ computations have not been triggered:*/
 	femmodel->parameters->FindParam(&isuq,QmuIsdakotaEnum);
-	_printf_("hello0\n");
 	if (!isuq && femmodel->inputs->GetInputObjectEnum(SolidearthExternalDisplacementEastRateEnum)==DatasetInputEnum){
-		_printf_("hello1\n");
-		
 		femmodel->parameters->FindParam(&modelid,SolidearthExternalModelidEnum);
 
 		/*replace dataset of forcings with only one, the modelid'th:*/
-		MmeToInputFromIdx(femmodel->inputs,femmodel->elements,modelid-1,SolidearthExternalDisplacementNorthRateEnum, P1Enum);
-		MmeToInputFromIdx(femmodel->inputs,femmodel->elements,modelid-1,SolidearthExternalDisplacementEastRateEnum, P1Enum);
-		MmeToInputFromIdx(femmodel->inputs,femmodel->elements,modelid-1,SolidearthExternalDisplacementUpRateEnum, P1Enum);
-		MmeToInputFromIdx(femmodel->inputs,femmodel->elements,modelid-1,SolidearthExternalGeoidRateEnum, P1Enum);
+		MmeToInputFromIdx(femmodel->inputs,femmodel->elements,femmodel->parameters, modelid-1,SolidearthExternalDisplacementNorthRateEnum, P1Enum);
+		MmeToInputFromIdx(femmodel->inputs,femmodel->elements,femmodel->parameters,modelid-1,SolidearthExternalDisplacementEastRateEnum, P1Enum);
+		MmeToInputFromIdx(femmodel->inputs,femmodel->elements,femmodel->parameters,modelid-1,SolidearthExternalDisplacementUpRateEnum, P1Enum);
+		MmeToInputFromIdx(femmodel->inputs,femmodel->elements,femmodel->parameters,modelid-1,SolidearthExternalGeoidRateEnum, P1Enum);
 
 	}		
Index: /issm/branches/trunk-larour-SLPS2022/src/c/classes/FemModel.cpp
===================================================================
--- /issm/branches/trunk-larour-SLPS2022/src/c/classes/FemModel.cpp	(revision 27074)
+++ /issm/branches/trunk-larour-SLPS2022/src/c/classes/FemModel.cpp	(revision 27075)
@@ -2112,5 +2112,5 @@
 void FemModel::MmeToInputFromId(int id, int rootenum, int interpolationenum){ /*{{{*/
 
-	MmeToInputFromIdx(this->inputs,this->elements,id,rootenum,interpolationenum);
+	MmeToInputFromIdx(this->inputs,this->elements,this->parameters,id,rootenum,interpolationenum);
 
 }	//}}}
Index: /issm/branches/trunk-larour-SLPS2022/src/c/cores/sealevelchange_core.cpp
===================================================================
--- /issm/branches/trunk-larour-SLPS2022/src/c/cores/sealevelchange_core.cpp	(revision 27074)
+++ /issm/branches/trunk-larour-SLPS2022/src/c/cores/sealevelchange_core.cpp	(revision 27075)
@@ -205,11 +205,4 @@
 		TransferForcing(femmodel,DeltaTwsEnum);
 
-		/*transfer external forcings back to Earth:*/
-		/*TransferSealevel(femmodel,BedEnum);
-		TransferSealevel(femmodel,SealevelEnum);
-		if(horiz){
-			TransferSealevel(femmodel,BedEastEnum);
-			TransferSealevel(femmodel,BedNorthEnum);
-		}*/
 	}
 
@@ -272,5 +265,5 @@
 	/*only run if grd was requested, if we are the earth, and we have reached
 	 * the necessary number of time steps dictated by :*/
-	if(!grd)            return;
+	if(!grd)return; 
 	if(count!=frequency)return;
 	femmodel->parameters->FindParam(&iscoupling,IsSlcCouplingEnum);
@@ -606,4 +599,34 @@
 	femmodel->results->AddResult(new GenericExternalResult<IssmDouble*>(femmodel->results->Size()+1,AreaeEnum,areae,nel,1,1,1));
 	#endif
+
+	//Specialty code: {{{
+	bool isuq=0;
+	IssmDouble modelid=0;
+		 
+	/*Resolve Mmes using the modelid, if necessary: meaning if we are running a transient model and that UQ computations have not been triggered:*/
+	femmodel->parameters->FindParam(&isuq,QmuIsdakotaEnum);
+	if(!isuq){
+
+		/*Deal with solid earth: */
+		if (femmodel->inputs->Exist(SolidearthExternalDisplacementEastRateEnum) && femmodel->inputs->GetInputObjectEnum(SolidearthExternalDisplacementEastRateEnum)==DatasetInputEnum){
+			femmodel->parameters->FindParam(&modelid,SolidearthExternalModelidEnum);
+
+			/*replace dataset of forcings with only one, the modelid'th:*/
+			MmeToInputFromIdx(femmodel->inputs,femmodel->elements,femmodel->parameters,modelid-1,SolidearthExternalDisplacementNorthRateEnum, P1Enum);
+			MmeToInputFromIdx(femmodel->inputs,femmodel->elements,femmodel->parameters,modelid-1,SolidearthExternalDisplacementEastRateEnum, P1Enum);
+			MmeToInputFromIdx(femmodel->inputs,femmodel->elements,femmodel->parameters,modelid-1,SolidearthExternalDisplacementUpRateEnum, P1Enum);
+			MmeToInputFromIdx(femmodel->inputs,femmodel->elements,femmodel->parameters,modelid-1,SolidearthExternalGeoidRateEnum, P1Enum);
+
+		}
+		/*Deal with ocean: */
+		if (femmodel->inputs->Exist(OceantransportSpcbottompressureEnum) && femmodel->inputs->GetInputObjectEnum(OceantransportSpcbottompressureEnum)==DatasetInputEnum){
+			femmodel->parameters->FindParam(&modelid,DslModelidEnum);
+
+			/*replace dataset of forcings with only one, the modelid'th:*/
+			MmeToInputFromIdx(femmodel->inputs,femmodel->elements,femmodel->parameters,modelid-1,OceantransportSpcbottompressureEnum, P1Enum);
+			MmeToInputFromIdx(femmodel->inputs,femmodel->elements,femmodel->parameters,modelid-1,OceantransportSpcdslEnum, P1Enum);
+			MmeToInputFromIdx(femmodel->inputs,femmodel->elements,femmodel->parameters,modelid-1,OceantransportSpcstrEnum, P0Enum);
+		}
+	} /*}}}*/
 
 	return;
Index: /issm/branches/trunk-larour-SLPS2022/src/c/modules/MmeToInputFromIdx/MmeToInputFromIdx.cpp
===================================================================
--- /issm/branches/trunk-larour-SLPS2022/src/c/modules/MmeToInputFromIdx/MmeToInputFromIdx.cpp	(revision 27074)
+++ /issm/branches/trunk-larour-SLPS2022/src/c/modules/MmeToInputFromIdx/MmeToInputFromIdx.cpp	(revision 27075)
@@ -16,5 +16,5 @@
 #include "./MmeToInputFromIdx.h"
 
-void MmeToInputFromIdx(Inputs* inputs, Elements* elements, int id, int rootenum, int interpolationenum){
+void MmeToInputFromIdx(Inputs* inputs, Elements* elements, Parameters* parameters, int id, int rootenum, int interpolationenum){
 
 	TransientInput* transientinput  = NULL;
@@ -33,5 +33,5 @@
 	transientinput->GetAllTimes(&times,&N);
 	inputs->SetTransientInput(DummyEnum,times,N);
-	transientinput2 = inputs->GetTransientInput(DummyEnum);
+	transientinput2 = inputs->GetTransientInput(DummyEnum); transientinput2->Configure(parameters);
 
 	for(Object* & object : elements->objects){
Index: /issm/branches/trunk-larour-SLPS2022/src/c/modules/MmeToInputFromIdx/MmeToInputFromIdx.h
===================================================================
--- /issm/branches/trunk-larour-SLPS2022/src/c/modules/MmeToInputFromIdx/MmeToInputFromIdx.h	(revision 27074)
+++ /issm/branches/trunk-larour-SLPS2022/src/c/modules/MmeToInputFromIdx/MmeToInputFromIdx.h	(revision 27075)
@@ -7,5 +7,5 @@
 
 /* local prototypes: */
-void  MmeToInputFromIdx(Inputs* inputs, Elements* elements, int id, int rootenum, int interpolationenum);
+void  MmeToInputFromIdx(Inputs* inputs, Elements* elements, Parameters* parameters, int id, int rootenum, int interpolationenum);
 
 #endif