Changeset 20645
- Timestamp:
- 05/25/16 10:52:44 (9 years ago)
- Location:
- issm/trunk-jpl/src/c
- Files:
-
- 36 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/BalancethicknessAnalysis.cpp
r20155 r20645 10 10 /*Fetch parameters: */ 11 11 int stabilization; 12 iomodel-> Constant(&stabilization,BalancethicknessStabilizationEnum);12 iomodel->FindConstant(&stabilization,BalancethicknessStabilizationEnum); 13 13 14 14 /*Do not add constraints in DG*/ … … 25 25 26 26 /*Fetch parameters: */ 27 iomodel-> Constant(&stabilization,BalancethicknessStabilizationEnum);27 iomodel->FindConstant(&stabilization,BalancethicknessStabilizationEnum); 28 28 29 29 /*Loads only in DG*/ … … 54 54 55 55 int stabilization; 56 iomodel-> Constant(&stabilization,BalancethicknessStabilizationEnum);56 iomodel->FindConstant(&stabilization,BalancethicknessStabilizationEnum); 57 57 58 58 /*Check in 3d*/ … … 77 77 78 78 /*Fetch data needed: */ 79 iomodel-> Constant(&stabilization,BalancethicknessStabilizationEnum);79 iomodel->FindConstant(&stabilization,BalancethicknessStabilizationEnum); 80 80 81 81 /*Finite element type*/ -
issm/trunk-jpl/src/c/analyses/DamageEvolutionAnalysis.cpp
r20459 r20645 9 9 10 10 int finiteelement; 11 iomodel-> Constant(&finiteelement,DamageElementinterpEnum);11 iomodel->FindConstant(&finiteelement,DamageElementinterpEnum); 12 12 13 13 /*Fetch parameters: */ 14 14 int stabilization; 15 iomodel-> Constant(&stabilization,DamageStabilizationEnum);15 iomodel->FindConstant(&stabilization,DamageStabilizationEnum); 16 16 17 17 /*Do not add constraints in DG, they are weakly imposed*/ … … 34 34 int finiteelement; 35 35 36 iomodel-> Constant(&finiteelement,DamageElementinterpEnum);36 iomodel->FindConstant(&finiteelement,DamageElementinterpEnum); 37 37 ::CreateNodes(nodes,iomodel,DamageEvolutionAnalysisEnum,finiteelement); 38 38 }/*}}}*/ … … 45 45 bool ismovingfront; 46 46 47 iomodel-> Constant(&finiteelement,DamageElementinterpEnum);48 iomodel-> Constant(&ismovingfront,TransientIsmovingfrontEnum);47 iomodel->FindConstant(&finiteelement,DamageElementinterpEnum); 48 iomodel->FindConstant(&ismovingfront,TransientIsmovingfrontEnum); 49 49 50 50 /*Update elements: */ … … 88 88 /*Retrieve law dependent parameters: */ 89 89 int law; 90 iomodel-> Constant(&law,DamageLawEnum);90 iomodel->FindConstant(&law,DamageLawEnum); 91 91 if (law==0){ 92 92 parameters->AddObject(iomodel->CopyConstantObject(DamageStressThresholdEnum)); -
issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
r20625 r20645 22 22 23 23 /*Fetch parameters: */ 24 iomodel-> Constant(&heatcapacity,MaterialsHeatcapacityEnum);25 iomodel-> Constant(&referencetemperature,ConstantsReferencetemperatureEnum);24 iomodel->FindConstant(&heatcapacity,MaterialsHeatcapacityEnum); 25 iomodel->FindConstant(&referencetemperature,ConstantsReferencetemperatureEnum); 26 26 27 27 /*return if 2d mesh*/ … … 115 115 116 116 /*Is enthalpy requested?*/ 117 iomodel-> Constant(&isenthalpy,ThermalIsenthalpyEnum);117 iomodel->FindConstant(&isenthalpy,ThermalIsenthalpyEnum); 118 118 if(!isenthalpy) return; 119 119 … … 131 131 } 132 132 133 iomodel-> Constant(&dakota_analysis,QmuIsdakotaEnum);134 iomodel-> Constant(&ismovingfront,TransientIsmovingfrontEnum);135 iomodel-> Constant(&frictionlaw,FrictionLawEnum);133 iomodel->FindConstant(&dakota_analysis,QmuIsdakotaEnum); 134 iomodel->FindConstant(&ismovingfront,TransientIsmovingfrontEnum); 135 iomodel->FindConstant(&frictionlaw,FrictionLawEnum); 136 136 137 137 iomodel->FetchDataToInput(elements,ThicknessEnum); … … 164 164 165 165 /*Basal forcings variables*/ 166 iomodel-> Constant(&basalforcing_model,BasalforcingsEnum);166 iomodel->FindConstant(&basalforcing_model,BasalforcingsEnum); 167 167 switch(basalforcing_model){ 168 168 case MantlePlumeGeothermalFluxEnum: … … 185 185 break; 186 186 case 3: 187 iomodel-> Constant(&FrictionCoupling,FrictionCouplingEnum);187 iomodel->FindConstant(&FrictionCoupling,FrictionCouplingEnum); 188 188 iomodel->FetchDataToInput(elements,FrictionCEnum); 189 189 iomodel->FetchDataToInput(elements,FrictionAsEnum); … … 237 237 /*Deal with friction parameters*/ 238 238 int frictionlaw; 239 iomodel-> Constant(&frictionlaw,FrictionLawEnum);239 iomodel->FindConstant(&frictionlaw,FrictionLawEnum); 240 240 if(frictionlaw==4 || frictionlaw==6) parameters->AddObject(iomodel->CopyConstantObject(FrictionGammaEnum)); 241 241 if(frictionlaw==3) parameters->AddObject(iomodel->CopyConstantObject(FrictionCouplingEnum)); -
issm/trunk-jpl/src/c/analyses/FreeSurfaceTopAnalysis.cpp
r20155 r20645 67 67 68 68 /*Fetch data needed: */ 69 iomodel-> Constant(&smb_model,SmbEnum);69 iomodel->FindConstant(&smb_model,SmbEnum); 70 70 71 71 /*Update elements: */ -
issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
r19095 r20645 16 16 bool isefficientlayer; 17 17 /*retrieve some parameters: */ 18 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);18 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 19 19 20 20 /*Now, do we really want DC?*/ … … 22 22 23 23 /*Do we want an efficient layer*/ 24 iomodel-> Constant(&isefficientlayer,HydrologydcIsefficientlayerEnum);24 iomodel->FindConstant(&isefficientlayer,HydrologydcIsefficientlayerEnum); 25 25 26 26 /*If not return*/ … … 42 42 43 43 /*Now, do we really want DC?*/ 44 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);44 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 45 45 if(hydrology_model!=HydrologydcEnum) return; 46 46 47 47 /*Do we want an efficient layer*/ 48 iomodel-> Constant(&isefficientlayer,HydrologydcIsefficientlayerEnum);48 iomodel->FindConstant(&isefficientlayer,HydrologydcIsefficientlayerEnum); 49 49 if(!isefficientlayer) return; 50 50 … … 76 76 /*Now, do we really want DC?*/ 77 77 int hydrology_model; 78 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);78 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 79 79 if(hydrology_model!=HydrologydcEnum) return; 80 80 81 81 /*Do we want an efficient layer*/ 82 82 bool isefficientlayer; 83 iomodel-> Constant(&isefficientlayer,HydrologydcIsefficientlayerEnum);83 iomodel->FindConstant(&isefficientlayer,HydrologydcIsefficientlayerEnum); 84 84 if(!isefficientlayer) return; 85 85 … … 93 93 /*Do we really want DC?*/ 94 94 int hydrology_model; 95 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);95 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 96 96 if(hydrology_model!=HydrologydcEnum) return; 97 97 98 98 /*Do we want an efficient layer*/ 99 99 bool isefficientlayer; 100 iomodel-> Constant(&isefficientlayer,HydrologydcIsefficientlayerEnum);100 iomodel->FindConstant(&isefficientlayer,HydrologydcIsefficientlayerEnum); 101 101 if(!isefficientlayer) return; 102 102 -
issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp
r19032 r20645 25 25 26 26 /*retrieve some parameters: */ 27 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);27 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 28 28 29 29 /*Now, do we really want DC?*/ … … 65 65 66 66 /*Fetch data needed: */ 67 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);67 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 68 68 69 69 /*Now, do we really want DC?*/ … … 71 71 72 72 /*Fetch data needed: */ 73 iomodel-> Constant(&isefficientlayer,HydrologydcIsefficientlayerEnum);73 iomodel->FindConstant(&isefficientlayer,HydrologydcIsefficientlayerEnum); 74 74 75 75 /*Update elements: */ … … 104 104 /*Fetch parameters: */ 105 105 int hydrology_model; 106 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);106 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 107 107 108 108 /*Now, do we really want DC?*/ … … 118 118 /*retrieve some parameters: */ 119 119 int hydrology_model; 120 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);120 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 121 121 if(hydrology_model!=HydrologydcEnum) return; 122 122 … … 128 128 /*Fetch parameters: */ 129 129 int hydrology_model; 130 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);130 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 131 131 if(hydrology_model!=HydrologydcEnum) return; 132 132 -
issm/trunk-jpl/src/c/analyses/HydrologyShreveAnalysis.cpp
r20155 r20645 10 10 /*retrieve some parameters: */ 11 11 int hydrology_model; 12 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);12 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 13 13 14 14 if(hydrology_model!=HydrologyshreveEnum) return; … … 24 24 /*Fetch parameters: */ 25 25 int hydrology_model; 26 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);26 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 27 27 28 28 /*Now, do we really want Shreve?*/ … … 40 40 /*Fetch data needed: */ 41 41 int hydrology_model; 42 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);42 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 43 43 44 44 /*Now, do we really want Shreve?*/ … … 74 74 /*retrieve some parameters: */ 75 75 int hydrology_model; 76 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);76 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 77 77 78 78 /*Now, do we really want Shreve?*/ -
issm/trunk-jpl/src/c/analyses/HydrologySommersAnalysis.cpp
r19931 r20645 14 14 /*retrieve some parameters: */ 15 15 int hydrology_model; 16 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);16 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 17 17 18 18 if(hydrology_model!=HydrologysommersEnum) return; … … 25 25 /*Fetch parameters: */ 26 26 int hydrology_model; 27 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);27 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 28 28 29 29 /*Now, do we really want Sommers?*/ … … 67 67 /*Fetch parameters: */ 68 68 int hydrology_model; 69 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);69 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 70 70 71 71 /*Now, do we really want Sommers?*/ … … 83 83 /*Fetch data needed: */ 84 84 int hydrology_model,frictionlaw; 85 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);85 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 86 86 87 87 /*Now, do we really want Sommers?*/ … … 119 119 iomodel->FetchDataToInput(elements,VyEnum); 120 120 121 iomodel-> Constant(&frictionlaw,FrictionLawEnum);121 iomodel->FindConstant(&frictionlaw,FrictionLawEnum); 122 122 /*Friction law variables*/ 123 123 switch(frictionlaw){ … … 138 138 /*retrieve some parameters: */ 139 139 int hydrology_model; 140 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);140 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 141 141 142 142 /*Now, do we really want Sommers?*/ -
issm/trunk-jpl/src/c/analyses/L2ProjectionEPLAnalysis.cpp
r18930 r20645 17 17 /*Now, do we really want DC?*/ 18 18 int hydrology_model; 19 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);19 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 20 20 if(hydrology_model!=HydrologydcEnum) return; 21 21 22 22 /*Do we want an efficient layer*/ 23 23 bool isefficientlayer; 24 iomodel-> Constant(&isefficientlayer,HydrologydcIsefficientlayerEnum);24 iomodel->FindConstant(&isefficientlayer,HydrologydcIsefficientlayerEnum); 25 25 if(!isefficientlayer) return; 26 26 … … 43 43 44 44 /*Now, do we really want DC?*/ 45 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);45 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 46 46 if(hydrology_model!=HydrologydcEnum) return; 47 47 48 48 /*Do we want an efficient layer*/ 49 iomodel-> Constant(&isefficientlayer,HydrologydcIsefficientlayerEnum);49 iomodel->FindConstant(&isefficientlayer,HydrologydcIsefficientlayerEnum); 50 50 if(!isefficientlayer) return; 51 51 -
issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
r20517 r20645 50 50 /*Get moving front parameters*/ 51 51 int calvinglaw; 52 iomodel-> Constant(&calvinglaw,CalvingLawEnum);52 iomodel->FindConstant(&calvinglaw,CalvingLawEnum); 53 53 switch(calvinglaw){ 54 54 case DefaultCalvingEnum: -
issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp
r20459 r20645 10 10 /*Fetch parameters: */ 11 11 int stabilization; 12 iomodel-> Constant(&stabilization,MasstransportStabilizationEnum);12 iomodel->FindConstant(&stabilization,MasstransportStabilizationEnum); 13 13 14 14 /*Do not add constraints in DG, they are weakly imposed*/ … … 32 32 33 33 /*Fetch parameters: */ 34 iomodel-> Constant(&stabilization,MasstransportStabilizationEnum);34 iomodel->FindConstant(&stabilization,MasstransportStabilizationEnum); 35 35 36 36 /*Loads only in DG*/ … … 98 98 /*Fetch parameters: */ 99 99 int stabilization; 100 iomodel-> Constant(&stabilization,MasstransportStabilizationEnum);100 iomodel->FindConstant(&stabilization,MasstransportStabilizationEnum); 101 101 102 102 /*Check in 3d*/ … … 124 124 125 125 /*Fetch data needed: */ 126 iomodel-> Constant(&stabilization,MasstransportStabilizationEnum);127 iomodel-> Constant(&dakota_analysis,QmuIsdakotaEnum);128 iomodel-> Constant(&isgroundingline,TransientIsgroundinglineEnum);129 iomodel-> Constant(&ismovingfront,TransientIsmovingfrontEnum);126 iomodel->FindConstant(&stabilization,MasstransportStabilizationEnum); 127 iomodel->FindConstant(&dakota_analysis,QmuIsdakotaEnum); 128 iomodel->FindConstant(&isgroundingline,TransientIsgroundinglineEnum); 129 iomodel->FindConstant(&ismovingfront,TransientIsmovingfrontEnum); 130 130 131 131 /*Finite element type*/ -
issm/trunk-jpl/src/c/analyses/SmbAnalysis.cpp
r19644 r20645 34 34 35 35 /*Figure out smb model: */ 36 iomodel-> Constant(&smb_model,SmbEnum);36 iomodel->FindConstant(&smb_model,SmbEnum); 37 37 38 38 … … 61 61 break; 62 62 case SMBpddEnum: 63 iomodel-> Constant(&isdelta18o,SmbIsdelta18oEnum);64 iomodel-> Constant(&ismungsm,SmbIsmungsmEnum);63 iomodel->FindConstant(&isdelta18o,SmbIsdelta18oEnum); 64 iomodel->FindConstant(&ismungsm,SmbIsmungsmEnum); 65 65 iomodel->FetchDataToInput(elements,ThermalSpctemperatureEnum); 66 66 iomodel->FetchDataToInput(elements,SmbS0pEnum); … … 78 78 break; 79 79 case SMBd18opddEnum: 80 iomodel-> Constant(&ismungsm,SmbIsmungsmEnum);81 iomodel-> Constant(&isd18opd,SmbIsd18opdEnum);80 iomodel->FindConstant(&ismungsm,SmbIsmungsmEnum); 81 iomodel->FindConstant(&isd18opd,SmbIsd18opdEnum); 82 82 iomodel->FetchDataToInput(elements,ThermalSpctemperatureEnum); 83 83 iomodel->FetchDataToInput(elements,SmbS0pEnum); … … 127 127 parameters->AddObject(iomodel->CopyConstantObject(SmbEnum)); 128 128 129 iomodel-> Constant(&smb_model,SmbEnum);130 iomodel-> Constant(&interp,TimesteppingInterpForcingsEnum);129 iomodel->FindConstant(&smb_model,SmbEnum); 130 iomodel->FindConstant(&interp,TimesteppingInterpForcingsEnum); 131 131 132 132 switch(smb_model){ … … 162 162 parameters->AddObject(iomodel->CopyConstantObject(SmbRlapsEnum)); 163 163 parameters->AddObject(iomodel->CopyConstantObject(SmbRlapslgmEnum)); 164 iomodel-> Constant(&isdelta18o,SmbIsdelta18oEnum);165 iomodel-> Constant(&ismungsm,SmbIsmungsmEnum);164 iomodel->FindConstant(&isdelta18o,SmbIsdelta18oEnum); 165 iomodel->FindConstant(&ismungsm,SmbIsmungsmEnum); 166 166 167 167 if(ismungsm){ … … 194 194 parameters->AddObject(iomodel->CopyConstantObject(SmbRlapsEnum)); 195 195 parameters->AddObject(iomodel->CopyConstantObject(SmbRlapslgmEnum)); 196 iomodel-> Constant(&ismungsm,SmbIsmungsmEnum);197 iomodel-> Constant(&isd18opd,SmbIsd18opdEnum);196 iomodel->FindConstant(&ismungsm,SmbIsmungsmEnum); 197 iomodel->FindConstant(&isd18opd,SmbIsd18opdEnum); 198 198 if(isd18opd){ 199 199 iomodel->FetchData(&temp,&N,&M,SmbDelta18oEnum); _assert_(N==2); -
issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
r20644 r20645 41 41 42 42 /*Fetch parameters: */ 43 iomodel-> Constant(&g,ConstantsGEnum);44 iomodel-> Constant(&rho_ice,MaterialsRhoIceEnum);45 iomodel-> Constant(&FSreconditioning,StressbalanceFSreconditioningEnum);46 iomodel-> Constant(&isSIA,FlowequationIsSIAEnum);47 iomodel-> Constant(&isSSA,FlowequationIsSSAEnum);48 iomodel-> Constant(&isL1L2,FlowequationIsL1L2Enum);49 iomodel-> Constant(&isHO,FlowequationIsHOEnum);50 iomodel-> Constant(&isFS,FlowequationIsFSEnum);43 iomodel->FindConstant(&g,ConstantsGEnum); 44 iomodel->FindConstant(&rho_ice,MaterialsRhoIceEnum); 45 iomodel->FindConstant(&FSreconditioning,StressbalanceFSreconditioningEnum); 46 iomodel->FindConstant(&isSIA,FlowequationIsSIAEnum); 47 iomodel->FindConstant(&isSSA,FlowequationIsSSAEnum); 48 iomodel->FindConstant(&isL1L2,FlowequationIsL1L2Enum); 49 iomodel->FindConstant(&isHO,FlowequationIsHOEnum); 50 iomodel->FindConstant(&isFS,FlowequationIsFSEnum); 51 51 52 52 /*Now, is the flag macayaealHO on? otherwise, do nothing: */ … … 63 63 64 64 /*Get finite element type*/ 65 if(isSSA) iomodel-> Constant(&finiteelement,FlowequationFeSSAEnum);65 if(isSSA) iomodel->FindConstant(&finiteelement,FlowequationFeSSAEnum); 66 66 else if(isL1L2) finiteelement = P1Enum; 67 else if(isHO) iomodel-> Constant(&finiteelement,FlowequationFeHOEnum);68 else if(isFS){ iomodel-> Constant(&finiteelement,FlowequationFeFSEnum);67 else if(isHO) iomodel->FindConstant(&finiteelement,FlowequationFeHOEnum); 68 else if(isFS){ iomodel->FindConstant(&finiteelement,FlowequationFeFSEnum); 69 69 /*Deduce velocity interpolation from finite element*/ 70 70 switch(finiteelement){ … … 451 451 452 452 /*Fetch parameters: */ 453 iomodel-> Constant(&isL1L2,FlowequationIsL1L2Enum);454 iomodel-> Constant(&isFS,FlowequationIsFSEnum);455 iomodel-> Constant(&isSSA,FlowequationIsSSAEnum);456 iomodel-> Constant(&isHO,FlowequationIsHOEnum);457 iomodel-> Constant(&numrifts,RiftsNumriftsEnum);453 iomodel->FindConstant(&isL1L2,FlowequationIsL1L2Enum); 454 iomodel->FindConstant(&isFS,FlowequationIsFSEnum); 455 iomodel->FindConstant(&isSSA,FlowequationIsSSAEnum); 456 iomodel->FindConstant(&isHO,FlowequationIsHOEnum); 457 iomodel->FindConstant(&numrifts,RiftsNumriftsEnum); 458 458 459 459 /*Now, is the flag macayaealHO on? otherwise, do nothing: */ … … 507 507 508 508 /*Fetch parameters: */ 509 iomodel-> Constant(&isSSA,FlowequationIsSSAEnum);510 iomodel-> Constant(&isL1L2,FlowequationIsL1L2Enum);511 iomodel-> Constant(&isHO,FlowequationIsHOEnum);512 iomodel-> Constant(&isFS,FlowequationIsFSEnum);509 iomodel->FindConstant(&isSSA,FlowequationIsSSAEnum); 510 iomodel->FindConstant(&isL1L2,FlowequationIsL1L2Enum); 511 iomodel->FindConstant(&isHO,FlowequationIsHOEnum); 512 iomodel->FindConstant(&isFS,FlowequationIsFSEnum); 513 513 514 514 /*Now, check that we have non SIA elements */ … … 527 527 if(isSSA){ 528 528 approximation=SSAApproximationEnum; 529 iomodel-> Constant(&finiteelement,FlowequationFeSSAEnum);529 iomodel->FindConstant(&finiteelement,FlowequationFeSSAEnum); 530 530 } 531 531 else if(isL1L2){ … … 535 535 else if(isHO){ 536 536 approximation = HOApproximationEnum; 537 iomodel-> Constant(&finiteelement,FlowequationFeHOEnum);537 iomodel->FindConstant(&finiteelement,FlowequationFeHOEnum); 538 538 } 539 539 else if(isFS){ 540 540 approximation = FSApproximationEnum; 541 iomodel-> Constant(&finiteelement,FlowequationFeFSEnum);541 iomodel->FindConstant(&finiteelement,FlowequationFeFSEnum); 542 542 } 543 543 iomodel->FetchData(3,FlowequationBorderSSAEnum,FlowequationVertexEquationEnum,StressbalanceReferentialEnum); … … 681 681 682 682 /*Fetch constants needed: */ 683 iomodel-> Constant(&isSSA,FlowequationIsSSAEnum);684 iomodel-> Constant(&isL1L2,FlowequationIsL1L2Enum);685 iomodel-> Constant(&isHO,FlowequationIsHOEnum);686 iomodel-> Constant(&isFS,FlowequationIsFSEnum);687 iomodel-> Constant(&control_analysis,InversionIscontrolEnum);688 iomodel-> Constant(&dakota_analysis,QmuIsdakotaEnum);689 iomodel-> Constant(&materials_type,MaterialsEnum);690 iomodel-> Constant(&ismovingfront,TransientIsmovingfrontEnum);691 iomodel-> Constant(&frictionlaw,FrictionLawEnum);683 iomodel->FindConstant(&isSSA,FlowequationIsSSAEnum); 684 iomodel->FindConstant(&isL1L2,FlowequationIsL1L2Enum); 685 iomodel->FindConstant(&isHO,FlowequationIsHOEnum); 686 iomodel->FindConstant(&isFS,FlowequationIsFSEnum); 687 iomodel->FindConstant(&control_analysis,InversionIscontrolEnum); 688 iomodel->FindConstant(&dakota_analysis,QmuIsdakotaEnum); 689 iomodel->FindConstant(&materials_type,MaterialsEnum); 690 iomodel->FindConstant(&ismovingfront,TransientIsmovingfrontEnum); 691 iomodel->FindConstant(&frictionlaw,FrictionLawEnum); 692 692 693 693 /*return if no processing required*/ … … 706 706 /*Get finite element type*/ 707 707 if(!iscoupling){ 708 if(isSSA) iomodel-> Constant(&finiteelement,FlowequationFeSSAEnum);708 if(isSSA) iomodel->FindConstant(&finiteelement,FlowequationFeSSAEnum); 709 709 else if(isL1L2) finiteelement = P1Enum; 710 else if(isHO) iomodel-> Constant(&finiteelement,FlowequationFeHOEnum);711 else if(isFS) iomodel-> Constant(&finiteelement,FlowequationFeFSEnum);710 else if(isHO) iomodel->FindConstant(&finiteelement,FlowequationFeHOEnum); 711 else if(isFS) iomodel->FindConstant(&finiteelement,FlowequationFeFSEnum); 712 712 for(int i=0;i<iomodel->numberofelements;i++){ 713 713 finiteelement_list[i]=finiteelement; … … 773 773 } 774 774 /*LATH parameters*/ 775 iomodel-> Constant(&fe_FS,FlowequationFeFSEnum);775 iomodel->FindConstant(&fe_FS,FlowequationFeFSEnum); 776 776 if(fe_FS==LATaylorHoodEnum || fe_FS==LACrouzeixRaviartEnum){ 777 777 InputUpdateFromConstantx(elements,0.,SigmaNNEnum); … … 790 790 break; 791 791 case 3: 792 iomodel-> Constant(&FrictionCoupling,FrictionCouplingEnum);792 iomodel->FindConstant(&FrictionCoupling,FrictionCouplingEnum); 793 793 iomodel->FetchDataToInput(elements,FrictionCEnum); 794 794 iomodel->FetchDataToInput(elements,FrictionAsEnum); … … 862 862 863 863 /*XTH LATH parameters*/ 864 iomodel-> Constant(&fe_FS,FlowequationFeFSEnum);864 iomodel->FindConstant(&fe_FS,FlowequationFeFSEnum); 865 865 if(fe_FS==XTaylorHoodEnum || fe_FS==LATaylorHoodEnum || fe_FS==LACrouzeixRaviartEnum){ 866 866 parameters->AddObject(iomodel->CopyConstantObject(AugmentedLagrangianREnum)); … … 869 869 } 870 870 871 iomodel-> Constant(&materials_type,MaterialsEnum);871 iomodel->FindConstant(&materials_type,MaterialsEnum); 872 872 if(materials_type==MatdamageiceEnum){ 873 873 parameters->AddObject(iomodel->CopyConstantObject(DamageLawEnum)); … … 884 884 /*Deal with friction parameters*/ 885 885 int frictionlaw; 886 iomodel-> Constant(&frictionlaw,FrictionLawEnum);886 iomodel->FindConstant(&frictionlaw,FrictionLawEnum); 887 887 if(frictionlaw==4 || frictionlaw==6) parameters->AddObject(iomodel->CopyConstantObject(FrictionGammaEnum)); 888 888 if(frictionlaw==3) parameters->AddObject(iomodel->CopyConstantObject(FrictionCouplingEnum)); -
issm/trunk-jpl/src/c/analyses/StressbalanceSIAAnalysis.cpp
r20459 r20645 13 13 14 14 /*Fetch parameters: */ 15 iomodel-> Constant(&isSIA,FlowequationIsSIAEnum);16 iomodel-> Constant(&isSSA,FlowequationIsSSAEnum);17 iomodel-> Constant(&isL1L2,FlowequationIsL1L2Enum);18 iomodel-> Constant(&isHO,FlowequationIsHOEnum);19 iomodel-> Constant(&isFS,FlowequationIsFSEnum);15 iomodel->FindConstant(&isSIA,FlowequationIsSIAEnum); 16 iomodel->FindConstant(&isSSA,FlowequationIsSSAEnum); 17 iomodel->FindConstant(&isL1L2,FlowequationIsL1L2Enum); 18 iomodel->FindConstant(&isHO,FlowequationIsHOEnum); 19 iomodel->FindConstant(&isFS,FlowequationIsFSEnum); 20 20 21 21 /*Now, is the flag isSIA on? otherwise, do nothing: */ … … 83 83 84 84 /*Fetch parameters: */ 85 iomodel-> Constant(&isSIA,FlowequationIsSIAEnum);85 iomodel->FindConstant(&isSIA,FlowequationIsSIAEnum); 86 86 87 87 /*Now, is the flag isSIA on? otherwise, do nothing: */ … … 124 124 bool ismovingfront; 125 125 int frictionlaw; 126 iomodel-> Constant(&isSIA,FlowequationIsSIAEnum);127 iomodel-> Constant(&ismovingfront,TransientIsmovingfrontEnum);128 iomodel-> Constant(&frictionlaw,FrictionLawEnum);126 iomodel->FindConstant(&isSIA,FlowequationIsSIAEnum); 127 iomodel->FindConstant(&ismovingfront,TransientIsmovingfrontEnum); 128 iomodel->FindConstant(&frictionlaw,FrictionLawEnum); 129 129 130 130 /*Now, is the flag SIA on? otherwise, do nothing: */ -
issm/trunk-jpl/src/c/analyses/StressbalanceVerticalAnalysis.cpp
r20155 r20645 18 18 19 19 /*Fetch parameters: */ 20 iomodel-> Constant(&isSIA,FlowequationIsSIAEnum);21 iomodel-> Constant(&isSSA,FlowequationIsSSAEnum);22 iomodel-> Constant(&isL1L2,FlowequationIsL1L2Enum);23 iomodel-> Constant(&isHO,FlowequationIsHOEnum);24 iomodel-> Constant(&isFS,FlowequationIsFSEnum);20 iomodel->FindConstant(&isSIA,FlowequationIsSIAEnum); 21 iomodel->FindConstant(&isSSA,FlowequationIsSSAEnum); 22 iomodel->FindConstant(&isL1L2,FlowequationIsL1L2Enum); 23 iomodel->FindConstant(&isHO,FlowequationIsHOEnum); 24 iomodel->FindConstant(&isFS,FlowequationIsFSEnum); 25 25 26 26 /*Do we have coupling*/ -
issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp
r20459 r20645 24 24 else if (iomodel->solution_enum==TransientSolutionEnum){ 25 25 int smb_model; 26 iomodel-> Constant(&smb_model,SmbEnum);26 iomodel->FindConstant(&smb_model,SmbEnum); 27 27 if(smb_model==SMBpddEnum) isdynamic=true; 28 28 if(smb_model==SMBd18opddEnum) isdynamic=true; … … 90 90 91 91 bool dakota_analysis, ismovingfront; 92 iomodel-> Constant(&dakota_analysis,QmuIsdakotaEnum);93 iomodel-> Constant(&ismovingfront,TransientIsmovingfrontEnum);94 iomodel-> Constant(&frictionlaw,FrictionLawEnum);92 iomodel->FindConstant(&dakota_analysis,QmuIsdakotaEnum); 93 iomodel->FindConstant(&ismovingfront,TransientIsmovingfrontEnum); 94 iomodel->FindConstant(&frictionlaw,FrictionLawEnum); 95 95 96 96 iomodel->FetchDataToInput(elements,ThicknessEnum); … … 120 120 } 121 121 /*Basal forcings variables*/ 122 iomodel-> Constant(&basalforcing_model,BasalforcingsEnum);122 iomodel->FindConstant(&basalforcing_model,BasalforcingsEnum); 123 123 switch(basalforcing_model){ 124 124 case MantlePlumeGeothermalFluxEnum: … … 140 140 break; 141 141 case 3: 142 iomodel-> Constant(&FrictionCoupling,FrictionCouplingEnum);142 iomodel->FindConstant(&FrictionCoupling,FrictionCouplingEnum); 143 143 iomodel->FetchDataToInput(elements,FrictionCEnum); 144 144 iomodel->FetchDataToInput(elements,FrictionAsEnum); … … 192 192 /*Deal with friction parameters*/ 193 193 int frictionlaw; 194 iomodel-> Constant(&frictionlaw,FrictionLawEnum);194 iomodel->FindConstant(&frictionlaw,FrictionLawEnum); 195 195 if(frictionlaw==4 || frictionlaw==6) parameters->AddObject(iomodel->CopyConstantObject(FrictionGammaEnum)); 196 196 if(frictionlaw==3) parameters->AddObject(iomodel->CopyConstantObject(FrictionCouplingEnum)); -
issm/trunk-jpl/src/c/analyses/UzawaPressureAnalysis.cpp
r18930 r20645 17 17 int fe_FS; 18 18 19 iomodel-> Constant(&fe_FS,FlowequationFeFSEnum);19 iomodel->FindConstant(&fe_FS,FlowequationFeFSEnum); 20 20 if(fe_FS==LATaylorHoodEnum) finiteelement = P1Enum; 21 21 else if(fe_FS==LACrouzeixRaviartEnum) finiteelement = P1DGEnum; … … 34 34 int fe_FS; 35 35 36 iomodel-> Constant(&fe_FS,FlowequationFeFSEnum);36 iomodel->FindConstant(&fe_FS,FlowequationFeFSEnum); 37 37 if(fe_FS==LATaylorHoodEnum) finiteelement = P1Enum; 38 38 else if(fe_FS==LACrouzeixRaviartEnum) finiteelement = P1DGEnum; -
issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
r20644 r20645 1399 1399 1400 1400 /*Fetch parameters: */ 1401 iomodel-> Constant(&yts,ConstantsYtsEnum);1402 iomodel-> Constant(&control_analysis,InversionIscontrolEnum);1403 if(control_analysis) iomodel-> Constant(&num_control_type,InversionNumControlParametersEnum);1404 if(control_analysis) iomodel-> Constant(&num_responses,InversionNumCostFunctionsEnum);1401 iomodel->FindConstant(&yts,ConstantsYtsEnum); 1402 iomodel->FindConstant(&control_analysis,InversionIscontrolEnum); 1403 if(control_analysis) iomodel->FindConstant(&num_control_type,InversionNumControlParametersEnum); 1404 if(control_analysis) iomodel->FindConstant(&num_responses,InversionNumCostFunctionsEnum); 1405 1405 1406 1406 /*Recover vertices ids needed to initialize inputs*/ … … 2696 2696 2697 2697 /*Fetch parameters: */ 2698 iomodel-> Constant(&yts,ConstantsYtsEnum);2699 iomodel-> Constant(&dakota_analysis,QmuIsdakotaEnum);2698 iomodel->FindConstant(&yts,ConstantsYtsEnum); 2699 iomodel->FindConstant(&dakota_analysis,QmuIsdakotaEnum); 2700 2700 2701 2701 /*Checks if debuging*/ -
issm/trunk-jpl/src/c/classes/Elements/Tetra.cpp
r20608 r20645 370 370 371 371 /*Fetch parameters: */ 372 iomodel-> Constant(&yts,ConstantsYtsEnum);373 iomodel-> Constant(&control_analysis,InversionIscontrolEnum);374 if(control_analysis) iomodel-> Constant(&num_control_type,InversionNumControlParametersEnum);375 if(control_analysis) iomodel-> Constant(&num_responses,InversionNumCostFunctionsEnum);372 iomodel->FindConstant(&yts,ConstantsYtsEnum); 373 iomodel->FindConstant(&control_analysis,InversionIscontrolEnum); 374 if(control_analysis) iomodel->FindConstant(&num_control_type,InversionNumControlParametersEnum); 375 if(control_analysis) iomodel->FindConstant(&num_responses,InversionNumCostFunctionsEnum); 376 376 377 377 /*Recover vertices ids needed to initialize inputs*/ … … 874 874 875 875 /*Fetch parameters: */ 876 iomodel-> Constant(&yts,ConstantsYtsEnum);877 iomodel-> Constant(&dakota_analysis,QmuIsdakotaEnum);878 iomodel-> Constant(&isFS,FlowequationIsFSEnum);876 iomodel->FindConstant(&yts,ConstantsYtsEnum); 877 iomodel->FindConstant(&dakota_analysis,QmuIsdakotaEnum); 878 iomodel->FindConstant(&isFS,FlowequationIsFSEnum); 879 879 880 880 /*Checks if debuging*/ -
issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
r20608 r20645 1811 1811 1812 1812 /*Get parameters: */ 1813 iomodel-> Constant(&yts,ConstantsYtsEnum);1814 iomodel-> Constant(&control_analysis,InversionIscontrolEnum);1815 if(control_analysis) iomodel-> Constant(&num_control_type,InversionNumControlParametersEnum);1816 if(control_analysis) iomodel-> Constant(&num_responses,InversionNumCostFunctionsEnum);1813 iomodel->FindConstant(&yts,ConstantsYtsEnum); 1814 iomodel->FindConstant(&control_analysis,InversionIscontrolEnum); 1815 if(control_analysis) iomodel->FindConstant(&num_control_type,InversionNumControlParametersEnum); 1816 if(control_analysis) iomodel->FindConstant(&num_responses,InversionNumCostFunctionsEnum); 1817 1817 1818 1818 /*Recover vertices ids needed to initialize inputs*/ -
issm/trunk-jpl/src/c/classes/FemModel.cpp
r20637 r20645 438 438 bool isSIA,isFS; 439 439 int fe_FS; 440 iomodel-> Constant(&fe_FS,FlowequationFeFSEnum);441 iomodel-> Constant(&isSIA,FlowequationIsSIAEnum);442 iomodel-> Constant(&isFS,FlowequationIsFSEnum);440 iomodel->FindConstant(&fe_FS,FlowequationFeFSEnum); 441 iomodel->FindConstant(&isSIA,FlowequationIsSIAEnum); 442 iomodel->FindConstant(&isFS,FlowequationIsFSEnum); 443 443 analyses_temp[numanalyses++]=StressbalanceAnalysisEnum; 444 444 analyses_temp[numanalyses++]=StressbalanceVerticalAnalysisEnum; … … 457 457 case SteadystateSolutionEnum:{ 458 458 bool isSIA,isenthalpy; 459 iomodel-> Constant(&isSIA,FlowequationIsSIAEnum);460 iomodel-> Constant(&isenthalpy,ThermalIsenthalpyEnum);459 iomodel->FindConstant(&isSIA,FlowequationIsSIAEnum); 460 iomodel->FindConstant(&isenthalpy,ThermalIsenthalpyEnum); 461 461 analyses_temp[numanalyses++]=StressbalanceAnalysisEnum; 462 462 analyses_temp[numanalyses++]=StressbalanceVerticalAnalysisEnum; … … 477 477 case ThermalSolutionEnum:{ 478 478 bool isenthalpy; 479 iomodel-> Constant(&isenthalpy,ThermalIsenthalpyEnum);479 iomodel->FindConstant(&isenthalpy,ThermalIsenthalpyEnum); 480 480 if(isenthalpy){ 481 481 analyses_temp[numanalyses++]=EnthalpyAnalysisEnum; … … 546 546 case TransientSolutionEnum:{ 547 547 bool isSIA,isFS,isthermal,isenthalpy,ismasstransport,isgroundingline,isstressbalance,ismovingfront,ishydrology,isdamage,issmb,isslr,isgia; 548 iomodel-> Constant(&isSIA,FlowequationIsSIAEnum);549 iomodel-> Constant(&isFS,FlowequationIsFSEnum);550 iomodel-> Constant(&isthermal,TransientIsthermalEnum);551 iomodel-> Constant(&isenthalpy,ThermalIsenthalpyEnum);552 iomodel-> Constant(&ismovingfront,TransientIsmovingfrontEnum);553 iomodel-> Constant(&ismasstransport,TransientIsmasstransportEnum);554 iomodel-> Constant(&isstressbalance,TransientIsstressbalanceEnum);555 iomodel-> Constant(&isgroundingline,TransientIsgroundinglineEnum);556 iomodel-> Constant(&isdamage,TransientIsdamageevolutionEnum);557 iomodel-> Constant(&ishydrology,TransientIshydrologyEnum);558 iomodel-> Constant(&issmb,TransientIssmbEnum);559 iomodel-> Constant(&isslr,TransientIsslrEnum);560 iomodel-> Constant(&isgia,TransientIsgiaEnum);548 iomodel->FindConstant(&isSIA,FlowequationIsSIAEnum); 549 iomodel->FindConstant(&isFS,FlowequationIsFSEnum); 550 iomodel->FindConstant(&isthermal,TransientIsthermalEnum); 551 iomodel->FindConstant(&isenthalpy,ThermalIsenthalpyEnum); 552 iomodel->FindConstant(&ismovingfront,TransientIsmovingfrontEnum); 553 iomodel->FindConstant(&ismasstransport,TransientIsmasstransportEnum); 554 iomodel->FindConstant(&isstressbalance,TransientIsstressbalanceEnum); 555 iomodel->FindConstant(&isgroundingline,TransientIsgroundinglineEnum); 556 iomodel->FindConstant(&isdamage,TransientIsdamageevolutionEnum); 557 iomodel->FindConstant(&ishydrology,TransientIshydrologyEnum); 558 iomodel->FindConstant(&issmb,TransientIssmbEnum); 559 iomodel->FindConstant(&isslr,TransientIsslrEnum); 560 iomodel->FindConstant(&isgia,TransientIsgiaEnum); 561 561 if(isstressbalance){ 562 562 int fe_FS; 563 iomodel-> Constant(&fe_FS,FlowequationFeFSEnum);563 iomodel->FindConstant(&fe_FS,FlowequationFeFSEnum); 564 564 analyses_temp[numanalyses++]=StressbalanceAnalysisEnum; 565 565 analyses_temp[numanalyses++]=StressbalanceVerticalAnalysisEnum; -
issm/trunk-jpl/src/c/classes/IoModel.cpp
r20635 r20645 244 244 } 245 245 /*}}}*/ 246 void IoModel:: Constant(bool* poutput,int constant_enum){/*{{{*/246 void IoModel::FindConstant(bool* poutput,int constant_enum){/*{{{*/ 247 247 248 248 _assert_(constant_enum>=0); … … 252 252 } 253 253 /*}}}*/ 254 void IoModel:: Constant(int* poutput,int constant_enum){/*{{{*/254 void IoModel::FindConstant(int* poutput,int constant_enum){/*{{{*/ 255 255 256 256 _assert_(constant_enum>=0); … … 260 260 } 261 261 /*}}}*/ 262 void IoModel:: Constant(IssmDouble* poutput,int constant_enum){/*{{{*/262 void IoModel::FindConstant(IssmDouble* poutput,int constant_enum){/*{{{*/ 263 263 264 264 _assert_(constant_enum>=0); … … 268 268 } 269 269 /*}}}*/ 270 void IoModel:: Constant(char** poutput,int constant_enum){/*{{{*/270 void IoModel::FindConstant(char** poutput,int constant_enum){/*{{{*/ 271 271 272 272 _assert_(constant_enum>=0); -
issm/trunk-jpl/src/c/classes/IoModel.h
r19138 r20645 66 66 /*Input/Output*/ 67 67 void CheckEnumSync(void); 68 void Constant(bool *poutput,int constant_enum);69 void Constant(int *poutput,int constant_enum);70 void Constant(IssmDouble *poutput,int constant_enum);71 void Constant(char **poutput,int constant_enum);68 void FindConstant(bool *poutput,int constant_enum); 69 void FindConstant(int *poutput,int constant_enum); 70 void FindConstant(IssmDouble *poutput,int constant_enum); 71 void FindConstant(char **poutput,int constant_enum); 72 72 Param *CopyConstantObject(int constant_enum); 73 73 IssmDouble *Data(int dataenum); -
issm/trunk-jpl/src/c/classes/Loads/Riftfront.cpp
r20462 r20645 46 46 47 47 /*Fetch parameters: */ 48 iomodel-> Constant(&penalty_lock,StressbalanceRiftPenaltyLockEnum);48 iomodel->FindConstant(&penalty_lock,StressbalanceRiftPenaltyLockEnum); 49 49 50 50 /*Ok, retrieve all the data needed to add a penalty between the two nodes: */ -
issm/trunk-jpl/src/c/classes/Materials/Matice.cpp
r20608 r20645 49 49 /*Other perporties*/ 50 50 int materialtype; 51 iomodel-> Constant(&materialtype,MaterialsEnum);51 iomodel->FindConstant(&materialtype,MaterialsEnum); 52 52 if(materialtype==MatdamageiceEnum) this->isdamaged = true; 53 53 else if(materialtype==MaticeEnum) this->isdamaged = false; -
issm/trunk-jpl/src/c/classes/Materials/Matpar.cpp
r20644 r20645 73 73 bool isefficientlayer; 74 74 int hydrology_model,smb_model,materials_type; 75 iomodel-> Constant(&hydrology_model,HydrologyModelEnum);76 iomodel-> Constant(&smb_model,SmbEnum);77 iomodel-> Constant(&materials_type,MaterialsEnum);75 iomodel->FindConstant(&hydrology_model,HydrologyModelEnum); 76 iomodel->FindConstant(&smb_model,SmbEnum); 77 iomodel->FindConstant(&materials_type,MaterialsEnum); 78 78 79 79 this->mid = matpar_mid; … … 83 83 case MatdamageiceEnum: 84 84 case MatbaflEnum: 85 iomodel-> Constant(&this->rho_ice,MaterialsRhoIceEnum);86 iomodel-> Constant(&this->rho_water,MaterialsRhoSeawaterEnum);87 iomodel-> Constant(&this->rho_freshwater,MaterialsRhoFreshwaterEnum);88 iomodel-> Constant(&this->mu_water,MaterialsMuWaterEnum);89 iomodel-> Constant(&this->heatcapacity,MaterialsHeatcapacityEnum);90 iomodel-> Constant(&this->thermalconductivity,MaterialsThermalconductivityEnum);91 iomodel-> Constant(&this->temperateiceconductivity,MaterialsTemperateiceconductivityEnum);92 iomodel-> Constant(&this->latentheat,MaterialsLatentheatEnum);93 iomodel-> Constant(&this->beta,MaterialsBetaEnum);94 iomodel-> Constant(&this->meltingpoint,MaterialsMeltingpointEnum);95 iomodel-> Constant(&this->referencetemperature,ConstantsReferencetemperatureEnum);96 iomodel-> Constant(&this->mixed_layer_capacity,MaterialsMixedLayerCapacityEnum);97 iomodel-> Constant(&this->thermal_exchange_velocity,MaterialsThermalExchangeVelocityEnum);98 iomodel-> Constant(&this->g,ConstantsGEnum);85 iomodel->FindConstant(&this->rho_ice,MaterialsRhoIceEnum); 86 iomodel->FindConstant(&this->rho_water,MaterialsRhoSeawaterEnum); 87 iomodel->FindConstant(&this->rho_freshwater,MaterialsRhoFreshwaterEnum); 88 iomodel->FindConstant(&this->mu_water,MaterialsMuWaterEnum); 89 iomodel->FindConstant(&this->heatcapacity,MaterialsHeatcapacityEnum); 90 iomodel->FindConstant(&this->thermalconductivity,MaterialsThermalconductivityEnum); 91 iomodel->FindConstant(&this->temperateiceconductivity,MaterialsTemperateiceconductivityEnum); 92 iomodel->FindConstant(&this->latentheat,MaterialsLatentheatEnum); 93 iomodel->FindConstant(&this->beta,MaterialsBetaEnum); 94 iomodel->FindConstant(&this->meltingpoint,MaterialsMeltingpointEnum); 95 iomodel->FindConstant(&this->referencetemperature,ConstantsReferencetemperatureEnum); 96 iomodel->FindConstant(&this->mixed_layer_capacity,MaterialsMixedLayerCapacityEnum); 97 iomodel->FindConstant(&this->thermal_exchange_velocity,MaterialsThermalExchangeVelocityEnum); 98 iomodel->FindConstant(&this->g,ConstantsGEnum); 99 99 100 100 switch(smb_model){ … … 103 103 break; 104 104 case SMBgembEnum: 105 iomodel-> Constant(&this->albedo_ice,SmbAIceEnum);106 iomodel-> Constant(&this->albedo_snow,SmbASnowEnum);105 iomodel->FindConstant(&this->albedo_ice,SmbAIceEnum); 106 iomodel->FindConstant(&this->albedo_snow,SmbASnowEnum); 107 107 break; 108 108 case SMBpddEnum: 109 iomodel-> Constant(&this->desfac,SmbDesfacEnum);110 iomodel-> Constant(&this->rlaps,SmbRlapsEnum);111 iomodel-> Constant(&this->rlapslgm,SmbRlapslgmEnum);109 iomodel->FindConstant(&this->desfac,SmbDesfacEnum); 110 iomodel->FindConstant(&this->rlaps,SmbRlapsEnum); 111 iomodel->FindConstant(&this->rlapslgm,SmbRlapslgmEnum); 112 112 break; 113 113 case SMBd18opddEnum: 114 iomodel-> Constant(&this->desfac,SmbDesfacEnum);115 iomodel-> Constant(&this->rlaps,SmbRlapsEnum);116 iomodel-> Constant(&this->rlapslgm,SmbRlapslgmEnum);117 iomodel-> Constant(&this->dpermil,SmbDpermilEnum);114 iomodel->FindConstant(&this->desfac,SmbDesfacEnum); 115 iomodel->FindConstant(&this->rlaps,SmbRlapsEnum); 116 iomodel->FindConstant(&this->rlapslgm,SmbRlapslgmEnum); 117 iomodel->FindConstant(&this->dpermil,SmbDpermilEnum); 118 118 case SMBgradientsEnum: 119 119 break; … … 133 133 } 134 134 if(hydrology_model==HydrologydcEnum){ 135 iomodel-> Constant(&this->sediment_compressibility,HydrologydcSedimentCompressibilityEnum);136 iomodel-> Constant(&this->sediment_porosity,HydrologydcSedimentPorosityEnum);137 iomodel-> Constant(&this->sediment_thickness,HydrologydcSedimentThicknessEnum);138 iomodel-> Constant(&this->water_compressibility,HydrologydcWaterCompressibilityEnum);139 iomodel-> Constant(&isefficientlayer,HydrologydcIsefficientlayerEnum);135 iomodel->FindConstant(&this->sediment_compressibility,HydrologydcSedimentCompressibilityEnum); 136 iomodel->FindConstant(&this->sediment_porosity,HydrologydcSedimentPorosityEnum); 137 iomodel->FindConstant(&this->sediment_thickness,HydrologydcSedimentThicknessEnum); 138 iomodel->FindConstant(&this->water_compressibility,HydrologydcWaterCompressibilityEnum); 139 iomodel->FindConstant(&isefficientlayer,HydrologydcIsefficientlayerEnum); 140 140 141 141 if(isefficientlayer){ 142 iomodel-> Constant(&this->epl_compressibility,HydrologydcEplCompressibilityEnum);143 iomodel-> Constant(&this->epl_porosity,HydrologydcEplPorosityEnum);144 iomodel-> Constant(&this->epl_init_thickness,HydrologydcEplInitialThicknessEnum);145 iomodel-> Constant(&this->epl_colapse_thickness,HydrologydcEplColapseThicknessEnum);146 iomodel-> Constant(&this->epl_max_thickness,HydrologydcEplMaxThicknessEnum);147 iomodel-> Constant(&this->epl_conductivity,HydrologydcEplConductivityEnum);142 iomodel->FindConstant(&this->epl_compressibility,HydrologydcEplCompressibilityEnum); 143 iomodel->FindConstant(&this->epl_porosity,HydrologydcEplPorosityEnum); 144 iomodel->FindConstant(&this->epl_init_thickness,HydrologydcEplInitialThicknessEnum); 145 iomodel->FindConstant(&this->epl_colapse_thickness,HydrologydcEplColapseThicknessEnum); 146 iomodel->FindConstant(&this->epl_max_thickness,HydrologydcEplMaxThicknessEnum); 147 iomodel->FindConstant(&this->epl_conductivity,HydrologydcEplConductivityEnum); 148 148 } 149 149 } … … 159 159 160 160 /*gia: */ 161 iomodel-> Constant(&this->lithosphere_shear_modulus,MaterialsLithosphereShearModulusEnum);162 iomodel-> Constant(&this->lithosphere_density,MaterialsLithosphereDensityEnum);163 iomodel-> Constant(&this->mantle_shear_modulus,MaterialsMantleShearModulusEnum);164 iomodel-> Constant(&this->mantle_density,MaterialsMantleDensityEnum);161 iomodel->FindConstant(&this->lithosphere_shear_modulus,MaterialsLithosphereShearModulusEnum); 162 iomodel->FindConstant(&this->lithosphere_density,MaterialsLithosphereDensityEnum); 163 iomodel->FindConstant(&this->mantle_shear_modulus,MaterialsMantleShearModulusEnum); 164 iomodel->FindConstant(&this->mantle_density,MaterialsMantleDensityEnum); 165 165 166 166 /*slr:*/ 167 iomodel-> Constant(&this->earth_density,MaterialsEarthDensityEnum);167 iomodel->FindConstant(&this->earth_density,MaterialsEarthDensityEnum); 168 168 169 169 break; -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Autodiff/CreateParametersAutodiff.cpp
r19497 r20645 28 28 29 29 /*retrieve some parameters: */ 30 iomodel-> Constant(&autodiff_analysis,AutodiffIsautodiffEnum);30 iomodel->FindConstant(&autodiff_analysis,AutodiffIsautodiffEnum); 31 31 32 32 #ifdef _HAVE_ADOLC_ … … 83 83 /*}}}*/ 84 84 /*retrieve driver: {{{*/ 85 iomodel-> Constant(&autodiff_driver,AutodiffDriverEnum);85 iomodel->FindConstant(&autodiff_driver,AutodiffDriverEnum); 86 86 parameters->AddObject(iomodel->CopyConstantObject(AutodiffDriverEnum)); 87 87 … … 101 101 /*}}}*/ 102 102 /*Deal with dependents first: {{{*/ 103 iomodel-> Constant(&num_dependent_objects,AutodiffNumDependentObjectsEnum);103 iomodel->FindConstant(&num_dependent_objects,AutodiffNumDependentObjectsEnum); 104 104 dependent_objects=new DataSet(); 105 105 num_dep=0; -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Control/CreateParametersControl.cpp
r18996 r20645 23 23 24 24 /*retrieve some parameters: */ 25 iomodel-> Constant(&control_analysis,InversionIscontrolEnum);26 iomodel-> Constant(&inversiontype,InversionTypeEnum);25 iomodel->FindConstant(&control_analysis,InversionIscontrolEnum); 26 iomodel->FindConstant(&inversiontype,InversionTypeEnum); 27 27 28 28 if(control_analysis){ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Control/UpdateElementsAndMaterialsControl.cpp
r20445 r20645 19 19 20 20 /*Fetch parameters: */ 21 iomodel-> Constant(&control_analysis,InversionIscontrolEnum);22 if(control_analysis) iomodel-> Constant(&num_control_type,InversionNumControlParametersEnum);21 iomodel->FindConstant(&control_analysis,InversionIscontrolEnum); 22 if(control_analysis) iomodel->FindConstant(&num_control_type,InversionNumControlParametersEnum); 23 23 24 24 /*Now, return if no control*/ … … 28 28 29 29 /*Fetch Observations */ 30 iomodel-> Constant(&num_cost_functions,InversionNumCostFunctionsEnum);31 iomodel-> Constant(&domaintype,DomainTypeEnum);30 iomodel->FindConstant(&num_cost_functions,InversionNumCostFunctionsEnum); 31 iomodel->FindConstant(&domaintype,DomainTypeEnum); 32 32 for(int i=0;i<num_cost_functions;i++){ 33 33 cost_function= reCast<int,IssmDouble>(iomodel->Data(InversionCostFunctionsEnum)[i]); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp
r20644 r20645 17 17 18 18 /*Fetch parameters: */ 19 iomodel-> Constant(&control_analysis,InversionIscontrolEnum);20 iomodel-> Constant(&dakota_analysis,QmuIsdakotaEnum);21 iomodel-> Constant(&materials_type,MaterialsEnum);19 iomodel->FindConstant(&control_analysis,InversionIscontrolEnum); 20 iomodel->FindConstant(&dakota_analysis,QmuIsdakotaEnum); 21 iomodel->FindConstant(&materials_type,MaterialsEnum); 22 22 23 23 /*Did we already create the elements? : */ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp
r20635 r20645 99 99 /*Basal forcing parameters*/ 100 100 parameters->AddObject(iomodel->CopyConstantObject(BasalforcingsEnum)); 101 iomodel-> Constant(&basalforcing_model,BasalforcingsEnum);101 iomodel->FindConstant(&basalforcing_model,BasalforcingsEnum); 102 102 switch(basalforcing_model){ 103 103 case FloatingMeltRateEnum: … … 135 135 parameters->AddObject(new IntParam(SolutionTypeEnum,solution_type)); 136 136 137 iomodel-> Constant(&time,TimesteppingStartTimeEnum);137 iomodel->FindConstant(&time,TimesteppingStartTimeEnum); 138 138 parameters->AddObject(new DoubleParam(TimeEnum,time)); 139 139 parameters->AddObject(new IntParam(StepEnum,0)); … … 153 153 iomodel->DeleteData(&requestedoutputs,numoutputs,SteadystateRequestedOutputsEnum); 154 154 155 iomodel-> Constant(&materialtype,MaterialsEnum);155 iomodel->FindConstant(&materialtype,MaterialsEnum); 156 156 if(materialtype==MatdamageiceEnum){ 157 157 iomodel->FetchData(&requestedoutputs,&numoutputs,DamageEvolutionRequestedOutputsEnum); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Dakota/CreateParametersDakota.cpp
r16543 r20645 36 36 37 37 /*recover parameters: */ 38 iomodel-> Constant(&dakota_analysis,QmuIsdakotaEnum);38 iomodel->FindConstant(&dakota_analysis,QmuIsdakotaEnum); 39 39 40 40 if(dakota_analysis){ 41 41 42 iomodel-> Constant(&name,MiscellaneousNameEnum);43 iomodel-> Constant(&numberofresponses,QmuNumberofresponsesEnum);42 iomodel->FindConstant(&name,MiscellaneousNameEnum); 43 iomodel->FindConstant(&numberofresponses,QmuNumberofresponsesEnum); 44 44 45 45 /*name of qmu input, error and output files*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Dakota/UpdateElementsAndMaterialsDakota.cpp
r16236 r20645 13 13 /*recover parameters: */ 14 14 bool dakota_analysis; 15 iomodel-> Constant(&dakota_analysis,QmuIsdakotaEnum);15 iomodel->FindConstant(&dakota_analysis,QmuIsdakotaEnum); 16 16 17 17 if(dakota_analysis){ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/ElementsAndVerticesPartitioning.cpp
r17692 r20645 44 44 45 45 /*Fetch parameters: */ 46 iomodel-> Constant(&numrifts,RiftsNumriftsEnum);46 iomodel->FindConstant(&numrifts,RiftsNumriftsEnum); 47 47 48 48 /*First, check that partitioning has not yet been carryed out. Just check whether my_elements pointers is not already assigned a value: */ … … 67 67 elements_width=6; 68 68 iomodel->FetchData(&elements2d,NULL,NULL,MeshElements2dEnum); 69 iomodel-> Constant(&numberofelements2d,MeshNumberofelements2dEnum);70 iomodel-> Constant(&numberofvertices2d,MeshNumberofvertices2dEnum);71 iomodel-> Constant(&numlayers,MeshNumberoflayersEnum);69 iomodel->FindConstant(&numberofelements2d,MeshNumberofelements2dEnum); 70 iomodel->FindConstant(&numberofvertices2d,MeshNumberofvertices2dEnum); 71 iomodel->FindConstant(&numlayers,MeshNumberoflayersEnum); 72 72 break; 73 73 default: -
issm/trunk-jpl/src/c/modules/ModelProcessorx/ModelProcessorx.cpp
r20635 r20645 27 27 28 28 /*Fetch parameters: */ 29 iomodel-> Constant(&verbose,VerboseEnum);29 iomodel->FindConstant(&verbose,VerboseEnum); 30 30 SetVerbosityLevel(verbose); 31 31 -
issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp
r19527 r20645 203 203 /* Get constants */ 204 204 femmodel->parameters->FindParam(&yts,ConstantsYtsEnum); 205 /*iomodel-> Constant(&yts,ConstantsYtsEnum);*/205 /*iomodel->FindConstant(&yts,ConstantsYtsEnum);*/ 206 206 /*this->parameters->FindParam(&yts,ConstantsYtsEnum);*/ 207 207 /*Mathieu original*/
Note:
See TracChangeset
for help on using the changeset viewer.