Changeset 22817


Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/cores/sealevelrise_core.cpp

    r21908 r22817  
    1515        Vector<IssmDouble> *Sg_absolute  = NULL;
    1616        Vector<IssmDouble> *Sg_eustatic  = NULL;
     17        Vector<IssmDouble> *slr_initial  = NULL;
    1718        Vector<IssmDouble> *steric_rate_g  = NULL;
    1819        Vector<IssmDouble> *U_radial  = NULL;
     
    8081        if(isslr){
    8182                Sg_eustatic=sealevelrise_core_eustatic(femmodel); //generalized eustatic (Farrel and Clark, Eq 4, 1st, 3rd and 4rd terms on the RHS.
     83       
     84                /* The following is for reproducing Farrell&Clark1976 Fig. 1, and aimed for the workshop!
     85                 * md.slr.sealevel is considered as the "initial sea-level" where 1 m slr is distributed uniformly over the ocean
     86                 * similar strategy may work for computin SAL due to "internal mass distribution" associated with ocean dynamics
     87                 * if it breaks the code, it will be reverted. And, we will strategize how we want to accomodate */
     88                GetVectorFromInputsx(&slr_initial,femmodel,SealevelEnum,VertexSIdEnum);
     89                Sg_eustatic->AXPY(slr_initial,1);
    8290
    8391                Sg=sealevelrise_core_noneustatic(femmodel,Sg_eustatic); //ocean loading tems  (2nd and 5th terms on the RHS of Farrel and Clark)
Note: See TracChangeset for help on using the changeset viewer.