[17802] | 1 | Index: ../trunk-jpl/src/c/classes/Elements/Tria.cpp
|
---|
| 2 | ===================================================================
|
---|
| 3 | --- ../trunk-jpl/src/c/classes/Elements/Tria.cpp (revision 16867)
|
---|
| 4 | +++ ../trunk-jpl/src/c/classes/Elements/Tria.cpp (revision 16868)
|
---|
| 5 | @@ -7027,6 +7027,7 @@
|
---|
| 6 | const int numdof = NDOF1 *NUMVERTICES;
|
---|
| 7 | int *doflist = NULL;
|
---|
| 8 | bool isefficientlayer;
|
---|
| 9 | + bool active_element;
|
---|
| 10 | int transfermethod;
|
---|
| 11 | IssmDouble sed_trans,sed_thick;
|
---|
| 12 | IssmDouble leakage,h_max;
|
---|
| 13 | @@ -7035,7 +7036,7 @@
|
---|
| 14 | IssmDouble epl_specificstoring[numdof],sedstoring[numdof];
|
---|
| 15 | IssmDouble epl_head[numdof],sed_head[numdof];
|
---|
| 16 |
|
---|
| 17 | - Input* active_element=NULL;
|
---|
| 18 | + Input* active_element_input=NULL;
|
---|
| 19 |
|
---|
| 20 | GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
|
---|
| 21 |
|
---|
| 22 | @@ -7053,9 +7054,8 @@
|
---|
| 23 | break;
|
---|
| 24 | case 1:
|
---|
| 25 |
|
---|
| 26 | - int analysis_type;
|
---|
| 27 | - parameters->FindParam(&analysis_type,AnalysisTypeEnum);
|
---|
| 28 | - active_element=inputs->GetInput(HydrologydcMaskEplactiveEltEnum); _assert_(active_element);
|
---|
| 29 | + active_element_input=inputs->GetInput(HydrologydcMaskEplactiveEltEnum); _assert_(active_element_input);
|
---|
| 30 | + active_element_input->GetInputValue(&active_element);
|
---|
| 31 |
|
---|
| 32 | GetInputListOnVertices(&activeEpl[0],HydrologydcMaskEplactiveNodeEnum);
|
---|
| 33 | GetInputListOnVertices(&sed_head[0],SedimentHeadEnum);
|
---|
| 34 | @@ -7067,12 +7067,8 @@
|
---|
| 35 | sed_trans = matpar->GetSedimentTransmitivity();
|
---|
| 36 | sed_thick = matpar->GetSedimentThickness();
|
---|
| 37 |
|
---|
| 38 | - // printf("Which analyse %i \n",analysis_type);
|
---|
| 39 | -
|
---|
| 40 | if(!active_element){
|
---|
| 41 |
|
---|
| 42 | - //printf("NoneActive\n");
|
---|
| 43 | -
|
---|
| 44 | /*No transfer if the EPL is not active*/
|
---|
| 45 | for(int i=0;i<numdof;i++){
|
---|
| 46 | wh_trans=0.0;
|
---|
| 47 | @@ -7100,7 +7096,7 @@
|
---|
| 48 | }
|
---|
| 49 | /*Assign output pointer*/
|
---|
| 50 | transfer->SetValue(doflist[i],wh_trans,INS_VAL);
|
---|
| 51 | - //printf("Computing, %e \n",wh_trans);
|
---|
| 52 | +
|
---|
| 53 | }
|
---|
| 54 | }
|
---|
| 55 | break;
|
---|