source: issm/trunk/src/c/modules/ModelProcessorx/Balancevelocities/CreateNodesBalancedvelocities.cpp@ 8287

Last change on this file since 8287 was 8287, checked in by Mathieu Morlighem, 14 years ago

moving Balanced* to Balance* TBC
File size: 1.9 KB
Line 
1/*
2 * CreateNodesBalancevelocities.c:
3 */
4
5#include "../../../Container/Container.h"
6#include "../../../toolkits/toolkits.h"
7#include "../../../io/io.h"
8#include "../../../EnumDefinitions/EnumDefinitions.h"
9#include "../../../objects/objects.h"
10#include "../../../shared/shared.h"
11#include "../../MeshPartitionx/MeshPartitionx.h"
12#include "../../../include/include.h"
13#include "../ModelProcessorx.h"
14
15void CreateNodesBalancevelocities(Nodes** pnodes, IoModel* iomodel,ConstDataHandle iomodel_handle){
16
17 /*Intermediary*/
18 int i;
19 bool continuous_galerkin=true;
20
21 /*DataSets: */
22 Nodes* nodes = NULL;
23
24 /*Recover pointer: */
25 nodes=*pnodes;
26
27 /*Create nodes if they do not exist yet*/
28 if(!nodes) nodes = new Nodes(NodesEnum);
29
30 /*Continuous Galerkin partition of nodes: */
31 NodesPartitioning(&iomodel->my_nodes,iomodel->my_elements,iomodel->my_vertices,iomodel,iomodel_handle,continuous_galerkin);
32
33 /*First fetch data: */
34 IoModelFetchData(&iomodel->gridonbed,NULL,NULL,iomodel_handle,"gridonbed");
35 IoModelFetchData(&iomodel->gridonsurface,NULL,NULL,iomodel_handle,"gridonsurface");
36 IoModelFetchData(&iomodel->gridonicesheet,NULL,NULL,iomodel_handle,"gridonicesheet");
37 IoModelFetchData(&iomodel->gridoniceshelf,NULL,NULL,iomodel_handle,"gridoniceshelf");
38 IoModelFetchData(&iomodel->vertices_type,NULL,NULL,iomodel_handle,"vertices_type");
39 IoModelFetchData(&iomodel->gridonwater,NULL,NULL,iomodel_handle,"gridonwater");
40 for (i=0;i<iomodel->numberofvertices;i++){
41
42 if(iomodel->my_vertices[i]){
43
44 /*Add node to nodes dataset: */
45 nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,i+1,i,iomodel,BalancevelocitiesAnalysisEnum));
46
47 }
48 }
49
50 /*Clean fetched data: */
51 xfree((void**)&iomodel->gridonbed);
52 xfree((void**)&iomodel->gridonsurface);
53 xfree((void**)&iomodel->gridonicesheet);
54 xfree((void**)&iomodel->gridoniceshelf);
55 xfree((void**)&iomodel->gridonwater);
56 xfree((void**)&iomodel->vertices_type);
57
58 /*Assign output pointer: */
59 *pnodes=nodes;
60}
Note: See TracBrowser for help on using the repository browser.