Changeset 19062

Timestamp:

02/02/15 14:46:12 (10 years ago)

Author:

Eric.Larour

Message:

CHG: added IceMass computation.

Location:

issm/trunk-jpl/src/c/classes

Files:

• FemModel.cpp (modified) (3 diffs)
• FemModel.h (modified) (1 diff)

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

@@ -694 +694 @@
                 case DivergenceEnum:               this->Divergencex(responses); break;
                 case MaxDivergenceEnum:            this->MaxDivergencex(responses); break;
+                case IceMassEnum:                  this->IceMassx(responses); break;
                 case IceVolumeEnum:                this->IceVolumex(responses); break;
                 case IceVolumeAboveFloatationEnum: this->IceVolumeAboveFloatationx(responses); break;
@@ -777 +778 @@
                                         case DivergenceEnum:               this->Divergencex(&double_result);               break;
                                         case MaxDivergenceEnum:            this->MaxDivergencex(&double_result);            break;
+                                        case IceMassEnum:                this->IceMassx(&double_result);                break;
                                         case IceVolumeEnum:                this->IceVolumex(&double_result);                break;
                                         case IceVolumeAboveFloatationEnum: this->IceVolumeAboveFloatationx(&double_result); break;
@@ -1394 +1396 @@
 
 }/*}}}*/
+void FemModel::IceMassx(IssmDouble* pM){/*{{{*/
+
+        IssmDouble local_ice_mass = 0;
+        IssmDouble total_ice_mass;
+
+        for(int i=0;i<this->elements->Size();i++){
+                Element* element=xDynamicCast<Element*>(this->elements->GetObjectByOffset(i));
+                local_ice_mass+=element->IceMass();
+        }
+        ISSM_MPI_Reduce(&local_ice_mass,&total_ice_mass,1,ISSM_MPI_DOUBLE,ISSM_MPI_SUM,0,IssmComm::GetComm() );
+        ISSM_MPI_Bcast(&total_ice_mass,1,ISSM_MPI_DOUBLE,0,IssmComm::GetComm());
+
+        /*Assign output pointers: */
+        *pM=total_ice_mass;
+
+}/*}}}*/
 void FemModel::IceVolumex(IssmDouble* pV){/*{{{*/
 

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

@@ -79 +79 @@
                 void Divergencex(IssmDouble* pdiv);
                 void MaxDivergencex(IssmDouble* pdiv);
+                void IceMassx(IssmDouble* pV);
                 void IceVolumex(IssmDouble* pV);
                 void IceVolumeAboveFloatationx(IssmDouble* pV);