| 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 | #include "../Update.h"
|
|---|
| 15 | class DataSet;
|
|---|
| 16 | class Parameters;
|
|---|
| 17 | class Parameter;
|
|---|
| 18 | class Elements;
|
|---|
| 19 | class Loads;
|
|---|
| 20 | class Nodes;
|
|---|
| 21 | class Node;
|
|---|
| 22 | class Vertices;
|
|---|
| 23 | class Vertex;
|
|---|
| 24 | class Materials;
|
|---|
| 25 | class Material;
|
|---|
| 26 | class Matpar;
|
|---|
| 27 | class Inputs;
|
|---|
| 28 | class Input;
|
|---|
| 29 | class Gauss;
|
|---|
| 30 | class ElementVector;
|
|---|
| 31 | template <class doublematrix> class Matrix;
|
|---|
| 32 | template <class doubletype> class Vector;
|
|---|
| 33 | class ElementMatrix;
|
|---|
| 34 | class ElementVector;
|
|---|
| 35 | /*}}}*/
|
|---|
| 36 |
|
|---|
| 37 | class Element: public Object,public Update{
|
|---|
| 38 |
|
|---|
| 39 | public:
|
|---|
| 40 | int id;
|
|---|
| 41 | int sid;
|
|---|
| 42 | Inputs *inputs;
|
|---|
| 43 | Node **nodes;
|
|---|
| 44 | Vertex **vertices;
|
|---|
| 45 | Material *material;
|
|---|
| 46 | Matpar *matpar;
|
|---|
| 47 | Parameters *parameters;
|
|---|
| 48 |
|
|---|
| 49 | int* element_type_list;
|
|---|
| 50 | int element_type;
|
|---|
| 51 |
|
|---|
| 52 | public:
|
|---|
| 53 | /*Constructors/Destructores*/
|
|---|
| 54 | Element();
|
|---|
| 55 | ~Element();
|
|---|
| 56 |
|
|---|
| 57 | /*Functions*/
|
|---|
| 58 | void AddInput(Input* input_in);
|
|---|
| 59 | /*bool AllActive(void);*/
|
|---|
| 60 | /*bool AnyActive(void);*/
|
|---|
| 61 | void ComputeLambdaS(void);
|
|---|
| 62 | void ComputeNewDamage();
|
|---|
| 63 | void ComputeStrainRate();
|
|---|
| 64 | void CoordinateSystemTransform(IssmDouble** ptransform,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 65 | void DeepEcho();
|
|---|
| 66 | void DeleteInput(int input_enum);
|
|---|
| 67 | void DeleteMaterials(void);
|
|---|
| 68 | void Delta18oParameterization(void);
|
|---|
| 69 | void Delta18opdParameterization(void);
|
|---|
| 70 | IssmDouble Divergence(void);
|
|---|
| 71 | void dViscositydBFS(IssmDouble* pdmudB,int dim,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input);
|
|---|
| 72 | void dViscositydBHO(IssmDouble* pdmudB,int dim,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 73 | void dViscositydBSSA(IssmDouble* pdmudB,int dim,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 74 | void dViscositydDSSA(IssmDouble* pdmudB,int dim,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 75 | void Echo();
|
|---|
| 76 | IssmDouble EnthalpyDiffusionParameter(IssmDouble enthalpy,IssmDouble pressure);
|
|---|
| 77 | IssmDouble EnthalpyDiffusionParameterVolume(int numvertices,IssmDouble* enthalpy,IssmDouble* pressure);
|
|---|
| 78 | void EnthalpyToThermal(IssmDouble* ptemperature,IssmDouble* pwaterfraction,IssmDouble enthalpy,IssmDouble pressure);
|
|---|
| 79 | void FindParam(bool* pvalue,int paramenum);
|
|---|
| 80 | void FindParam(int* pvalue,int paramenum);
|
|---|
| 81 | void FindParam(IssmDouble* pvalue,int paramenum);
|
|---|
| 82 | void FindParam(int** pvalues,int* psize,int paramenum);
|
|---|
| 83 | IssmDouble FloatingArea(IssmDouble* mask);
|
|---|
| 84 | void GetDofList(int** pdoflist,int approximation_enum,int setenum);
|
|---|
| 85 | void GetDofListPressure(int** pdoflist,int setenum);
|
|---|
| 86 | void GetDofListVelocity(int** pdoflist,int setenum);
|
|---|
| 87 | Input* GetInput(int inputenum);
|
|---|
| 88 | void GetInputListOnNodes(IssmDouble* pvalue,int enumtype);
|
|---|
| 89 | void GetInputListOnNodes(IssmDouble* pvalue,int enumtype,IssmDouble defaultvalue);
|
|---|
| 90 | void GetInputListOnNodesVelocity(IssmDouble* pvalue,int enumtype);
|
|---|
| 91 | void GetInputListOnVertices(IssmDouble* pvalue,int enumtype);
|
|---|
| 92 | void GetInputListOnVertices(IssmDouble* pvalue,int enumtype,IssmDouble defaultvalue);
|
|---|
| 93 | void GetInputLocalMinMaxOnNodes(IssmDouble* min,IssmDouble* max,IssmDouble* ug);
|
|---|
| 94 | void GetInputValue(bool* pvalue,int enum_type);
|
|---|
| 95 | void GetInputValue(int* pvalue,int enum_type);
|
|---|
| 96 | void GetInputValue(IssmDouble* pvalue,int enum_type);
|
|---|
| 97 | void GetInputValue(IssmDouble* pvalue,Gauss* gauss,int enum_type);
|
|---|
| 98 | void GetInputsInterpolations(Vector<IssmDouble>* interps);
|
|---|
| 99 | IssmDouble GetMaterialParameter(int enum_in);
|
|---|
| 100 | void GetNodesLidList(int* lidlist);
|
|---|
| 101 | void GetNodesSidList(int* sidlist);
|
|---|
| 102 | void GetPhi(IssmDouble* phi, IssmDouble* epsilon, IssmDouble viscosity);
|
|---|
| 103 | void GetVectorFromInputs(Vector<IssmDouble>* vector, int name_enum, int type);
|
|---|
| 104 | void GetVertexPidList(int* pidlist);
|
|---|
| 105 | void GetVerticesConnectivityList(int* connectivitylist);
|
|---|
| 106 | void GetVerticesCoordinates(IssmDouble** xyz_list);
|
|---|
| 107 | void GetVerticesSidList(int* sidlist);
|
|---|
| 108 | IssmDouble GetXcoord(IssmDouble* xyz_list,Gauss* gauss);
|
|---|
| 109 | IssmDouble GetYcoord(IssmDouble* xyz_list,Gauss* gauss);
|
|---|
| 110 | IssmDouble GetZcoord(IssmDouble* xyz_list,Gauss* gauss);
|
|---|
| 111 | void GradientIndexing(int* indexing,int control_index,bool onsid=false);
|
|---|
| 112 | IssmDouble GroundedArea(IssmDouble* mask);
|
|---|
| 113 | bool HasNodeOnBase();
|
|---|
| 114 | bool HasNodeOnSurface();
|
|---|
| 115 | IssmDouble IceMass();
|
|---|
| 116 | IssmDouble IceMass(IssmDouble* mask);
|
|---|
| 117 | IssmDouble IceVolume(IssmDouble* mask);
|
|---|
| 118 | IssmDouble IceVolumeAboveFloatation(IssmDouble* mask);
|
|---|
| 119 | int Id();
|
|---|
| 120 | void InputChangeName(int enum_type,int enum_type_old);
|
|---|
| 121 | void InputCreate(IssmDouble* vector,IoModel* iomodel,int M,int N,int vector_type,int vector_enum,int code);
|
|---|
| 122 | void InputDuplicate(int original_enum,int new_enum);
|
|---|
| 123 | void InputUpdateFromConstant(IssmDouble constant, int name);
|
|---|
| 124 | void InputUpdateFromConstant(int constant, int name);
|
|---|
| 125 | void InputUpdateFromConstant(bool constant, int name);
|
|---|
| 126 | bool IsFloating();
|
|---|
| 127 | bool IsIceInElement();
|
|---|
| 128 | bool IsInput(int name);
|
|---|
| 129 | bool IsLandInElement();
|
|---|
| 130 | bool IsWaterInElement();
|
|---|
| 131 | void LinearFloatingiceMeltingRate();
|
|---|
| 132 | void MantlePlumeGeothermalFlux();
|
|---|
| 133 | void MarshallElement(char** pmarshalled_data,int* pmarshalled_data_size, int marshall_direction,int numanalyses);
|
|---|
| 134 | void MigrateGroundingLine(IssmDouble* sheet_ungrounding);
|
|---|
| 135 | void MismipFloatingiceMeltingRate();
|
|---|
| 136 | void MungsmtpParameterization(void);
|
|---|
| 137 | ElementMatrix* NewElementMatrix(int approximation_enum=NoneApproximationEnum);
|
|---|
| 138 | ElementMatrix* NewElementMatrixCoupling(int number_nodes,int approximation_enum=NoneApproximationEnum);
|
|---|
| 139 | ElementVector* NewElementVector(int approximation_enum=NoneApproximationEnum);
|
|---|
| 140 | void PositiveDegreeDay(IssmDouble* pdds,IssmDouble* pds,IssmDouble signorm,bool ismungsm);
|
|---|
| 141 | IssmDouble PureIceEnthalpy(IssmDouble pressure);
|
|---|
| 142 | void ResultInterpolation(int* pinterpolation,int*nodesperelement,int* parray_size, int output_enum);
|
|---|
| 143 | void ResultToPatch(IssmDouble* values,int nodesperelement,int output_enum);
|
|---|
| 144 | void ResultToMatrix(IssmDouble* values,int ncols,int output_enum);
|
|---|
| 145 | void ResultToVector(Vector<IssmDouble>* vector,int output_enum);
|
|---|
| 146 | void SetwiseNodeConnectivity(int* d_nz,int* o_nz,Node* node,bool* flags,int* flagsindices,int set1_enum,int set2_enum);
|
|---|
| 147 | int Sid();
|
|---|
| 148 | void SmbGemb();
|
|---|
| 149 | void StrainRateESA(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 150 | void StrainRateFS(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input);
|
|---|
| 151 | void StrainRateHO(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 152 | void StrainRateHO2dvertical(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 153 | void StrainRateSSA(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 154 | void StrainRateSSA1d(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input);
|
|---|
| 155 | void StressMaxPrincipalCreateInput(void);
|
|---|
| 156 | void ThermalToEnthalpy(IssmDouble* penthalpy,IssmDouble temperature,IssmDouble waterfraction,IssmDouble pressure);
|
|---|
| 157 | IssmDouble TMeltingPoint(IssmDouble pressure);
|
|---|
| 158 | IssmDouble TotalFloatingBmb(IssmDouble* mask);
|
|---|
| 159 | IssmDouble TotalGroundedBmb(IssmDouble* mask);
|
|---|
| 160 | IssmDouble TotalSmb(IssmDouble* mask);
|
|---|
| 161 | void TransformInvStiffnessMatrixCoord(ElementMatrix* Ke,int cs_enum);
|
|---|
| 162 | void TransformInvStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int cs_enum);
|
|---|
| 163 | void TransformInvStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 164 | void TransformLoadVectorCoord(ElementVector* pe,int cs_enum);
|
|---|
| 165 | void TransformLoadVectorCoord(ElementVector* pe,int* cs_array);
|
|---|
| 166 | void TransformLoadVectorCoord(ElementVector* pe,Node** nodes,int numnodes,int cs_enum);
|
|---|
| 167 | void TransformLoadVectorCoord(ElementVector* pe,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 168 | void TransformLoadVectorCoord(ElementVector* pe,int numnodes,int transformenum){_error_("not implemented yet");};/*Tiling only*/
|
|---|
| 169 | void TransformLoadVectorCoord(ElementVector* pe,int numnodes,int* transformenum_list){_error_("not implemented yet");};/*Tiling only*/
|
|---|
| 170 | void TransformSolutionCoord(IssmDouble* solution,int cs_enum);
|
|---|
| 171 | void TransformSolutionCoord(IssmDouble* solution,int* cs_array);
|
|---|
| 172 | void TransformSolutionCoord(IssmDouble* solution,int numnodes,int cs_enum);
|
|---|
| 173 | void TransformSolutionCoord(IssmDouble* solution,int numnodes,int* cs_array);
|
|---|
| 174 | void TransformSolutionCoord(IssmDouble* solution,Node** nodes,int numnodes,int cs_enum);
|
|---|
| 175 | void TransformSolutionCoord(IssmDouble* solution,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 176 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,int cs_enum);
|
|---|
| 177 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,int* cs_array);
|
|---|
| 178 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int cs_enum);
|
|---|
| 179 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 180 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,int numnodes,int* transformenum_list){_error_("not implemented yet");};/*Tiling only*/
|
|---|
| 181 | void ViscousHeatingCreateInput(void);
|
|---|
| 182 |
|
|---|
| 183 | /*Virtual functions*/
|
|---|
| 184 | virtual void AddBasalInput(int input_enum, IssmDouble* values, int interpolation_enum)=0;
|
|---|
| 185 | virtual void AddInput(int input_enum, IssmDouble* values, int interpolation_enum)=0;
|
|---|
| 186 | virtual void AverageOntoPartition(Vector<IssmDouble>* partition_contributions,Vector<IssmDouble>* partition_areas,IssmDouble* vertex_response,IssmDouble* qmu_part)=0;
|
|---|
| 187 | virtual void BasalNodeIndices(int* pnumindices,int** pindices,int finiteelement){_error_("not implemented yet");};
|
|---|
| 188 | virtual void CalvingRateDev(void){_error_("not implemented yet");};
|
|---|
| 189 | virtual void CalvingRateLevermann(void)=0;
|
|---|
| 190 | virtual IssmDouble CharacteristicLength(void)=0;
|
|---|
| 191 | virtual void ComputeBasalStress(Vector<IssmDouble>* sigma_b)=0;
|
|---|
| 192 | virtual void ComputeDeviatoricStressTensor(void)=0;
|
|---|
| 193 | virtual void ComputeSigmaNN(void)=0;
|
|---|
| 194 | virtual void ComputeStressTensor(void)=0;
|
|---|
| 195 | virtual void ComputeEsaStrainAndVorticity(void)=0;
|
|---|
| 196 | virtual void Configure(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters)=0;
|
|---|
| 197 | virtual void ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index)=0;
|
|---|
| 198 | virtual void ControlToVectors(Vector<IssmPDouble>* vector_control, Vector<IssmPDouble>* vector_gradient,int control_enum)=0;
|
|---|
| 199 | virtual void ElementResponse(IssmDouble* presponse,int response_enum)=0;
|
|---|
| 200 | virtual void ElementSizes(IssmDouble* phx,IssmDouble* phy,IssmDouble* phz)=0;
|
|---|
| 201 | virtual int FiniteElement(void)=0;
|
|---|
| 202 | virtual IssmDouble FloatingArea(void)=0;
|
|---|
| 203 | virtual void FSContactMigration(Vector<IssmDouble>* vertexgrounded,Vector<IssmDouble>* vertexfloating)=0;
|
|---|
| 204 | virtual Element* GetBasalElement(void)=0;
|
|---|
| 205 | virtual int GetElementType(void)=0;
|
|---|
| 206 | virtual void GetGroundedPart(int* point1,IssmDouble* fraction1,IssmDouble* fraction2, bool* mainlyfloating)=0;
|
|---|
| 207 | virtual IssmDouble GetGroundedPortion(IssmDouble* xyz_list)=0;
|
|---|
| 208 | virtual void GetIcefrontCoordinates(IssmDouble** pxyz_front,IssmDouble* xyz_list,int levelsetenum)=0;
|
|---|
| 209 | virtual void GetInputValue(IssmDouble* pvalue,Node* node,int enumtype)=0;
|
|---|
| 210 | virtual void GetLevelCoordinates(IssmDouble** pxyz_front,IssmDouble* xyz_list,int levelsetenum,IssmDouble level)=0;
|
|---|
| 211 | virtual void GetLevelsetPositivePart(int* point1,IssmDouble* fraction1,IssmDouble* fraction2, bool* mainlynegative,IssmDouble* levelsetvalues)=0;
|
|---|
| 212 | virtual Node* GetNode(int node_number)=0;
|
|---|
| 213 | virtual int GetNodeIndex(Node* node)=0;
|
|---|
| 214 | virtual int GetNumberOfNodes(void)=0;
|
|---|
| 215 | virtual int GetNumberOfNodes(int enum_type)=0;
|
|---|
| 216 | virtual int GetNumberOfVertices(void)=0;
|
|---|
| 217 | virtual void GetSolutionFromInputsOneDof(Vector<IssmDouble>* solution,int solutionenum)=0;
|
|---|
| 218 | virtual void GetVectorFromControlInputs(Vector<IssmDouble>* gradient,int control_enum,int control_index,const char* data,bool onsid)=0;
|
|---|
| 219 | virtual void GetVerticesCoordinatesBase(IssmDouble** xyz_list)=0;
|
|---|
| 220 | virtual void GetVerticesCoordinatesTop(IssmDouble** xyz_list)=0;
|
|---|
| 221 | virtual IssmDouble GroundedArea(void)=0;
|
|---|
| 222 | virtual IssmDouble IceVolume(void)=0;
|
|---|
| 223 | virtual IssmDouble IceVolumeAboveFloatation(void)=0;
|
|---|
| 224 | virtual void InputControlUpdate(IssmDouble scalar,bool save_parameter)=0;
|
|---|
| 225 | virtual void InputDepthAverageAtBase(int enum_type,int average_enum_type)=0;
|
|---|
| 226 | virtual void InputExtrude(int input_enum,int start)=0;
|
|---|
| 227 | virtual void InputScale(int enum_type,IssmDouble scale_factor)=0;
|
|---|
| 228 | virtual void InputUpdateFromSolutionOneDofCollapsed(IssmDouble* solution,int inputenum)=0;
|
|---|
| 229 | virtual void InputUpdateFromSolutionOneDof(IssmDouble* solution,int inputenum)=0;
|
|---|
| 230 | virtual bool IsFaceOnBoundary(void)=0;
|
|---|
| 231 | virtual bool IsIcefront(void)=0;
|
|---|
| 232 | virtual bool IsNodeOnShelfFromFlags(IssmDouble* flags)=0;
|
|---|
| 233 | virtual bool IsOnBase()=0;
|
|---|
| 234 | virtual bool IsOnSurface()=0;
|
|---|
| 235 | virtual bool IsZeroLevelset(int levelset_enum)=0;
|
|---|
| 236 | virtual void JacobianDeterminant(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 237 | virtual void JacobianDeterminantBase(IssmDouble* Jdet,IssmDouble* xyz_list_base,Gauss* gauss)=0;
|
|---|
| 238 | virtual void JacobianDeterminantLine(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 239 | virtual void JacobianDeterminantSurface(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 240 | virtual void JacobianDeterminantTop(IssmDouble* Jdet,IssmDouble* xyz_list_base,Gauss* gauss)=0;
|
|---|
| 241 | virtual void Marshall(char** pmarshalled_data,int* pmarshalled_data_size, int marshall_direction)=0;
|
|---|
| 242 | virtual IssmDouble Masscon(IssmDouble* levelset)=0;
|
|---|
| 243 | virtual IssmDouble MassFlux(IssmDouble* segment)=0;
|
|---|
| 244 | virtual IssmDouble MassFlux(IssmDouble x1,IssmDouble y1, IssmDouble x2, IssmDouble y2,int segment_id)=0;
|
|---|
| 245 | virtual IssmDouble MinEdgeLength(IssmDouble* xyz_list)=0;
|
|---|
| 246 | virtual IssmDouble Misfit(int modelenum,int observationenum,int weightsenum)=0;
|
|---|
| 247 | virtual IssmDouble MisfitArea(int weightsenum)=0;
|
|---|
| 248 | virtual Gauss* NewGauss(void)=0;
|
|---|
| 249 | virtual Gauss* NewGauss(int order)=0;
|
|---|
| 250 | virtual Gauss* NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order)=0;
|
|---|
| 251 | virtual Gauss* NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order_horiz,int order_vert)=0;
|
|---|
| 252 | virtual Gauss* NewGauss(int point1,IssmDouble fraction1,IssmDouble fraction2,bool mainlyfloating,int order)=0;
|
|---|
| 253 | virtual Gauss* NewGaussBase(int order)=0;
|
|---|
| 254 | virtual Gauss* NewGaussLine(int vertex1,int vertex2,int order)=0;
|
|---|
| 255 | virtual Gauss* NewGaussTop(int order)=0;
|
|---|
| 256 | virtual void NodalFunctions(IssmDouble* basis,Gauss* gauss)=0;
|
|---|
| 257 | virtual void NodalFunctionsDerivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 258 | virtual void NodalFunctionsDerivativesVelocity(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 259 | virtual void NodalFunctionsMINIDerivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 260 | virtual void NodalFunctionsPressure(IssmDouble* basis, Gauss* gauss)=0;
|
|---|
| 261 | virtual void NodalFunctionsP1(IssmDouble* basis,Gauss* gauss)=0;
|
|---|
| 262 | virtual void NodalFunctionsP1Derivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 263 | virtual void NodalFunctionsP2(IssmDouble* basis,Gauss* gauss)=0;
|
|---|
| 264 | virtual void NodalFunctionsVelocity(IssmDouble* basis, Gauss* gauss)=0;
|
|---|
| 265 | virtual void NodalFunctionsTensor(IssmDouble* basis, Gauss* gauss)=0;
|
|---|
| 266 | virtual int NodalValue(IssmDouble* pvalue, int index, int natureofdataenum)=0;
|
|---|
| 267 | virtual void NormalBase(IssmDouble* normal,IssmDouble* xyz_list)=0;
|
|---|
| 268 | virtual void NormalSection(IssmDouble* normal,IssmDouble* xyz_list)=0;
|
|---|
| 269 | virtual void NormalTop(IssmDouble* normal,IssmDouble* xyz_list)=0;
|
|---|
| 270 | virtual int NumberofNodesPressure(void)=0;
|
|---|
| 271 | virtual int NumberofNodesVelocity(void)=0;
|
|---|
| 272 | virtual void PotentialUngrounding(Vector<IssmDouble>* potential_sheet_ungrounding)=0;
|
|---|
| 273 | virtual int PressureInterpolation()=0;
|
|---|
| 274 | virtual void ReduceMatrices(ElementMatrix* Ke,ElementVector* pe)=0;
|
|---|
| 275 | virtual void ResetFSBasalBoundaryCondition()=0;
|
|---|
| 276 | virtual void ResetHooks()=0;
|
|---|
| 277 | virtual void ResetLevelsetFromSegmentlist(IssmDouble* segments,int numsegments){_error_("not implemented yet");};
|
|---|
| 278 | virtual void SetControlInputsFromVector(IssmDouble* vector,int control_enum,int control_index)=0;
|
|---|
| 279 | virtual void SetCurrentConfiguration(Elements* elements,Loads* loads,Nodes* nodes,Materials* materials,Parameters* parameters)=0;
|
|---|
| 280 | virtual void SetTemporaryElementType(int element_type_in)=0;
|
|---|
| 281 | virtual Element* SpawnBasalElement(void)=0;
|
|---|
| 282 | virtual Element* SpawnTopElement(void)=0;
|
|---|
| 283 | virtual IssmDouble StabilizationParameter(IssmDouble u, IssmDouble v, IssmDouble w, IssmDouble diameter, IssmDouble kappa)=0;
|
|---|
| 284 | virtual void StrainRateparallel(void)=0;
|
|---|
| 285 | virtual void StrainRateperpendicular(void)=0;
|
|---|
| 286 | virtual void StressIntensityFactor(void)=0;
|
|---|
| 287 | virtual IssmDouble SurfaceArea(void)=0;
|
|---|
| 288 | virtual int TensorInterpolation()=0;
|
|---|
| 289 | virtual IssmDouble TimeAdapt()=0;
|
|---|
| 290 | virtual IssmDouble TotalFloatingBmb(void)=0;
|
|---|
| 291 | virtual IssmDouble TotalGroundedBmb(void)=0;
|
|---|
| 292 | virtual IssmDouble TotalSmb(void)=0;
|
|---|
| 293 | virtual void Update(int index, IoModel* iomodel,int analysis_counter,int analysis_type,int finite_element)=0;
|
|---|
| 294 | virtual void UpdateConstraintsExtrudeFromBase(void)=0;
|
|---|
| 295 | virtual void UpdateConstraintsExtrudeFromTop(void)=0;
|
|---|
| 296 | virtual int UpdatePotentialUngrounding(IssmDouble* potential_sheet_ungrounding,Vector<IssmDouble>* vec_nodes_on_iceshelf,IssmDouble* nodes_on_iceshelf)=0;
|
|---|
| 297 | virtual void ValueP1DerivativesOnGauss(IssmDouble* dvalue,IssmDouble* values,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 298 | virtual void ValueP1OnGauss(IssmDouble* pvalue,IssmDouble* values,Gauss* gauss)=0;
|
|---|
| 299 | virtual int VelocityInterpolation()=0;
|
|---|
| 300 | virtual int VertexConnectivity(int vertexindex)=0;
|
|---|
| 301 | virtual void VerticalSegmentIndices(int** pindices,int* pnumseg)=0;
|
|---|
| 302 | virtual void VerticalSegmentIndicesBase(int** pindices,int* pnumseg)=0;
|
|---|
| 303 | virtual void ViscousHeating(IssmDouble* pphi,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input)=0;
|
|---|
| 304 | virtual void WriteLevelsetSegment(DataSet* segments){_error_("not implemented yet");};
|
|---|
| 305 | virtual void ZeroLevelsetCoordinates(IssmDouble** pxyz_zero,IssmDouble* xyz_list,int levelsetenum)=0;
|
|---|
| 306 |
|
|---|
| 307 | #ifdef _HAVE_GIAIVINS_
|
|---|
| 308 | virtual void GiaDeflection(Vector<IssmDouble>* wg,Vector<IssmDouble>* dwgdt,IssmDouble* x,IssmDouble* y)=0;
|
|---|
| 309 | #endif
|
|---|
| 310 | #ifdef _HAVE_ESA_
|
|---|
| 311 | virtual void EsaGeodetic2D(Vector<IssmDouble>* pUp,Vector<IssmDouble>* pNorth,Vector<IssmDouble>* pEast,IssmDouble* xx,IssmDouble* yy)=0;
|
|---|
| 312 | virtual void EsaGeodetic3D(Vector<IssmDouble>* pUp,Vector<IssmDouble>* pNorth,Vector<IssmDouble>* pEast,IssmDouble* latitude,IssmDouble* longitude,IssmDouble* radius,IssmDouble* xx,IssmDouble* yy,IssmDouble* zz,IssmDouble eartharea)=0;
|
|---|
| 313 | #endif
|
|---|
| 314 | #ifdef _HAVE_SEALEVELRISE_
|
|---|
| 315 | virtual IssmDouble GetArea3D(void)=0;
|
|---|
| 316 | virtual IssmDouble GetAreaSpherical(void)=0;
|
|---|
| 317 | virtual IssmDouble OceanAverage(IssmDouble* Sg)=0;
|
|---|
| 318 | virtual IssmDouble OceanArea(void)=0;
|
|---|
| 319 | virtual void SealevelriseMomentOfInertia(IssmDouble* dI_list,IssmDouble* Sg_old,IssmDouble eartharea)=0;
|
|---|
| 320 | virtual void SealevelriseEustatic(Vector<IssmDouble>* pSgi,IssmDouble* peustatic,IssmDouble* latitude,IssmDouble* longitude,IssmDouble* radius,IssmDouble oceanarea,IssmDouble eartharea)=0;
|
|---|
| 321 | virtual void SealevelriseNonEustatic(Vector<IssmDouble>* pSgo,IssmDouble* Sg_old,IssmDouble* latitude,IssmDouble* longitude,IssmDouble* radius,IssmDouble eartharea)=0;
|
|---|
| 322 | virtual void SealevelriseGeodetic(Vector<IssmDouble>* pUp,Vector<IssmDouble>* pNorth,Vector<IssmDouble>* pEast,IssmDouble* Sg,IssmDouble* latitude,IssmDouble* longitude,IssmDouble* radius,IssmDouble* xx,IssmDouble* yy,IssmDouble* zz,IssmDouble eartharea)=0;
|
|---|
| 323 | #endif
|
|---|
| 324 |
|
|---|
| 325 | };
|
|---|
| 326 | #endif
|
|---|
![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
