Ignore:
Timestamp:
06/17/13 16:23:06 (12 years ago)
Author:
bdef
Message:

CHG: So far the last working version of the hydrological Code (Still some differencies with ELmer/Ice results)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

    r15273 r15276  
    65266526                        GetInputListOnVertices(&sed_head[0],SedimentHeadEnum);
    65276527                        GetInputListOnVertices(&epl_head[0],EplHeadEnum);
    6528 
     6528                       
    65296529                        this->parameters->FindParam(&leakage,HydrologydcLeakageFactorEnum);
    65306530
     
    65346534                        for(int i=0;i<numdof;i++){
    65356535                                if(activeEpl[i]==0.0)continue;
    6536 
     6536                               
    65376537                                eplstoring[i]=matpar->GetEplStoring();         
    65386538                                sedstoring[i]=matpar->GetSedimentStoring();
    6539 
     6539                               
    65406540                                if(sed_head[i]>epl_head[i]){
    65416541                                        wh_trans=sed_trans*sedstoring[i]*(epl_head[i]-sed_head[i])/(leakage*sed_thick);                         
     
    65476547                                                wh_trans=0.0;
    65486548                                        }
    6549                                         //                                      if((sed_head[i]+(wh_trans/sedstoring[i]))>h_max){
    6550                                         //(h_max-sed_head[i])*eplstoring[i];
    6551                                         //}
    65526549                                }
    65536550                                /*Assign output pointer*/
Note: See TracChangeset for help on using the changeset viewer.