source: issm/trunk/src/c/classes/Elements/Tria.h@ 23189

Last change on this file since 23189 was 23189, checked in by Mathieu Morlighem, 7 years ago

merged trunk-jpl and trunk for revision 23187
File size: 12.9 KB
Line 
1/*! \file Tria.h
2 * \brief: header file for tria object
3 */
4
5#ifndef _TRIA_H_
6#define _TRIA_H_
7
8/*Headers:*/
9/*{{{*/
10#include "./Element.h"
11#include "./ElementHook.h"
12#include "./TriaRef.h"
13class Parameters;
14class Inputs;
15class IoModel;
16class Results;
17class Node;
18class Material;
19class Matpar;
20class Seg;
21class ElementMatrix;
22class ElementVector;
23class Vertex;
24class GaussTria;
25
26#include "../../shared/Exceptions/exceptions.h"
27#include "../../shared/Enum/Enum.h"
28/*}}}*/
29
30class Tria: public Element,public ElementHook,public TriaRef{
31
32 public:
33
34 /*Tria constructors, destructors {{{*/
35 Tria(){};
36 Tria(int tria_id,int tria_sid,IoModel* iomodel,int nummodels);
37 ~Tria();
38 /*}}}*/
39 /*Object virtual functions definitions:{{{ */
40 Object *copy();
41 void Marshall(char** pmarshalled_data,int* pmarshalled_data_size, int marshall_direction);
42 int ObjectEnum();
43 /*}}}*/
44 /*Update virtual functions resolution: {{{*/
45 #ifdef _HAVE_DAKOTA_
46 void InputUpdateFromMatrixDakota(IssmDouble* matrix, int nows, int ncols, int name, int type);
47 void InputUpdateFromVectorDakota(IssmDouble* vector, int name, int type);
48 #endif
49 void InputUpdateFromIoModel(int index, IoModel* iomodel);
50 void InputUpdateFromVector(IssmDouble* vector, int name, int type);
51 /*}}}*/
52 /*Element virtual functions definitions: {{{*/
53 void AverageOntoPartition(Vector<IssmDouble>* partition_contributions,Vector<IssmDouble>* partition_areas,IssmDouble* vertex_response,IssmDouble* qmu_part);
54 void CalvingRateVonmises();
55 void CalvingCrevasseDepth();
56 void CalvingRateLevermann();
57 IssmDouble CharacteristicLength(void);
58 void ComputeBasalStress(Vector<IssmDouble>* sigma_b);
59 void ComputeDeviatoricStressTensor();
60 void ComputeEsaStrainAndVorticity();
61 void ComputeSigmaNN();
62 void ComputeStressTensor();
63 void ComputeSurfaceNormalVelocity();
64 void Configure(Elements* elements,Loads* loads,Nodes* nodesin,Vertices* verticesin,Materials* materials,Parameters* parameters);
65 void ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int N, int M);
66 void ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index);
67 void ControlToVectors(Vector<IssmPDouble>* vector_control, Vector<IssmPDouble>* vector_gradient,int control_enum);
68 void CreateDistanceInputFromSegmentlist(IssmDouble* distances,int distanceenum);
69 int EdgeOnBaseIndex();
70 void EdgeOnBaseIndices(int* pindex1,int* pindex);
71 int EdgeOnSurfaceIndex();
72 void EdgeOnSurfaceIndices(int* pindex1,int* pindex);
73 void ElementResponse(IssmDouble* presponse,int response_enum);
74 void ElementSizes(IssmDouble* hx,IssmDouble* hy,IssmDouble* hz);
75 int FiniteElement(void);
76 IssmDouble FloatingArea(bool scaled);
77 void FSContactMigration(Vector<IssmDouble>* vertex_sigmann,Vector<IssmDouble>* vertex_waterpressure);
78 Element* GetBasalElement(void){_error_("not implemented yet");};
79 void GetLevelsetPositivePart(int* point1,IssmDouble* fraction1,IssmDouble* fraction2, bool* mainlynegative,IssmDouble* levelsetvalues);
80 void GetGroundedPart(int* point1,IssmDouble* fraction1, IssmDouble* fraction2,bool* mainlyfloating);
81 IssmDouble GetGroundedPortion(IssmDouble* xyz_list);
82 void GetIcefrontCoordinates(IssmDouble** pxyz_front,IssmDouble* xyz_list,int levelsetenum);
83 void GetLevelCoordinates(IssmDouble** pxyz_front,IssmDouble* xyz_list,int levelsetenum,IssmDouble level);
84 int GetNodeIndex(Node* node);
85 int GetNumberOfNodes(void);
86 int GetNumberOfNodes(int enum_type);
87 int GetNumberOfVertices(void);
88 void GetVectorFromControlInputs(Vector<IssmDouble>* gradient,int control_enum,int control_index,const char* data,int offset);
89 void GetVectorFromControlInputs(Vector<IssmDouble>* gradient,int control_enum,int control_index,const char* data);
90 void GetVerticesCoordinatesBase(IssmDouble** pxyz_list);
91 void GetVerticesCoordinatesTop(IssmDouble** pxyz_list);
92 IssmDouble GroundedArea(bool scaled);
93 bool HasEdgeOnBase();
94 bool HasEdgeOnSurface();
95 IssmDouble IceVolume(bool scaled);
96 IssmDouble IceVolumeAboveFloatation(bool scaled);
97 void InputDepthAverageAtBase(int enum_type,int average_enum_type);
98 void InputExtrude(int enum_type,int start){_error_("not implemented"); /*For penta only*/};
99 bool IsFaceOnBoundary(void);
100 bool IsIcefront(void);
101 bool IsNodeOnShelfFromFlags(IssmDouble* flags);
102 bool IsOnBase();
103 bool IsOnSurface();
104 bool IsZeroLevelset(int levelset_enum);
105 IssmDouble Masscon(IssmDouble* levelset);
106 IssmDouble MassFlux(IssmDouble* segment);
107 IssmDouble MassFlux(IssmDouble x1,IssmDouble y1, IssmDouble x2, IssmDouble y2,int segment_id);
108 void MaterialUpdateFromTemperature(void){_error_("not implemented yet");};
109 IssmDouble Misfit(int modelenum,int observationenum,int weightsenum);
110 IssmDouble MisfitArea(int weightsenum);
111 int NodalValue(IssmDouble* pvalue, int index, int natureofdataenum);
112 int NumberofNodesPressure(void);
113 int NumberofNodesVelocity(void);
114 void PicoUpdateBoxid(int* pmax_boxid_basin);
115 void PicoUpdateBox(int loopboxid);
116 void PicoComputeBasalMelt(void);
117 void PotentialUngrounding(Vector<IssmDouble>* potential_sheet_ungrounding);
118 int PressureInterpolation();
119 void ReduceMatrices(ElementMatrix* Ke,ElementVector* pe);
120 void ResetFSBasalBoundaryCondition(void);
121 void ResetHooks();
122 void ResetLevelsetFromSegmentlist(IssmDouble* segments,int numsegments);
123 void SetControlInputsFromVector(IssmDouble* vector,int control_enum,int control_index,int offset,int N,int M);
124 void SetControlInputsFromVector(IssmDouble* vector,int control_enum,int control_index);
125 void SetCurrentConfiguration(Elements* elements,Loads* loads,Nodes* nodes,Materials* materials,Parameters* parameters);
126 Element* SpawnBasalElement(void);
127 Element* SpawnTopElement(void);
128 void StrainRateparallel();
129 void StrainRateperpendicular();
130 void StressIntensityFactor(void){_error_("not implemented yet");};
131 IssmDouble SurfaceArea(void);
132 int TensorInterpolation();
133 IssmDouble TimeAdapt();
134 IssmDouble TotalFloatingBmb(bool scaled);
135 IssmDouble TotalGroundedBmb(bool scaled);
136 IssmDouble TotalSmb(bool scaled);
137 void Update(int index, IoModel* iomodel,int analysis_counter,int analysis_type,int finitelement);
138 int UpdatePotentialUngrounding(IssmDouble* vertices_potentially_ungrounding,Vector<IssmDouble>* vec_nodes_on_iceshelf,IssmDouble* nodes_on_iceshelf);
139 void ValueP1DerivativesOnGauss(IssmDouble* dvalue,IssmDouble* values,IssmDouble* xyz_list,Gauss* gauss);
140 void ValueP1OnGauss(IssmDouble* pvalue,IssmDouble* values,Gauss* gauss);
141 int VelocityInterpolation();
142 int VertexConnectivity(int vertexindex);
143 void VerticalSegmentIndices(int** pindices,int* pnumseg){_error_("not implemented yet");};
144 void VerticalSegmentIndicesBase(int** pindices,int* pnumseg){_error_("not implemented yet");};
145 void WriteFieldIsovalueSegment(DataSet* segments,int fieldenum,IssmDouble fieldvalue);
146
147 #ifdef _HAVE_GIAIVINS_
148 void GiaDeflection(Vector<IssmDouble>* wg,Vector<IssmDouble>* dwgdt,IssmDouble* x,IssmDouble* y);
149 #endif
150 #ifdef _HAVE_ESA_
151 void EsaGeodetic2D(Vector<IssmDouble>* pUp,Vector<IssmDouble>* pNorth,Vector<IssmDouble>* pEast, Vector<IssmDouble>* pX,Vector<IssmDouble>* pY,IssmDouble* xx,IssmDouble* yy);
152 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);
153 #endif
154 #ifdef _HAVE_SEALEVELRISE_
155 IssmDouble OceanArea(void);
156 IssmDouble OceanAverage(IssmDouble* Sg);
157 void SealevelriseMomentOfInertia(IssmDouble* dI_list,IssmDouble* Sg_old,IssmDouble eartharea);
158 void SealevelriseEustatic(Vector<IssmDouble>* pSgi,IssmDouble* peustatic,IssmDouble* latitude,IssmDouble* longitude,IssmDouble* radius,IssmDouble oceanarea,IssmDouble eartharea);
159 void SealevelriseNonEustatic(Vector<IssmDouble>* pSgo,IssmDouble* Sg_old,IssmDouble* latitude,IssmDouble* longitude,IssmDouble* radius,IssmDouble eartharea);
160 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,int horiz);
161 #endif
162 /*}}}*/
163 /*Tria specific routines:{{{*/
164 void AddBasalInput(int input_enum, IssmDouble* values, int interpolation_enum);
165 void AddInput(int input_enum, IssmDouble* values, int interpolation_enum);
166 void AddControlInput(int input_enum, IssmDouble* values,IssmDouble* values_min,IssmDouble* values_max, int interpolation_enum,int id);
167 void DatasetInputCreate(IssmDouble* array,int M,int N,int* individual_enums,int num_inputs,IoModel* iomodel,int input_enum);
168 IssmDouble GetArea(void);
169 IssmDouble GetHorizontalSurfaceArea(void);
170 IssmDouble GetArea3D(void);
171 IssmDouble GetAreaIce(void);
172 IssmDouble GetAreaSpherical(void);
173 void GetAreaCoordinates(IssmDouble *area_coordinates,IssmDouble* xyz_zero,IssmDouble* xyz_list,int numpoints);
174 int GetElementType(void);
175 void GetInputValue(IssmDouble* pvalue,Node* node,int enumtype);
176 void GetLevelsetIntersection(int** pindices, int* pnumiceverts, IssmDouble* fraction, int levelset_enum, IssmDouble level);
177 void GetMaterialInputValue(IssmDouble* pvalue,Node* node,int enumtype);
178 Node* GetNode(int node_number);
179 void InputUpdateFromSolutionOneDof(IssmDouble* solution,int enum_type);
180 void InputUpdateFromSolutionOneDofCollapsed(IssmDouble* solution,int enum_type){_error_("not implemented yet");};
181 void JacobianDeterminant(IssmDouble* pJdet, IssmDouble* xyz_list,Gauss* gauss);
182 void JacobianDeterminantBase(IssmDouble* pJdet,IssmDouble* xyz_list_base,Gauss* gauss);
183 void JacobianDeterminantLine(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss){_error_("not implemented yet");};
184 void JacobianDeterminantSurface(IssmDouble* pJdet, IssmDouble* xyz_list,Gauss* gauss);
185 void JacobianDeterminantTop(IssmDouble* pJdet,IssmDouble* xyz_list_base,Gauss* gauss);
186 IssmDouble MinEdgeLength(IssmDouble* xyz_list){_error_("not implemented yet");};
187 Gauss* NewGauss(void);
188 Gauss* NewGauss(int order);
189 Gauss* NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order);
190 Gauss* NewGauss(int point1,IssmDouble fraction1,IssmDouble fraction2,bool mainlyfloating,int order);
191 Gauss* NewGauss(int point1,IssmDouble fraction1,IssmDouble fraction2,int order);
192 Gauss* NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order_horiz,int order_vert);
193 Gauss* NewGaussBase(int order);
194 Gauss* NewGaussLine(int vertex1,int vertex2,int order){_error_("not implemented yet");};
195 Gauss* NewGaussTop(int order);
196 void NodalFunctions(IssmDouble* basis,Gauss* gauss);
197 void NodalFunctionsDerivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss);
198 void NodalFunctionsDerivativesVelocity(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss);
199 void NodalFunctionsMINIDerivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss){_error_("not implemented yet");};
200 void NodalFunctionsPressure(IssmDouble* basis,Gauss* gauss);
201 void NodalFunctionsP1(IssmDouble* basis,Gauss* gauss);
202 void NodalFunctionsP1Derivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss);
203 void NodalFunctionsP2(IssmDouble* basis,Gauss* gauss);
204 void NodalFunctionsTensor(IssmDouble* basis,Gauss* gauss);
205 void NodalFunctionsVelocity(IssmDouble* basis,Gauss* gauss);
206 void NormalBase(IssmDouble* normal,IssmDouble* xyz_list);
207 void NormalSection(IssmDouble* normal,IssmDouble* xyz_list);
208 void NormalTop(IssmDouble* normal,IssmDouble* xyz_list);
209 void SetClone(int* minranks);
210 void SetTemporaryElementType(int element_type_in){_error_("not implemented yet");};
211 Seg* SpawnSeg(int index1,int index2);
212 IssmDouble StabilizationParameter(IssmDouble u, IssmDouble v, IssmDouble w, IssmDouble diameter, IssmDouble kappa){_error_("not implemented yet");};
213 void UpdateConstraintsExtrudeFromBase(void);
214 void UpdateConstraintsExtrudeFromTop(void);
215 void ViscousHeating(IssmDouble* pphi,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input){_error_("not implemented yet");};
216 /*}}}*/
217
218};
219#endif /* _TRIA_H */
Note: See TracBrowser for help on using the repository browser.