Index: /issm/trunk-jpl/src/c/cores/sealevelrise_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/sealevelrise_core.cpp	(revision 22816)
+++ /issm/trunk-jpl/src/c/cores/sealevelrise_core.cpp	(revision 22817)
@@ -15,4 +15,5 @@
 	Vector<IssmDouble> *Sg_absolute  = NULL; 
 	Vector<IssmDouble> *Sg_eustatic  = NULL; 
+	Vector<IssmDouble> *slr_initial  = NULL; 
 	Vector<IssmDouble> *steric_rate_g  = NULL; 
 	Vector<IssmDouble> *U_radial  = NULL; 
@@ -80,4 +81,11 @@
 	if(isslr){
 		Sg_eustatic=sealevelrise_core_eustatic(femmodel); //generalized eustatic (Farrel and Clark, Eq 4, 1st, 3rd and 4rd terms on the RHS.
+	
+		/* The following is for reproducing Farrell&Clark1976 Fig. 1, and aimed for the workshop! 
+		 * md.slr.sealevel is considered as the "initial sea-level" where 1 m slr is distributed uniformly over the ocean 
+		 * similar strategy may work for computin SAL due to "internal mass distribution" associated with ocean dynamics 
+		 * if it breaks the code, it will be reverted. And, we will strategize how we want to accomodate */
+		GetVectorFromInputsx(&slr_initial,femmodel,SealevelEnum,VertexSIdEnum);
+		Sg_eustatic->AXPY(slr_initial,1); 
 
 		Sg=sealevelrise_core_noneustatic(femmodel,Sg_eustatic); //ocean loading tems  (2nd and 5th terms on the RHS of Farrel and Clark)