SurfaceAreax.cpp File Reference

#include "./SurfaceAreax.h"
#include "../../shared/shared.h"
#include "../../toolkits/toolkits.h"
#include "../InputUpdateFromConstantx/InputUpdateFromConstantx.h"

Go to the source code of this file.

Functions
void	SurfaceAreax (IssmDouble *pS, FemModel *femmodel)

Function Documentation

SurfaceAreax()

void SurfaceAreax	(	IssmDouble *	pS,
		FemModel *	femmodel
	)

Definition at line 11 of file SurfaceAreax.cpp.

                                                    {
 
   /*Intermediary*/
   Element* element=NULL;
 
   /*output: */
   IssmDouble S = 0.;
   IssmDouble S_sum;
 
   /*Compute gradients: */
   for(int i=0;i<femmodel->elements->Size();i++){
      element=xDynamicCast<Element*>(femmodel->elements->GetObjectByOffset(i));
      S+=element->SurfaceArea();
   }
 
   /*Sum all J from all cpus of the cluster:*/
   ISSM_MPI_Reduce (&S,&S_sum,1,ISSM_MPI_DOUBLE,ISSM_MPI_SUM,0,IssmComm::GetComm() );
   ISSM_MPI_Bcast(&S_sum,1,ISSM_MPI_DOUBLE,0,IssmComm::GetComm()); 
   S=S_sum;
 
   /*add surface area to element inputs:*/
   InputUpdateFromConstantx(femmodel,S,SurfaceAreaEnum);
 
   /*Assign output pointers: */
   if(pS) *pS=S;
}