Changeset 19161


Ignore:
Timestamp:
02/26/15 09:36:07 (10 years ago)
Author:
bdef
Message:

BUG: fix for the friction coupling

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

Legend:

Unmodified
Added
Removed

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

    r19002 r19161  
    20772077        /*Intermediaries*/
    20782078        int      domaintype,dim;
     2079
    20792080        Element* basalelement;
    20802081
    … …  
    21412142        Input* cinput = basalelement->GetInput(FrictionCEnum);
    21422143        Input* Asinput = basalelement->GetInput(FrictionAsEnum);
    2143         Input* Ninput = basalelement->GetInput(FrictionEffectivePressureEnum);
    21442144        Input* nInput =basalelement->GetInput(MaterialsRheologyNEnum);
    2145        
     2145        Input* Ninput = basalelement->GetInput(FrictionEffectivePressureEnum); 
    21462146        /* Start  looping on the number of gaussian points: */
    21472147        Gauss* gauss=basalelement->NewGauss(4);

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

    r18930 r19161  
    109109        bool dakota_analysis,islevelset,isenthalpy;
    110110        int frictionlaw;
    111 
     111        int FrictionCoupling;
     112       
    112113        /*Now, is the model 3d? otherwise, do nothing: */
    113114        if(iomodel->domaintype==Domain2DhorizontalEnum)return;
    … …  
    175176                        break;
    176177                case 3:
     178                        iomodel->Constant(&FrictionCoupling,FrictionCouplingEnum);
    177179                        iomodel->FetchDataToInput(elements,FrictionCEnum);
    178180                        iomodel->FetchDataToInput(elements,FrictionAsEnum);
    179181                        iomodel->FetchDataToInput(elements,FrictionQEnum);
    180                         iomodel->FetchDataToInput(elements,FrictionEffectivePressureEnum);
     182                        if (FrictionCoupling==0){
     183                                iomodel->FetchDataToInput(elements,FrictionEffectivePressureEnum);
     184                        }
    181185                        break;
    182186                case 4:

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

    r19033 r19161  
    673673        int    materials_type,finiteelement,fe_FS;
    674674        int    approximation,frictionlaw;
     675        int    FrictionCoupling;
    675676        int*   finiteelement_list=NULL;
    676677        bool   isSSA,isL1L2,isHO,isFS,iscoupling;
    … …  
    800801                        break;
    801802                case 3:
     803                        iomodel->Constant(&FrictionCoupling,FrictionCouplingEnum);
    802804                        iomodel->FetchDataToInput(elements,FrictionCEnum);
    803805                        iomodel->FetchDataToInput(elements,FrictionAsEnum);
    804806                        iomodel->FetchDataToInput(elements,FrictionQEnum);
    805                         iomodel->FetchDataToInput(elements,FrictionEffectivePressureEnum);
     807                        if (FrictionCoupling==0){
     808                                iomodel->FetchDataToInput(elements,FrictionEffectivePressureEnum);
     809                        }
    806810                        break;
    807811                case 4:

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

    r19143 r19161  
    6262
    6363        int finiteelement = P1Enum;
     64       
    6465        if(iomodel->domaintype==Domain3DEnum) iomodel->FetchData(2,MeshVertexonbaseEnum,MeshVertexonsurfaceEnum);
    6566        ::CreateNodes(nodes,iomodel,ThermalAnalysisEnum,finiteelement);
    … …  
    7273
    7374        int frictionlaw;
    74 
     75        int FrictionCoupling;
    7576        /*Now, is the model 3d? otherwise, do nothing: */
    7677        if(iomodel->domaintype==Domain2DhorizontalEnum)return;
    … …  
    130131                        break;
    131132                case 3:
     133                        iomodel->Constant(&FrictionCoupling,FrictionCouplingEnum);
    132134                        iomodel->FetchDataToInput(elements,FrictionCEnum);
    133135                        iomodel->FetchDataToInput(elements,FrictionAsEnum);
    134136                        iomodel->FetchDataToInput(elements,FrictionQEnum);
    135                         iomodel->FetchDataToInput(elements,FrictionEffectivePressureEnum);
     137                        if (FrictionCoupling==0){
     138                                iomodel->FetchDataToInput(elements,FrictionEffectivePressureEnum);
     139                        }
    136140                        break;
    137141                case 4:
Note: See TracChangeset for help on using the changeset viewer.