Changeset 821


Ignore:
Timestamp:
06/05/09 15:07:39 (16 years ago)
Author:
seroussi
Message:

no parameters needed

Location:
issm/trunk/src/c/UpdateGeometryx
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/UpdateGeometryx/UpdateGeometryx.cpp

    r816 r821  
    1515
    1616void 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,
    1818                Vec vec_newthickness,Vec vec_bed,Vec vec_surface){
    1919
     
    6363                }
    6464        }
     65
    6566        /*recover material parameters: */
    6667        rho_ice=matpar->GetRhoIce();
     
    9192                //For grids on ice shelt, we have hydrostatic equilibrium (for now)
    9293                if (node->IsOnShelf()){
    93 
    9494                        b=-rho_ice/rho_water*h;
    9595                        s=(1-rho_ice/rho_water)*h;
    96 
    9796                }
    9897
     
    101100                VecSetValues(outbed,1,&dof,&b,INSERT_VALUES);
    102101                VecSetValues(outsurface,1,&dof,&s,INSERT_VALUES);
    103 
    104102        }
    105103
     
    118116        xfree((void**)&bed);
    119117        xfree((void**)&surface);
    120 
    121118       
    122119        /*Assign output pointers: */
  • issm/trunk/src/c/UpdateGeometryx/UpdateGeometryx.h

    r816 r821  
    1111/* local prototypes: */
    1212void 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,,
    1414                Vec newthickness,Vec bed,Vec surface);
    1515
Note: See TracChangeset for help on using the changeset viewer.