Changeset 16834


Ignore:
Timestamp:
11/19/13 15:40:54 (11 years ago)
Author:
Mathieu Morlighem
Message:

CHG: some simplifications

Location:
issm/trunk-jpl/src/c/classes/Elements
Files:
2 edited

Legend:

Unmodified
Added
Removed

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

    r16833 r16834  
    87488748ElementVector* Penta::CreatePVectorStressbalanceHOFront(void){
    87498749
    8750         /*Intermediaries */
    8751         IssmDouble  xyz_list[NUMVERTICES][3];
    8752         bool        isfront;
    8753 
    8754         isfront=IsZeroLevelset(MaskIceLevelsetEnum);
    8755 
    87568750        /*If no front, return NULL*/
    8757         if(!isfront) return NULL;
     8751        if(!IsZeroLevelset(MaskIceLevelsetEnum)) return NULL;
    87588752
    87598753        /*Intermediaries*/
    … …  
    87618755        IssmDouble  Jdet,surface,z_g,water_pressure,ice_pressure,air_pressure;
    87628756        IssmDouble  surface_under_water,base_under_water,pressure;
     8757        IssmDouble  xyz_list[NUMVERTICES][3];
    87638758        IssmDouble  xyz_list_front[4][3];
    87648759        IssmDouble  area_coordinates[4][3];
    … …  
    88428837ElementVector* Penta::CreatePVectorStressbalanceFSFront(void){
    88438838
    8844         /*Intermediaries */
     8839        /*If no front, return NULL*/
     8840        if(!IsZeroLevelset(MaskIceLevelsetEnum)) return NULL;
     8841
     8842        /*Fetch number of nodes and dof for this finite element*/
    88458843        int         i;
    8846         IssmDouble  xyz_list[NUMVERTICES][3];
    8847         bool        isfront;
    8848 
    8849         /*Retrieve all inputs and parameters*/
    8850         isfront=IsZeroLevelset(MaskIceLevelsetEnum);
    8851 
    8852         /*If no front, return NULL*/
    8853         if(!isfront) return NULL;
    8854 
    8855         /*Fetch number of nodes and dof for this finite element*/
    88568844        IssmDouble  rho_ice,rho_water,gravity;
    88578845        IssmDouble  Jdet,z_g,water_pressure,air_pressure;
    88588846        IssmDouble  surface_under_water,base_under_water,pressure;
     8847        IssmDouble  xyz_list[NUMVERTICES][3];
    88598848        IssmDouble  xyz_list_front[4][3];
    88608849        IssmDouble  area_coordinates[4][3];

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

    r16833 r16834  
    39223922ElementVector* Tria::CreatePVectorStressbalanceFSFront(void){
    39233923
    3924         /*Intermediaries */
     3924        /*If no front, return NULL*/
     3925        if(!IsZeroLevelset(MaskIceLevelsetEnum)) return NULL;
     3926
     3927        /*Intermediaries*/
    39253928        int         i;
    3926         IssmDouble  xyz_list[NUMVERTICES][3];
    3927         bool        isfront;
    3928 
    3929         /*If the level set is awlays <=0, there is no ice front here*/
    3930         isfront=IsZeroLevelset(MaskIceLevelsetEnum);
    3931 
    3932         /*If no front, return NULL*/
    3933         if(!isfront) return NULL;
    3934 
    3935         /*Intermediaries*/
    39363929        IssmDouble  rho_ice,rho_water,gravity,y_g;
    39373930        IssmDouble  Jdet,water_pressure,air_pressure,pressure;
     3931        IssmDouble  xyz_list[NUMVERTICES][3];
    39383932        IssmDouble  xyz_list_front[2][3];
    39393933        IssmDouble  area_coordinates[2][3];
    … …  
    41994193ElementVector* Tria::CreatePVectorStressbalanceSSAFront(){
    42004194
    4201         /*Intermediaries */
    4202         IssmDouble  xyz_list[NUMVERTICES][3];
    4203         bool        isfront;
    4204 
    4205         /*Determine if there is an icefront in this element*/
    4206         isfront=IsZeroLevelset(MaskIceLevelsetEnum);
    4207 
    42084195        /*If no front, return NULL*/
    4209         if(!isfront) return NULL;
     4196        if(!IsZeroLevelset(MaskIceLevelsetEnum)) return NULL;
    42104197
    42114198        /*Intermediaries*/
    … …  
    42134200        IssmDouble  Jdet,thickness,bed,water_pressure,ice_pressure,air_pressure;
    42144201        IssmDouble  surface_under_water,base_under_water,pressure;
     4202        IssmDouble  xyz_list[NUMVERTICES][3];
    42154203        IssmDouble  xyz_list_front[2][3];
    42164204        IssmDouble  area_coordinates[2][3];
Note: See TracChangeset for help on using the changeset viewer.