Changeset 17964


Ignore:
Timestamp:
05/08/14 13:09:13 (11 years ago)
Author:
bdef
Message:

CHG: change in transfer computation

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

Legend:

Unmodified
Added
Removed

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

    r17950 r17964  
    288288                        thickness_input ->GetInputValue(&epl_thickness,gauss);
    289289                       
    290                         /*Dealing with the sediment part of the transfer term*/
     290                        /*Dealing with the epl part of the transfer term*/
    291291                        transfer=GetHydrologyPVectorTransfer(basalelement,gauss,thickness_input,sed_head_input,epl_head_input,sed_trans_input,residual_input);
    292292                        scalar = Jdet*gauss->weight*((water_head*epl_specificstoring*epl_thickness)+(transfer*dt));
    … …  
    410410                residual_input->GetInputValue(&residual,gauss);
    411411                element->FindParam(&leakage,HydrologydcLeakageFactorEnum);
    412 
    413                 if(epl_head>sed_head){
    414                         if(residual>0.0){       
    415                                 transfer=0.0;
    416                         }
    417                         else{
    418                                 transfer=(sediment_transmitivity)/(sediment_thickness*leakage);
    419                                 //                              transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity)/(sediment_thickness*leakage);
    420                         }
    421                 }
    422                 else{
    423                         transfer=(sediment_transmitivity)/(sediment_thickness*leakage);
    424                         //                      transfer=(sediment_storing*sediment_transmitivity)/(sediment_thickness*leakage);
    425                 }
     412                transfer=(sediment_transmitivity)/(sediment_thickness*leakage);
     413                /* if(epl_head>sed_head){ */
     414                /*      if(residual>0.0){        */
     415                /*              transfer=0.0; */
     416                /*      } */
     417                /*      else{ */
     418                /*              transfer=(sediment_transmitivity)/(sediment_thickness*leakage); */
     419                /*              //transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity)/(sediment_thickness*leakage); */
     420                /*      } */
     421                /* } */
     422                /* else{ */
     423                /*      transfer=(sediment_transmitivity)/(sediment_thickness*leakage); */
     424                /*      //transfer=(sediment_storing*sediment_transmitivity)/(sediment_thickness*leakage); */
     425                /* } */
    426426                break;
    427427        default:
    … …  
    463463                residual_input->GetInputValue(&residual,gauss);
    464464                element->FindParam(&leakage,HydrologydcLeakageFactorEnum);
    465                
    466                 if(epl_head>sediment_head){
    467                         if(residual>0.0){       
    468                                 transfer=0.0;
    469                         }
    470                         else{
    471                                 transfer=(sediment_transmitivity*sediment_head)/(sediment_thickness*leakage);
    472                                 //                              transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity*sediment_head)/(sediment_thickness*leakage);
    473                         }
    474                 }
    475                 else{
    476                         transfer=(sediment_transmitivity*sediment_head)/(sediment_thickness*leakage);
    477                         //                      transfer=(sediment_storing*sediment_transmitivity*sediment_head)/(sediment_thickness*leakage);
    478                 }
     465                transfer=(sediment_transmitivity*sediment_head)/(sediment_thickness*leakage);
     466                /* if(epl_head>sediment_head){ */
     467                /*      if(residual>0.0){        */
     468                /*              transfer=0.0; */
     469                /*      } */
     470                /*      else{ */
     471                /*              transfer=(sediment_transmitivity*sediment_head)/(sediment_thickness*leakage); */
     472                /*              //transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity*sediment_head)/(sediment_thickness*leakage); */
     473                /*      } */
     474                /* } */
     475                /* else{ */
     476                /*      transfer=(sediment_transmitivity*sediment_head)/(sediment_thickness*leakage); */
     477                /*      //transfer=(sediment_storing*sediment_transmitivity*sediment_head)/(sediment_thickness*leakage); */
     478                /* } */
    479479                break;
    480480        default:

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

    r17950 r17964  
    329329
    330330                scalar = Jdet*gauss->weight*(water_load);
     331
    331332                if(dt!=0.) scalar = scalar*dt;
    332333                for(int i=0;i<numnodes;i++){
    … …  
    601602                sed_trans_input->GetInputValue(&sediment_transmitivity,gauss);
    602603                element->FindParam(&leakage,HydrologydcLeakageFactorEnum);
    603        
    604                 if(epl_head>sed_head){
    605                         h_max=GetHydrologyDCInefficientHmax(element,gauss,thickness_input,base_input);
    606                         if(sed_head>=h_max){
    607                                 transfer=0.0;
    608                         }
    609                         else{
    610                                 transfer=(sediment_transmitivity)/(sediment_thickness*leakage);
    611                                 //                              transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity)/(sediment_thickness*leakage);
    612                         }
    613                 }
    614                 else{
    615                         transfer=(sediment_transmitivity)/(sediment_thickness*leakage);
    616                         //                      transfer=(sediment_storing*sediment_transmitivity)/(sediment_thickness*leakage);
    617                 }
     604                transfer=(sediment_transmitivity)/(sediment_thickness*leakage);
     605                /* if(epl_head>sed_head){ */
     606                /*      h_max=GetHydrologyDCInefficientHmax(element,gauss,thickness_input,base_input); */
     607                /*      if(sed_head>=h_max){ */
     608                /*              transfer=0.0; */
     609                /*      } */
     610                /*      else{ */
     611                /*              transfer=(sediment_transmitivity)/(sediment_thickness*leakage); */
     612                /*              //transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity)/(sediment_thickness*leakage); */
     613                /*      } */
     614                /* } */
     615                /* else{ */
     616                /*      transfer=(sediment_transmitivity)/(sediment_thickness*leakage); */
     617                /*      //transfer=(sediment_storing*sediment_transmitivity)/(sediment_thickness*leakage); */
     618                /* } */
    618619                break;
    619620        default:
    … …  
    658659
    659660                element->FindParam(&leakage,HydrologydcLeakageFactorEnum);
    660 
    661                 if(epl_head>sediment_head){
    662                         h_max=GetHydrologyDCInefficientHmax(element,gauss,thickness_input,base_input);
    663                         if(sediment_head>=h_max){
    664                                 transfer=0.0;
    665                         }
    666                         else{
    667                                 transfer=(sediment_transmitivity*epl_head)/(sediment_thickness*leakage);
    668                                 //                              transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity*epl_head)/(sediment_thickness*leakage);
    669                         }
    670                 }
    671                 else{
    672                         transfer=(sediment_transmitivity*epl_head)/(sediment_thickness*leakage);
    673                         //                      transfer=(sediment_storing*sediment_transmitivity*epl_head)/(sediment_thickness*leakage);
    674                 }
     661                transfer=(sediment_transmitivity*epl_head)/(sediment_thickness*leakage);
     662                /* if(epl_head>sediment_head){ */
     663                /*      h_max=GetHydrologyDCInefficientHmax(element,gauss,thickness_input,base_input); */
     664                /*      if(sediment_head>=h_max){ */
     665                /*              transfer=0.0; */
     666                /*      } */
     667                /*      else{ */
     668                /*              transfer=(sediment_transmitivity*epl_head)/(sediment_thickness*leakage); */
     669                /*              //transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity*epl_head)/(sediment_thickness*leakage); */
     670                /*      } */
     671                /* } */
     672                /* else{ */
     673                /*      transfer=(sediment_transmitivity*epl_head)/(sediment_thickness*leakage); */
     674                /*      //transfer=(sediment_storing*sediment_transmitivity*epl_head)/(sediment_thickness*leakage); */
     675                /* } */
    675676                break;
    676677        default:

  • issm/trunk-jpl/src/c/classes/Loads/Pengrid.cpp

    r17585 r17964  
    646646                if(zigzag_counter>penalty_lock){
    647647                        unstable=0;
    648                         active=0;
     648                        active=1;
    649649                }
    650650        }
Note: See TracChangeset for help on using the changeset viewer.