| 1 | /*!\file: Element.h
|
|---|
| 2 | * \brief abstract class for Element object
|
|---|
| 3 | * This class is a place holder for the Tria and the Penta elements.
|
|---|
| 4 | * It is derived from Element, so DataSets can contain them.
|
|---|
| 5 | */
|
|---|
| 6 |
|
|---|
| 7 | #ifndef _ELEMENT_H_
|
|---|
| 8 | #define _ELEMENT_H_
|
|---|
| 9 |
|
|---|
| 10 | /*Headers:*/
|
|---|
| 11 | /*{{{*/
|
|---|
| 12 | #include "../../datastructures/datastructures.h"
|
|---|
| 13 | #include "../../toolkits/toolkits.h"
|
|---|
| 14 | class DataSet;
|
|---|
| 15 | class Parameters;
|
|---|
| 16 | class Parameter;
|
|---|
| 17 | class Elements;
|
|---|
| 18 | class Loads;
|
|---|
| 19 | class Nodes;
|
|---|
| 20 | class Node;
|
|---|
| 21 | class Vertices;
|
|---|
| 22 | class Vertex;
|
|---|
| 23 | class Materials;
|
|---|
| 24 | class Material;
|
|---|
| 25 | class Matlitho;
|
|---|
| 26 | class Inputs;
|
|---|
| 27 | class Inputs;
|
|---|
| 28 | class Input;
|
|---|
| 29 | class Input;
|
|---|
| 30 | class ElementInput;
|
|---|
| 31 | class DatasetInput;
|
|---|
| 32 | class IoModel;
|
|---|
| 33 | class SealevelGeometry;
|
|---|
| 34 | class Gauss;
|
|---|
| 35 | class GrdLoads;
|
|---|
| 36 | class ElementVector;
|
|---|
| 37 | template <class doublematrix> class Matrix;
|
|---|
| 38 | template <class doubletype> class Vector;
|
|---|
| 39 | class ElementMatrix;
|
|---|
| 40 | class ElementVector;
|
|---|
| 41 | class BarystaticContributions;
|
|---|
| 42 | /*}}}*/
|
|---|
| 43 |
|
|---|
| 44 | class Element: public Object{
|
|---|
| 45 |
|
|---|
| 46 | public:
|
|---|
| 47 | int id;
|
|---|
| 48 | int sid;
|
|---|
| 49 | int lid;
|
|---|
| 50 | Inputs *inputs;
|
|---|
| 51 | Node **nodes;
|
|---|
| 52 | Vertex **vertices;
|
|---|
| 53 | Material *material;
|
|---|
| 54 | Parameters *parameters;
|
|---|
| 55 | bool isonsurface;
|
|---|
| 56 | bool isonbase;
|
|---|
| 57 |
|
|---|
| 58 | int* element_type_list;
|
|---|
| 59 | int element_type;
|
|---|
| 60 |
|
|---|
| 61 | public:
|
|---|
| 62 | /*Constructors/Destructores*/
|
|---|
| 63 | Element();
|
|---|
| 64 | ~Element();
|
|---|
| 65 |
|
|---|
| 66 | /*Functions*/
|
|---|
| 67 | /*bool AllActive(void);*/
|
|---|
| 68 | /*bool AnyActive(void);*/
|
|---|
| 69 | bool AnyFSet(void);
|
|---|
| 70 | void ArmaProcess_pre18Oct2022(bool isstepforarma,int arorder,int maorder,IssmDouble telapsed_arma,IssmDouble tstep_arma,IssmDouble* termconstant,IssmDouble* trend,IssmDouble* arlagcoefs,IssmDouble* malagcoefs,bool isfieldstochastic,int enum_type);
|
|---|
| 71 | void ArmaProcess(bool isstepforarma,int arorder,int maorder,int numparams,int numbreaks,IssmDouble tstep_arma,IssmDouble* polyparams,IssmDouble* arlagcoefs,IssmDouble* malagcoefs,IssmDouble* datebreaks,bool isfieldstochastic,int enum_type);
|
|---|
| 72 | void BasinLinearFloatingiceMeltingRate(IssmDouble* deepwaterel,IssmDouble* upperwatermelt,IssmDouble* upperwaterel,IssmDouble* perturbation);
|
|---|
| 73 | void CalvingSetZeroRate(void);
|
|---|
| 74 | void CalvingRateToVector(void);
|
|---|
| 75 | void ComputeLambdaS(void);
|
|---|
| 76 | void ComputeNewDamage();
|
|---|
| 77 | void ComputeStrainRate();
|
|---|
| 78 | void CoordinateSystemTransform(IssmDouble** ptransform,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 79 | void DeepEcho();
|
|---|
| 80 | void DeleteMaterials(void);
|
|---|
| 81 | void Delta18oParameterization(void);
|
|---|
| 82 | void Delta18opdParameterization(void);
|
|---|
| 83 | void SmbGradCompParameterization(void);
|
|---|
| 84 | IssmDouble Divergence(void);
|
|---|
| 85 | void dViscositydBFS(IssmDouble* pdmudB,int dim,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input);
|
|---|
| 86 | void dViscositydBHO(IssmDouble* pdmudB,int dim,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 87 | void dViscositydBSSA(IssmDouble* pdmudB,int dim,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 88 | void dViscositydBMOLHO(IssmDouble* pdmudB,int dim,IssmDouble* xyz_list,Gauss* gauss,Input* vxbase_input,Input* vybase_input, Input* vxshear_input ,Input* vyshear_input,Input* thickness_input,Input* n_input, IssmDouble zeta);
|
|---|
| 89 | void dViscositydDSSA(IssmDouble* pdmudB,int dim,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 90 | void Echo();
|
|---|
| 91 | void FindParam(bool* pvalue,int paramenum);
|
|---|
| 92 | void FindParam(int* pvalue,int paramenum);
|
|---|
| 93 | void FindParam(IssmDouble* pvalue,int paramenum);
|
|---|
| 94 | IssmDouble FindParam(int paramenum);
|
|---|
| 95 | void FindParam(int** pvalues,int* psize,int paramenum);
|
|---|
| 96 | IssmDouble FloatingArea(IssmDouble* mask, bool scaled);
|
|---|
| 97 | void GetDofList(int** pdoflist,int approximation_enum,int setenum);
|
|---|
| 98 | void GetDofListPressure(int** pdoflist,int setenum);
|
|---|
| 99 | void GetDofListVelocity(int** pdoflist,int setenum);
|
|---|
| 100 | void GetDofListLocal(int** pdoflist,int approximation_enum,int setenum);
|
|---|
| 101 | void GetDofListLocalPressure(int** pdoflist,int setenum);
|
|---|
| 102 | void GetDofListLocalVelocity(int** pdoflist,int setenum);
|
|---|
| 103 | void GetInputListOnNodes(IssmDouble* pvalue,int enumtype);
|
|---|
| 104 | void GetInputListOnNodes(IssmDouble* pvalue,int enumtype,IssmDouble defaultvalue);
|
|---|
| 105 | void GetInputListOnNodesVelocity(IssmDouble* pvalue,int enumtype);
|
|---|
| 106 | void GetInputListOnVertices(IssmDouble* pvalue,int enumtype);
|
|---|
| 107 | void GetInputListOnVerticesAtTime(IssmDouble* pvalue,int enumtype,IssmDouble time);
|
|---|
| 108 | void GetInputListOnVertices(IssmDouble* pvalue,int enumtype,IssmDouble defaultvalue);
|
|---|
| 109 | void GetInputLocalMinMaxOnNodes(IssmDouble* min,IssmDouble* max,IssmDouble* ug);
|
|---|
| 110 | void GetInputValue(bool* pvalue,int enum_type);
|
|---|
| 111 | void GetInputValue(int* pvalue,int enum_type);
|
|---|
| 112 | void GetInputValue(IssmDouble* pvalue,int enum_type);
|
|---|
| 113 | void GetInputValue(IssmDouble* pvalue,Gauss* gauss,int enum_type);
|
|---|
| 114 | Node* GetNode(int nodeindex);
|
|---|
| 115 | int GetNodeIndex(Node* node);
|
|---|
| 116 | void GetNodesLidList(int* lidlist);
|
|---|
| 117 | void GetNodesSidList(int* sidlist);
|
|---|
| 118 | void GetPhi(IssmDouble* phi, IssmDouble* epsilon, IssmDouble viscosity);
|
|---|
| 119 | void GetSolutionFromInputsOneDof(Vector<IssmDouble>* solution,int solutionenum);
|
|---|
| 120 | void GetVectorFromInputs(Vector<IssmDouble>* vector, int name_enum, int type);
|
|---|
| 121 | void GetVectorFromInputs(Vector<IssmDouble>* vector, int name_enum, int type,IssmDouble time);
|
|---|
| 122 | void GetVerticesLidList(int* lidlist);
|
|---|
| 123 | void GetVerticesPidList(int* pidlist);
|
|---|
| 124 | void GetVerticesConnectivityList(int* connectivitylist);
|
|---|
| 125 | void GetVerticesCoordinates(IssmDouble** xyz_list);
|
|---|
| 126 | void GetVerticesSidList(int* sidlist);
|
|---|
| 127 | IssmDouble GetXcoord(IssmDouble* xyz_list,Gauss* gauss);
|
|---|
| 128 | IssmDouble GetYcoord(IssmDouble* xyz_list,Gauss* gauss);
|
|---|
| 129 | IssmDouble GetZcoord(IssmDouble* xyz_list,Gauss* gauss);
|
|---|
| 130 | void GradientIndexing(int* indexing,int control_index);
|
|---|
| 131 | IssmDouble GroundedArea(IssmDouble* mask, bool scaled);
|
|---|
| 132 | bool HasNodeOnBase();
|
|---|
| 133 | bool HasNodeOnSurface();
|
|---|
| 134 | IssmDouble IceMass(bool scaled);
|
|---|
| 135 | IssmDouble IceMass(IssmDouble* mask, bool scaled);
|
|---|
| 136 | IssmDouble IceVolume(IssmDouble* mask, bool scaled);
|
|---|
| 137 | IssmDouble IceVolumeAboveFloatation(IssmDouble* mask, bool scaled);
|
|---|
| 138 | int Id();
|
|---|
| 139 | void InputCreate(IssmDouble* vector,Inputs* inputs,IoModel* iomodel,int M,int N,int vector_type,int vector_enum,int code);
|
|---|
| 140 | void InputCreateP1FromConstant(Inputs* inputs,IoModel* iomodel,IssmDouble value,int vector_enum);
|
|---|
| 141 | void ControlInputCreate(IssmDouble* doublearray,IssmDouble* independents_min,IssmDouble* independents_max,Inputs*inputs,IoModel* iomodel,int M,int N,IssmDouble scale,int input_enum,int id);
|
|---|
| 142 | void DatasetInputAdd(int enum_type,IssmDouble* vector,Inputs* inputs,IoModel* iomodel,int M,int N,int vector_type,int vector_enum,int input_enum);
|
|---|
| 143 | void InputUpdateFromConstant(IssmDouble constant, int name);
|
|---|
| 144 | void InputUpdateFromConstant(int constant, int name);
|
|---|
| 145 | void InputUpdateFromConstant(bool constant, int name);
|
|---|
| 146 | void InputUpdateFromConstant(IssmDouble constant, int name, int type);
|
|---|
| 147 |
|
|---|
| 148 | bool IsAllFloating();
|
|---|
| 149 | bool IsAllGrounded();
|
|---|
| 150 | bool IsGrounded();
|
|---|
| 151 | bool IsOnBase();
|
|---|
| 152 | bool IsOnSurface();
|
|---|
| 153 | bool IsIceInElement();
|
|---|
| 154 | bool IsIceOnlyInElement();
|
|---|
| 155 | bool IsOceanInElement();
|
|---|
| 156 | bool IsOceanOnlyInElement();
|
|---|
| 157 | bool IsAllMinThicknessInElement();
|
|---|
| 158 | bool IsLandInElement();
|
|---|
| 159 | void Ismip6FloatingiceMeltingRate();
|
|---|
| 160 | void LapseRateBasinSMB(int numelevbins, IssmDouble* lapserates, IssmDouble* elevbins,IssmDouble* refelevation);
|
|---|
| 161 | void LinearFloatingiceMeltingRate();
|
|---|
| 162 | void SpatialLinearFloatingiceMeltingRate();
|
|---|
| 163 | void MantlePlumeGeothermalFlux();
|
|---|
| 164 | void MarshallElement2(MarshallHandle* marshallhandle,int numanalyses);
|
|---|
| 165 | void MigrateGroundingLine(IssmDouble* sheet_ungrounding);
|
|---|
| 166 | void MismipFloatingiceMeltingRate();
|
|---|
| 167 | void MonthlyFactorBasin(IssmDouble* monthlyfac,int enum_type);
|
|---|
| 168 | void MonthlyPiecewiseLinearEffectBasin(int nummonthbreaks,IssmDouble* monthlyintercepts,IssmDouble* monthlytrends,IssmDouble* monthlydatebreaks,int enum_type);
|
|---|
| 169 | void BeckmannGoosseFloatingiceMeltingRate();
|
|---|
| 170 | void MungsmtpParameterization(void);
|
|---|
| 171 | ElementMatrix* NewElementMatrix(int approximation_enum=NoneApproximationEnum);
|
|---|
| 172 | ElementMatrix* NewElementMatrixCoupling(int number_nodes,int approximation_enum=NoneApproximationEnum);
|
|---|
| 173 | ElementVector* NewElementVector(int approximation_enum=NoneApproximationEnum);
|
|---|
| 174 | void PicoUpdateBoxid(int* pmax_boxid_basin);
|
|---|
| 175 | void PicoUpdateBox(int loopboxid);
|
|---|
| 176 | void PicoComputeBasalMelt();
|
|---|
| 177 | void PositiveDegreeDay(IssmDouble* pdds,IssmDouble* pds,IssmDouble signorm,bool ismungsm,bool issetpddfac);
|
|---|
| 178 | void PositiveDegreeDaySicopolis(bool isfirnwarming);
|
|---|
| 179 | void RignotMeltParameterization();
|
|---|
| 180 | void ResultInterpolation(int* pinterpolation,int*nodesperelement,int* parray_size, int output_enum);
|
|---|
| 181 | void ResultToPatch(IssmDouble* values,int nodesperelement,int output_enum);
|
|---|
| 182 | void ResultToMatrix(IssmDouble* values,int ncols,int output_enum);
|
|---|
| 183 | void ResultToVector(Vector<IssmDouble>* vector,int output_enum);
|
|---|
| 184 | void SetwiseNodeConnectivity(int* d_nz,int* o_nz,Node* node,bool* flags,int* flagsindices,int* flagsindices_counter,int set1_enum,int set2_enum, int analysis_type);
|
|---|
| 185 | void SetBoolInput(Inputs* inputs,int enum_in,bool value);
|
|---|
| 186 |
|
|---|
| 187 | void SetIntInput(Inputs* inputs,int enum_in,int value);
|
|---|
| 188 | void SmbSemic();
|
|---|
| 189 | void SmbSemicTransient();
|
|---|
| 190 | int Sid();
|
|---|
| 191 | void SmbGemb(IssmDouble timeinputs, int count);
|
|---|
| 192 | void StrainRateESA(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 193 | void StrainRateFS(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input);
|
|---|
| 194 | void StrainRateHO(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 195 | void StrainRateHO2dvertical(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 196 | void StrainRateMOLHO(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vxshear_input,Input* vyshear_input,Input* thickness_input,Input* n_input,IssmDouble zeta);
|
|---|
| 197 | void StrainRateSSA(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 198 | void StrainRateSSA1d(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input);
|
|---|
| 199 | void StressMaxPrincipalCreateInput(void);
|
|---|
| 200 | void SubglacialWaterPressure(int output_enum);
|
|---|
| 201 | IssmDouble TotalFloatingBmb(IssmDouble* mask, bool scaled);
|
|---|
| 202 | IssmDouble TotalGroundedBmb(IssmDouble* mask, bool scaled);
|
|---|
| 203 | IssmDouble TotalSmb(IssmDouble* mask, bool scaled);
|
|---|
| 204 | void TransformInvStiffnessMatrixCoord(ElementMatrix* Ke,int cs_enum);
|
|---|
| 205 | void TransformInvStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int cs_enum);
|
|---|
| 206 | void TransformInvStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 207 | void TransformLoadVectorCoord(ElementVector* pe,int cs_enum);
|
|---|
| 208 | void TransformLoadVectorCoord(ElementVector* pe,int* cs_array);
|
|---|
| 209 | void TransformLoadVectorCoord(ElementVector* pe,Node** nodes,int numnodes,int cs_enum);
|
|---|
| 210 | void TransformLoadVectorCoord(ElementVector* pe,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 211 | void TransformLoadVectorCoord(ElementVector* pe,int numnodes,int transformenum){_error_("not implemented yet");};/*Tiling only*/
|
|---|
| 212 | void TransformLoadVectorCoord(ElementVector* pe,int numnodes,int* transformenum_list){_error_("not implemented yet");};/*Tiling only*/
|
|---|
| 213 | void TransformSolutionCoord(IssmDouble* solution,int cs_enum);
|
|---|
| 214 | void TransformSolutionCoord(IssmDouble* solution,int* cs_array);
|
|---|
| 215 | void TransformSolutionCoord(IssmDouble* solution,int numnodes,int cs_enum);
|
|---|
| 216 | void TransformSolutionCoord(IssmDouble* solution,int numnodes,int* cs_array);
|
|---|
| 217 | void TransformSolutionCoord(IssmDouble* solution,Node** nodes,int numnodes,int cs_enum);
|
|---|
| 218 | void TransformSolutionCoord(IssmDouble* solution,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 219 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,int cs_enum);
|
|---|
| 220 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,int* cs_array);
|
|---|
| 221 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int cs_enum);
|
|---|
| 222 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 223 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,int numnodes,int* transformenum_list){_error_("not implemented yet");};/*Tiling only*/
|
|---|
| 224 | void ViscousHeatingCreateInput(void);
|
|---|
| 225 | void ThermalToEnthalpy(IssmDouble * penthalpy,IssmDouble temperature,IssmDouble waterfraction,IssmDouble pressure);
|
|---|
| 226 | IssmDouble TMeltingPoint(IssmDouble pressure);
|
|---|
| 227 | void EnthalpyToThermal(IssmDouble* ptemperature,IssmDouble* pwaterfraction,IssmDouble enthalpy,IssmDouble pressure);
|
|---|
| 228 | IssmDouble EnthalpyDiffusionParameter(IssmDouble enthalpy,IssmDouble pressure);
|
|---|
| 229 | IssmDouble EnthalpyDiffusionParameterVolume(int numvertices,IssmDouble* enthalpy,IssmDouble* pressure);
|
|---|
| 230 | IssmDouble PureIceEnthalpy(IssmDouble pressure);
|
|---|
| 231 |
|
|---|
| 232 | /*Virtual functions*/
|
|---|
| 233 | virtual void AddBasalInput(int input_enum, IssmDouble* values, int interpolation_enum){_error_("not implemented");};
|
|---|
| 234 | virtual void AddInput(int input_enum, IssmDouble* values, int interpolation_enum){_error_("not implemented");};
|
|---|
| 235 | virtual void AddControlInput(int input_enum,Inputs* inputs,IoModel* iomodel,IssmDouble* values,IssmDouble* values_min,IssmDouble* values_max, int interpolation_enum,int id){_error_("not supported yet");};
|
|---|
| 236 | virtual void DatasetInputCreate(IssmDouble* array,int M,int N,int* individual_enums,int num_inputs,Inputs* inputs,IoModel* iomodel,int input_enum){_error_("not supported");};
|
|---|
| 237 | virtual void AverageOntoPartition(Vector<IssmDouble>* partition_contributions,Vector<IssmDouble>* partition_areas,IssmDouble* vertex_response,IssmDouble* qmu_part)=0;
|
|---|
| 238 | virtual void BasalNodeIndices(int* pnumindices,int** pindices,int finiteelement){_error_("not implemented yet");};
|
|---|
| 239 | virtual void CalvingRateParameterization(void){_error_("not implemented yet");};
|
|---|
| 240 | virtual void CalvingRateCalvingMIP(void){_error_("not implemented yet");};
|
|---|
| 241 | virtual void CalvingRateVonmises(void){_error_("not implemented yet");};
|
|---|
| 242 | virtual void CalvingRateVonmisesAD(void){_error_("not implemented yet");};
|
|---|
| 243 | virtual void CalvingRateTest(void){_error_("not implemented yet");};
|
|---|
| 244 | virtual void CalvingCrevasseDepth(void){_error_("not implemented yet");};
|
|---|
| 245 | virtual void CalvingRateLevermann(void)=0;
|
|---|
| 246 | virtual void CalvingPollard(void){_error_("not implemented yet");};
|
|---|
| 247 | virtual void CalvingFluxLevelset(void){_error_("not implemented yet");};
|
|---|
| 248 | virtual void CalvingMeltingFluxLevelset(void){_error_("not implemented yet");};
|
|---|
| 249 | virtual IssmDouble CharacteristicLength(void)=0;
|
|---|
| 250 | virtual void ComputeBasalStress(void){_error_("not implemented yet");};
|
|---|
| 251 | virtual void ComputeDeviatoricStressTensor(void)=0;
|
|---|
| 252 | virtual void ComputeSigmaNN(void)=0;
|
|---|
| 253 | virtual void ComputeSigmaVM(void){_error_("not implemented yet");};
|
|---|
| 254 | virtual void ComputeStressTensor(void)=0;
|
|---|
| 255 | virtual void ComputeEsaStrainAndVorticity(void)=0;
|
|---|
| 256 | virtual void Configure(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters,Inputs* inputsin)=0;
|
|---|
| 257 | virtual void ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int M,int N,int interp)=0;
|
|---|
| 258 | virtual void ControlToVectors(Vector<IssmPDouble>* vector_control, Vector<IssmPDouble>* vector_gradient,int control_enum,int control_interp)=0;
|
|---|
| 259 | virtual void CreateDistanceInputFromSegmentlist(IssmDouble* distances,int distanceenum){_error_("not implemented yet");};
|
|---|
| 260 | virtual void CreateInputTimeAverage(int transientinput_enum,int averagedinput_enum,IssmDouble init_time,IssmDouble end_time,int averaging_method){_error_("not implemented yet "<<this->ObjectEnum());};
|
|---|
| 261 | virtual void ElementResponse(IssmDouble* presponse,int response_enum)=0;
|
|---|
| 262 | virtual void ElementSizes(IssmDouble* phx,IssmDouble* phy,IssmDouble* phz)=0;
|
|---|
| 263 | virtual void ElementCoordinates(Vector<IssmDouble>* vxe,Vector<IssmDouble>* vye,Vector<IssmDouble>* vze, Vector<IssmDouble>* vareae, bool spherical=false)=0;
|
|---|
| 264 | virtual void ElementCoordinates(Vector<IssmDouble>* vlonge,Vector<IssmDouble>* vlate,Vector<IssmDouble>* vareae)=0;
|
|---|
| 265 | virtual int FiniteElement(void)=0;
|
|---|
| 266 | virtual IssmDouble FloatingArea(bool scaled)=0;
|
|---|
| 267 | virtual void FSContactMigration(Vector<IssmDouble>* vertex_sigmann,Vector<IssmDouble>* vertex_waterpressure)=0;
|
|---|
| 268 | virtual Element* GetBasalElement(void)=0;
|
|---|
| 269 | virtual int GetElementType(void)=0;
|
|---|
| 270 | virtual IssmDouble GetHorizontalSurfaceArea(void){_error_("not implemented");};
|
|---|
| 271 | virtual void GetGroundedPart(int* point1,IssmDouble* fraction1,IssmDouble* fraction2, bool* mainlyfloating)=0;
|
|---|
| 272 | virtual IssmDouble GetGroundedPortion(IssmDouble* xyz_list)=0;
|
|---|
| 273 | virtual void GetFractionGeometry(IssmDouble* weights, IssmDouble* pphi, int* ppoint1,IssmDouble* pfraction1,IssmDouble* pfraction2, bool* ptrapezeisnegative, IssmDouble* gl)=0;
|
|---|
| 274 | virtual void GetNodalWeightsAndAreaAndCentroidsFromLeveset(IssmDouble* loadweights, IssmDouble* ploadarea, IssmDouble* platbar, IssmDouble* plongbar, IssmDouble late, IssmDouble longe, IssmDouble area, int levelsetenum)=0;
|
|---|
| 275 | virtual void GetNodalWeightsAndAreaAndCentroidsFromLeveset(IssmDouble* loadweights, IssmDouble* ploadarea, IssmDouble* platbar, IssmDouble* plongbar, IssmDouble late, IssmDouble longe, IssmDouble area, int levelset1enum, int levelset2enum)=0;
|
|---|
| 276 | virtual void GetBarycenterFromLevelset(IssmDouble* platbar, IssmDouble* plongbar,IssmDouble phi,IssmDouble fraction1,IssmDouble fraction2,IssmDouble late, IssmDouble longe, int point1,int istrapeze1, IssmDouble planetradius)=0;
|
|---|
| 277 |
|
|---|
| 278 | virtual IssmDouble GetIcefrontArea(){_error_("not implemented");};
|
|---|
| 279 | virtual void GetIcefrontCoordinates(IssmDouble** pxyz_front,IssmDouble* xyz_list,int levelsetenum)=0;
|
|---|
| 280 | virtual DatasetInput* GetDatasetInput(int inputenum){_error_("not implemented");};
|
|---|
| 281 | virtual Input* GetInput(int inputenum)=0;
|
|---|
| 282 | virtual Input* GetInput(int inputenum,IssmDouble time)=0;
|
|---|
| 283 | virtual Input* GetInput(int inputenum,IssmDouble start_time,IssmDouble end_time,int averaging_method)=0;
|
|---|
| 284 | virtual void GetInputValue(IssmDouble* pvalue,Vertex* vertex,int enumtype){_error_("not implemented yet");};
|
|---|
| 285 | virtual void GetInputValue(IssmDouble* pvalue,Node* node,int enumtype){_error_("not implemented yet");};
|
|---|
| 286 | virtual void GetInputListOnVertices(IssmDouble* pvalue,Input* input,IssmDouble default_value)=0;
|
|---|
| 287 | virtual void GetInputListOnNodes(IssmDouble* pvalue,Input* input,IssmDouble default_value)=0;
|
|---|
| 288 | virtual void GetLevelCoordinates(IssmDouble** pxyz_front,IssmDouble* xyz_list,int levelsetenum,IssmDouble level)=0;
|
|---|
| 289 | virtual void GetLevelsetPositivePart(int* point1,IssmDouble* fraction1,IssmDouble* fraction2, bool* mainlynegative,IssmDouble* levelsetvalues)=0;
|
|---|
| 290 | virtual int GetVertexIndex(Vertex* vertex){_error_("not implemented");};;
|
|---|
| 291 | virtual int GetNumberOfNodes(void)=0;
|
|---|
| 292 | virtual int GetNumberOfNodes(int enum_type)=0;
|
|---|
| 293 | virtual int GetNumberOfVertices(void)=0;
|
|---|
| 294 | virtual void GetVectorFromControlInputs(Vector<IssmDouble>* gradient,int control_enum,int control_index,int N,const char* data,int offset)=0;
|
|---|
| 295 | virtual void GetVerticesCoordinatesBase(IssmDouble** xyz_list)=0;
|
|---|
| 296 | virtual void GetVerticesCoordinatesTop(IssmDouble** xyz_list)=0;
|
|---|
| 297 | virtual IssmDouble GroundedArea(bool scaled)=0;
|
|---|
| 298 | virtual IssmDouble IceVolume(bool scaled)=0;
|
|---|
| 299 | virtual IssmDouble IceVolumeAboveFloatation(bool scaled)=0;
|
|---|
| 300 | virtual IssmDouble IcefrontMassFlux(bool scaled){_error_("not implemented");};
|
|---|
| 301 | virtual IssmDouble IcefrontMassFluxLevelset(bool scaled){_error_("not implemented");};
|
|---|
| 302 | virtual IssmDouble GroundinglineMassFlux(bool scaled){_error_("not implemented");};
|
|---|
| 303 | virtual void InputDepthAverageAtBase(int enum_type,int average_enum_type)=0;
|
|---|
| 304 | virtual void DatasetInputExtrude(int input_enum,int start){_error_("not implemented yet");};
|
|---|
| 305 | virtual void InputExtrude(int input_enum,int start)=0;
|
|---|
| 306 | virtual void InputUpdateFromSolutionOneDofCollapsed(IssmDouble* solution,int inputenum)=0;
|
|---|
| 307 | virtual void InputUpdateFromSolutionOneDof(IssmDouble* solution,int inputenum)=0;
|
|---|
| 308 | #ifdef _HAVE_DAKOTA_
|
|---|
| 309 | virtual void InputUpdateFromMatrixDakota(IssmDouble* matrix, int rows, int ncols, int name, int type)=0;
|
|---|
| 310 | virtual void InputUpdateFromVectorDakota(IssmDouble* vector, int name, int type)=0;
|
|---|
| 311 | virtual void InputScaleFromDakota(IssmDouble* distributed_values, IssmDouble* partition, int npart, int nt, int name)=0;
|
|---|
| 312 | #endif
|
|---|
| 313 | virtual void InputUpdateFromIoModel(int index, IoModel* iomodel)=0;
|
|---|
| 314 | virtual void InputUpdateFromVector(IssmDouble* vector, int name, int type)=0;
|
|---|
| 315 | virtual bool IsFaceOnBoundary(void)=0;
|
|---|
| 316 | virtual bool IsIcefront(void)=0;
|
|---|
| 317 | virtual bool IsNodeOnShelfFromFlags(IssmDouble* flags)=0;
|
|---|
| 318 |
|
|---|
| 319 | virtual bool IsZeroLevelset(int levelset_enum)=0;
|
|---|
| 320 | virtual void JacobianDeterminant(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 321 | virtual void JacobianDeterminantBase(IssmDouble* Jdet,IssmDouble* xyz_list_base,Gauss* gauss)=0;
|
|---|
| 322 | virtual void JacobianDeterminantLine(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 323 | virtual void JacobianDeterminantSurface(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 324 | virtual void JacobianDeterminantTop(IssmDouble* Jdet,IssmDouble* xyz_list_base,Gauss* gauss)=0;
|
|---|
| 325 | virtual void Marshall(MarshallHandle* marshallhandle)=0;
|
|---|
| 326 | virtual IssmDouble Masscon(IssmDouble* levelset)=0;
|
|---|
| 327 | virtual IssmDouble MassFlux(IssmDouble* segment)=0;
|
|---|
| 328 | virtual IssmDouble MassFlux(IssmDouble x1,IssmDouble y1, IssmDouble x2, IssmDouble y2,int segment_id)=0;
|
|---|
| 329 | virtual IssmDouble MinEdgeLength(IssmDouble* xyz_list)=0;
|
|---|
| 330 | virtual IssmDouble Misfit(int modelenum,int observationenum,int weightsenum)=0;
|
|---|
| 331 | virtual IssmDouble MisfitArea(int weightsenum)=0;
|
|---|
| 332 | virtual void MovingFrontalVelocity(void){_error_("not implemented yet");};
|
|---|
| 333 | virtual Gauss* NewGauss(void)=0;
|
|---|
| 334 | virtual Gauss* NewGauss(int order)=0;
|
|---|
| 335 | virtual Gauss* NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order)=0;
|
|---|
| 336 | virtual Gauss* NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order_horiz,int order_vert)=0;
|
|---|
| 337 | virtual Gauss* NewGauss(int point1,IssmDouble fraction1,IssmDouble fraction2,bool mainlyfloating,int order)=0;
|
|---|
| 338 | virtual Gauss* NewGauss(int point1,IssmDouble fraction1,IssmDouble fraction2,int order)=0;
|
|---|
| 339 | virtual Gauss* NewGauss(IssmDouble fraction1,IssmDouble fraction2,int order)=0;
|
|---|
| 340 | virtual Gauss* NewGaussBase(int order)=0;
|
|---|
| 341 | virtual Gauss* NewGaussLine(int vertex1,int vertex2,int order)=0;
|
|---|
| 342 | virtual Gauss* NewGaussTop(int order)=0;
|
|---|
| 343 | virtual void NodalFunctions(IssmDouble* basis,Gauss* gauss)=0;
|
|---|
| 344 | virtual void NodalFunctionsDerivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 345 | virtual void NodalFunctionsDerivativesVelocity(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 346 | virtual void NodalFunctionsMINIDerivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 347 | virtual void NodalFunctionsPressure(IssmDouble* basis, Gauss* gauss)=0;
|
|---|
| 348 | virtual void NodalFunctionsP1(IssmDouble* basis,Gauss* gauss)=0;
|
|---|
| 349 | virtual void NodalFunctionsP1Derivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 350 | virtual void NodalFunctionsP2(IssmDouble* basis,Gauss* gauss)=0;
|
|---|
| 351 | virtual void NodalFunctionsVelocity(IssmDouble* basis, Gauss* gauss)=0;
|
|---|
| 352 | virtual void NodalFunctionsTensor(IssmDouble* basis, Gauss* gauss)=0;
|
|---|
| 353 | virtual int NodalValue(IssmDouble* pvalue, int index, int natureofdataenum)=0;
|
|---|
| 354 | virtual void NormalBase(IssmDouble* normal,IssmDouble* xyz_list)=0;
|
|---|
| 355 | virtual void NormalSection(IssmDouble* normal,IssmDouble* xyz_list)=0;
|
|---|
| 356 | virtual void NormalTop(IssmDouble* normal,IssmDouble* xyz_list)=0;
|
|---|
| 357 | virtual int NumberofNodesPressure(void)=0;
|
|---|
| 358 | virtual int NumberofNodesVelocity(void)=0;
|
|---|
| 359 | virtual void PotentialUngrounding(Vector<IssmDouble>* potential_sheet_ungrounding)=0;
|
|---|
| 360 | virtual int PressureInterpolation()=0;
|
|---|
| 361 | virtual void Recover3DMOLHOInput(int targetVel_enum, int numnodes, IssmDouble* vb, IssmDouble* vsh, IssmDouble* n, IssmDouble* H, IssmDouble* s){_error_("not implemented yet");};
|
|---|
| 362 | virtual void ReduceMatrices(ElementMatrix* Ke,ElementVector* pe)=0;
|
|---|
| 363 | virtual void ResetFSBasalBoundaryCondition()=0;
|
|---|
| 364 | virtual void ResetHooks()=0;
|
|---|
| 365 | virtual void SetElementInput(int enum_in,IssmDouble value){_error_("not implemented yet");};
|
|---|
| 366 | virtual void SetElementInput(int enum_in,IssmDouble value, int type)=0;
|
|---|
| 367 | virtual void SetElementInput(Inputs* inputs,int enum_in,IssmDouble values){_error_("not implemented yet");};
|
|---|
| 368 | virtual void SetElementInput(Inputs* inputs,int numindices,int* indices,IssmDouble* values,int enum_in){_error_("not implemented yet");};
|
|---|
| 369 | virtual void SetControlInputsFromVector(IssmDouble* vector,int control_enum,int control_index,int offset,int M,int N)=0;
|
|---|
| 370 | virtual void SetCurrentConfiguration(Elements* elements,Loads* loads,Nodes* nodes,Materials* materials,Parameters* parameters)=0;
|
|---|
| 371 | virtual void SetTemporaryElementType(int element_type_in)=0;
|
|---|
| 372 | virtual Element* SpawnBasalElement(bool depthaverage_materials=false)=0;
|
|---|
| 373 | virtual Element* SpawnTopElement(void)=0;
|
|---|
| 374 | virtual bool IsSpawnedElement(void)=0;
|
|---|
| 375 | virtual IssmDouble StabilizationParameter(IssmDouble u, IssmDouble v, IssmDouble w, IssmDouble diameter, IssmDouble kappa)=0;
|
|---|
| 376 | virtual void StabilizationParameterAnisotropic(IssmDouble* tau_parameter_anisotropic, IssmDouble u, IssmDouble v, IssmDouble w, IssmDouble hx, IssmDouble hy, IssmDouble hz, IssmDouble kappa)=0;
|
|---|
| 377 | virtual void StrainRateparallel(void)=0;
|
|---|
| 378 | virtual void StrainRateperpendicular(void)=0;
|
|---|
| 379 | virtual void StressIntensityFactor(void)=0;
|
|---|
| 380 | virtual IssmDouble SurfaceArea(void)=0;
|
|---|
| 381 | virtual void TangentBase(IssmDouble* bed_tangent,IssmDouble* bed_normal){_error_("not implemented yet");};
|
|---|
| 382 | virtual int TensorInterpolation()=0;
|
|---|
| 383 | virtual IssmDouble TimeAdapt()=0;
|
|---|
| 384 | virtual IssmDouble TotalCalvingFluxLevelset(bool scaled){_error_("not implemented");};
|
|---|
| 385 | virtual IssmDouble TotalCalvingMeltingFluxLevelset(bool scaled){_error_("not implemented");};
|
|---|
| 386 | virtual IssmDouble TotalFloatingBmb(bool scaled)=0;
|
|---|
| 387 | virtual IssmDouble TotalGroundedBmb(bool scaled)=0;
|
|---|
| 388 | virtual IssmDouble TotalSmb(bool scaled)=0;
|
|---|
| 389 | virtual void Update(Inputs* inputs,int index, IoModel* iomodel,int analysis_counter,int analysis_type,int finite_element)=0;
|
|---|
| 390 | virtual void UpdateConstraintsExtrudeFromBase(void)=0;
|
|---|
| 391 | virtual void UpdateConstraintsExtrudeFromTop(void)=0;
|
|---|
| 392 | virtual int UpdatePotentialUngrounding(IssmDouble* potential_sheet_ungrounding,Vector<IssmDouble>* vec_nodes_on_iceshelf,IssmDouble* nodes_on_iceshelf)=0;
|
|---|
| 393 | virtual void ValueP1DerivativesOnGauss(IssmDouble* dvalue,IssmDouble* values,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 394 | virtual void ValueP1OnGauss(IssmDouble* pvalue,IssmDouble* values,Gauss* gauss)=0;
|
|---|
| 395 | virtual int VelocityInterpolation()=0;
|
|---|
| 396 | virtual int VertexConnectivity(int vertexindex)=0;
|
|---|
| 397 | virtual void VerticalSegmentIndices(int** pindices,int* pnumseg)=0;
|
|---|
| 398 | virtual void VerticalSegmentIndicesBase(int** pindices,int* pnumseg)=0;
|
|---|
| 399 | virtual void ViscousHeating(IssmDouble* pphi,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input){_error_("not implemented yet");};
|
|---|
| 400 | virtual void WriteFieldIsovalueSegment(DataSet* segments,int fieldenum,IssmDouble fieldvalue){_error_("not implemented yet");};
|
|---|
| 401 |
|
|---|
| 402 | #ifdef _HAVE_ESA_
|
|---|
| 403 | virtual void EsaGeodetic2D(Vector<IssmDouble>* pUp,Vector<IssmDouble>* pNorth,Vector<IssmDouble>* pEast, Vector<IssmDouble>* pX, Vector<IssmDouble>* pY,IssmDouble* xx,IssmDouble* yy)=0;
|
|---|
| 404 | virtual void EsaGeodetic3D(Vector<IssmDouble>* pUp,Vector<IssmDouble>* pNorth,Vector<IssmDouble>* pEast,IssmDouble* latitude,IssmDouble* longitude,IssmDouble* radius,IssmDouble* xx,IssmDouble* yy,IssmDouble* zz)=0;
|
|---|
| 405 | #endif
|
|---|
| 406 | #ifdef _HAVE_SEALEVELCHANGE_
|
|---|
| 407 | virtual IssmDouble GetArea3D(void)=0;
|
|---|
| 408 | virtual IssmDouble GetAreaSpherical(void)=0;
|
|---|
| 409 | virtual IssmDouble GetTriangleAreaSpherical(IssmDouble xyz_list[3][3])=0;
|
|---|
| 410 | virtual void GiaDeflection(Vector<IssmDouble>* wg,Vector<IssmDouble>* dwgdt,Matlitho* litho, IssmDouble* x,IssmDouble* y)=0;
|
|---|
| 411 |
|
|---|
| 412 | virtual void SealevelchangeGeometrySubElementKernel(SealevelGeometry* slgeom)=0;
|
|---|
| 413 | virtual void SealevelchangeShift(GrdLoads* loads, IssmDouble offset, SealevelGeometry* slgeom)=0;
|
|---|
| 414 | virtual void SealevelchangeGeometryInitial(IssmDouble* xxe, IssmDouble* yye, IssmDouble* zze, IssmDouble* areae, int* lids, int* vcount)=0;
|
|---|
| 415 | virtual void SealevelchangeGeometryCentroidLoads(SealevelGeometry* slgeom, IssmDouble* xxe, IssmDouble* yye, IssmDouble* zze, IssmDouble* areae)=0;
|
|---|
| 416 | virtual void SealevelchangeGeometrySubElementLoads(SealevelGeometry* slgeom, IssmDouble* areae)=0;
|
|---|
| 417 | virtual void SealevelchangeBarystaticLoads(GrdLoads* loads, BarystaticContributions* barycontrib, SealevelGeometry* slgeom)=0;
|
|---|
| 418 | virtual void SealevelchangeConvolution(IssmDouble* sealevelpercpu, GrdLoads* loads, IssmDouble* rotationvector,SealevelGeometry* slgeom)=0;
|
|---|
| 419 | virtual void SealevelchangeOceanAverage(GrdLoads* loads, Vector<IssmDouble>* oceanareas, Vector<IssmDouble>* subelementoceanareas, IssmDouble* sealevelpercpu, SealevelGeometry* slgeom)=0;
|
|---|
| 420 | virtual void SealevelchangeDeformationConvolution(IssmDouble* sealevelpercpu, GrdLoads* loads, IssmDouble* rotationvector,SealevelGeometry* slgeom)=0;
|
|---|
| 421 | virtual void SealevelchangeUpdateViscousFields(IssmDouble lincoeff, int newindex, int offset)=0;
|
|---|
| 422 | #endif
|
|---|
| 423 |
|
|---|
| 424 | };
|
|---|
| 425 | #endif
|
|---|
![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
