- Timestamp:
- 02/09/15 13:55:09 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
r19094 r19095 667 667 default: _error_("mesh "<<EnumToStringx(domaintype)<<" not supported yet"); 668 668 } 669 670 669 /*Intermediaries*/ 671 672 int eplflip_lock;673 670 int numnodes =basalelement->GetNumberOfNodes(); 674 671 IssmDouble* epl_thickness =xNew<IssmDouble>(numnodes); … … 684 681 active_element_input->GetInputValue(&active_element); 685 682 686 basalelement->parameters->FindParam(&eplflip_lock,HydrologydcEplflipLockEnum);687 688 683 basalelement-> GetInputListOnVertices(&old_active[0],HydrologydcMaskEplactiveNodeEnum); 689 684 basalelement-> GetInputListOnVertices(&epl_thickness[0],HydrologydcEplThicknessEnum); … … 695 690 sedheadmin=sedhead[0]; 696 691 for(i=1;i<numnodes;i++) if(sedhead[i]<=sedheadmin)sedheadmin=sedhead[i]; 697 698 692 for(i=0;i<numnodes;i++){ 699 700 693 /*If node is now closed bring its thickness back to initial*/ 701 694 if (old_active[i]==0.){
Note:
See TracChangeset
for help on using the changeset viewer.