Changeset 25287

Timestamp:

07/17/20 13:42:53 (5 years ago)

Author:

Eric.Larour

Message:

CHG: missing bottom pressure in geodetic routine.

File:

• issm/trunk-jpl/src/c/classes/Elements/Tria.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -5945 +5945 @@
         /*diverse:*/
         int gsize;
-        IssmDouble I, S;                //change in relative ice thickness and sea level
+        IssmDouble I, S, BP;            //change in relative ice thickness and sea level
         IssmDouble rho_ice,rho_water;
         int horiz;
@@ -5987 +5987 @@
                 /*convert to kg/m^2:*/
                 S=rho_water*S;
+                
+                /*If bottom pressures are available, retrieve them to load the bedrock:*/
+                if(bp_compute_fingerprints){
+                        Input2* bottompressure_change_input=this->GetInput2(DslSeaWaterPressureChangeAtSeaFloor);
+                        if (!bottompressure_change_input)_error_("bottom pressure input needed to compute sea level rise fingerprint!");
+                        bottompressure_change_input->GetInputAverage(&BP);
+
+                        /*convert from m to kg/m^2:*/
+                        BP=BP*rho_water;
+
+                        S+=BP;
+                }
 
                 for(int i=0;i<gsize;i++){