Changeset 3597


Ignore:
Timestamp:
04/21/10 16:12:34 (15 years ago)
Author:
Mathieu Morlighem
Message:

removed useless ISSM_DEBUG outputs

Location:
issm/trunk/src/c/objects
Files:
2 edited

Legend:

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

    r3570 r3597  
    426426        SegmentPressureLoad(pe_g,matpar->GetRhoWater(),matpar->GetRhoIce(),matpar->GetG(),thickness_list,bed_list,normal,length);
    427427               
    428         #ifdef _DEBUGELEMENTS_
    429         if(my_rank==RANK && this->icefront.eid==ELID){
    430                 printf("\nicefront load\n");
    431                 printf("grids %i %i\n",grid1,grid2);
    432                 printf("rho_water %g\n",rho_water);
    433                 printf("rho_ice %g\n",rho_ice);
    434                 printf("gravity %g\n",gravity);
    435                 printf("thickness_list (%g,%g)\n",thickness_list[0],thickness_list[1]);
    436                 printf("bed_list (%g,%g)\n",bed_list[0],bed_list[1]);
    437                 printf("normal (%g,%g)\n",normal[0],normal[1]);
    438                 printf("length %g\n",length);
    439                 printf("pe_g->terms\n");
    440                 for(i=0;i<numgrids*NDOF2;i++){
    441                         printf("%g ",*(pe_g->terms+i));
    442                 }
    443         }
    444         #endif
    445 
    446428        /*Plug pe_g into vector: */
    447429        VecSetValues(pg,numdofs,doflist,pe_g,ADD_VALUES);
     
    580562        QuadPressureLoad(pe_g,matpar->GetRhoWater(),matpar->GetRhoIce(),matpar->GetG(),thickness_list_quad,bed_list_quad,normal1,normal2,normal3,normal4,&xyz_list_quad[0][0]);
    581563
    582 
    583         #ifdef _DEBUGELEMENTS_
    584         if(my_rank==RANK && eid==ELID){
    585                 printf("\nicefront load\n");
    586                 printf("grids %i %i %i %i\n",grid1,grid2,grid3,grid4);
    587                 printf("rho_water %g\n",rho_water);
    588                 printf("rho_ice %g\n",rho_ice);
    589                 printf("gravity %g\n",gravity);
    590                 printf("thickness_list (%g,%g,%g,%g)\n",thickness_list[0],thickness_list[1],thickness_list[2],thickness_list[3]);
    591                 printf("bed_list (%g,%g,%g,%g)\n",bed_list[0],bed_list[1],bed_list[2],bed_list[3]);
    592                 printf("normal1 (%g,%g,%g)\n",normal1[0],normal1[1],normal1[2]);
    593                 printf("normal2 (%g,%g,%g)\n",normal2[0],normal2[1],normal2[2]);
    594                 printf("normal3 (%g,%g,%g)\n",normal3[0],normal3[1],normal3[2]);
    595                 printf("normal4 (%g,%g,%g)\n",normal4[0],normal4[1],normal4[2]);
    596                 printf("pe_g->terms\n");
    597                 for(i=0;i<numgridsload*NDOF2;i++){
    598                         printf("%g ",*(pe_g->terms+i));
    599                 }
    600         }
    601         #endif
    602 
    603564        /*Plug pe_g into vector: */
    604565        VecSetValues(pg,numdofs,doflist,pe_g,ADD_VALUES);
     
    736697        //Compute load contribution for this quad:
    737698        QuadPressureLoadStokes(pe_g,matpar->GetRhoWater(),matpar->GetRhoIce(),matpar->GetG(),thickness_list_quad,bed_list_quad,normal1,normal2,normal3,normal4,&xyz_list_quad[0][0]);
    738 
    739         #ifdef _DEBUGELEMENTS_
    740         if(my_rank==RANK && eid==ELID){
    741                 printf("\nicefront load\n");
    742                 printf("grids %i %i %i %i\n",grid1,grid2,grid3,grid4);
    743                 printf("rho_water %g\n",rho_water);
    744                 printf("rho_ice %g\n",rho_ice);
    745                 printf("gravity %g\n",gravity);
    746                 printf("thickness_list (%g,%g,%g,%g)\n",thickness_list[0],thickness_list[1],thickness_list[2],thickness_list[3]);
    747                 printf("bed_list (%g,%g,%g,%g)\n",bed_list[0],bed_list[1],bed_list[2],bed_list[3]);
    748                 printf("normal1 (%g,%g,%g)\n",normal1[0],normal1[1],normal1[2]);
    749                 printf("normal2 (%g,%g,%g)\n",normal2[0],normal2[1],normal2[2]);
    750                 printf("normal3 (%g,%g,%g)\n",normal3[0],normal3[1],normal3[2]);
    751                 printf("normal4 (%g,%g,%g)\n",normal4[0],normal4[1],normal4[2]);
    752                 printf("pe_g->terms\n");
    753                 for(i=0;i<numdofloads*NDOF4;i++){
    754                         printf("%g ",*(pe_g->terms+i));
    755                 }
    756         }
    757         #endif
    758699
    759700        /*Plug pe_g into vector: */
  • issm/trunk/src/c/objects/Penta.cpp

    r3595 r3597  
    738738                }
    739739
    740 #ifdef _DEBUGELEMENTS_
    741                 if(my_rank==RANK && id==ELID){
    742                         printf("El id %i Rank %i PentaElement input list before gaussian loop: \n",ELID,RANK);
    743                         printf("   rho_ice: %g\n",matice->GetRhoIce());
    744                         printf("   rho_water:%g \n",matice->GetRhoWater());
    745                         printf("   gravity: %g\n",matpar->GetG());
    746                         printf("   Velocity: \n");
    747                         for (i=0;i<numgrids;i++){
    748                                 printf("      grid %i  [%g,%g,%g]\n",i,vxvy_list[i][0],vxvy_list[i][1],vxvy_list[i][2]);
    749                         }
    750                         printf("   B [%g %g %g %g %g %g]\n",B_list[0],B_list[1],B_list[2],B_list[3],B_list[4],B_list[5]);
    751                         printf("   K [%g %g %g %g %g %g]\n",K_list[0],K_list[1],K_list[2],K_list[3],K_list[4],K_list[5]);
    752                         printf("   thickness [%g %g %g %g %g %g]\n",thickness_list[0],thickness_list[1],thickness_list[2],thickness_list[3],thickness_list[4],thickness_list[5]);
    753                         printf("   surface [%g %g %g %g %g %g]\n",surface_list[0],surface_list[1],surface_list[2],surface_list[3],surface_list[4],surface_list[5]);
    754                         printf("   bed [%g %g %g %g %g %g]\n",bed_list[0],bed_list[1],bed_list[2],bed_list[3],bed_list[4],bed_list[5]);
    755                         printf("   temperature_average [%g %g %g %g %g %g]\n",temperature_average_list[0],temperature_average_list[1],temperature_average_list[2],temperature_average_list[3],temperature_average_list[4],temperature_average_list[5]);
    756                 }
    757 #endif
    758 
    759740                /*Get gaussian points and weights. Penta is an extrusion of a Tria, we therefore
    760741                  get tria gaussian points as well as segment gaussian points. For tria gaussian
     
    768749                GaussPenta( &num_area_gauss, &first_gauss_area_coord, &second_gauss_area_coord, &third_gauss_area_coord, &area_gauss_weights, &fourth_gauss_vert_coord,&vert_gauss_weights,order_area_gauss,num_vert_gauss);
    769750
    770 #ifdef _DEBUGGAUSS_
    771                 if(my_rank==RANK && id==ELID){
    772                         printf("El id %i Rank %i PentaElement gauss points\n",ELID,RANK);
    773                         for (i=0;i<num_area_gauss;i++){
    774                                 printf("   Area Gauss coord %i: %lf %lf %lf Weight: %lf\n",i,*(first_gauss_area_coord+i),*(second_gauss_area_coord+i),*(third_gauss_area_coord+i),*(area_gauss_weights+i));
    775                         }
    776                         for (i=0;i<num_vert_gauss;i++){
    777                                 printf("   Vert Gauss coord %i: %lf Weight: %lf\n",i,*(fourth_gauss_vert_coord+i),*(vert_gauss_weights+i));
    778                         }       
    779                 }
    780 #endif
    781751                /* Start  looping on the number of gaussian points: */
    782752                for (ig1=0; ig1<num_area_gauss; ig1++){
     
    32253195        l1l2l3[1]=gauss_coord_tria[1];
    32263196        l1l2l3[2]=gauss_coord_tria[2];
    3227 
    3228 
    3229 #ifdef _DELUG_
    3230         for (i=0;i<3;i++){
    3231                 printf("Node %i  h=%lf \n",i,l1l2l3[i]);
    3232         }
    3233 #endif
    32343197
    32353198        /*Build LStokes: */
     
    33343297        GetNodalFunctionsDerivatives(&dh1dh6[0][0],xyz_list,gauss_coord);
    33353298
    3336 #ifdef _DELUG_
    3337         for (i=0;i<3;i++){
    3338                 printf("Node %i  h=%lf \n",i,l1l2l3[i]);
    3339         }
    3340 #endif
    3341 
    33423299        /*Build LprimeStokes: */
    33433300        for (i=0;i<3;i++){
Note: See TracChangeset for help on using the changeset viewer.