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

BUG: fixing some free surface problems

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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){
Note: See TracChangeset for help on using the changeset viewer.