Regionaloutput.cpp

Go to the documentation of this file.

 
/*Include files: {{{*/
#ifdef HAVE_CONFIG_H
   #include <config.h>
#else
#error "Cannot compile with HAVE_CONFIG_H symbol! run configure first!"
#endif
 
/*Headers:*/
#include "./classes.h"
#include "./Definition.h"
#include "./Elements/Element.h"
#include "./Elements/Elements.h"
#include "./FemModel.h"
#include "../classes/Params/Parameters.h"
 
/*}}}*/
 
Regionaloutput::Regionaloutput(char* in_name, int in_definitionenum, char* in_outputname, IssmDouble* maskin, int Min){ /*{{{*/
 
   this->definitionenum=in_definitionenum;
   this->outputname = xNew<char>(strlen(in_outputname)+1);
   xMemCpy<char>(this->outputname,in_outputname,strlen(in_outputname)+1);
   this->name = xNew<char>(strlen(in_name)+1);
   xMemCpy<char>(this->name,in_name,strlen(in_name)+1);
 
   this->mask   = xNew<IssmDouble>(Min);
   xMemCpy<IssmDouble>(this->mask, maskin, Min);
 
   this->M=Min;
 
}
/*}}}*/
Regionaloutput::~Regionaloutput(){/*{{{*/
   if(this->name)xDelete(this->name); 
   if(this->outputname)xDelete(this->outputname);
   if(this->mask)xDelete(this->mask);
}
/*}}}*/
 
/*Object virtual function resolutoin: */
Object* Regionaloutput::copy() {/*{{{*/
   Regionaloutput* mf = new Regionaloutput(this->name,this->definitionenum,this->outputname,this->mask,this->M);
   return (Object*) mf;
}
/*}}}*/
void Regionaloutput::DeepEcho(void){/*{{{*/
   this->Echo();
}
/*}}}*/
void Regionaloutput::Echo(void){/*{{{*/
   _printf_(" Regionaloutput: " << this->name << " " << this->definitionenum << "\n");
   _printf_("    outputname enum: " << this->outputname << "Enum\n");
   _printf_("    mask: " << this->mask << "\n");
   _printf_("    M: " << this->M << "\n");
}
/*}}}*/
int Regionaloutput::Id(void){/*{{{*/
   return -1;
}
/*}}}*/
void Regionaloutput::Marshall(char** pmarshalled_data,int* pmarshalled_data_size, int marshall_direction){/*{{{*/
   _error_("not implemented yet!"); 
} 
/*}}}*/
int Regionaloutput::ObjectEnum(void){/*{{{*/
   return RegionaloutputEnum;
}
/*}}}*/
 
/*Definition virtual function resolutoin: */
int Regionaloutput::DefinitionEnum(){/*{{{*/
 
   return this->definitionenum;
}
/*}}}*/
char* Regionaloutput::Name(){/*{{{*/
 
   char* name2=xNew<char>(strlen(this->name)+1);
   xMemCpy(name2,this->name,strlen(this->name)+1);
 
   return name2;
}
/*}}}*/
IssmDouble Regionaloutput::Response(FemModel* femmodel){/*{{{*/
 
   int i;
   IssmDouble val_t=0.;
   IssmDouble all_val_t=0.;
   int outputenum = StringToEnumx(this->outputname);
 
   for(i=0;i<femmodel->elements->Size();i++){
      Element* element=(Element*)femmodel->elements->GetObjectByOffset(i);
      switch(outputenum){
         case GroundedAreaEnum:
            val_t+=element->GroundedArea(this->mask,false);
            break;
         case GroundedAreaScaledEnum:
            val_t+=element->GroundedArea(this->mask,true);
            break;
         case FloatingAreaEnum:
            val_t+=element->FloatingArea(this->mask,false);
            break;
         case FloatingAreaScaledEnum:
            val_t+=element->FloatingArea(this->mask,true);
            break;
         case IceMassEnum:
            val_t+=element->IceMass(this->mask,false);
            break;
         case IceMassScaledEnum:
            val_t+=element->IceMass(this->mask,true);
            break;
         case IceVolumeEnum:
            val_t+=element->IceVolume(this->mask,false);
            break;
         case IceVolumeScaledEnum:
            val_t+=element->IceVolume(this->mask,true);
            break;
         case IceVolumeAboveFloatationEnum:
            val_t+=element->IceVolumeAboveFloatation(this->mask,false);
            break;
         case IceVolumeAboveFloatationScaledEnum:
            val_t+=element->IceVolumeAboveFloatation(this->mask,true);
            break;
         case TotalFloatingBmbEnum:
            val_t+=element->TotalFloatingBmb(this->mask,false);
            break;
         case TotalFloatingBmbScaledEnum:
            val_t+=element->TotalFloatingBmb(this->mask,true);
            break;
         case TotalGroundedBmbEnum:
            val_t+=element->TotalGroundedBmb(this->mask,false);
            break;
         case TotalGroundedBmbScaledEnum:
            val_t+=element->TotalGroundedBmb(this->mask,true);
            break;
         case TotalSmbEnum:
            val_t+=element->TotalSmb(this->mask,false);
            break;
         case TotalSmbScaledEnum:
            val_t+=element->TotalSmb(this->mask,true);
            break;
         default:
            _error_("Regional output type " << this->outputname << " not supported yet!");
      }
   }
 
   ISSM_MPI_Allreduce ( (void*)&val_t,(void*)&all_val_t,1,ISSM_MPI_DOUBLE,ISSM_MPI_SUM,IssmComm::GetComm());
   val_t=all_val_t;
 
   return val_t;
}
/*}}}*/