Changeset 27829 for issm/trunk-jpl/src/c/classes/Elements/Element.cpp

Timestamp:

07/16/23 18:22:25 (20 months ago)

Author:

inwoo

Message:

CHG: add calculation of TotalSmbMelt and TotalSmbRefreeze for SEMIC.

File:

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

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

@@ -4318 +4318 @@
                 _printf0_("smb core: dt       : " << dt <<"\n");
         }
-        /* loop over vertices and days */ //FIXME account for leap years (365 -> 366)
+        /* loop over vertices and days */
         Gauss* gauss=this->NewGauss();
         /* Retrieve inputs: */
@@ -4347 +4347 @@
                 enum_qmr        =SmbSemicQmrInitEnum;
         } 
-        if(isverbose && this->Sid()==0)_printf0_("smb core: assign temp.\n");
+        //if(isverbose && this->Sid()==0)_printf0_("smb core: assign temp.\n");
         Input* tsurf_input       = this->GetInput(enum_temp); _assert_(tsurf_in);
-        if(isverbose && this->Sid()==0)_printf0_("smb core: assign mask.\n");
+        //if(isverbose && this->Sid()==0)_printf0_("smb core: assign mask.\n");
         Input* mask_input        = this->GetInput(SmbMaskEnum); _assert_(mask_input);
-        if(isverbose && this->Sid()==0)_printf0_("smb core: assign Tamp.\n");
+        //if(isverbose && this->Sid()==0)_printf0_("smb core: assign Tamp.\n");
         Input* Tamp_input        = this->GetInput(SmbTampEnum); _assert_(Tamp_input);
-        if(isverbose && this->Sid()==0)_printf0_("smb core: assign albedo.\n");
+        //if(isverbose && this->Sid()==0)_printf0_("smb core: assign albedo.\n");
         Input* albedo_input      = this->GetInput(enum_albedo); _assert_(albedo_input);
         Input* albedo_snow_input = this->GetInput(enum_albedo_snow); _assert_(albedo_snow_input);
         Input* hice_input        = this->GetInput(enum_hice); _assert_(hice_input);
         Input* hsnow_input       = this->GetInput(enum_hsnow); _assert_(hsnow_input);
-        if(isverbose && this->Sid()==0)_printf0_("smb core: assign qmr.\n");
         Input* qmr_input         = this->GetInput(enum_qmr); _assert_(qmr_input);
 
@@ -4438 +4437 @@
                  */
                 smb_out[iv]  = smb_out[iv]*yts;  // w.e. m/sec -> m/yr
-                smbi_out[iv] = smbi_out[iv]*rho_water/rho_ice;
+                smbi_out[iv] = smbi_out[iv]*rho_water/rho_ice; // w.e. m/sec -> ice m/yr
                 smbs_out[iv] = smbs_out[iv]*yts; // w.e. m/sec -> m/yr
                 saccu_out[iv] = saccu_out[iv]*yts; // w.e. m/sec -> m/yr
-        }
-        /*
-         * unit conversion is not required for smelt and refr_out.
-         */
+                smelt_out[iv] = smelt_out[iv]*rho_water/rho_ice; // w.e. m/sec -> ice m/yr
+                refr_out[iv]  = refr_out[iv]*rho_water/rho_ice; // w.e. m/sec -> ice m/yr
+        }
 
         if(isverbose && this->Sid()==0){
@@ -5406 +5404 @@
 }
 /*}}}*/
+IssmDouble Element::TotalSmbMelt(IssmDouble* mask, bool scaled){/*{{{*/
+
+        /*Retrieve values of the mask defining the element: */
+        for(int i=0;i<this->GetNumberOfVertices();i++){
+                if(mask[this->vertices[i]->Sid()]<=0.){
+                        return 0.;
+                }
+        }
+
+        /*Return: */
+        return this->TotalSmbMelt(scaled);
+}
+/*}}}*/
+IssmDouble Element::TotalSmbRefreeze(IssmDouble* mask, bool scaled){/*{{{*/
+
+        /*Retrieve values of the mask defining the element: */
+        for(int i=0;i<this->GetNumberOfVertices();i++){
+                if(mask[this->vertices[i]->Sid()]<=0.){
+                        return 0.;
+                }
+        }
+
+        /*Return: */
+        return this->TotalSmbRefreeze(scaled);
+}
+/*}}}*/
 void       Element::TransformInvStiffnessMatrixCoord(ElementMatrix* Ke,int transformenum){/*{{{*/