Changeset 21489
- Timestamp:
- 01/13/17 04:35:51 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
r21488 r21489 674 674 recurence->SetValue(basalelement->nodes[i]->Sid(),1.,INS_VAL); 675 675 } 676 /* If epl head gets under base elevation, close the layer */676 /* //If epl head gets under base elevation, close the layer */ 677 677 /* else if(eplhead[i]<(base[i]-1.0e-8)){ */ 678 678 /* vec_mask->SetValue(basalelement->nodes[i]->Sid(),0.,INS_VAL); */ … … 683 683 if (old_active[i]==0.){ 684 684 epl_thickness[i]=init_thick; 685 }686 /*Activate if we have a residual from sediment*/687 if(residual[i]>0.){688 vec_mask->SetValue(basalelement->nodes[i]->Sid(),1.,INS_VAL);689 if(old_active[i]==0.){690 recurence->SetValue(basalelement->nodes[i]->Sid(),1.,INS_VAL);685 /*Activate if we have a residual from sediment*/ 686 if(residual[i]>0.){ 687 vec_mask->SetValue(basalelement->nodes[i]->Sid(),1.,INS_VAL); 688 if(old_active[i]==0.){ 689 recurence->SetValue(basalelement->nodes[i]->Sid(),1.,INS_VAL); 690 } 691 691 } 692 692 } … … 703 703 } 704 704 } 705 706 705 } 707 706 basalelement->AddInput(HydrologydcEplThicknessEnum,epl_thickness,basalelement->GetElementType()); … … 745 744 /* epl_storing=storing; */ 746 745 /* } */ 746 /* return epl_storing; */ 747 747 return storing; 748 //return epl_storing;749 748 }/*}}}*/ 750 749 … … 761 760 762 761 epl_transmitivity=epl_conductivity*epl_thickness; 763 / * epl_transmitivity=max(1.0e-6,(epl_conductivity*min(water_sheet,epl_thickness))); */762 //epl_transmitivity=max(1.0e-6,(epl_conductivity*min(water_sheet,epl_thickness))); 764 763 return epl_transmitivity; 765 764 }/*}}}*/
Note:
See TracChangeset
for help on using the changeset viewer.