Ignore:
Timestamp:
07/28/10 15:07:43 (15 years ago)
Author:
seroussi
Message:

fixed problem with pressure steadystate

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/objects/Elements/Tria.cpp

    r4840 r4846  
    57925792
    57935793        int          dummy;
    5794         Input*       VzInput=NULL;
    5795         double*      VzPtr=NULL;
     5794        Input*       vz_input=NULL;
     5795        double*      vz_ptr=NULL;
    57965796       
    57975797        /*Get dof list: */
     
    58105810
    58115811        /*Get Vz*/
    5812         VzInput=inputs->GetInput(VzEnum);
    5813         if (VzInput){
    5814                 if (VzInput->Enum()!=TriaVertexInputEnum){
    5815                         ISSMERROR("Cannot compute Vel as Vz is of type %s",EnumAsString(VzInput->Enum()));
     5812        vz_input=inputs->GetInput(VzEnum);
     5813        if (vz_input){
     5814                if (vz_input->Enum()!=TriaVertexInputEnum){
     5815                        ISSMERROR("Cannot compute Vel as Vz is of type %s",EnumAsString(vz_input->Enum()));
    58165816                }
    5817                 VzInput->GetValuesPtr(&VzPtr,&dummy);
    5818                 for(i=0;i<numvertices;i++) vz[i]=VzPtr[i];
     5817                vz_input->GetValuesPtr(&vz_ptr,&dummy);
     5818                for(i=0;i<numvertices;i++) vz[i]=vz_ptr[i];
    58195819        }
    58205820        else{
     
    58705870
    58715871        int          dummy;
    5872         Input*       VzInput=NULL;
    5873         double*      VzPtr=NULL;
     5872        Input*       vz_input=NULL;
     5873        double*      vz_ptr=NULL;
    58745874       
    58755875        /*Get dof list: */
     
    58885888
    58895889        /*Get Vz*/
    5890         VzInput=inputs->GetInput(VzEnum);
    5891         if (VzInput){
    5892                 if (VzInput->Enum()!=TriaVertexInputEnum){
    5893                         ISSMERROR("Cannot compute Vel as Vz is of type %s",EnumAsString(VzInput->Enum()));
     5890        vz_input=inputs->GetInput(VzEnum);
     5891        if (vz_input){
     5892                if (vz_input->Enum()!=TriaVertexInputEnum){
     5893                        ISSMERROR("Cannot compute Vel as Vz is of type %s",EnumAsString(vz_input->Enum()));
    58945894                }
    5895                 VzInput->GetValuesPtr(&VzPtr,&dummy);
    5896                 for(i=0;i<numvertices;i++) vz[i]=VzPtr[i];
     5895                vz_input->GetValuesPtr(&vz_ptr,&dummy);
     5896                for(i=0;i<numvertices;i++) vz[i]=vz_ptr[i];
    58975897        }
    58985898        else{
Note: See TracChangeset for help on using the changeset viewer.