Changeset 16019


Ignore:
Timestamp:
08/29/13 15:00:14 (12 years ago)
Author:
Mathieu Morlighem
Message:

CHG: fixed balanced velocities now working

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r16018 r16019  
    55875587        /*Intermediaries */
    55885588        IssmDouble diffusivity;
    5589         IssmDouble Jdettria,D_scalar,dt,h;
     5589        IssmDouble Jdet,D_scalar,dt,h;
    55905590        IssmDouble vx,vy,vel,dvxdx,dvydy;
    55915591        IssmDouble dvx[2],dvy[2];
    … …  
    56225622                gauss->GaussPoint(ig);
    56235623
    5624                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     5624                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    56255625                GetNodalFunctions(basis,gauss);
    56265626
    … …  
    56305630                vy_input->GetInputDerivativeValue(&dvy[0],&xyz_list[0][0],gauss);
    56315631
    5632                 D_scalar=gauss->weight*Jdettria;
     5632                D_scalar=gauss->weight*Jdet;
    56335633
    56345634                TripleMultiply(basis,1,numnodes,1,
    … …  
    56425642                dvxdx=dvx[0];
    56435643                dvydy=dvy[1];
    5644                 D_scalar=dt*gauss->weight*Jdettria;
     5644                D_scalar=dt*gauss->weight*Jdet;
    56455645
    56465646                D[0][0]=D_scalar*dvxdx;
    … …  
    58125812        /*Intermediaries */
    58135813        int        i;
    5814         IssmDouble Jdettria,dt;
     5814        IssmDouble Jdet,dt;
    58155815        IssmDouble basal_melting_g;
    58165816        IssmDouble old_watercolumn_g;
    … …  
    58405840                gauss->GaussPoint(ig);
    58415841
    5842                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     5842                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    58435843                GetNodalFunctions(basis, gauss);
    58445844
    … …  
    58475847
    58485848                if(reCast<int,IssmDouble>(dt)){
    5849                         for(i=0;i<numnodes;i++) pe->values[i]+=Jdettria*gauss->weight*(old_watercolumn_g+dt*basal_melting_g)*basis[i];
     5849                        for(i=0;i<numnodes;i++) pe->values[i]+=Jdet*gauss->weight*(old_watercolumn_g+dt*basal_melting_g)*basis[i];
    58505850                }
    58515851                else{
    5852                         for(i=0;i<numnodes;i++) pe->values[i]+=Jdettria*gauss->weight*basal_melting_g*basis[i];
     5852                        for(i=0;i<numnodes;i++) pe->values[i]+=Jdet*gauss->weight*basal_melting_g*basis[i];
    58535853                }
    58545854        }
    … …  
    63116311        int        stabilization;
    63126312        int        dim;
    6313         IssmDouble Jdettria,D_scalar,dt,h;
     6313        IssmDouble Jdet,D_scalar,dt,h;
    63146314        IssmDouble vel,vx,vy,dvxdx,dvydy;
    63156315        IssmDouble dvx[2],dvy[2];
    … …  
    63496349                gauss->GaussPoint(ig);
    63506350
    6351                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     6351                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    63526352                GetNodalFunctions(basis,gauss);
    63536353
    … …  
    63576357                vyaverage_input->GetInputDerivativeValue(&dvy[0],&xyz_list[0][0],gauss);
    63586358
    6359                 D_scalar=gauss->weight*Jdettria;
     6359                D_scalar=gauss->weight*Jdet;
    63606360
    63616361                TripleMultiply(basis,1,numnodes,1,
    … …  
    63696369                dvxdx=dvx[0];
    63706370                dvydy=dvy[1];
    6371                 D_scalar=dt*gauss->weight*Jdettria;
     6371                D_scalar=dt*gauss->weight*Jdet;
    63726372
    63736373                D[0][0]=D_scalar*dvxdx;
    … …  
    64306430        int        dim;
    64316431        IssmDouble xyz_list[NUMVERTICES][3];
    6432         IssmDouble Jdettria,D_scalar,dt,vx,vy;
     6432        IssmDouble Jdet,D_scalar,dt,vx,vy;
    64336433
    64346434        /*Fetch number of nodes for this finite element*/
    … …  
    64666466                vyaverage_input->GetInputValue(&vy,gauss);
    64676467
    6468                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     6468                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    64696469                GetNodalFunctions(basis,gauss);
    64706470
    6471                 D_scalar=gauss->weight*Jdettria;
     6471                D_scalar=gauss->weight*Jdet;
    64726472
    64736473                TripleMultiply(basis,1,numnodes,1,
    … …  
    64806480                GetBprimeMasstransport(B, &xyz_list[0][0], gauss);
    64816481
    6482                 D_scalar=-dt*gauss->weight*Jdettria;
     6482                D_scalar=-dt*gauss->weight*Jdet;
    64836483                D[0][0]=D_scalar*vx;
    64846484                D[0][1]=0.;
    … …  
    65066506        int        stabilization;
    65076507        int        dim;
    6508         IssmDouble Jdettria,D_scalar,dt,h;
     6508        IssmDouble Jdet,D_scalar,dt,h;
    65096509        IssmDouble vel,vx,vy;
    65106510        IssmDouble xyz_list[NUMVERTICES][3];
    … …  
    65346534                gauss->GaussPoint(ig);
    65356535
    6536                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     6536                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    65376537                GetNodalFunctions(basis,gauss);
    65386538
    … …  
    65406540                vy_input->GetInputValue(&vy,gauss);
    65416541
    6542                 D_scalar=gauss->weight*Jdettria;
     6542                D_scalar=gauss->weight*Jdet;
    65436543
    65446544                TripleMultiply(basis,1,numnodes,1,
    … …  
    65506550                GetBprimeMasstransport(Bprime,&xyz_list[0][0],gauss);
    65516551
    6552                 D_scalar=dt*gauss->weight*Jdettria;
     6552                D_scalar=dt*gauss->weight*Jdet;
    65536553
    65546554                D[0][0]=D_scalar*vx;
    … …  
    66046604        int        stabilization;
    66056605        int        dim;
    6606         IssmDouble Jdettria,D_scalar,dt,h;
     6606        IssmDouble Jdet,D_scalar,dt,h;
    66076607        IssmDouble vel,vx,vy;
    66086608        IssmDouble xyz_list[NUMVERTICES][3];
    … …  
    66326632                gauss->GaussPoint(ig);
    66336633
    6634                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     6634                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    66356635                GetNodalFunctions(basis,gauss);
    66366636
    … …  
    66386638                vy_input->GetInputValue(&vy,gauss);
    66396639
    6640                 D_scalar=gauss->weight*Jdettria;
     6640                D_scalar=gauss->weight*Jdet;
    66416641
    66426642                TripleMultiply(basis,1,numnodes,1,
    … …  
    66486648                GetBprimeMasstransport(Bprime,&xyz_list[0][0],gauss);
    66496649
    6650                 D_scalar=dt*gauss->weight*Jdettria;
     6650                D_scalar=dt*gauss->weight*Jdet;
    66516651
    66526652                D[0][0]=D_scalar*vx;
    … …  
    67136713
    67146714        /*Intermediaries */
    6715         IssmDouble Jdettria,dt;
     6715        IssmDouble Jdet,dt;
    67166716        IssmDouble ms,mb,mb_correction,thickness;
    67176717        IssmDouble xyz_list[NUMVERTICES][3];
    … …  
    67396739                gauss->GaussPoint(ig);
    67406740
    6741                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     6741                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    67426742                GetNodalFunctions(basis,gauss);
    67436743
    … …  
    67506750                 mb_correction=0.;
    67516751
    6752                 for(int i=0;i<numnodes;i++) pe->values[i]+=Jdettria*gauss->weight*(thickness+dt*(ms-mb-mb_correction))*basis[i];
     6752                for(int i=0;i<numnodes;i++) pe->values[i]+=Jdet*gauss->weight*(thickness+dt*(ms-mb-mb_correction))*basis[i];
    67536753        }
    67546754
    … …  
    67636763
    67646764        /*Intermediaries */
    6765         IssmDouble Jdettria,dt;
     6765        IssmDouble Jdet,dt;
    67666766        IssmDouble ms,mb,thickness;
    67676767        IssmDouble xyz_list[NUMVERTICES][3];
    … …  
    67876787                gauss->GaussPoint(ig);
    67886788
    6789                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     6789                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    67906790                GetNodalFunctions(basis,gauss);
    67916791
    … …  
    67946794                thickness_input->GetInputValue(&thickness,gauss);
    67956795
    6796                 for(int i=0;i<numnodes;i++) pe->values[i]+=Jdettria*gauss->weight*(thickness+dt*(ms-mb))*basis[i];
     6796                for(int i=0;i<numnodes;i++) pe->values[i]+=Jdet*gauss->weight*(thickness+dt*(ms-mb))*basis[i];
    67976797        }
    67986798
    … …  
    68076807
    68086808        /*Intermediaries */
    6809         IssmDouble Jdettria,dt;
     6809        IssmDouble Jdet,dt;
    68106810        IssmDouble ms,surface,vz;
    68116811        IssmDouble xyz_list[NUMVERTICES][3];
    … …  
    68326832                gauss->GaussPoint(ig);
    68336833
    6834                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     6834                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    68356835                GetNodalFunctions(basis,gauss);
    68366836
    … …  
    68396839                surface_input->GetInputValue(&surface,gauss);
    68406840
    6841                 for(int i=0;i<numnodes;i++) pe->values[i]+=Jdettria*gauss->weight*(surface + dt*ms + dt*vz)*basis[i];
     6841                for(int i=0;i<numnodes;i++) pe->values[i]+=Jdet*gauss->weight*(surface + dt*ms + dt*vz)*basis[i];
    68426842        }
    68436843
    … …  
    68526852
    68536853        /*Intermediaries */
    6854         IssmDouble Jdettria,dt;
     6854        IssmDouble Jdet,dt;
    68556855        IssmDouble mb,mb_correction,bed,vz;
    68566856        IssmDouble xyz_list[NUMVERTICES][3];
    … …  
    68786878                gauss->GaussPoint(ig);
    68796879
    6880                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     6880                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    68816881                GetNodalFunctions(basis,gauss);
    68826882
    … …  
    68896889                 mb_correction=0.;
    68906890
    6891                 for(int i=0;i<numnodes;i++) pe->values[i]+=Jdettria*gauss->weight*(bed+dt*(mb-mb_correction) + dt*vz)*basis[i];
     6891                for(int i=0;i<numnodes;i++) pe->values[i]+=Jdet*gauss->weight*(bed+dt*(mb-mb_correction) + dt*vz)*basis[i];
    68926892        }
    68936893
    … …  
    70627062        int        stabilization;
    70637063        int        i,j,dim;
    7064         IssmDouble Jdettria,vx,vy,dvxdx,dvydy,vel,h;
     7064        IssmDouble Jdet,vx,vy,dvxdx,dvydy,vel,h;
    70657065        IssmDouble D_scalar;
    70667066        IssmDouble dvx[2],dvy[2];
    … …  
    70997099                gauss->GaussPoint(ig);
    71007100
    7101                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     7101                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    71027102                GetBMasstransport(B,&xyz_list[0][0],gauss);
    71037103                GetBprimeMasstransport(Bprime,&xyz_list[0][0],gauss);
    … …  
    71107110                dvxdx=dvx[0];
    71117111                dvydy=dvy[1];
    7112                 D_scalar=gauss->weight*Jdettria;
     7112                D_scalar=gauss->weight*Jdet;
    71137113
    71147114                D[0][0]=D_scalar*dvxdx;
    … …  
    71707170        /*Intermediaries*/
    71717171        int        dim;
    7172         IssmDouble vx,vy,D_scalar,Jdettria;
     7172        IssmDouble vx,vy,D_scalar,Jdet;
    71737173        IssmDouble xyz_list[NUMVERTICES][3];
    71747174
    … …  
    71957195                gauss->GaussPoint(ig);
    71967196
    7197                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     7197                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    71987198                /*WARNING: B and Bprime are inverted compared to usual masstransport!!!!*/
    71997199                GetBMasstransport(Bprime,&xyz_list[0][0],gauss);
    … …  
    72037203                vy_input->GetInputValue(&vy,gauss);
    72047204
    7205                 D_scalar=-gauss->weight*Jdettria;
     7205                D_scalar=-gauss->weight*Jdet;
    72067206                D[0][0]=D_scalar*vx;
    72077207                D[0][1]=0.;
    … …  
    72287228        /*Intermediaries */
    72297229        IssmDouble xyz_list[NUMVERTICES][3];
    7230         IssmDouble dhdt_g,mb_g,ms_g,Jdettria;
     7230        IssmDouble dhdt_g,mb_g,ms_g,Jdet;
    72317231        IssmDouble h,gamma,thickness;
    72327232        IssmDouble hnx,hny,dhnx[2],dhny[2];
    72337233        IssmDouble vel,vx,vy,dvxdx,dvydy;
    72347234        IssmDouble D[2][2];
    7235         int        stabilization=1;
     7235        int        stabilization=0;
    72367236
    72377237        /*Fetch number of nodes and dof for this finite element*/
    … …  
    72867286                gamma=h/(2.*thickness+1.e-10);
    72877287
    7288                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     7288                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    72897289                GetNodalFunctions(basis,gauss);
    72907290                GetNodalFunctionsDerivatives(dbasis,&xyz_list[0][0],gauss);
    … …  
    72927292                for(int i=0;i<numnodes;i++){
    72937293                        for(int j=0;j<numnodes;j++){
    7294                                 Ke->values[i*numnodes+j] += gauss->weight*Jdettria*(
     7294                                Ke->values[i*numnodes+j] += gauss->weight*Jdet*(
    72957295                                                        (basis[i]+gamma*(basis[i]*(dhnx[0]+dhny[1]) + dbasis[0*numnodes+i]*hnx + dbasis[1*numnodes+i]*hny))*
    72967296                                                        (basis[j]*(dhnx[0]+dhny[1])  + dbasis[0*numnodes+j]*hnx + dbasis[1*numnodes+j]*hny)
    … …  
    73187318                if(stabilization==1 || stabilization==2){
    73197319                        GetBprimeMasstransport(B,&xyz_list[0][0],gauss);
    7320                         D[0][0]=gauss->weight*Jdettria*D[0][0];
    7321                         D[1][0]=gauss->weight*Jdettria*D[1][0];
    7322                         D[0][1]=gauss->weight*Jdettria*D[0][1];
    7323                         D[1][1]=gauss->weight*Jdettria*D[1][1];
     7320                        D[0][0]=gauss->weight*Jdet*D[0][0];
     7321                        D[1][0]=gauss->weight*Jdet*D[1][0];
     7322                        D[0][1]=gauss->weight*Jdet*D[0][1];
     7323                        D[1][1]=gauss->weight*Jdet*D[1][1];
    73247324                        TripleMultiply(B,2,numnodes,1,
    73257325                                                &D[0][0],2,2,0,
    … …  
    73497349        IssmDouble xyz_list[NUMVERTICES][3];
    73507350        IssmDouble D[2][2];
    7351         IssmDouble l=15.;
     7351        IssmDouble l=12.;
    73527352
    73537353        /*Fetch number of nodes and dof for this finite element*/
    … …  
    73557355
    73567356        /*Initialize Element matrix and vectors*/
    7357         ElementMatrix* Ke    = CreateMassMatrix();
    7358         IssmDouble*    B     = xNew<IssmDouble>(2*numnodes);
     7357        ElementMatrix* Ke     = new ElementMatrix(nodes,numnodes,this->parameters);
     7358        IssmDouble*    dbasis = xNew<IssmDouble>(2*numnodes);
     7359        IssmDouble*    basis  = xNew<IssmDouble>(numnodes);
    73597360
    73607361        /*Retrieve all inputs and parameters*/
    … …  
    73687369                gauss->GaussPoint(ig);
    73697370
    7370                 GetBMasstransport(B,&xyz_list[0][0], gauss);
    73717371                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    73727372                thickness_input->GetInputValue(&thickness,gauss);
    73737373                if(thickness<50.) thickness=50.;
    73747374
    7375                 D_scalar=gauss->weight*Jdet*(l*thickness)*(l*thickness);
    7376                 D[0][0]=D_scalar;
    7377            D[1][0]=0.;
    7378                 D[0][1]=0.;
    7379                 D[1][1]=D_scalar;
    7380 
    7381                 TripleMultiply(B,2,numnodes,1,
    7382                                         &D[0][0],2,2,0,
    7383                                         B,2,numnodes,0,
    7384                                         &Ke->values[0],1);
     7375                GetNodalFunctions(basis,gauss);
     7376                GetNodalFunctionsDerivatives(dbasis,&xyz_list[0][0],gauss);
     7377
     7378                for(int i=0;i<numnodes;i++){
     7379                        for(int j=0;j<numnodes;j++){
     7380                                Ke->values[i*numnodes+j] += gauss->weight*Jdet*(
     7381                                                        basis[i]*basis[j]
     7382                                                        +(l*thickness)*(l*thickness)*(dbasis[0*numnodes+i]*dbasis[0*numnodes+j] + dbasis[1*numnodes+i]*dbasis[1*numnodes+j])
     7383                                                        );
     7384                        }
     7385                }
    73857386        }
    73867387
    73877388        /*Clean up and return*/
    73887389        delete gauss;
    7389         xDelete<IssmDouble>(B);
     7390        xDelete<IssmDouble>(dbasis);
     7391        xDelete<IssmDouble>(basis);
    73907392        return Ke;
    73917393}
    … …  
    74107412        /*Intermediaries */
    74117413        IssmDouble xyz_list[NUMVERTICES][3];
    7412         IssmDouble dhdt_g,mb_g,ms_g,Jdettria;
     7414        IssmDouble dhdt_g,mb_g,ms_g,Jdet;
    74137415
    74147416        /*Fetch number of nodes and dof for this finite element*/
    … …  
    74357437                dhdt_input->GetInputValue(&dhdt_g,gauss);
    74367438
    7437                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     7439                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    74387440                GetNodalFunctions(basis,gauss);
    74397441
    7440                 for(int i=0;i<numnodes;i++) pe->values[i]+=Jdettria*gauss->weight*(ms_g-mb_g-dhdt_g)*basis[i];
     7442                for(int i=0;i<numnodes;i++) pe->values[i]+=Jdet*gauss->weight*(ms_g-mb_g-dhdt_g)*basis[i];
    74417443        }
    74427444
    … …  
    74527454        /*Intermediaries */
    74537455        IssmDouble xyz_list[NUMVERTICES][3];
    7454         IssmDouble mb_g,ms_g,dhdt_g,Jdettria;
     7456        IssmDouble mb_g,ms_g,dhdt_g,Jdet;
    74557457
    74567458        /*Fetch number of nodes and dof for this finite element*/
    … …  
    74777479                dhdt_input->GetInputValue(&dhdt_g,gauss);
    74787480
    7479                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     7481                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    74807482                GetNodalFunctions(basis,gauss);
    74817483
    7482                 for(int i=0;i<numnodes;i++) pe->values[i]+=Jdettria*gauss->weight*(ms_g-mb_g-dhdt_g)*basis[i];
     7484                for(int i=0;i<numnodes;i++) pe->values[i]+=Jdet*gauss->weight*(ms_g-mb_g-dhdt_g)*basis[i];
    74837485        }
    74847486
    … …  
    74947496        /*Intermediaries */
    74957497        IssmDouble xyz_list[NUMVERTICES][3];
    7496         IssmDouble dhdt_g,mb_g,ms_g,Jdettria;
     7498        IssmDouble dhdt_g,mb_g,ms_g,Jdet;
    74977499        IssmDouble h,gamma,thickness;
    74987500        IssmDouble hnx,hny,dhnx[2],dhny[2];
    … …  
    75487550                if(thickness<50.) thickness=50.;
    75497551
    7550                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     7552                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    75517553                GetNodalFunctions(basis,gauss);
    75527554                GetNodalFunctionsDerivatives(dbasis,&xyz_list[0][0],gauss);
    … …  
    75607562
    75617563                for(int i=0;i<numnodes;i++){
    7562                         pe->values[i]+=Jdettria*gauss->weight*(ms_g-mb_g-dhdt_g)*(
     7564                        pe->values[i]+=Jdet*gauss->weight*(ms_g-mb_g-dhdt_g)*(
    75637565                                                basis[i] + gamma*(basis[i]*(dhnx[0]+dhny[1])+hnx*dbasis[0*numnodes+i] + hny*dbasis[1*numnodes+i])
    75647566                                                );
    … …  
    75837585        /*Intermediaries */
    75847586        IssmDouble xyz_list[NUMVERTICES][3];
    7585         IssmDouble Jdettria;
     7587        IssmDouble Jdet;
    75867588        IssmDouble thickness,slope[2];
    75877589        IssmDouble taud_x;
    … …  
    76087610                surface_input->GetInputDerivativeValue(&slope[0],&xyz_list[0][0],gauss);
    76097611                taud_x = matpar->GetRhoIce()*matpar->GetG()*thickness*slope[0];
    7610 
    7611                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     7612                //taud_x = slope[0];
     7613
     7614                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    76127615                GetNodalFunctions(basis,gauss);
    76137616
    7614                 for(int i=0;i<numnodes;i++) pe->values[i]+=Jdettria*gauss->weight*taud_x*basis[i];
     7617                for(int i=0;i<numnodes;i++) pe->values[i]+=Jdet*gauss->weight*taud_x*basis[i];
    76157618        }
    76167619
    … …  
    76267629        /*Intermediaries */
    76277630        IssmDouble xyz_list[NUMVERTICES][3];
    7628         IssmDouble Jdettria;
     7631        IssmDouble Jdet;
    76297632        IssmDouble thickness,slope[2];
    76307633        IssmDouble taud_y;
    … …  
    76517654                surface_input->GetInputDerivativeValue(&slope[0],&xyz_list[0][0],gauss);
    76527655                taud_y = matpar->GetRhoIce()*matpar->GetG()*thickness*slope[1];
    7653 
    7654                 GetJacobianDeterminant(&Jdettria, &xyz_list[0][0],gauss);
     7656                //taud_y = slope[1];
     7657
     7658                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    76557659                GetNodalFunctions(basis,gauss);
    76567660
    7657                 for(int i=0;i<numnodes;i++) pe->values[i]+=Jdettria*gauss->weight*taud_y*basis[i];
     7661                for(int i=0;i<numnodes;i++) pe->values[i]+=Jdet*gauss->weight*taud_y*basis[i];
    76587662        }
    76597663
Note: See TracChangeset for help on using the changeset viewer.