Changeset 22817

Timestamp:

05/30/18 22:33:25 (7 years ago)

Author:

adhikari

Message:

CHG: minor adjustment to let intial slr relax according to SAL

File:

• issm/trunk-jpl/src/c/cores/sealevelrise_core.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/cores/sealevelrise_core.cpp

@@ -15 +15 @@
         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 +81 @@
         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)