Changeset 1764
- Timestamp:
- 08/18/09 15:55:14 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/ModelProcessorx/DiagnosticStokes/CreateLoadsDiagnosticStokes.cpp
r483 r1764 62 62 /*Create pressure loads as boundary conditions. Pay attention to the partitioning if we are running in parallel (the grids 63 63 * referenced by a certain load must belong to the cluster node): */ 64 ModelFetchData((void**)&model-> segmentonneumann_diag_stokes,&numberofsegs_diag_stokes,NULL,model_handle,"segmentonneumann_diag_stokes","Matrix","Mat");64 ModelFetchData((void**)&model->pressureloads_stokes,&numberofpressureloads_stokes,NULL,model_handle,"pressureload_stokes","Matrix","Mat"); 65 65 ModelFetchData((void**)&model->elements_type,NULL,NULL,model_handle,"elements_type","Matrix","Mat"); 66 66 ModelFetchData((void**)&model->thickness,NULL,NULL,model_handle,"thickness","Matrix","Mat"); … … 70 70 71 71 /*First load data:*/ 72 for (i=0;i<numberof segs_diag_stokes;i++){72 for (i=0;i<numberofpressureloads_stokes;i++){ 73 73 74 74 segment_width=5; 75 75 76 element=(int)(*(model-> segmentonneumann_diag_stokes+segment_width*i+segment_width-1)-1); //element is in the last column76 element=(int)(*(model->pressureload_stokes+segment_width*i+segment_width-1)-1); //element is in the last column 77 77 78 78 #ifdef _PARALLEL_ … … 86 86 icefront_mparid=model->numberofelements+1; //matlab indexing 87 87 icefront_sid=count+1; //matlab indexing 88 icefront_eid=(int)*(model-> segmentonneumann_diag_stokes+segment_width*i+segment_width-1); //matlab indexing88 icefront_eid=(int)*(model->pressureload_stokes+segment_width*i+segment_width-1); //matlab indexing 89 89 icefront_element_type=PentaEnum(); 90 90 91 i1=(int)*(model-> segmentonneumann_diag_stokes+segment_width*i+0);92 i2=(int)*(model-> segmentonneumann_diag_stokes+segment_width*i+1);93 i3=(int)*(model-> segmentonneumann_diag_stokes+segment_width*i+2);94 i4=(int)*(model-> segmentonneumann_diag_stokes+segment_width*i+3);91 i1=(int)*(model->pressureload_stokes+segment_width*i+0); 92 i2=(int)*(model->pressureload_stokes+segment_width*i+1); 93 i3=(int)*(model->pressureload_stokes+segment_width*i+2); 94 i4=(int)*(model->pressureload_stokes+segment_width*i+3); 95 95 96 96 icefront_node_ids[0]=i1; … … 116 116 } 117 117 /*Free data: */ 118 xfree((void**)&model-> segmentonneumann_diag_stokes);118 xfree((void**)&model->pressureload_stokes); 119 119 xfree((void**)&model->elements_type); 120 120 xfree((void**)&model->thickness);
Note:
See TracChangeset
for help on using the changeset viewer.