Changeset 15968
- Timestamp:
- 08/27/13 09:39:47 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/modules/GroundinglineMigrationx/GroundinglineMigrationx.cpp
r15967 r15968 12 12 int migration_style,analysis_type; 13 13 IssmDouble *vertices_potentially_ungrounding = NULL; 14 IssmDouble * vertices_ungrounding= NULL;14 IssmDouble *phi_ungrounding = NULL; 15 15 Element *element = NULL; 16 16 … … 37 37 38 38 /*propagate ice shelf into connex areas of the ice sheet that potentially want to unground: */ 39 vertices_ungrounding=PropagateFloatingiceToGroundedNeighbors(elements,nodes,vertices,parameters,vertices_potentially_ungrounding);39 phi_ungrounding=PropagateFloatingiceToGroundedNeighbors(elements,nodes,vertices,parameters,vertices_potentially_ungrounding); 40 40 } 41 41 … … 43 43 for(int i=0;i<elements->Size();i++){ 44 44 element=dynamic_cast<Element*>(elements->GetObjectByOffset(i)); 45 element->MigrateGroundingLine( vertices_ungrounding);45 element->MigrateGroundingLine(phi_ungrounding); 46 46 } 47 47 48 48 /*free ressouces: */ 49 49 xDelete<IssmDouble>(vertices_potentially_ungrounding); 50 xDelete<IssmDouble>( vertices_ungrounding);50 xDelete<IssmDouble>(phi_ungrounding); 51 51 } 52 52 … … 54 54 IssmDouble* PotentialUngrounding(Elements* elements,Vertices* vertices,Parameters* parameters){ 55 55 56 int i,numberofvertices;57 IssmDouble* vertices_potentially_ungrounding = NULL;58 Vector<IssmDouble>* 59 Element* element = NULL;56 int i,numberofvertices; 57 IssmDouble* vertices_potentially_ungrounding = NULL; 58 Vector<IssmDouble>* vec_vertices_potentially_ungrounding = NULL; 59 Element* element = NULL; 60 60 61 61 /*Initialize vector with number of vertices*/ … … 80 80 /*FUNCTION PropagateFloatingiceToGroundedNeighbors {{{*/ 81 81 IssmDouble* PropagateFloatingiceToGroundedNeighbors(Elements* elements,Nodes* nodes,Vertices* vertices,Parameters* parameters,IssmDouble* vertices_potentially_ungrounding){ 82 83 int i,analysis_type; 84 int numberofvertices; 85 int nflipped,local_nflipped; 86 IssmDouble* phi = NULL; 87 IssmDouble* elements_neighboring_floatingce = NULL; 88 Vector<IssmDouble>* vec_elements_neighboring_floatingice = NULL; 89 Vector<IssmDouble>* vec_phi = NULL; 90 Node* node = NULL; 91 Element* element = NULL; 82 int i,analysis_type,numberofvertices; 83 int nflipped,local_nflipped; 84 IssmDouble* phi = NULL; 85 IssmDouble* elements_neighboring_floatingce = NULL; 86 Vector<IssmDouble>* vec_elements_neighboring_floatingice = NULL; 87 Vector<IssmDouble>* vec_phi = NULL; 88 Node* node = NULL; 89 Element* element = NULL; 92 90 93 91 /*recover parameters: */
Note:
See TracChangeset
for help on using the changeset viewer.