Changeset 18641


Ignore:
Timestamp:
10/14/14 15:09:05 (10 years ago)
Author:
bdef
Message:

CHG: cleaning print

Location:
issm/trunk-jpl/src/c/analyses
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp

    r18633 r18641  
    352352        active_element_input->GetInputValue(&active_element);
    353353
    354         /* if(!active_element){ */
    355         /*      /\*Keeping thickness to initial value if EPL is not active*\/ */
    356         /*      for(i=0;i<numnodes;i++){ */
    357         /*              eplHeads[i]=sedhead[i]; */
    358         /*      } */
    359         /* } */
    360         /* else{ */
    361                 /*Use the dof list to index into the solution vector: */
     354        /*Use the dof list to index into the solution vector: */
    362355        for(i=0;i<numnodes;i++){
    363356                eplHeads[i]=solution[doflist[i]];
    364357                if(xIsNan<IssmDouble>(eplHeads[i])) _error_("NaN found in solution vector");
    365358        }
    366                 //      }
    367359        /*Add input to the element: */
    368360        element->AddBasalInput(EplHeadEnum,eplHeads,P1Enum);
    … …  
    439431                _error_("no case higher than 1 for the Transfer method");
    440432        }
    441         if(element->Id()==42){
    442                 printf("Transferefficient Kmat %e, %e, %e, %i \n",transfer,sediment_head,epl_head,element->Id());
    443         }
     433        /* if(element->Id()==42){ */
     434        /*      printf("Transferefficient Kmat %e, %e, %e, %e \n",transfer,sediment_head,epl_head,sed_residual); */
     435        /* } */
    444436        return transfer;
    445437}/*}}}*/
    … …  
    485477                _error_("no case higher than 1 for the Transfer method");
    486478        }
    487         if(element->Id()==42){
    488                 printf("Transferefficient Pvec %e, %e, %e\n",transfer,sediment_head,epl_head);
    489         }
     479        /* if(element->Id()==42){ */
     480        /*      printf("Transferefficient Pvec %e, %e, %e, %e\n",transfer,sediment_head,epl_head,sed_residual); */
     481        /* } */
    490482        return transfer;
    491483}/*}}}*/

  • issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp

    r18633 r18641  
    197197        basalelement ->FindParam(&dt,TimesteppingTimeStepEnum);
    198198        basalelement ->FindParam(&isefficientlayer,HydrologydcIsefficientlayerEnum);
    199         Input* SedTrans_input    = basalelement->GetInput(HydrologydcSedimentTransmitivityEnum); _assert_(SedTrans_input);
    200         Input* sed_head_input    = basalelement->GetInput(SedimentHeadEnum);
    201         Input* epl_head_input    = basalelement->GetInput(EplHeadEnum);
    202         Input* residual_input  = basalelement->GetInput(SedimentHeadResidualEnum);
     199        Input* SedTrans_input = basalelement->GetInput(HydrologydcSedimentTransmitivityEnum); _assert_(SedTrans_input);
     200        Input* sed_head_input = basalelement->GetInput(SedimentHeadEnum);
     201        Input* epl_head_input = basalelement->GetInput(EplHeadEnum);
     202        Input* base_input     = basalelement->GetInput(BaseEnum);
     203        Input* thick_input    = basalelement->GetInput(ThicknessEnum);
    203204
    204205        IssmDouble sediment_storing = SedimentStoring(basalelement);
    … …  
    238239                                active_element_input->GetInputValue(&active_element);
    239240                                if(active_element){
    240                                         transfer=GetHydrologyKMatrixTransfer(basalelement,gauss,sed_head_input,epl_head_input,residual_input);
     241                                        transfer=GetHydrologyKMatrixTransfer(basalelement,gauss,sed_head_input,epl_head_input,thick_input,base_input);
    241242                                        basalelement->NodalFunctions(&basis[0],gauss);
    242243                                        D_scalar=transfer*gauss->weight*Jdet*dt;
    … …  
    300301        basalelement->FindParam(&isefficientlayer,HydrologydcIsefficientlayerEnum);
    301302
    302         Input* sed_head_input    = basalelement->GetInput(SedimentHeadEnum);
    303         Input* epl_head_input    = basalelement->GetInput(EplHeadEnum);
    304         Input* water_input       = basalelement->GetInput(BasalforcingsGroundediceMeltingRateEnum); _assert_(water_input);
    305         Input* residual_input    = basalelement->GetInput(SedimentHeadResidualEnum);
     303        Input* sed_head_input = basalelement->GetInput(SedimentHeadEnum);
     304        Input* epl_head_input = basalelement->GetInput(EplHeadEnum);
     305        Input* water_input    = basalelement->GetInput(BasalforcingsGroundediceMeltingRateEnum); _assert_(water_input);
     306        Input* base_input     = basalelement->GetInput(BaseEnum);
     307        Input* thick_input    = basalelement->GetInput(ThicknessEnum);
    306308        if(dt!= 0.){old_wh_input = basalelement->GetInput(SedimentHeadOldEnum);                     _assert_(old_wh_input);}
    307309
    … …  
    338340                                active_element_input->GetInputValue(&active_element);
    339341                                if(active_element){
    340                                         transfer=GetHydrologyPVectorTransfer(basalelement,gauss,sed_head_input,epl_head_input,residual_input);
     342                                        transfer=GetHydrologyPVectorTransfer(basalelement,gauss,sed_head_input,epl_head_input,thick_input,base_input);
    341343                                }
    342344                                else{
    … …  
    568570/*}}}*/
    569571
    570 IssmDouble HydrologyDCInefficientAnalysis::GetHydrologyKMatrixTransfer(Element* element, Gauss* gauss, Input* sed_head_input, Input* epl_head_input, Input* residual_input){/*{{{*/
     572IssmDouble HydrologyDCInefficientAnalysis::GetHydrologyKMatrixTransfer(Element* element, Gauss* gauss, Input* sed_head_input, Input* epl_head_input, Input* thick_input, Input* base_input){/*{{{*/
    571573
    572574        int transfermethod;
    573         IssmDouble sed_residual;
     575        IssmDouble hmax;
    574576        IssmDouble epl_head,sediment_head;
    575577        IssmDouble leakage,transfer;
    … …  
    586588                _assert_(sed_head_input);
    587589                _assert_(epl_head_input);
    588                 _assert_(residual_input);
    589590               
    590591                sed_head_input->GetInputValue(&sediment_head,gauss);
    591592                epl_head_input->GetInputValue(&epl_head,gauss);
    592                 residual_input->GetInputValue(&sed_residual,gauss);
    593593                element->FindParam(&leakage,HydrologydcLeakageFactorEnum);
     594
     595                hmax=GetHydrologyDCInefficientHmax(element, gauss, thick_input, base_input);
    594596               
    595597                if(epl_head>sediment_head){ 
    596                         if(sediment_head>=400.0){
     598                        if(sediment_head>=hmax){
    597599                                transfer=0.0;
     600                        }
     601                        else if{
     602                                transfer=(leakage);
    598603                        }
    599604                        else{
    … …  
    608613                _error_("no case higher than 1 for the Transfer method");
    609614        }
    610         if(element->Id()==42){
    611                 printf("TransferInfficient Kmat %e, %e, %e, %i\n",transfer,sediment_head,epl_head,element->Id());
    612         }
     615        /* if(element->Id()==42){ */
     616        /*      printf("TransferInfficient Kmat %e, %e, %e\n",transfer,sediment_head,epl_head); */
     617        /* } */
    613618        return transfer;
    614619}/*}}}*/
    615620
    616 IssmDouble HydrologyDCInefficientAnalysis::GetHydrologyPVectorTransfer(Element* element, Gauss* gauss, Input* sed_head_input, Input* epl_head_input, Input* residual_input){/*{{{*/
     621IssmDouble HydrologyDCInefficientAnalysis::GetHydrologyPVectorTransfer(Element* element, Gauss* gauss, Input* sed_head_input, Input* epl_head_input, Input* thick_input, Input* base_input){/*{{{*/
    617622
    618623        int transfermethod;
    619         IssmDouble sed_residual;
     624        IssmDouble hmax;
    620625        IssmDouble epl_head,sediment_head;
    621626        IssmDouble leakage,transfer;
    … …  
    632637                _assert_(sed_head_input);
    633638                _assert_(epl_head_input);
    634                 _assert_(residual_input);
    635639               
    636640                sed_head_input->GetInputValue(&sediment_head,gauss);
    637641                epl_head_input->GetInputValue(&epl_head,gauss);
    638                 residual_input->GetInputValue(&sed_residual,gauss);
    639642                element->FindParam(&leakage,HydrologydcLeakageFactorEnum);
    640643
    641                 if(epl_head>sediment_head){ 
    642                         if(sediment_head>=400.0){
     644                hmax=GetHydrologyDCInefficientHmax(element, gauss, thick_input, base_input);
     645               
     646                if(epl_head>sediment_head){ 
     647                        if(sediment_head>=hmax){
    643648                                transfer=0.0;
    644649                        }
    … …  
    654659                _error_("no case higher than 1 for the Transfer method");
    655660        }
    656         if(element->Id()==42){
    657                 printf("TransferInfficient Pvec %e, %e, %e, %i\n",transfer,sediment_head,epl_head,element->Id());
    658         }
     661        /* if(element->Id()==42){ */
     662        /*      printf("TransferInfficient Pvec %e, %e, %e\n",transfer,sediment_head,epl_head); */
     663        /* } */
    659664        return transfer;
    660665}/*}}}*/

  • issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.h

    r18633 r18641  
    3838                IssmDouble GetHydrologyDCInefficientHmax(Element* element, Gauss* gauss, Input* thickness_input, Input* base_input);
    3939                void GetHydrologyDCInefficientHmax(IssmDouble* ph_max,Element* element, Node* innode);
    40                 IssmDouble GetHydrologyKMatrixTransfer(Element* element, Gauss* gauss, Input* sed_head_input, Input* epl_head_input, Input* residual_input);
    41                 IssmDouble GetHydrologyPVectorTransfer(Element* element, Gauss* gauss, Input* sed_head_input, Input* epl_head_input, Input* residual_input);
     40                IssmDouble GetHydrologyKMatrixTransfer(Element* element, Gauss* gauss, Input* sed_head_input, Input* epl_head_input, Input* thick_input, Input* base_input);
     41                IssmDouble GetHydrologyPVectorTransfer(Element* element, Gauss* gauss, Input* sed_head_input, Input* epl_head_input, Input* thick_input, Input* base_input);
    4242                void ElementizeEplMask(FemModel* femmodel);
    4343};
Note: See TracChangeset for help on using the changeset viewer.