Changeset 19731


Ignore:
Timestamp:
11/17/15 11:07:24 (9 years ago)
Author:
Mathieu Morlighem
Message:

CHG: changing the sign of conductivity

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/analyses/HydrologySommersAnalysis.cpp

    r19725 r19731  
    118118
    119119        /* Start  looping on the number of gaussian points: */
    120         Gauss* gauss=element->NewGauss(2);
     120        Gauss* gauss=element->NewGauss(1);
    121121        for(int ig=gauss->begin();ig<gauss->end();ig++){
    122122                gauss->GaussPoint(ig);
     
    127127                for(int i=0;i<numnodes;i++){
    128128                        for(int j=0;j<numnodes;j++){
    129                                 Ke->values[i*numnodes+j] += -conductivity*gauss->weight*Jdet*(dbasis[0*numnodes+i]*dbasis[0*numnodes+j] + dbasis[1*numnodes+i]*dbasis[1*numnodes+j]);
     129                                Ke->values[i*numnodes+j] += conductivity*gauss->weight*Jdet*(dbasis[0*numnodes+i]*dbasis[0*numnodes+j] + dbasis[1*numnodes+i]*dbasis[1*numnodes+j]);
    130130                        }
    131131                }
Note: See TracChangeset for help on using the changeset viewer.