Changeset 3597
- Timestamp:
- 04/21/10 16:12:34 (15 years ago)
- Location:
- issm/trunk/src/c/objects
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/objects/Icefront.cpp
r3570 r3597 426 426 SegmentPressureLoad(pe_g,matpar->GetRhoWater(),matpar->GetRhoIce(),matpar->GetG(),thickness_list,bed_list,normal,length); 427 427 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 #endif445 446 428 /*Plug pe_g into vector: */ 447 429 VecSetValues(pg,numdofs,doflist,pe_g,ADD_VALUES); … … 580 562 QuadPressureLoad(pe_g,matpar->GetRhoWater(),matpar->GetRhoIce(),matpar->GetG(),thickness_list_quad,bed_list_quad,normal1,normal2,normal3,normal4,&xyz_list_quad[0][0]); 581 563 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 #endif602 603 564 /*Plug pe_g into vector: */ 604 565 VecSetValues(pg,numdofs,doflist,pe_g,ADD_VALUES); … … 736 697 //Compute load contribution for this quad: 737 698 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 #endif758 699 759 700 /*Plug pe_g into vector: */ -
issm/trunk/src/c/objects/Penta.cpp
r3595 r3597 738 738 } 739 739 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 #endif758 759 740 /*Get gaussian points and weights. Penta is an extrusion of a Tria, we therefore 760 741 get tria gaussian points as well as segment gaussian points. For tria gaussian … … 768 749 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); 769 750 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 #endif781 751 /* Start looping on the number of gaussian points: */ 782 752 for (ig1=0; ig1<num_area_gauss; ig1++){ … … 3225 3195 l1l2l3[1]=gauss_coord_tria[1]; 3226 3196 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 #endif3234 3197 3235 3198 /*Build LStokes: */ … … 3334 3297 GetNodalFunctionsDerivatives(&dh1dh6[0][0],xyz_list,gauss_coord); 3335 3298 3336 #ifdef _DELUG_3337 for (i=0;i<3;i++){3338 printf("Node %i h=%lf \n",i,l1l2l3[i]);3339 }3340 #endif3341 3342 3299 /*Build LprimeStokes: */ 3343 3300 for (i=0;i<3;i++){
Note:
See TracChangeset
for help on using the changeset viewer.