Index: /issm/branches/trunk-larour-NatClimateChange2016/src/c/cores/sealevelrise_core_noneustatic.cpp
===================================================================
--- /issm/branches/trunk-larour-NatClimateChange2016/src/c/cores/sealevelrise_core_noneustatic.cpp	(revision 21395)
+++ /issm/branches/trunk-larour-NatClimateChange2016/src/c/cores/sealevelrise_core_noneustatic.cpp	(revision 21396)
@@ -18,4 +18,5 @@
 
 	Vector<IssmDouble> *Sgo    = NULL; //ocean convolution of the perturbation to gravity potential.
+	Vector<IssmDouble> *Sgo_rot= NULL; // rotational feedback 
 	IssmDouble          Sgo_oceanaverage = 0;  //average of Sgo over the ocean.
 
@@ -27,4 +28,5 @@
 	bool spherical=true;
 	bool converged=true;
+	bool rotation=true;
 	bool verboseconvolution=true;
 	int max_nonlinear_iterations;
@@ -41,4 +43,7 @@
 	femmodel->parameters->FindParam(&eps_abs,SealevelriseAbstolEnum);
 	femmodel->parameters->FindParam(&configuration_type,ConfigurationTypeEnum);
+	
+	/*computational flag: */
+	femmodel->parameters->FindParam(&rotation,SealevelriseRotationEnum);
 
 	/*first, recover lat,long and radius vectors from vertices: */
@@ -75,4 +80,13 @@
 		Sgo->Assemble(); 
 
+		if(rotation){
+			/*call rotational feedback  module: */
+			Sgo_rot = new Vector<IssmDouble>(gsize); Sgo_rot->Assemble();
+			femmodel->SealevelriseRotationalFeedback(Sgo_rot,Sg_old,latitude,longitude,radius); 
+			Sgo_rot->Assemble(); 
+
+			Sgo->AXPY(Sgo_rot,1); 
+		}
+		
 		/*we need to average Sgo over the ocean: RHS term  5 in Eq.4 of Farrel and clarke. Only the elements can do that: */
 		Sgo_oceanaverage=femmodel->SealevelriseOceanAverage(Sgo);