Changeset 13104


Ignore:
Timestamp:
08/21/12 10:59:17 (13 years ago)
Author:
cborstad
Message:

BUG: re-added depth averaging of Z

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r13101 r13104  
    15551555        /*Go through all the input objects, and find the one corresponding to enum_type, if it exists: */
    15561556        if (enum_type==MaterialsRheologyBbarEnum) input=this->matice->inputs->GetInput(MaterialsRheologyBEnum);
     1557        else if (enum_type==MaterialsRheologyZbarEnum) input=this->matice->inputs->GetInput(MaterialsRheologyZEnum);
    15571558        else input=this->inputs->GetInput(enum_type);
    15581559        //if (!input) _error_("Input " << EnumToStringx(enum_type) << " not found in penta->inputs"); why error out? if the requested input does not exist, we should still
    … …  
    16671668                                        break;
    16681669                                case MaterialsRheologyBbarEnum:
     1670                                case MaterialsRheologyZbarEnum:
    16691671                                        /*Matice will take care of it*/ break;
    16701672                                default:
    … …  
    32303232                        *presponse=this->matice->GetBbar();
    32313233                        break;
     3234                case MaterialsRheologyZbarEnum:
     3235                        *presponse=this->matice->GetZbar();
     3236                        break;
    32323237                case VelEnum:
    32333238                        {
    … …  
    43474352                                for(i=0;i<numdof;i++) B[i]=B_average;
    43484353                                this->matice->inputs->AddInput(new PentaP1Input(MaterialsRheologyBEnum,B));
     4354                                this->matice->inputs->AddInput(new PentaP1Input(MaterialsRheologyZEnum,B));
    43494355                                break;
    43504356                        case ArrheniusEnum:
    … …  
    43554361                                for(i=0;i<numdof;i++) B[i]=B_average;
    43564362                                this->matice->inputs->AddInput(new PentaP1Input(MaterialsRheologyBEnum,B));
     4363                                this->matice->inputs->AddInput(new PentaP1Input(MaterialsRheologyZEnum,B));
    43574364                                break;
    43584365                        default:
    … …  
    44254432                                for(i=0;i<numdof;i++) B[i]=B_average;
    44264433                                this->matice->inputs->AddInput(new PentaP1Input(MaterialsRheologyBEnum,B));
     4434                                this->matice->inputs->AddInput(new PentaP1Input(MaterialsRheologyZEnum,B));
    44274435                                break;
    44284436                        case ArrheniusEnum:
    … …  
    44334441                                for(i=0;i<numdof;i++) B[i]=B_average;
    44344442                                this->matice->inputs->AddInput(new PentaP1Input(MaterialsRheologyBEnum,B));
     4443                                this->matice->inputs->AddInput(new PentaP1Input(MaterialsRheologyZEnum,B));
    44354444                                break;
    44364445                        default:
    … …  
    45404549        this->InputDepthAverageAtBase(MaterialsRheologyBEnum,MaterialsRheologyBbarEnum,MaterialsEnum);
    45414550
     4551        /*Depth Averaging Z*/
     4552        this->InputDepthAverageAtBase(MaterialsRheologyZEnum,MaterialsRheologyZbarEnum,MaterialsEnum);
     4553
    45424554        /*Call Tria function*/
    45434555        Tria* tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria.
    … …  
    45474559        /*Delete B averaged*/
    45484560        this->matice->inputs->DeleteInput(MaterialsRheologyBbarEnum);
     4561
     4562        /*Delete Z averaged*/
     4563        this->matice->inputs->DeleteInput(MaterialsRheologyZbarEnum);
    45494564
    45504565        /*clean up and return*/
    … …  
    63046319        this->InputDepthAverageAtBase(MaterialsRheologyBEnum,MaterialsRheologyBbarEnum,MaterialsEnum);
    63056320
     6321        /*Depth Averaging Z*/
     6322        this->InputDepthAverageAtBase(MaterialsRheologyZEnum,MaterialsRheologyZbarEnum,MaterialsEnum);
     6323
    63066324        /*Call Tria function*/
    63076325        Tria* tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria.
    … …  
    63116329        /*Delete B averaged*/
    63126330        this->matice->inputs->DeleteInput(MaterialsRheologyBbarEnum);
     6331
     6332        /*Delete Z averaged*/
     6333        this->matice->inputs->DeleteInput(MaterialsRheologyZbarEnum);
    63136334
    63146335        /*clean up and return*/
    … …  
    78607881        this->InputDepthAverageAtBase(MaterialsRheologyBEnum,MaterialsRheologyBbarEnum,MaterialsEnum);
    78617882
     7883        /*Depth Averaging Z*/
     7884        this->InputDepthAverageAtBase(MaterialsRheologyZEnum,MaterialsRheologyZbarEnum,MaterialsEnum);
     7885
    78627886        /*Call Tria function*/
    78637887        Tria* tria=(Tria*)SpawnTria(0,1,2); //nodes 0, 1 and 2 make the new tria.
    … …  
    78677891        /*Delete B averaged*/
    78687892        this->matice->inputs->DeleteInput(MaterialsRheologyBbarEnum);
     7893
     7894        /*Delete Z averaged*/
     7895        this->matice->inputs->DeleteInput(MaterialsRheologyZbarEnum);
    78697896
    78707897        /*clean up and return*/
Note: See TracChangeset for help on using the changeset viewer.