Changeset 9600


Ignore:
Timestamp:
09/02/11 17:33:03 (14 years ago)
Author:
Mathieu Morlighem
Message:

Final bug fix: meanvel was not properly initialized

File:
1 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/c/objects/Elements/Tria.cpp

    r9597 r9600  
    15121512        double     dux,duy;
    15131513        double     epsvel=2.220446049250313e-16;
    1514         double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
     1514        double     meanvel=3.170979198376458e-05; /*1000 m/yr*/
    15151515        double     scalex=0,scaley=0,scale=0,S=0;
    15161516        double     vx,vy,vxobs,vyobs,weight;
    … …  
    16141614                                                velocity_mag    =sqrt(pow(vx,   2.)+pow(vy,   2.))+epsvel;
    16151615                                                obs_velocity_mag=sqrt(pow(vxobs,2.)+pow(vyobs,2.))+epsvel;
     1616                                                printf("-------------- file: Tria.cpp line: %i\n",__LINE__);
     1617                                                printf("%g %g\n",velocity_mag,obs_velocity_mag);
    16161618                                                scale=-8*pow(meanvel,2.)/pow(velocity_mag,2.)*log(velocity_mag/obs_velocity_mag);
     1619                                                printf("%g\n",meanvel);
    16171620                                                dux=scale*vx;
    16181621                                                duy=scale*vy;
     1622                                                printf("%g %g %g\n",vx,vy,scale);
    16191623                                                pe->values[i*NDOF2+0]+=dux*weight*Jdet*gauss->weight*basis[i];
    16201624                                                pe->values[i*NDOF2+1]+=duy*weight*Jdet*gauss->weight*basis[i];
    … …  
    16881692        double     dux,duy;
    16891693        double     epsvel=2.220446049250313e-16;
    1690         double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
     1694        double     meanvel=3.170979198376458e-05; /*1000 m/yr*/
    16911695        double     scalex=0,scaley=0,scale=0,S=0;
    16921696        double     vx,vy,vxobs,vyobs,weight;
    … …  
    48584862        double     misfit,Jdet;
    48594863        double     epsvel=2.220446049250313e-16;
    4860         double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
     4864        double     meanvel=3.170979198376458e-05; /*1000 m/yr*/
    48614865        double     velocity_mag,obs_velocity_mag;
    48624866        double     xyz_list[NUMVERTICES][3];
    … …  
    49234927        double     Jelem=0, S=0;
    49244928        double     epsvel=2.220446049250313e-16;
    4925         double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
     4929        double     meanvel=3.170979198376458e-05; /*1000 m/yr*/
    49264930        double     misfit, Jdet;
    49274931        double     vx,vy,vxobs,vyobs,weight;
    … …  
    50135017        double     misfit,Jdet;
    50145018        double     epsvel=2.220446049250313e-16;
    5015         double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
     5019        double     meanvel=3.170979198376458e-05; /*1000 m/yr*/
    50165020        double     vx,vy,vxobs,vyobs,weight;
    50175021        double     xyz_list[NUMVERTICES][3];
Note: See TracChangeset for help on using the changeset viewer.