Changeset 7107 for issm/trunk/src/c/objects/Elements/Tria.cpp
- Timestamp:
- 01/14/11 14:11:33 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/objects/Elements/Tria.cpp
r7097 r7107 4402 4402 void Tria::MigrateGroundingLine(void){ 4403 4403 4404 4405 double *values = NULL; 4406 double h[3],s[3],b[3],ba[3]; 4407 double bed_hydro; 4408 double rho_water,rho_ice,density; 4409 int isonshelf[3]; 4404 4410 int i; 4405 4406 double h[3];4407 double s[3];4408 double b[3];4409 double ba[3];4410 double bed_hydro;4411 int isonshelf[3];4412 Input *surface_input = NULL;4413 Input *bed_input = NULL;4414 double rho_water;4415 double rho_ice;4416 double density;4417 double *values = NULL;4418 4411 bool elementonshelf = false; 4419 4412 4420 4413 4421 4414 /*Recover info at the vertices: */ 4422 surface_input =inputs->GetInput(SurfaceEnum);_assert_(surface_input);4423 bed_input =inputs->GetInput(BedEnum);_assert_(bed_input);4415 Input* surface_input =inputs->GetInput(SurfaceEnum); _assert_(surface_input); 4416 Input* bed_input =inputs->GetInput(BedEnum); _assert_(bed_input); 4424 4417 if((surface_input->Enum()!=TriaVertexInputEnum) | (bed_input->Enum()!=TriaVertexInputEnum))_error_(" not supported yet for bed and surface interpolations not P1!"); 4425 4418 … … 4465 4458 /*Surface and bed are updated. Update inputs:*/ 4466 4459 surface_input->GetValuesPtr(&values,NULL); for(i=0;i<3;i++)values[i]=s[i]; 4467 bed_input->GetValuesPtr(&values,NULL); for(i=0;i<3;i++)values[i]=b[i];4460 bed_input->GetValuesPtr(&values,NULL); for(i=0;i<3;i++)values[i]=b[i]; 4468 4461 4469 4462 }
Note:
See TracChangeset
for help on using the changeset viewer.