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

Timestamp:

01/23/19 13:41:28 (6 years ago)

Author:

youngmc3

Message:

NEW: add new class md.frontalforcings

File:

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

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

@@ -1496 +1496 @@
 
 }/*}}}*/
+void FemModel::IcefrontAreax(){/*{{{*/
+
+        int numvertices      = this->GetElementsWidth();
+        int numbasins;
+        IssmDouble* BasinId   = xNew<IssmDouble>(numvertices);
+        this->parameters->FindParam(&numbasins,FrontalForcingsNumberofBasinsEnum);
+        IssmDouble* basin_icefront_area           = xNewZeroInit<IssmDouble>(numbasins);
+
+        for(int basin=1;basin<numbasins+1;basin++){
+                IssmDouble local_icefront_area = 0;
+                IssmDouble total_icefront_area;
+
+                for(int i=0;i<this->elements->Size();i++){
+                        Element* element=xDynamicCast<Element*>(this->elements->GetObjectByOffset(i));
+                        element->GetInputListOnVertices(BasinId,FrontalForcingsBasinIdEnum);
+                        for(int j=0;j<numvertices;j++){
+                                if(BasinId[j]==basin){
+                                        local_icefront_area+=element->GetIcefrontArea();
+                                        break;
+                                }
+                        }
+                }
+                ISSM_MPI_Reduce(&local_icefront_area,&total_icefront_area,1,ISSM_MPI_DOUBLE,ISSM_MPI_SUM,0,IssmComm::GetComm());
+                ISSM_MPI_Bcast(&total_icefront_area,1,ISSM_MPI_DOUBLE,0,IssmComm::GetComm());
+        
+                basin_icefront_area[basin-1]=total_icefront_area;
+        }
+        
+        this->parameters->AddObject(new DoubleVecParam(FrontalForcingsBasinIcefrontAreaEnum,basin_icefront_area,numbasins));
+        
+        xDelete<IssmDouble>(basin_icefront_area);
+}/*}}}*/
 void FemModel::IceMassx(IssmDouble* pM, bool scaled){/*{{{*/
 
@@ -2390 +2422 @@
 }
 /*}}}*/
+void FemModel::RignotMeltParameterizationx(){/*{{{*/
+
+        for(int i=0;i<elements->Size();i++){
+                Element* element=dynamic_cast<Element*>(this->elements->GetObjectByOffset(i));
+                element->RignotMeltParameterization();
+        }
+}
+/*}}}*/
 void FemModel::StrainRateparallelx(){/*{{{*/