Changeset 23888 for issm/trunk-jpl/src/c/classes/FemModel.cpp

Timestamp:

04/23/19 09:55:09 (6 years ago)

Author:

Mathieu Morlighem

Message:

CHG: adding more requested output, ice front flux seems to be working, grounding line flux still fails

File:

• issm/trunk-jpl/src/c/classes/FemModel.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/classes/FemModel.cpp

@@ -1546 +1546 @@
 
 }/*}}}*/
+void FemModel::GroundinglineMassFluxx(IssmDouble* pM, bool scaled){/*{{{*/
+
+        IssmDouble local_mass_flux = 0;
+        IssmDouble total_mass_flux;
+
+        for(int i=0;i<this->elements->Size();i++){
+                Element* element=xDynamicCast<Element*>(this->elements->GetObjectByOffset(i));
+                local_mass_flux+=element->GroundinglineMassFlux(scaled);
+        }
+        ISSM_MPI_Reduce(&local_mass_flux,&total_mass_flux,1,ISSM_MPI_DOUBLE,ISSM_MPI_SUM,0,IssmComm::GetComm() );
+        ISSM_MPI_Bcast(&total_mass_flux,1,ISSM_MPI_DOUBLE,0,IssmComm::GetComm());
+
+        /*Assign output pointers: */
+        *pM=total_mass_flux;
+
+}/*}}}*/
 void FemModel::IceMassx(IssmDouble* pM, bool scaled){/*{{{*/
 
@@ -2206 +2222 @@
                                         case GroundedAreaEnum:                   this->GroundedAreax(&double_result,false);             break;
                                         case GroundedAreaScaledEnum:             this->GroundedAreax(&double_result,true);              break;
+                                        case GroundinglineMassFluxEnum:          this->GroundinglineMassFluxx(&double_result,false);    break;
                                         case FloatingAreaEnum:                   this->FloatingAreax(&double_result,false);             break;
                                         case FloatingAreaScaledEnum:             this->FloatingAreax(&double_result,true);              break;
@@ -2403 +2420 @@
                 case IceVolumeAboveFloatationEnum:       this->IceVolumeAboveFloatationx(responses, false); break;
                 case IceVolumeAboveFloatationScaledEnum: this->IceVolumeAboveFloatationx(responses, true); break;
-                case IcefrontMassFluxEnum:               this->IcefrontMassFluxx(responses, true); break;
+                case IcefrontMassFluxEnum:               this->IcefrontMassFluxx(responses, false); break;
                 case GroundedAreaEnum:                   this->GroundedAreax(responses, false); break;
                 case GroundedAreaScaledEnum:             this->GroundedAreax(responses, true); break;