Changeset 17177


Ignore:
Timestamp:
01/28/14 07:44:35 (11 years ago)
Author:
Mathieu Morlighem
Message:

BUG: fixing some free surface problems

Location:
issm/trunk-jpl/src/c
Files:
5 edited

Legend:

Unmodified
Added
Removed

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

    r17172 r17177  
    145145        Input* vx_input=basalelement->GetInput(VxEnum); _assert_(vx_input);
    146146        Input* vy_input=NULL;
    147         if(dim>1) basalelement->GetInput(VyEnum); _assert_(vy_input);
     147        if(dim>1){basalelement->GetInput(VyEnum); _assert_(vy_input);}
    148148        h = basalelement->CharacteristicLength();
    149149
    … …  
    165165                                        &Ke->values[0],1);
    166166
    167                 GetB(B,element,dim,xyz_list,gauss);
    168                 GetBprime(Bprime,element,dim,xyz_list,gauss);
     167                GetB(B,basalelement,dim,xyz_list,gauss);
     168                GetBprime(Bprime,basalelement,dim,xyz_list,gauss);
    169169
    170170                D_scalar=dt*gauss->weight*Jdet;

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

    r17172 r17177  
    174174                                        &Ke->values[0],1);
    175175
    176                 GetB(B,element,dim,xyz_list,gauss);
    177                 GetBprime(Bprime,element,dim,xyz_list,gauss);
     176                GetB(B,topelement,dim,xyz_list,gauss);
     177                GetBprime(Bprime,topelement,dim,xyz_list,gauss);
    178178
    179179                D_scalar=dt*gauss->weight*Jdet;

  • issm/trunk-jpl/src/c/classes/Elements/Seg.cpp

    r16993 r17177  
    5151/*}}}*/
    5252
     53/*FUNCTION Seg::CharacteristicLength{{{*/
     54IssmDouble Seg::CharacteristicLength(void){
     55
     56        IssmDouble xyz_list[NUMVERTICES][3];
     57        IssmDouble x1,y1,x2,y2;
     58
     59        /*Get xyz list: */
     60        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     61        x1=xyz_list[0][0]; y1=xyz_list[0][1];
     62        x2=xyz_list[1][0]; y2=xyz_list[1][1];
     63
     64        return sqrt((x2-x1)*(x2-x1) + (y2-y1)*(y2-y1));
     65}
     66/*}}}*/
    5367/*FUNCTION Seg::Echo{{{*/
    5468void Seg::Echo(void){
    … …  
    124138/*}}}*/
    125139
    126 /*FUNCTION Seg::GetSize{{{*/
    127 IssmDouble Seg::GetSize(void){
    128 
    129         IssmDouble xyz_list[NUMVERTICES][3];
    130         IssmDouble x1,y1,x2,y2;
    131 
    132         /*Get xyz list: */
    133         ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    134         x1=xyz_list[0][0]; y1=xyz_list[0][1];
    135         x2=xyz_list[1][0]; y2=xyz_list[1][1];
    136 
    137         return sqrt((x2-x1)*(x2-x1) + (y2-y1)*(y2-y1));
    138 }
    139 /*}}}*/
    140140/*FUNCTION Seg::GetNumberOfNodes;{{{*/
    141141int Seg::GetNumberOfNodes(void){

  • issm/trunk-jpl/src/c/classes/Elements/Seg.h

    r17172 r17177  
    6262                void        AddInput(int input_enum, IssmDouble* values, int interpolation_enum){_error_("not implemented yet");};
    6363                void        AddMaterialInput(int input_enum, IssmDouble* values, int interpolation_enum){_error_("not implemented yet");};
    64                 IssmDouble  CharacteristicLength(void){_error_("not implemented yet");};
     64                IssmDouble  CharacteristicLength(void);
    6565                void        ComputeBasalStress(Vector<IssmDouble>* sigma_b){_error_("not implemented yet");};
    6666                void        ComputeStrainRate(Vector<IssmDouble>* eps){_error_("not implemented yet");};
    … …  
    256256#endif
    257257                /*}}}*/
    258                 /*Seg specific routines:*/
    259                 IssmDouble     GetSize(void);
    260258};
    261259#endif  /* _SEG_H */

  • issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

    r17176 r17177  
    185185}
    186186/*}}}*/
    187 /*FUNCTION Tria::AddInput{{{*/
     187/*FUNCTION Tria::CharacteristicLength{{{*/
    188188IssmDouble Tria::CharacteristicLength(void){
    189189
Note: See TracChangeset for help on using the changeset viewer.