Changeset 15130 for issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
- Timestamp:
- 05/28/13 15:16:02 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
r15129 r15130 2429 2429 /*}}}*/ 2430 2430 /*FUNCTION Penta::NodalValue {{{*/ 2431 int Penta::NodalValue(IssmDouble* pvalue, int index, int natureofdataenum ,bool process_units){2431 int Penta::NodalValue(IssmDouble* pvalue, int index, int natureofdataenum){ 2432 2432 2433 2433 int i; … … 2984 2984 2985 2985 /*Get for Vx and Vy, the max of abs value: */ 2986 this->MaxAbsVx(&maxabsvx ,false);2987 this->MaxAbsVy(&maxabsvy ,false);2988 this->MaxAbsVz(&maxabsvz ,false);2986 this->MaxAbsVx(&maxabsvx); 2987 this->MaxAbsVy(&maxabsvy); 2988 this->MaxAbsVz(&maxabsvz); 2989 2989 2990 2990 /* Get node coordinates and dof list: */ … … 3266 3266 /*}}}*/ 3267 3267 /*FUNCTION Penta::MinVel{{{*/ 3268 void Penta::MinVel(IssmDouble* pminvel , bool process_units){3268 void Penta::MinVel(IssmDouble* pminvel){ 3269 3269 3270 3270 /*Get minimum:*/ … … 3276 3276 /*}}}*/ 3277 3277 /*FUNCTION Penta::MinVx{{{*/ 3278 void Penta::MinVx(IssmDouble* pminvx , bool process_units){3278 void Penta::MinVx(IssmDouble* pminvx){ 3279 3279 3280 3280 /*Get minimum:*/ … … 3286 3286 /*}}}*/ 3287 3287 /*FUNCTION Penta::MinVy{{{*/ 3288 void Penta::MinVy(IssmDouble* pminvy , bool process_units){3288 void Penta::MinVy(IssmDouble* pminvy){ 3289 3289 3290 3290 /*Get minimum:*/ … … 3296 3296 /*}}}*/ 3297 3297 /*FUNCTION Penta::MinVz{{{*/ 3298 void Penta::MinVz(IssmDouble* pminvz , bool process_units){3298 void Penta::MinVz(IssmDouble* pminvz){ 3299 3299 3300 3300 /*Get minimum:*/ … … 3306 3306 /*}}}*/ 3307 3307 /*FUNCTION Penta::MassFlux {{{*/ 3308 IssmDouble Penta::MassFlux( IssmDouble* segment ,bool process_units){3308 IssmDouble Penta::MassFlux( IssmDouble* segment){ 3309 3309 3310 3310 IssmDouble mass_flux=0; … … 3318 3318 /*Spawn Tria element from the base of the Penta: */ 3319 3319 Tria* tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria. 3320 mass_flux=tria->MassFlux(segment ,process_units);3320 mass_flux=tria->MassFlux(segment); 3321 3321 delete tria->material; delete tria; 3322 3322 … … 3330 3330 /*}}}*/ 3331 3331 /*FUNCTION Penta::MaxAbsVx{{{*/ 3332 void Penta::MaxAbsVx(IssmDouble* pmaxabsvx , bool process_units){3332 void Penta::MaxAbsVx(IssmDouble* pmaxabsvx){ 3333 3333 3334 3334 /*Get maximum:*/ … … 3340 3340 /*}}}*/ 3341 3341 /*FUNCTION Penta::MaxAbsVy{{{*/ 3342 void Penta::MaxAbsVy(IssmDouble* pmaxabsvy , bool process_units){3342 void Penta::MaxAbsVy(IssmDouble* pmaxabsvy){ 3343 3343 3344 3344 /*Get maximum:*/ … … 3350 3350 /*}}}*/ 3351 3351 /*FUNCTION Penta::MaxAbsVz{{{*/ 3352 void Penta::MaxAbsVz(IssmDouble* pmaxabsvz , bool process_units){3352 void Penta::MaxAbsVz(IssmDouble* pmaxabsvz){ 3353 3353 3354 3354 /*Get maximum:*/ … … 3360 3360 /*}}}*/ 3361 3361 /*FUNCTION Penta::MaxVel{{{*/ 3362 void Penta::MaxVel(IssmDouble* pmaxvel , bool process_units){3362 void Penta::MaxVel(IssmDouble* pmaxvel){ 3363 3363 3364 3364 /*Get maximum:*/ … … 3371 3371 /*}}}*/ 3372 3372 /*FUNCTION Penta::MaxVx{{{*/ 3373 void Penta::MaxVx(IssmDouble* pmaxvx , bool process_units){3373 void Penta::MaxVx(IssmDouble* pmaxvx){ 3374 3374 3375 3375 /*Get maximum:*/ … … 3381 3381 /*}}}*/ 3382 3382 /*FUNCTION Penta::MaxVy{{{*/ 3383 void Penta::MaxVy(IssmDouble* pmaxvy , bool process_units){3383 void Penta::MaxVy(IssmDouble* pmaxvy){ 3384 3384 3385 3385 /*Get maximum:*/ … … 3391 3391 /*}}}*/ 3392 3392 /*FUNCTION Penta::MaxVz{{{*/ 3393 void Penta::MaxVz(IssmDouble* pmaxvz , bool process_units){3393 void Penta::MaxVz(IssmDouble* pmaxvz){ 3394 3394 3395 3395 /*Get maximum:*/ … … 5415 5415 /*}}}*/ 5416 5416 /*FUNCTION Penta::SurfaceAverageVelMisfit {{{*/ 5417 IssmDouble Penta::SurfaceAverageVelMisfit( bool process_units,int weight_index){5417 IssmDouble Penta::SurfaceAverageVelMisfit(int weight_index){ 5418 5418 5419 5419 int approximation; … … 5438 5438 * and compute SurfaceAverageVelMisfit*/ 5439 5439 tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria (lower face). 5440 J=tria->SurfaceAverageVelMisfit( process_units,weight_index);5440 J=tria->SurfaceAverageVelMisfit(weight_index); 5441 5441 delete tria->material; delete tria; 5442 5442 return J; … … 5445 5445 5446 5446 tria=(Tria*)SpawnTria(3,4,5); //nodes 3, 4 and 5 make the new tria (upper face). 5447 J=tria->SurfaceAverageVelMisfit( process_units,weight_index);5447 J=tria->SurfaceAverageVelMisfit(weight_index); 5448 5448 delete tria->material; delete tria; 5449 5449 return J; … … 5452 5452 /*}}}*/ 5453 5453 /*FUNCTION Penta::SurfaceAbsVelMisfit {{{*/ 5454 IssmDouble Penta::SurfaceAbsVelMisfit( bool process_units,int weight_index){5454 IssmDouble Penta::SurfaceAbsVelMisfit(int weight_index){ 5455 5455 5456 5456 int approximation; … … 5475 5475 * and compute SurfaceAbsVelMisfit*/ 5476 5476 tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria (lower face). 5477 J=tria->SurfaceAbsVelMisfit( process_units,weight_index);5477 J=tria->SurfaceAbsVelMisfit(weight_index); 5478 5478 delete tria->material; delete tria; 5479 5479 return J; … … 5482 5482 5483 5483 tria=(Tria*)SpawnTria(3,4,5); //nodes 3, 4 and 5 make the new tria (upper face). 5484 J=tria->SurfaceAbsVelMisfit( process_units,weight_index);5484 J=tria->SurfaceAbsVelMisfit(weight_index); 5485 5485 delete tria->material; delete tria; 5486 5486 return J; … … 5489 5489 /*}}}*/ 5490 5490 /*FUNCTION Penta::SurfaceLogVelMisfit {{{*/ 5491 IssmDouble Penta::SurfaceLogVelMisfit( bool process_units,int weight_index){5491 IssmDouble Penta::SurfaceLogVelMisfit(int weight_index){ 5492 5492 5493 5493 int approximation; … … 5512 5512 * and compute SurfaceLogVelMisfit*/ 5513 5513 tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria (lower face). 5514 J=tria->SurfaceLogVelMisfit( process_units,weight_index);5514 J=tria->SurfaceLogVelMisfit(weight_index); 5515 5515 delete tria->material; delete tria; 5516 5516 return J; … … 5519 5519 5520 5520 tria=(Tria*)SpawnTria(3,4,5); //nodes 3, 4 and 5 make the new tria (upper face). 5521 J=tria->SurfaceLogVelMisfit( process_units,weight_index);5521 J=tria->SurfaceLogVelMisfit(weight_index); 5522 5522 delete tria->material; delete tria; 5523 5523 return J; … … 5526 5526 /*}}}*/ 5527 5527 /*FUNCTION Penta::SurfaceLogVxVyMisfit {{{*/ 5528 IssmDouble Penta::SurfaceLogVxVyMisfit( bool process_units,int weight_index){5528 IssmDouble Penta::SurfaceLogVxVyMisfit(int weight_index){ 5529 5529 5530 5530 IssmDouble J; … … 5551 5551 * and compute SurfaceLogVxVyMisfit*/ 5552 5552 tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria (lower face). 5553 J=tria->SurfaceLogVxVyMisfit( process_units,weight_index);5553 J=tria->SurfaceLogVxVyMisfit(weight_index); 5554 5554 delete tria->material; delete tria; 5555 5555 return J; … … 5558 5558 5559 5559 tria=(Tria*)SpawnTria(3,4,5); //nodes 3, 4 and 5 make the new tria (upper face). 5560 J=tria->SurfaceLogVxVyMisfit( process_units,weight_index);5560 J=tria->SurfaceLogVxVyMisfit(weight_index); 5561 5561 delete tria->material; delete tria; 5562 5562 return J; … … 5565 5565 /*}}}*/ 5566 5566 /*FUNCTION Penta::SurfaceRelVelMisfit {{{*/ 5567 IssmDouble Penta::SurfaceRelVelMisfit( bool process_units,int weight_index){5567 IssmDouble Penta::SurfaceRelVelMisfit(int weight_index){ 5568 5568 5569 5569 int approximation; … … 5588 5588 * and compute SurfaceRelVelMisfit*/ 5589 5589 tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria (lower face). 5590 J=tria->SurfaceRelVelMisfit( process_units,weight_index);5590 J=tria->SurfaceRelVelMisfit(weight_index); 5591 5591 delete tria->material; delete tria; 5592 5592 return J; … … 5595 5595 5596 5596 tria=(Tria*)SpawnTria(3,4,5); //nodes 3, 4 and 5 make the new tria (upper face). 5597 J=tria->SurfaceRelVelMisfit( process_units,weight_index);5597 J=tria->SurfaceRelVelMisfit(weight_index); 5598 5598 delete tria->material; delete tria; 5599 5599 return J; … … 5602 5602 /*}}}*/ 5603 5603 /*FUNCTION Penta::ThicknessAbsGradient{{{*/ 5604 IssmDouble Penta::ThicknessAbsGradient( bool process_units,int weight_index){5604 IssmDouble Penta::ThicknessAbsGradient(int weight_index){ 5605 5605 5606 5606 _error_("Not implemented yet"); … … 5608 5608 /*}}}*/ 5609 5609 /*FUNCTION Penta::ThicknessAbsMisfit {{{*/ 5610 IssmDouble Penta::ThicknessAbsMisfit( bool process_units,int weight_index){5610 IssmDouble Penta::ThicknessAbsMisfit(int weight_index){ 5611 5611 5612 5612 int approximation; … … 5622 5622 5623 5623 tria=(Tria*)SpawnTria(0,1,2); 5624 J=tria->ThicknessAbsMisfit( process_units,weight_index);5624 J=tria->ThicknessAbsMisfit(weight_index); 5625 5625 delete tria->material; delete tria; 5626 5626 return J; … … 5628 5628 /*}}}*/ 5629 5629 /*FUNCTION Penta::DragCoefficientAbsGradient{{{*/ 5630 IssmDouble Penta::DragCoefficientAbsGradient( bool process_units,int weight_index){5630 IssmDouble Penta::DragCoefficientAbsGradient(int weight_index){ 5631 5631 5632 5632 IssmDouble J; … … 5637 5637 5638 5638 tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria 5639 J=tria->DragCoefficientAbsGradient( process_units,weight_index);5639 J=tria->DragCoefficientAbsGradient(weight_index); 5640 5640 delete tria->material; delete tria; 5641 5641 return J; … … 5643 5643 /*}}}*/ 5644 5644 /*FUNCTION Penta::RheologyBbarAbsGradient{{{*/ 5645 IssmDouble Penta::RheologyBbarAbsGradient( bool process_units,int weight_index){5645 IssmDouble Penta::RheologyBbarAbsGradient(int weight_index){ 5646 5646 5647 5647 IssmDouble J; … … 5652 5652 5653 5653 tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria 5654 J=tria->RheologyBbarAbsGradient( process_units,weight_index);5654 J=tria->RheologyBbarAbsGradient(weight_index); 5655 5655 delete tria->material; delete tria; 5656 5656 return J;
Note:
See TracChangeset
for help on using the changeset viewer.