Changeset 976


Ignore:
Timestamp:
06/13/09 09:41:11 (15 years ago)
Author:
seroussi
Message:

added Bbar in matice

Location:
issm/trunk/src/c
Files:
11 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp

    r816 r976  
    7575        int    matice_mid;
    7676        double matice_B;
     77        double matice_Bbar;
    7778        double matice_n;
    7879       
    … …  
    312313                                B_avg=B_avg/3;
    313314                                matice_B=B_avg;
     315                                matice_Bbar=B_avg;
    314316                                matice_n=(double)*(model->n+i);
    315317                               
    316318                                /*Create matice using its constructor:*/
    317                                 matice= new Matice(matice_mid,matice_B,matice_n);
     319                                matice= new Matice(matice_mid,matice_B,matice_Bbar,matice_n);
    318320               
    319321                                /*Add matice element to materials dataset: */
    … …  
    439441                                B_avg=B_avg/6;
    440442                                matice_B= B_avg;
     443                                matice_Bbar= B_avg;
    441444                                matice_n=(double)*(model->n+i);
    442445               
    443446                                /*Create matice using its constructor:*/
    444                                 matice= new Matice(matice_mid,matice_B,matice_n);
     447                                matice= new Matice(matice_mid,matice_B,matice_Bbar,matice_n);
    445448               
    446449                                /*Add matice element to materials dataset: */

  • issm/trunk/src/c/ModelProcessorx/DiagnosticHutter/CreateElementsNodesAndMaterialsDiagnosticHutter.cpp

    r816 r976  
    5555        int    matice_mid;
    5656        double matice_B;
     57        double matice_Bbar;
    5758        double matice_n;
    5859
    … …  
    248249                                matice_mid=i+1; //same as the material id from the geom2 elements.
    249250                                matice_B=model->B[i];   
     251                                matice_Bbar=model->B[i];       
    250252                                matice_n=(double)model->n[1]; //n defined on elements not grids, so take the first value everywhere
    251253                       
    252254                                /*Create matice using its constructor:*/
    253                                 matice= new Matice(matice_mid,matice_B,matice_n);
     255                                matice= new Matice(matice_mid,matice_B,matice_Bbar,matice_n);
    254256       
    255257                                /*Add matice element to materials dataset: */
    … …  
    301303                                        matice_mid=i+1; //same as the material id from the geom2 elements.
    302304                                        matice_B=model->B[i];   
     305                                        matice_Bbar=model->B[i];       
    303306                                        matice_n=(double)model->n[1]; //n defined on elements not grids, so take the first value everywhere
    304307                               
    305308                                        /*Create matice ubeam its constructor:*/
    306                                         matice= new Matice(matice_mid,matice_B,matice_n);
     309                                        matice= new Matice(matice_mid,matice_B,matice_Bbar,matice_n);
    307310               
    308311                                        /*Add matice element to materials dataset: */

  • issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateElementsNodesAndMaterialsDiagnosticStokes.cpp

    r864 r976  
    5454        int    matice_mid;
    5555        double matice_B;
     56        double matice_Bbar;
    5657        double matice_n;
    5758       
    … …  
    285286                        B_avg=B_avg/6;
    286287                        matice_B= B_avg;
     288                        matice_Bbar= B_avg;
    287289                        matice_n=(double)*(model->n+i);
    288290                       
    … …  
    290292       
    291293                                /*Create matice using its constructor:*/
    292                                 matice= new Matice(matice_mid,matice_B,matice_n);
     294                                matice= new Matice(matice_mid,matice_B,matice_Bbar,matice_n);
    293295               
    294296                                /*Add matice element to materials dataset: */

  • issm/trunk/src/c/ModelProcessorx/DiagnosticVert/CreateElementsNodesAndMaterialsDiagnosticVert.cpp

    r816 r976  
    5353        int    matice_mid;
    5454        double matice_B;
     55        double matice_Bbar;
    5556        double matice_n;
    5657       
    … …  
    247248
    248249                /*Create matice using its constructor:*/
    249                 matice= new Matice(matice_mid,matice_B,matice_n);
     250                matice= new Matice(matice_mid,matice_B,matice_Bbar,matice_n);
    250251
    251252                /*Add matice element to materials dataset: */

  • issm/trunk/src/c/ModelProcessorx/Melting/CreateElementsNodesAndMaterialsMelting.cpp

    r816 r976  
    5353        int    matice_mid;
    5454        double matice_B;
     55        double matice_Bbar;
    5556        double matice_n;
    5657       
    … …  
    256257                B_avg=B_avg/6;
    257258                matice_B= B_avg;
     259                matice_Bbar= B_avg;
    258260                matice_n=(double)*(model->n+i);
    259261
    260262                /*Create matice using its constructor:*/
    261                 matice= new Matice(matice_mid,matice_B,matice_n);
     263                matice= new Matice(matice_mid,matice_B,matice_Bbar,matice_n);
    262264
    263265                /*Add matice element to materials dataset: */

  • issm/trunk/src/c/ModelProcessorx/Prognostic/CreateElementsNodesAndMaterialsPrognostic.cpp

    r820 r976  
    7575        int    matice_mid;
    7676        double matice_B;
     77        double matice_Bbar;
    7778        double matice_n;
    7879       

  • issm/trunk/src/c/ModelProcessorx/Thermal/CreateElementsNodesAndMaterialsThermal.cpp

    r816 r976  
    5454        int    matice_mid;
    5555        double matice_B;
     56        double matice_Bbar;
    5657        double matice_n;
    5758       
    … …  
    252253                B_avg=B_avg/6;
    253254                matice_B= B_avg;
     255                matice_Bbar= B_avg;
    254256                matice_n=(double)*(model->n+i);
    255257
    256258                /*Create matice using its constructor:*/
    257                 matice= new Matice(matice_mid,matice_B,matice_n);
     259                matice= new Matice(matice_mid,matice_B,matice_Bbar,matice_n);
    258260
    259261                /*Add matice element to materials dataset: */

  • issm/trunk/src/c/objects/Matice.cpp

    r945 r976  
    2020        return;
    2121}
    22 Matice::Matice(int matice_mid,double matice_B,double matice_n){
     22Matice::Matice(int matice_mid,double matice_B,double matice_Bbar,double matice_n){
    2323        mid=matice_mid;
    2424        B=matice_B;
     25        Bbar=matice_Bbar;
    2526        n=matice_n;
    2627        return;
    … …  
    3536        printf("   mid: %i\n",mid);
    3637        printf("   B: %g\n",B);
     38        printf("   Bbar: %g\n",Bbar);
    3739        printf("   n: %g\n",n);
    3840        return;
    … …  
    4345        printf("   mid: %i\n",mid);
    4446        printf("   B: %g\n",B);
     47        printf("   Bbar: %g\n",Bbar);
    4548        printf("   n: %g\n",n);
    4649        return;
    … …  
    6366        memcpy(marshalled_dataset,&mid,sizeof(mid));marshalled_dataset+=sizeof(mid);
    6467        memcpy(marshalled_dataset,&B,sizeof(B));marshalled_dataset+=sizeof(B);
     68        memcpy(marshalled_dataset,&Bbar,sizeof(Bbar));marshalled_dataset+=sizeof(Bbar);
    6569        memcpy(marshalled_dataset,&n,sizeof(n));marshalled_dataset+=sizeof(n);
    6670
    … …  
    7175int   Matice::MarshallSize(){
    7276
    73         return sizeof(mid)+sizeof(B)+sizeof(n)+sizeof(int); //sizeof(int) for enum type
     77        return sizeof(mid)+sizeof(B)+sizeof(Bbar)+sizeof(n)+sizeof(int); //sizeof(int) for enum type
    7478}
    7579
    … …  
    9195        memcpy(&mid,marshalled_dataset,sizeof(mid));marshalled_dataset+=sizeof(mid);
    9296        memcpy(&B,marshalled_dataset,sizeof(B));marshalled_dataset+=sizeof(B);
     97        memcpy(&Bbar,marshalled_dataset,sizeof(Bbar));marshalled_dataset+=sizeof(Bbar);
    9398        memcpy(&n,marshalled_dataset,sizeof(n));marshalled_dataset+=sizeof(n);
    9499
    … …  
    127132               
    128133
     134#undef __FUNCT__
     135#define __FUNCT__ "Matice::SetBbar"
     136void  Matice::SetBbar(double Bbar_param){
     137        Bbar=Bbar_param;
     138}
     139               
    129140
    130141#undef __FUNCT__
    … …  
    158169        if (n==1){
    159170                /*Viscous behaviour! viscosity=B: */
    160                 viscosity=B;
     171                viscosity=Bbar;
    161172        }
    162173        else{
    … …  
    178189                        else{
    179190                                e=(n-1)/2/n;
    180                                 viscosity=2*B/(2*pow(A,e));
     191                                viscosity=2*Bbar/(2*pow(A,e));
    181192                        }
    182193                }
    … …  
    387398}
    388399
     400double Matice::GetBbar(){
     401        return Bbar;
     402}
     403
    389404double Matice::GetN(){
    390405        return n;

  • issm/trunk/src/c/objects/Matice.h

    r803 r976  
    1313                int        mid;
    1414                double B;
     15                double Bbar;
    1516                double n;
    1617
    … …  
    1819
    1920                Matice();
    20                 Matice(int mid,double B,double n);
     21                Matice(int mid,double B,double Bbar,double n);
    2122                ~Matice();
    2223
    … …  
    3334                void  UpdateFromInputs(void* inputs);
    3435                void  SetB(double B_param);
     36                void  SetBbar(double Bbar_param);
    3537                void  GetViscosity2d(double* pviscosity, double* pepsilon);
    3638                void  GetViscosity2(double* pviscosity2, double* pepsilon);
    … …  
    3941                Object* copy();
    4042                double GetB();
     43                double GetBbar();
    4144                double GetN();
    4245

  • issm/trunk/src/c/objects/Penta.cpp

    r950 r976  
    11011101
    11021102        //Update material if necessary
    1103         if(inputs->Recover("temperature_average",&temperature_list[0],1,dofs,6,(void**)nodes)){
    1104                 if(matice){
    1105                         temperature_average=(temperature_list[0]+temperature_list[1]+temperature_list[2]+temperature_list[3]+temperature_list[4]+temperature_list[5])/6.0;
    1106                         B_average=Paterson(temperature_average);
    1107                         matice->SetB(B_average);
    1108                 }
    1109         }
    11101103        if(inputs->Recover("temperature",&temperature_list[0],1,dofs,6,(void**)nodes)){
    11111104                if(matice){
    … …  
    11131106                        B_average=Paterson(temperature_average);
    11141107                        matice->SetB(B_average);
     1108                }
     1109        }
     1110       
     1111        if(inputs->Recover("temperature_average",&temperature_list[0],1,dofs,6,(void**)nodes)){
     1112                if(matice){
     1113                        temperature_average=(temperature_list[0]+temperature_list[1]+temperature_list[2]+temperature_list[3]+temperature_list[4]+temperature_list[5])/6.0;
     1114                        B_average=Paterson(temperature_average);
     1115                        matice->SetBbar(B_average);
    11151116                }
    11161117        }

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

    r950 r976  
    14501450                temperature_average=(temperature_list[0]+temperature_list[1]+temperature_list[2])/3.0;
    14511451                B_average=Paterson(temperature_average);
    1452                 matice->SetB(B_average);
     1452                matice->SetBbar(B_average);
    14531453        }
    14541454       
    14551455        if(inputs->Recover("B",&B_list[0],1,dofs,3,(void**)nodes)){
    14561456                B_average=(B_list[0]+B_list[1]+B_list[2])/3.0;
    1457                 matice->SetB(B_average);
     1457                matice->SetBbar(B_average);
    14581458        }
    14591459
Note: See TracChangeset for help on using the changeset viewer.