Changeset 22817
- Timestamp:
- 05/30/18 22:33:25 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/cores/sealevelrise_core.cpp
r21908 r22817 15 15 Vector<IssmDouble> *Sg_absolute = NULL; 16 16 Vector<IssmDouble> *Sg_eustatic = NULL; 17 Vector<IssmDouble> *slr_initial = NULL; 17 18 Vector<IssmDouble> *steric_rate_g = NULL; 18 19 Vector<IssmDouble> *U_radial = NULL; … … 80 81 if(isslr){ 81 82 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); 82 90 83 91 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.