Changeset 21431
- Timestamp:
- 12/05/16 03:51:12 (8 years ago)
- Location:
- issm/trunk-jpl/src/c/analyses
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
r21430 r21431 405 405 bool active_element; 406 406 int domaintype; 407 Element* basalelement ;407 Element* basalelement=NULL; 408 408 409 409 /*Get basal element*/ … … 414 414 break; 415 415 case Domain3DEnum: 416 if(!element->IsOnBase()) return NULL;416 if(!element->IsOnBase()) return; 417 417 basalelement = element->SpawnBasalElement(); 418 418 break; … … 437 437 438 438 /*Use the dof list to index into the solution vector: */ 439 for(i =0;i<numnodes;i++){439 for(int i=0;i<numnodes;i++){ 440 440 eplHeads[i]=solution[doflist[i]]; 441 441 if(xIsNan<IssmDouble>(eplHeads[i])) _error_("NaN found in solution vector"); … … 653 653 654 654 bool active_element; 655 int i,j;656 655 int domaintype; 657 656 IssmDouble h_max; 658 657 IssmDouble sedheadmin; 659 Element* basalelement ;658 Element* basalelement=NULL; 660 659 661 660 /*Get basal element*/ … … 666 665 break; 667 666 case Domain3DEnum: 668 if(!element->IsOnBase()) return NULL;667 if(!element->IsOnBase()) return; 669 668 basalelement = element->SpawnBasalElement(); 670 669 break; … … 693 692 /*Get minimum sediment head of the element*/ 694 693 sedheadmin=sedhead[0]; 695 for(i =1;i<numnodes;i++) if(sedhead[i]<=sedheadmin)sedheadmin=sedhead[i];696 for(i =0;i<numnodes;i++){694 for(int i=1;i<numnodes;i++) if(sedhead[i]<=sedheadmin)sedheadmin=sedhead[i]; 695 for(int i=0;i<numnodes;i++){ 697 696 /*If node is now closed bring its thickness back to initial*/ 698 697 if (old_active[i]==0.){ … … 718 717 GetHydrologyDCInefficientHmax(&h_max,basalelement,basalelement->nodes[i]); 719 718 if(eplhead[i]>=h_max && active_element){ 720 for( j=0;j<numnodes;j++){719 for(int j=0;j<numnodes;j++){ 721 720 /*Increase of the domain is on the downstream node in term of sediment head*/ 722 721 if(sedhead[j] == sedheadmin){ … … 743 742 /*Constants*/ 744 743 int domaintype; 745 Element* basalelement ;744 Element* basalelement=NULL; 746 745 747 746 /*Get basal element*/ … … 752 751 break; 753 752 case Domain3DEnum: 754 if(!element->IsOnBase()) return NULL;753 if(!element->IsOnBase()) return; 755 754 basalelement = element->SpawnBasalElement(); 756 755 break; -
issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp
r21430 r21431 425 425 /*Intermediaries*/ 426 426 int domaintype; 427 Element* basalelement ;427 Element* basalelement=NULL; 428 428 bool converged; 429 429 int* doflist = NULL; … … 436 436 break; 437 437 case Domain3DEnum: 438 if(!element->IsOnBase()) return NULL;438 if(!element->IsOnBase()) return; 439 439 basalelement = element->SpawnBasalElement(); 440 440 break;
Note:
See TracChangeset
for help on using the changeset viewer.