Changeset 16725


Ignore:
Timestamp:
11/12/13 18:04:16 (11 years ago)
Author:
bdef
Message:

NEW: There is now an initial thickness for the EPL

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

Legend:

Unmodified
Added
Removed

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

    r16690 r16725  
    5757        iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateEnum);
    5858        iomodel->FetchDataToInput(elements,EplHeadEnum);
    59         iomodel->FetchDataToInput(elements,HydrologydcEplThicknessEnum);
     59        iomodel->FetchDataToInput(elements,HydrologydcEplInitialThicknessEnum);
    6060       
    6161}/*}}}*/

  • issm/trunk-jpl/src/c/classes/Elements/Penta.cpp

    r16721 r16725  
    27242724                                name==EplHeadOldEnum ||
    27252725                                name==HydrologydcEplThicknessOldEnum ||
     2726                                name==HydrologydcEplInitialThicknessEnum ||
    27262727                                name==HydrologydcEplThicknessEnum ||
    27272728                                name==HydrologydcMaskEplactiveEnum ||
    … …  
    1087310874        const int   numdof2d = NDOF1 *NUMVERTICES2D;
    1087410875        bool        isefficientlayer;
    10875         IssmDouble  n,A,dt;
     10876        IssmDouble  n,A,dt,init_thick;
    1087610877        IssmDouble  rho_water,rho_ice;
    1087710878        IssmDouble  gravity,latentheat,EPLgrad;
    … …  
    1089610897                latentheat       = matpar->GetLatentHeat();
    1089710898                epl_conductivity = matpar->GetEplConductivity();
     10899                init_thick       = matpar->GetEplInitialThickness();
    1089810900                n                = material->GetN();
    1089910901                A                = material->GetA();
    10900                
    1090110902               
    1090210903                GetInputListOnVertices(&activeEpl[0],HydrologydcMaskEplactiveEnum);
    … …  
    1091110912                        /*Keeping thickness to 1 if EPL is not active*/
    1091210913                        if(activeEpl[i]==0.0){
    10913                                 thickness[i]=1.0;
     10914                                thickness[i]=init_thick;
    1091410915                                thickness[i+numdof2d]=thickness[i];
    1091510916                        }

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

    r16721 r16725  
    21662166                                name==EplHeadEnum ||
    21672167                                name==HydrologydcEplThicknessOldEnum ||
     2168                                name==HydrologydcEplInitialThicknessEnum ||
    21682169                                name==HydrologydcEplThicknessEnum ||
    21692170                                name==HydrologydcMaskEplactiveEnum ||
    … …  
    73707371        const int   numdof         = NDOF1 *NUMVERTICES;
    73717372        bool        isefficientlayer;
    7372         IssmDouble  n,A,dt;
     7373        IssmDouble  n,A,dt,init_thick;
    73737374        IssmDouble  rho_water,rho_ice;
    73747375        IssmDouble  gravity,latentheat,EPLgrad;
    … …  
    73867387        if(isefficientlayer){
    73877388                /*For now, assuming just one way to compute EPL thickness*/
    7388                 rho_water  = matpar->GetRhoWater();
    7389                 rho_ice    = matpar->GetRhoIce();
    7390                 gravity    = matpar->GetG();
    7391                 latentheat = matpar->GetLatentHeat();
    7392                 n          = material->GetN();
    7393                 A          = material->GetAbar();
     7389                rho_water        = matpar->GetRhoWater();
     7390                rho_ice          = matpar->GetRhoIce();
     7391                gravity          = matpar->GetG();
     7392                latentheat       = matpar->GetLatentHeat();
     7393                epl_conductivity = matpar->GetEplConductivity();
     7394                init_thick       = matpar->GetEplInitialThickness();
     7395                n                = material->GetN();
     7396                A                = material->GetAbar();
    73947397               
    73957398                GetInputListOnVertices(&activeEpl[0],HydrologydcMaskEplactiveEnum);
    … …  
    74017404                GetInputListOnVertices(&bed[0],BedEnum);
    74027405               
    7403                 epl_conductivity = matpar->GetEplConductivity();
    7404                
    74057406                for(int i=0;i<numdof;i++){
    74067407                        /*Keeping thickness to 1 if EPL is not active*/
    74077408                        if(activeEpl[i]==0.0){
    7408                                 thickness[i]=1.0;
     7409                                thickness[i]=init_thick;
    74097410                        }
    74107411                        else{

  • issm/trunk-jpl/src/c/classes/Materials/Matpar.cpp

    r16697 r16725  
    6666                                iomodel->Constant(&this->epl_compressibility,HydrologydcEplCompressibilityEnum);
    6767                                iomodel->Constant(&this->epl_porosity,HydrologydcEplPorosityEnum);
     68                                iomodel->Constant(&this->epl_init_thickness,HydrologydcEplInitialThicknessEnum);
    6869                                iomodel->Constant(&this->epl_conductivity,HydrologydcEplConductivityEnum);
    6970                }
    … …  
    392393}               
    393394/*}}}*/
     395/*FUNCTION Matpar::GetEplInitialThickness {{{*/
     396IssmDouble Matpar::GetEplInitialThickness(){
     397        return epl_init_thickness;               
     398}               
     399/*}}}*/
    394400/*FUNCTION Matpar::GetEplConductivity {{{*/
    395401IssmDouble Matpar::GetEplConductivity(){

  • issm/trunk-jpl/src/c/classes/Materials/Matpar.h

    r16697 r16725  
    4848
    4949                IssmDouble  epl_compressibility;
    50                 IssmDouble  epl_porosity;       
     50                IssmDouble  epl_porosity;
     51                IssmDouble  epl_init_thickness;
    5152                IssmDouble  epl_conductivity;   
    5253
    … …  
    125126                IssmDouble GetSedimentThickness();
    126127                IssmDouble GetEplConductivity();
     128                IssmDouble GetEplInitialThickness();
    127129                IssmDouble TMeltingPoint(IssmDouble pressure);
    128130                IssmDouble PureIceEnthalpy(IssmDouble pressure);

  • issm/trunk-jpl/src/c/cores/hydrology_core.cpp

    r16638 r16725  
    8484                        if (isefficientlayer){
    8585                                InputDuplicatex(femmodel,EplHeadEnum,EplHeadOldEnum);
    86                                 InputDuplicatex(femmodel,HydrologydcEplThicknessEnum,HydrologydcEplThicknessOldEnum);
     86                                InputDuplicatex(femmodel,HydrologydcEplInitialThicknessEnum,HydrologydcEplThicknessEnum);
     87                                InputDuplicatex(femmodel,HydrologydcEplInitialThicknessEnum,HydrologydcEplThicknessOldEnum);
    8788                        }
    8889                       
Note: See TracChangeset for help on using the changeset viewer.