Changeset 22146


Ignore:
Timestamp:
10/06/17 11:47:24 (8 years ago)
Author:
seroussi
Message:

FIX: added missing checks on ice presence in coupling routines

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r22105 r22146  
    53665366ElementMatrix* StressbalanceAnalysis::CreateKMatrixCouplingHOFS(Element* element){/*{{{*/
    53675367
     5368        /* Check if ice in element */
     5369        if(!element->IsIceInElement()) return NULL;
     5370
    53685371        /*Constants*/
    53695372        int numnodes       = 3*6+1;
    … …  
    54435446}/*}}}*/
    54445447ElementMatrix* StressbalanceAnalysis::CreateKMatrixCouplingSSAFSFriction(Element* element){/*{{{*/
     5448
     5449        /* Check if ice in element */
     5450        if(!element->IsIceInElement()) return NULL;
    54455451
    54465452        /*Constants*/
    … …  
    55665572ElementMatrix* StressbalanceAnalysis::CreateKMatrixCouplingSSAFSViscous(Element* element){/*{{{*/
    55675573
     5574        /* Check if ice in element */
     5575        if(!element->IsIceInElement()) return NULL;
     5576
    55685577        /*Constants*/
    55695578        const int numdofm     = 2 *3;
    … …  
    56845693ElementMatrix* StressbalanceAnalysis::CreateKMatrixCouplingSSAHOFriction(Element* element){/*{{{*/
    56855694
     5695        /* Check if ice in element */
     5696        if(!element->IsIceInElement()) return NULL;
     5697
    56865698        if(element->IsFloating() || !element->IsOnBase()) return NULL;
    56875699
    … …  
    57605772}/*}}}*/
    57615773ElementMatrix* StressbalanceAnalysis::CreateKMatrixCouplingSSAHOViscous(Element* element){/*{{{*/
     5774
     5775        /* Check if ice in element */
     5776        if(!element->IsIceInElement()) return NULL;
    57625777
    57635778        /*Constants*/
    … …  
    59145929ElementMatrix* StressbalanceAnalysis::CreateKMatrixSSA3dFriction(Element* element){/*{{{*/
    59155930
     5931        /* Check if ice in element */
     5932        if(!element->IsIceInElement()) return NULL;
     5933
    59165934        /*Initialize Element matrix and return if necessary*/
    59175935        if(element->IsFloating() || !element->IsOnBase()) return NULL;
    … …  
    59285946}/*}}}*/
    59295947ElementMatrix* StressbalanceAnalysis::CreateKMatrixSSA3dViscous(Element* element){/*{{{*/
     5948
     5949        /* Check if ice in element */
     5950        if(!element->IsIceInElement()) return NULL;
    59305951
    59315952        /*Constants*/
    … …  
    60216042ElementVector* StressbalanceAnalysis::CreatePVectorCouplingHOFSFriction(Element* element){/*{{{*/
    60226043
     6044        /* Check if ice in element */
     6045        if(!element->IsIceInElement()) return NULL;
     6046
    60236047        /*Intermediaries*/
    60246048        int         i,approximation;
    … …  
    61046128ElementVector* StressbalanceAnalysis::CreatePVectorCouplingHOFSViscous(Element* element){/*{{{*/
    61056129
     6130        /* Check if ice in element */
     6131        if(!element->IsIceInElement()) return NULL;
     6132
    61066133        /*Intermediaries */
    61076134        int         i,approximation;
    … …  
    61846211}/*}}}*/
    61856212ElementVector* StressbalanceAnalysis::CreatePVectorCouplingSSAFSFriction(Element* element){/*{{{*/
     6213
     6214        /* Check if ice in element */
     6215        if(!element->IsIceInElement()) return NULL;
    61866216
    61876217        /*Intermediaries*/
    … …  
    62646294}/*}}}*/
    62656295ElementVector* StressbalanceAnalysis::CreatePVectorCouplingSSAFSViscous(Element* element){/*{{{*/
     6296
     6297        /* Check if ice in element */
     6298        if(!element->IsIceInElement()) return NULL;
    62666299
    62676300        /*Intermediaries */
Note: See TracChangeset for help on using the changeset viewer.