Changeset 26277
- Timestamp:
- 05/18/21 11:07:46 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp
r26234 r26277 196 196 } /*}}}*/ 197 197 break; 198 case MatlithoEnum: 199 { /*{{{*/ 198 case MatlithoEnum: { /*{{{*/ 200 199 iomodel->FetchData(13,"md.materials.radius","md.materials.viscosity","md.materials.lame_lambda","md.materials.lame_mu","md.materials.burgers_viscosity","md.materials.burgers_mu","md.materials.ebm_alpha","md.materials.ebm_delta","md.materials.ebm_taul","md.materials.ebm_tauh","md.materials.rheologymodel","md.materials.issolid","md.materials.density"); 201 200 materials->AddObject(new Matlitho(iomodel->numberofelements+1,iomodel)); … … 204 203 /*}}}*/ 205 204 break; 206 case MathydroEnum: /*{{{*/ 207 { 205 case MathydroEnum: {/*{{{*/ 208 206 /*If we don't have any materials pointed to by elements (meaning, if we are running only litho or hydro), 209 207 * then we need to zero out the hmaterial pointers inside the elements dataset so that it won't error out … … 239 237 iomodel->FetchDataToInput(inputs,elements,"md.materials.rheology_n",MaterialsRheologyNEnum); 240 238 iomodel->FetchDataToInput(inputs,elements,"md.materials.rheology_E",MaterialsRheologyEEnum); 241 for (int i=0;i<iomodel->numberofelements;i++) if(iomodel->my_elements[i]) materials->AddObject(new Matice(i+1,i,MatenhancediceEnum));239 for(int j=0;j<iomodel->numberofelements;j++) if(iomodel->my_elements[j]) materials->AddObject(new Matice(j+1,j,MatenhancediceEnum)); 242 240 switch(iomodel->domaindim){ 243 241 case 2: … … 256 254 iomodel->FetchDataToInput(inputs,elements,"md.materials.rheology_n",MaterialsRheologyNEnum); 257 255 iomodel->FetchDataToInput(inputs,elements,"md.damage.D",DamageDEnum); 258 for (int i=0;i<iomodel->numberofelements;i++) if(iomodel->my_elements[i]) materials->AddObject(new Matice(i+1,i,MatdamageiceEnum));256 for (int j=0;j<iomodel->numberofelements;j++) if(iomodel->my_elements[j]) materials->AddObject(new Matice(j+1,j,MatdamageiceEnum)); 259 257 switch(iomodel->domaindim){ 260 258 case 2: … … 273 271 iomodel->FetchDataToInput(inputs,elements,"md.materials.rheology_Ec",MaterialsRheologyEcEnum); 274 272 iomodel->FetchDataToInput(inputs,elements,"md.materials.rheology_Es",MaterialsRheologyEsEnum); 275 for(int i=0;i<iomodel->numberofelements;i++) if(iomodel->my_elements[i]) materials->AddObject(new Matestar(i+1,i,iomodel));273 for(int j=0;j<iomodel->numberofelements;j++) if(iomodel->my_elements[j]) materials->AddObject(new Matestar(j+1,j,iomodel)); 276 274 switch(iomodel->domaindim){ 277 275 case 2:
Note:
See TracChangeset
for help on using the changeset viewer.