Ignore:
Timestamp:
08/21/12 16:23:21 (13 years ago)
Author:
cborstad
Message:

BUG: debugging 3D NR on branch, rheology_Z was getting missed

File:
1 edited

Legend:

Unmodified
Added
Removed

  • TabularUnified issm/branches/trunk-jpl-damage/src/c/classes/objects/Materials/Matice.cpp

    r13101 r13117  
    314314        /*Intermediaries: */
    315315        IssmDouble A,e;
    316         IssmDouble B,n;
    317 
    318         /*Get B and n*/
    319         B=GetB();
     316        IssmDouble B,n,Z;
     317
     318        /*Get B, Z and n*/
    320319        n=GetN();
     320        Z=GetZ();
     321        B=Z*GetB();
    321322
    322323        if (n==1){
    … …  
    384385        /*Intermediaries: */
    385386        IssmDouble A,e;
    386         IssmDouble B,n;
     387        IssmDouble B,n,Z;
    387388        IssmDouble eps0;
    388389
    389390        /*Get B and n*/
    390391        eps0=pow((IssmDouble)10,(IssmDouble)-27);
    391         B=GetB();
    392392        n=GetN();
     393        Z=GetZ();
     394        B=Z*GetB();
    393395       
    394396        if (n==1){
Note: See TracChangeset for help on using the changeset viewer.