Changeset 22192
- Timestamp:
- 10/26/17 16:24:04 (7 years ago)
- Location:
- issm/trunk-jpl
- Files:
-
- 1 added
- 1 deleted
- 11 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
r22181 r22192 60 60 iomodel->FetchDataToInput(elements,"md.calving.meltingrate",CalvinglevermannMeltingrateEnum); 61 61 break; 62 case Calving DevEnum:62 case CalvingVonmisesEnum: 63 63 iomodel->FetchDataToInput(elements,"md.calving.stress_threshold_groundedice",CalvingStressThresholdGroundediceEnum); 64 64 iomodel->FetchDataToInput(elements,"md.calving.stress_threshold_floatingice",CalvingStressThresholdFloatingiceEnum); … … 89 89 case DefaultCalvingEnum: 90 90 case CalvingLevermannEnum: 91 case Calving DevEnum:91 case CalvingVonmisesEnum: 92 92 break; 93 93 case CalvingMinthicknessEnum: … … 220 220 switch(calvinglaw){ 221 221 case DefaultCalvingEnum: 222 case Calving DevEnum:222 case CalvingVonmisesEnum: 223 223 lsf_slopex_input = basalelement->GetInput(LevelsetfunctionSlopeXEnum); _assert_(lsf_slopex_input); 224 224 if(dim==2) lsf_slopey_input = basalelement->GetInput(LevelsetfunctionSlopeYEnum); _assert_(lsf_slopey_input); … … 291 291 switch(calvinglaw){ 292 292 case DefaultCalvingEnum: 293 case Calving DevEnum:293 case CalvingVonmisesEnum: 294 294 lsf_slopex_input->GetInputValue(&dlsf[0],gauss); 295 295 if(dim==2) lsf_slopey_input->GetInputValue(&dlsf[1],gauss); -
issm/trunk-jpl/src/c/classes/Elements/Element.cpp
r22181 r22192 2353 2353 this->CalvingRateLevermann(); 2354 2354 break; 2355 case Calving DevEnum:2356 this->CalvingRate Dev();2355 case CalvingVonmisesEnum: 2356 this->CalvingRateVonmises(); 2357 2357 break; 2358 2358 case CalvingCrevasseDepthEnum: … … 2366 2366 case StrainRateperpendicularEnum: this->StrainRateperpendicular(); break; 2367 2367 case SurfaceCrevasseEnum: this->CalvingCrevasseDepth(); break; 2368 case SigmaVMEnum: this->CalvingRate Dev(); break;2368 case SigmaVMEnum: this->CalvingRateVonmises(); break; 2369 2369 } 2370 2370 -
issm/trunk-jpl/src/c/classes/Elements/Element.h
r22181 r22192 188 188 virtual void AverageOntoPartition(Vector<IssmDouble>* partition_contributions,Vector<IssmDouble>* partition_areas,IssmDouble* vertex_response,IssmDouble* qmu_part)=0; 189 189 virtual void BasalNodeIndices(int* pnumindices,int** pindices,int finiteelement){_error_("not implemented yet");}; 190 virtual void CalvingRate Dev(void){_error_("not implemented yet");};190 virtual void CalvingRateVonmises(void){_error_("not implemented yet");}; 191 191 virtual void CalvingCrevasseDepth(void){_error_("not implemented yet");}; 192 192 virtual void CalvingRateLevermann(void)=0; -
issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
r22124 r22192 182 182 } 183 183 /*}}}*/ 184 void Penta::CalvingRate Dev(){/*{{{*/184 void Penta::CalvingRateVonmises(){/*{{{*/ 185 185 186 186 if(!this->IsOnBase()) return; -
issm/trunk-jpl/src/c/classes/Elements/Penta.h
r21931 r22192 49 49 void AverageOntoPartition(Vector<IssmDouble>* partition_contributions,Vector<IssmDouble>* partition_areas,IssmDouble* vertex_response,IssmDouble* qmu_part); 50 50 void BasalNodeIndices(int* pnumindices,int** pindices,int finiteelement); 51 void CalvingRate Dev();51 void CalvingRateVonmises(); 52 52 void CalvingRateLevermann(); 53 53 IssmDouble CharacteristicLength(void){_error_("not implemented yet");}; -
issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
r22183 r22192 208 208 } 209 209 /*}}}*/ 210 void Tria::CalvingRate Dev(){/*{{{*/210 void Tria::CalvingRateVonmises(){/*{{{*/ 211 211 212 212 IssmDouble xyz_list[NUMVERTICES][3]; -
issm/trunk-jpl/src/c/classes/Elements/Tria.h
r22181 r22192 52 52 /*Element virtual functions definitions: {{{*/ 53 53 void AverageOntoPartition(Vector<IssmDouble>* partition_contributions,Vector<IssmDouble>* partition_areas,IssmDouble* vertex_response,IssmDouble* qmu_part); 54 void CalvingRate Dev();54 void CalvingRateVonmises(); 55 55 void CalvingCrevasseDepth(); 56 56 void CalvingRateLevermann(); -
issm/trunk-jpl/src/c/classes/FemModel.cpp
r22181 r22192 920 920 921 921 }/*}}}*/ 922 void FemModel::CalvingRate Devx(){/*{{{*/922 void FemModel::CalvingRateVonmisesx(){/*{{{*/ 923 923 924 924 for(int i=0;i<elements->Size();i++){ 925 925 Element* element=dynamic_cast<Element*>(this->elements->GetObjectByOffset(i)); 926 element->CalvingRate Dev();926 element->CalvingRateVonmises(); 927 927 } 928 928 } -
issm/trunk-jpl/src/c/classes/FemModel.h
r22181 r22192 82 82 /*Modules*/ 83 83 void BalancethicknessMisfitx(IssmDouble* pV); 84 void CalvingRate Devx();84 void CalvingRateVonmisesx(); 85 85 void CalvingRateLevermannx(); 86 86 void DeviatoricStressx(); -
issm/trunk-jpl/src/c/modules/Calvingx/Calvingx.cpp
r22181 r22192 31 31 femmodel->CalvingRateLevermannx(); 32 32 break; 33 case Calving DevEnum:34 femmodel->CalvingRate Devx();35 femmodel->ElementOperationx(&Element::CalvingRate Dev);33 case CalvingVonmisesEnum: 34 femmodel->CalvingRateVonmisesx(); 35 femmodel->ElementOperationx(&Element::CalvingRateVonmises); 36 36 break; 37 37 default: -
issm/trunk-jpl/test/NightlyRun/test540.m
r22046 r22192 9 9 %calving parameters 10 10 md.mask.ice_levelset = 1e4*(md.mask.ice_levelset + 0.5); 11 md.calving=calving dev();11 md.calving=calvingvonmises(); 12 12 md.calving.meltingrate = zeros(md.mesh.numberofvertices,1); 13 13 md.transient.ismovingfront = 1;
Note:
See TracChangeset
for help on using the changeset viewer.