Changeset 821
- Timestamp:
- 06/05/09 15:07:39 (16 years ago)
- Location:
- issm/trunk/src/c/UpdateGeometryx
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/UpdateGeometryx/UpdateGeometryx.cpp
r816 r821 15 15 16 16 void UpdateGeometryx(Vec* poutthickness,Vec* poutbed,Vec* poutsurface, 17 DataSet* elements, DataSet* nodes,DataSet* loads, DataSet* materials, DataSet* parameters,17 DataSet* elements, DataSet* nodes,DataSet* loads, DataSet* materials, 18 18 Vec vec_newthickness,Vec vec_bed,Vec vec_surface){ 19 19 … … 63 63 } 64 64 } 65 65 66 /*recover material parameters: */ 66 67 rho_ice=matpar->GetRhoIce(); … … 91 92 //For grids on ice shelt, we have hydrostatic equilibrium (for now) 92 93 if (node->IsOnShelf()){ 93 94 94 b=-rho_ice/rho_water*h; 95 95 s=(1-rho_ice/rho_water)*h; 96 97 96 } 98 97 … … 101 100 VecSetValues(outbed,1,&dof,&b,INSERT_VALUES); 102 101 VecSetValues(outsurface,1,&dof,&s,INSERT_VALUES); 103 104 102 } 105 103 … … 118 116 xfree((void**)&bed); 119 117 xfree((void**)&surface); 120 121 118 122 119 /*Assign output pointers: */ -
issm/trunk/src/c/UpdateGeometryx/UpdateGeometryx.h
r816 r821 11 11 /* local prototypes: */ 12 12 void UpdateGeometryx(Vec* poutthickness,Vec* poutbed,Vec* poutsurface, 13 DataSet* elements, DataSet* nodes,DataSet* loads, DataSet* materials, DataSet* parameters,13 DataSet* elements, DataSet* nodes,DataSet* loads, DataSet* materials,, 14 14 Vec newthickness,Vec bed,Vec surface); 15 15
Note:
See TracChangeset
for help on using the changeset viewer.