Changeset 25278

Timestamp:

07/14/20 17:15:53 (5 years ago)

Author:

Eric.Larour

Message:

CHG: do not include BP into RSL, as it's present in the dynamic sea level part!

File:

• issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp (modified) (3 diffs)

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

@@ -449 +449 @@
         Vector<IssmDouble> *RSLg    = NULL;
         Vector<IssmDouble> *RSLg_old    = NULL;
-        Vector<IssmDouble> *BPg    = NULL;
 
         Vector<IssmDouble> *RSLgo    = NULL; //ocean convolution of the perturbation to gravity potential.
@@ -467 +466 @@
         IssmDouble           eustatic;
         IssmDouble                      Ixz, Iyz, Izz; 
-        int bp_compute_fingerprints= 0;
         
         if(VerboseSolution()) _printf0_("         converging on ocean components\n");
@@ -557 +555 @@
         
         
-        /*if we had bottom pressure loading, add dynamic sea level 
-         * to RSL:*/
-        femmodel->parameters->FindParam(&bp_compute_fingerprints,DslComputeFingerprintsEnum);
-        if(bp_compute_fingerprints){
-                GetVectorFromInputsx(&BPg,femmodel,DslSeaWaterPressureChangeAtSeaFloor,VertexSIdEnum);
-                RSLg->AXPY(BPg,1);
-        }
-
         delete RSLg_old;
-        delete BPg;
 
         return RSLg;