Changeset 4461


Ignore:
Timestamp:
07/08/10 11:20:58 (15 years ago)
Author:
seroussi
Message:

problem with element_type in friction

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/objects/Loads/Friction.cpp

    r4413 r4461  
    9595        if (Neff<0)Neff=0;
    9696
    97         if(strcmp(element_type,"2d")){
     97        if(strcmp(element_type,"2d")==0){
    9898                inputs->GetParameterValue(&vx, gauss,vxenum);
    9999                inputs->GetParameterValue(&vy, gauss,vyenum);
    100100                vmag=sqrt(pow(vx,2)+pow(vy,2));
    101101        }
    102         else{
     102        else if (strcmp(element_type,"3d")==0){
    103103                inputs->GetParameterValue(&vx, gauss,vxenum);
    104104                inputs->GetParameterValue(&vy, gauss,vyenum);
     
    106106                vmag=sqrt(pow(vx,2)+pow(vy,2)+pow(vz,2));
    107107        }
     108        else ISSMERROR("element_type %s not supported yet",element_type);
    108109       
    109110        alpha2=pow(drag_coefficient,2)*pow(Neff,r)*pow(vmag,(s-1));
Note: See TracChangeset for help on using the changeset viewer.