Changeset 21431
- Timestamp:
 - 12/05/16 03:51:12 (9 years ago)
 - Location:
 - issm/trunk-jpl/src/c/analyses
 - Files:
 - 
      
- 2 edited
 
- 
          
  HydrologyDCEfficientAnalysis.cpp (modified) (9 diffs)
 - 
          
  HydrologyDCInefficientAnalysis.cpp (modified) (2 diffs)
 
 
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.
  ![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)