Changeset 7673 for issm/trunk/src/c/objects/Elements/Tria.cpp
- Timestamp:
- 03/22/11 15:17:24 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/objects/Elements/Tria.cpp
r7672 r7673 4450 4450 // for(j=0;j<3;j++)thickness[j]=values[j]; 4451 4451 4452 printf("Tria::InputUpdateFromVectorDakota,id=%d\n",this->Id());4453 printf("nodes: %d,%d,%d\n",this->nodes[0]->Sid()+1,this->nodes[1]->Sid()+1,this->nodes[2]->Sid())+1;4454 printf("init thickness: %e,%e,%e\n",thickness_init[0],thickness_init[1],thickness_init[2]);4455 printf("thickness_coeff: %e,%e,%e\n",thickness_coeff[0],thickness_coeff[1],thickness_coeff[2]);4456 printf("values: %e,%e,%e\n",values[0],values[1],values[2]);4457 printf("old bed: %e,%e,%e\n",bed[0],bed[1],bed[2]);4458 printf("old surface: %e,%e,%e\n",surface[0],surface[1],surface[2]);4459 4460 4452 /*build new bed and surface: */ 4461 4453 if (this->IsOnShelf()){ 4462 printf("shelf\n");4463 4454 /*hydrostatic equilibrium: */ 4464 4455 double rho_ice,rho_water,di; … … 4489 4480 } 4490 4481 else{ 4491 printf("sheet\n");4492 4482 /*build new thickness: */ 4493 4483 for (j=0; j<3; j++) { … … 4504 4494 for(j=0;j<3;j++)bed[j]=surface[j]-thickness[j]; 4505 4495 } 4506 4507 printf("new thickness: %e,%e,%e\n",thickness[0],thickness[1],thickness[2]);4508 printf("new bed: %e,%e,%e\n",bed[0],bed[1],bed[2]);4509 printf("new surface: %e,%e,%e\n",surface[0],surface[1],surface[2]);4510 4496 4511 4497 /*Add new inputs: */
Note:
See TracChangeset
for help on using the changeset viewer.