Changeset 25118
- Timestamp:
- 06/22/20 22:29:05 (5 years ago)
- Location:
- issm/trunk-jpl/src
- Files:
-
- 6 added
- 24 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/Balancethickness2Analysis.cpp
r24335 r25118 32 32 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.surface",SurfaceEnum); 33 33 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.base",BaseEnum); 34 iomodel->FetchDataToInput(inputs2,elements,"md.s lr.sealevel",SealevelEnum,0);34 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.sealevel",SealevelEnum,0); 35 35 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); 36 36 iomodel->FetchDataToInput(inputs2,elements,"md.initialization.vx",VxEnum); -
issm/trunk-jpl/src/c/analyses/BalancethicknessAnalysis.cpp
r24936 r25118 99 99 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.surface",SurfaceEnum); 100 100 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.base",BaseEnum); 101 iomodel->FetchDataToInput(inputs2,elements,"md.s lr.sealevel",SealevelEnum,0);101 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.sealevel",SealevelEnum,0); 102 102 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); 103 103 iomodel->FetchDataToInput(inputs2,elements,"md.initialization.vx",VxEnum); -
issm/trunk-jpl/src/c/analyses/BalancevelocityAnalysis.cpp
r24673 r25118 41 41 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.surface",SurfaceEnum); 42 42 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.base",BaseEnum); 43 iomodel->FetchDataToInput(inputs2,elements,"md.s lr.sealevel",SealevelEnum,0);43 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.sealevel",SealevelEnum,0); 44 44 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); 45 45 iomodel->FetchDataToInput(inputs2,elements,"md.initialization.vx",VxEnum); -
issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
r24939 r25118 136 136 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.thickness",ThicknessEnum); 137 137 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.surface",SurfaceEnum); 138 iomodel->FetchDataToInput(inputs2,elements,"md.s lr.sealevel",SealevelEnum,0);138 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.sealevel",SealevelEnum,0); 139 139 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.base",BaseEnum); 140 140 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); -
issm/trunk-jpl/src/c/analyses/EsaAnalysis.cpp
r24959 r25118 61 61 /*some constant parameters: */ 62 62 parameters->AddObject(iomodel->CopyConstantObject("md.esa.hemisphere",EsaHemisphereEnum)); 63 parameters->AddObject(iomodel->CopyConstantObject("md.s lr.rigid",SealevelriseRigidEnum));64 parameters->AddObject(iomodel->CopyConstantObject("md.s lr.horiz",SealevelriseHorizEnum));65 parameters->AddObject(iomodel->CopyConstantObject("md.s lr.elastic",SealevelriseElasticEnum));66 parameters->AddObject(iomodel->CopyConstantObject("md.s lr.rotation",SealevelriseRotationEnum));63 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.settings.rigid",SolidearthSettingsRigidEnum)); 64 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.settings.horiz",SolidearthSettingsHorizEnum)); 65 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.settings.elastic",SolidearthSettingsElasticEnum)); 66 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.settings.rotation",SolidearthSettingsRotationEnum)); 67 67 68 68 /*deal with planet radius and area: */ 69 parameters->AddObject(iomodel->CopyConstantObject("md.s lr.planetradius",SealevelPlanetRadiusEnum));70 iomodel->FetchData(&planetradius,"md.s lr.planetradius");69 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.planetradius",SolidearthPlanetRadiusEnum)); 70 iomodel->FetchData(&planetradius,"md.solidearth.planetradius"); 71 71 planetarea=4*PI*planetradius*planetradius; 72 parameters->AddObject(new DoubleParam(S ealevelPlanetAreaEnum,planetarea));72 parameters->AddObject(new DoubleParam(SolidearthPlanetAreaEnum,planetarea)); 73 73 74 74 /*love numbers: */ -
issm/trunk-jpl/src/c/analyses/FreeSurfaceBaseAnalysis.cpp
r24933 r25118 74 74 75 75 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.surface",SurfaceEnum); 76 iomodel->FetchDataToInput(inputs2,elements,"md.s lr.sealevel",SealevelEnum,0);76 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.sealevel",SealevelEnum,0); 77 77 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); 78 78 iomodel->FetchDataToInput(inputs2,elements,"md.basalforcings.groundedice_melting_rate",BasalforcingsGroundediceMeltingRateEnum,0.); -
issm/trunk-jpl/src/c/analyses/FreeSurfaceTopAnalysis.cpp
r24933 r25118 77 77 78 78 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.surface",SurfaceEnum); 79 iomodel->FetchDataToInput(inputs2,elements,"md.s lr.sealevel",SealevelEnum,0);79 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.sealevel",SealevelEnum,0); 80 80 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); 81 81 iomodel->FetchDataToInput(inputs2,elements,"md.initialization.vx",VxEnum); -
issm/trunk-jpl/src/c/analyses/HydrologyShreveAnalysis.cpp
r24933 r25118 58 58 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.surface",SurfaceEnum); 59 59 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.base",BaseEnum); 60 iomodel->FetchDataToInput(inputs2,elements,"md.s lr.sealevel",SealevelEnum,0);60 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.sealevel",SealevelEnum,0); 61 61 if(iomodel->domaintype!=Domain2DhorizontalEnum){ 62 62 iomodel->FetchDataToInput(inputs2,elements,"md.mesh.vertexonbase",MeshVertexonbaseEnum); -
issm/trunk-jpl/src/c/analyses/L2ProjectionBaseAnalysis.cpp
r24335 r25118 42 42 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.surface",SurfaceEnum); 43 43 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.base",BaseEnum); 44 iomodel->FetchDataToInput(inputs2,elements,"md.s lr.sealevel",SealevelEnum,0);44 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.sealevel",SealevelEnum,0); 45 45 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); 46 46 if(iomodel->domaintype!=Domain2DhorizontalEnum & iomodel->domaintype!=Domain3DsurfaceEnum){ -
issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp
r24861 r25118 148 148 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.surface",SurfaceEnum); 149 149 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.base",BaseEnum); 150 iomodel->FetchDataToInput(inputs2,elements,"md.s lr.sealevel",SealevelEnum,0);150 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.sealevel",SealevelEnum,0); 151 151 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); 152 152 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ocean_levelset",MaskOceanLevelsetEnum); … … 907 907 element->AddBasalInput2(BaseEnum,newbase,P1Enum); 908 908 element->AddBasalInput2(SealevelriseCumDeltathicknessEnum,cumdeltathickness,P1Enum); 909 element->AddBasalInput2(S ealevelriseDeltathicknessEnum,deltathickness,P1Enum);909 element->AddBasalInput2(SurfaceloadIceThicknessChangeEnum,deltathickness,P1Enum); 910 910 911 911 /*Free ressources:*/ -
issm/trunk-jpl/src/c/analyses/MeltingAnalysis.cpp
r24335 r25118 56 56 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.surface",SurfaceEnum); 57 57 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.base",BaseEnum); 58 iomodel->FetchDataToInput(inputs2,elements,"md.s lr.sealevel",SealevelEnum,0);58 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.sealevel",SealevelEnum,0); 59 59 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); 60 60 if(iomodel->domaintype!=Domain2DhorizontalEnum){ -
issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp
r25066 r25118 38 38 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ocean_levelset",MaskOceanLevelsetEnum); 39 39 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); 40 iomodel->FetchData(&geodetic,"md.slr.geodetic"); 41 iomodel->FetchDataToInput(inputs2,elements,"md.slr.deltathickness",SealevelriseDeltathicknessEnum); 42 iomodel->FetchDataToInput(inputs2,elements,"md.slr.spcthickness",SealevelriseSpcthicknessEnum); 43 iomodel->FetchDataToInput(inputs2,elements,"md.slr.sealevel",SealevelEnum,0); 40 iomodel->FetchData(&geodetic,"md.solidearth.settings.computesealevelchange"); 41 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.surfaceload.icethicknesschange",SurfaceloadIceThicknessChangeEnum); 42 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.sealevel",SealevelEnum,0); 44 43 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.bed",BedEnum); 45 iomodel->FetchDataToInput(inputs2,elements,"md.s lr.hydro_rate",SealevelriseHydroRateEnum);44 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.surfaceload.waterheightchange",SurfaceloadWaterHeightChangeEnum); 46 45 47 46 /*dynamic sea level: */ … … 278 277 279 278 int nl; 279 int lovemodel=1; 280 280 IssmDouble* love_h=NULL; 281 281 IssmDouble* love_k=NULL; … … 307 307 /*some constant parameters: */ 308 308 parameters->AddObject(iomodel->CopyConstantObject("md.dsl.model",DslModelEnum)); 309 parameters->AddObject(iomodel->CopyConstantObject("md.slr.geodetic_run_frequency",SealevelriseGeodeticRunFrequencyEnum)); 310 parameters->AddObject(iomodel->CopyConstantObject("md.slr.reltol",SealevelriseReltolEnum)); 311 parameters->AddObject(iomodel->CopyConstantObject("md.slr.abstol",SealevelriseAbstolEnum)); 312 parameters->AddObject(iomodel->CopyConstantObject("md.slr.maxiter",SealevelriseMaxiterEnum)); 313 parameters->AddObject(iomodel->CopyConstantObject("md.slr.rigid",SealevelriseRigidEnum)); 314 parameters->AddObject(iomodel->CopyConstantObject("md.slr.horiz",SealevelriseHorizEnum)); 315 parameters->AddObject(iomodel->CopyConstantObject("md.slr.elastic",SealevelriseElasticEnum)); 316 parameters->AddObject(iomodel->CopyConstantObject("md.slr.rotation",SealevelriseRotationEnum)); 317 parameters->AddObject(iomodel->CopyConstantObject("md.slr.tide_love_h",SealevelriseTidalLoveHEnum)); 318 parameters->AddObject(iomodel->CopyConstantObject("md.slr.tide_love_k",SealevelriseTidalLoveKEnum)); 319 parameters->AddObject(iomodel->CopyConstantObject("md.slr.fluid_love",SealevelriseFluidLoveEnum)); 320 parameters->AddObject(iomodel->CopyConstantObject("md.slr.equatorial_moi",SealevelriseEquatorialMoiEnum)); 321 parameters->AddObject(iomodel->CopyConstantObject("md.slr.polar_moi",SealevelrisePolarMoiEnum)); 322 parameters->AddObject(iomodel->CopyConstantObject("md.slr.angular_velocity",SealevelriseAngularVelocityEnum)); 323 parameters->AddObject(iomodel->CopyConstantObject("md.slr.ocean_area_scaling",SealevelriseOceanAreaScalingEnum)); 324 parameters->AddObject(iomodel->CopyConstantObject("md.slr.geodetic",SealevelriseGeodeticEnum)); 325 parameters->AddObject(iomodel->CopyConstantObject("md.slr.planetradius",SealevelPlanetRadiusEnum)); 309 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.settings.runfrequency",SolidearthSettingsRunFrequencyEnum)); 310 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.settings.reltol",SolidearthSettingsReltolEnum)); 311 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.settings.abstol",SolidearthSettingsAbstolEnum)); 312 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.settings.maxiter",SolidearthSettingsMaxiterEnum)); 313 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.settings.rigid",SolidearthSettingsRigidEnum)); 314 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.settings.horiz",SolidearthSettingsHorizEnum)); 315 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.settings.elastic",SolidearthSettingsElasticEnum)); 316 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.settings.rotation",SolidearthSettingsRotationEnum)); 317 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.rotational.equatorialmoi",RotationalEquatorialMoiEnum)); 318 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.rotational.polarmoi",RotationalPolarMoiEnum)); 319 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.rotational.angularvelocity",RotationalAngularVelocityEnum)); 320 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.settings.ocean_area_scaling",SolidearthSettingsOceanAreaScalingEnum)); 321 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.settings.computesealevelchange",SolidearthSettingsComputesealevelchangeEnum)); 322 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.planetradius",SolidearthPlanetRadiusEnum)); 326 323 327 324 /*compute planet area and plug into parameters:*/ 328 iomodel->FetchData(&planetradius,"md.s lr.planetradius");325 iomodel->FetchData(&planetradius,"md.solidearth.planetradius"); 329 326 planetarea=4*PI*planetradius*planetradius; 330 parameters->AddObject(new DoubleParam(S ealevelPlanetAreaEnum,planetarea));327 parameters->AddObject(new DoubleParam(SolidearthPlanetAreaEnum,planetarea)); 331 328 332 329 /*Deal with dsl multi-model ensembles: {{{*/ … … 348 345 /*Deal with elasticity {{{*/ 349 346 /*love numbers: */ 350 iomodel->FetchData(&love_h,&nl,NULL,"md.slr.love_h"); 351 iomodel->FetchData(&love_k,&nl,NULL,"md.slr.love_k"); 352 iomodel->FetchData(&love_l,&nl,NULL,"md.slr.love_l"); 347 iomodel->FetchData(&lovemodel,"md.solidearth.love.model"); 348 if(lovemodel==1){ 349 iomodel->FetchData(&love_h,&nl,NULL,"md.solidearth.love.h"); 350 iomodel->FetchData(&love_k,&nl,NULL,"md.solidearth.love.k"); 351 iomodel->FetchData(&love_l,&nl,NULL,"md.solidearth.love.l"); 352 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.love.th",TidalLoveHEnum)); 353 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.love.tk",TidalLoveKEnum)); 354 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.love.tk2secular",TidalLoveK2SecularEnum)); 355 } 356 else{ 357 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.love.h",TidalLoveHEnum)); 358 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.love.k",TidalLoveKEnum)); 359 parameters->AddObject(iomodel->CopyConstantObject("md.solidearth.love.k2secular",TidalLoveK2SecularEnum)); 360 } 353 361 354 362 /*compute elastic green function for a range of angles*/ 355 iomodel->FetchData(°acc,"md.s lr.degacc");363 iomodel->FetchData(°acc,"md.solidearth.settings.degacc"); 356 364 M=reCast<int,IssmDouble>(180./degacc+1.); 357 365 … … 480 488 /*}}}*/ 481 489 /*Transitions:{{{ */ 482 iomodel->FetchData(&transitions,&transitions_M,&transitions_N,&ntransitions,"md.s lr.transitions");490 iomodel->FetchData(&transitions,&transitions_M,&transitions_N,&ntransitions,"md.solidearth.transitions"); 483 491 if(transitions){ 484 492 parameters->AddObject(new DoubleMatArrayParam(SealevelriseTransitionsEnum,transitions,ntransitions,transitions_M,transitions_N)); … … 493 501 } /*}}}*/ 494 502 /*Requested outputs {{{*/ 495 iomodel->FindConstant(&requestedoutputs,&numoutputs,"md.s lr.requested_outputs");503 iomodel->FindConstant(&requestedoutputs,&numoutputs,"md.solidearth.requested_outputs"); 496 504 if(numoutputs)parameters->AddObject(new StringArrayParam(SealevelriseRequestedOutputsEnum,requestedoutputs,numoutputs)); 497 iomodel->DeleteData(&requestedoutputs,numoutputs,"md.s lr.requested_outputs");505 iomodel->DeleteData(&requestedoutputs,numoutputs,"md.solidearth.requested_outputs"); 498 506 /*}}}*/ 499 507 -
issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
r24933 r25118 756 756 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.surface",SurfaceEnum); 757 757 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.base",BaseEnum); 758 iomodel->FetchDataToInput(inputs2,elements,"md.s lr.sealevel",SealevelEnum,0);758 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.sealevel",SealevelEnum,0); 759 759 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); 760 760 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ocean_levelset",MaskOceanLevelsetEnum); -
issm/trunk-jpl/src/c/analyses/StressbalanceVerticalAnalysis.cpp
r24936 r25118 104 104 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.surface",SurfaceEnum); 105 105 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.base",BaseEnum); 106 iomodel->FetchDataToInput(inputs2,elements,"md.s lr.sealevel",SealevelEnum,0);106 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.sealevel",SealevelEnum,0); 107 107 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); 108 108 if(iomodel->domaintype!=Domain2DhorizontalEnum){ -
issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp
r24933 r25118 133 133 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.surface",SurfaceEnum); 134 134 iomodel->FetchDataToInput(inputs2,elements,"md.geometry.base",BaseEnum); 135 iomodel->FetchDataToInput(inputs2,elements,"md.s lr.sealevel",SealevelEnum,0);135 iomodel->FetchDataToInput(inputs2,elements,"md.solidearth.sealevel",SealevelEnum,0); 136 136 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); 137 137 iomodel->FetchDataToInput(inputs2,elements,"md.mask.ocean_levelset",MaskOceanLevelsetEnum); -
issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
r25066 r25118 5233 5233 5234 5234 /*recover earth area: */ 5235 this->parameters->FindParam(&planetarea,S ealevelPlanetAreaEnum);5235 this->parameters->FindParam(&planetarea,SolidearthPlanetAreaEnum); 5236 5236 5237 5237 /*how many dofs are we working with here? */ … … 5387 5387 5388 5388 /*recover earth area: */ 5389 this->parameters->FindParam(&planetarea,S ealevelPlanetAreaEnum);5389 this->parameters->FindParam(&planetarea,SolidearthPlanetAreaEnum); 5390 5390 5391 5391 /*Compute lat,long,radius of elemental centroid: */ … … 5456 5456 5457 5457 /*Compute ice thickness change: */ 5458 Input2* deltathickness_input=this->GetInput2(S ealevelriseDeltathicknessEnum);5458 Input2* deltathickness_input=this->GetInput2(SurfaceloadIceThicknessChangeEnum); 5459 5459 if (!deltathickness_input)_error_("delta thickness input needed to compute sea level rise!"); 5460 5460 deltathickness_input->GetInputAverage(&I); … … 5530 5530 /*recover parameters: */ 5531 5531 rho_earth=FindParam(MaterialsEarthDensityEnum); 5532 this->parameters->FindParam(&computerigid,S ealevelriseRigidEnum);5533 this->parameters->FindParam(&computeelastic,S ealevelriseElasticEnum);5532 this->parameters->FindParam(&computerigid,SolidearthSettingsRigidEnum); 5533 this->parameters->FindParam(&computeelastic,SolidearthSettingsElasticEnum); 5534 5534 this->parameters->FindParam(&gsize,MeshNumberofverticesEnum); 5535 this->parameters->FindParam(&planetarea,S ealevelPlanetAreaEnum);5536 this->parameters->FindParam(&planetradius,S ealevelPlanetRadiusEnum);5537 this->parameters->FindParam(&horiz,S ealevelriseHorizEnum);5535 this->parameters->FindParam(&planetarea,SolidearthPlanetAreaEnum); 5536 this->parameters->FindParam(&planetradius,SolidearthPlanetRadiusEnum); 5537 this->parameters->FindParam(&horiz,SolidearthSettingsHorizEnum); 5538 5538 5539 5539 /*recover precomputed green function kernels:*/ … … 5719 5719 5720 5720 /*recover ocean area scaling: */ 5721 this->parameters->FindParam(&scaleoceanarea,S ealevelriseOceanAreaScalingEnum);5721 this->parameters->FindParam(&scaleoceanarea,SolidearthSettingsOceanAreaScalingEnum); 5722 5722 5723 5723 /*retrieve precomputed G:*/ … … 5737 5737 5738 5738 /*Retrieve ice thickness at vertices: */ 5739 Input2* deltathickness_input=this->GetInput2(S ealevelriseDeltathicknessEnum);5739 Input2* deltathickness_input=this->GetInput2(SurfaceloadIceThicknessChangeEnum); 5740 5740 if (!deltathickness_input)_error_("delta thickness input needed to compute sea level rise!"); 5741 5741 … … 5894 5894 5895 5895 /*recover parameters:*/ 5896 this->parameters->FindParam(&computeelastic,S ealevelriseElasticEnum);5897 this->parameters->FindParam(&horiz,S ealevelriseHorizEnum);5896 this->parameters->FindParam(&computeelastic,SolidearthSettingsElasticEnum); 5897 this->parameters->FindParam(&horiz,SolidearthSettingsHorizEnum); 5898 5898 this->parameters->FindParam(&bp_compute_fingerprints,DslComputeFingerprintsEnum); 5899 5899 … … 5931 5931 5932 5932 /*Compute ice thickness change: */ 5933 Input2* deltathickness_input=this->GetInput2(S ealevelriseDeltathicknessEnum);5933 Input2* deltathickness_input=this->GetInput2(SurfaceloadIceThicknessChangeEnum); 5934 5934 if (!deltathickness_input)_error_("delta thickness input needed to compute sea level rise!"); 5935 5935 deltathickness_input->GetInputAverage(&I); -
issm/trunk-jpl/src/c/classes/FemModel.cpp
r25066 r25118 4760 4760 4761 4761 /*recover computational flags: */ 4762 this->parameters->FindParam(&computerigid,S ealevelriseRigidEnum);4763 this->parameters->FindParam(&computeelastic,S ealevelriseElasticEnum);4762 this->parameters->FindParam(&computerigid,SolidearthSettingsRigidEnum); 4763 this->parameters->FindParam(&computeelastic,SolidearthSettingsElasticEnum); 4764 4764 4765 4765 /*Initialize temporary vector that will be used to sum eustatic components on all local elements, prior … … 4824 4824 4825 4825 /*pull out some useful parameters: */ 4826 parameters->FindParam(&tide_love_h, SealevelriseTidalLoveHEnum);4827 parameters->FindParam(&tide_love_k, SealevelriseTidalLoveKEnum);4828 parameters->FindParam(&fluid_love, SealevelriseFluidLoveEnum);4829 parameters->FindParam(&moi_e, SealevelriseEquatorialMoiEnum);4830 parameters->FindParam(&moi_p, SealevelrisePolarMoiEnum);4831 parameters->FindParam(&omega, SealevelriseAngularVelocityEnum);4826 parameters->FindParam(&tide_love_h,TidalLoveHEnum); 4827 parameters->FindParam(&tide_love_k,TidalLoveKEnum); 4828 parameters->FindParam(&fluid_love,TidalLoveK2SecularEnum); 4829 parameters->FindParam(&moi_e,RotationalEquatorialMoiEnum); 4830 parameters->FindParam(&moi_p,RotationalPolarMoiEnum); 4831 parameters->FindParam(&omega,RotationalAngularVelocityEnum); 4832 4832 4833 4833 /*compute perturbation terms for angular velocity vector: */ … … 4888 4888 4889 4889 /*retrieve parameters:*/ 4890 this->parameters->FindParam(&horiz,S ealevelriseHorizEnum);4890 this->parameters->FindParam(&horiz,SolidearthSettingsHorizEnum); 4891 4891 4892 4892 /*Serialize vectors from previous iteration:*/ -
issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp
r25069 r25118 110 110 111 111 /*Should we even be here?:*/ 112 femmodel->parameters->FindParam(&geodetic,S ealevelriseGeodeticEnum); if(!geodetic)return;112 femmodel->parameters->FindParam(&geodetic,SolidearthSettingsComputesealevelchangeEnum); if(!geodetic)return; 113 113 114 114 /*Verbose: */ … … 117 117 /*retrieve more parameters:*/ 118 118 femmodel->parameters->FindParam(&iscoupler,TransientIscouplerEnum); 119 femmodel->parameters->FindParam(&frequency,S ealevelriseGeodeticRunFrequencyEnum);119 femmodel->parameters->FindParam(&frequency,SolidearthSettingsRunFrequencyEnum); 120 120 femmodel->parameters->FindParam(&count,SealevelriseRunCountEnum); 121 femmodel->parameters->FindParam(&horiz,S ealevelriseHorizEnum);121 femmodel->parameters->FindParam(&horiz,SolidearthSettingsHorizEnum); 122 122 femmodel->parameters->FindParam(&dt,TimesteppingTimeStepEnum); 123 123 femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum); … … 154 154 155 155 /*we have accumulated thicknesses, dump them in deltathcikness: */ 156 if(modelid==earthid)InputDuplicatex(femmodel,SealevelriseCumDeltathicknessEnum,S ealevelriseDeltathicknessEnum);156 if(modelid==earthid)InputDuplicatex(femmodel,SealevelriseCumDeltathicknessEnum,SurfaceloadIceThicknessChangeEnum); 157 157 } 158 158 … … 274 274 femmodel->parameters->FindParam(&dt,TimesteppingTimeStepEnum); 275 275 femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum); 276 femmodel->parameters->FindParam(&geodetic,S ealevelriseGeodeticEnum);276 femmodel->parameters->FindParam(&geodetic,SolidearthSettingsComputesealevelchangeEnum); 277 277 278 278 /*in case we are running SealevelriseSolutionEnum, then bypass transient settings:*/ … … 290 290 GetStericRate(&steric_rate_g,femmodel); 291 291 GetDynamicRate(&dynamic_rate_g,femmodel); 292 GetVectorFromInputsx(&hydro_rate_g,femmodel,S ealevelriseHydroRateEnum,VertexSIdEnum);292 GetVectorFromInputsx(&hydro_rate_g,femmodel,SurfaceloadWaterHeightChangeEnum,VertexSIdEnum); 293 293 if(geodetic){ 294 294 GetVectorFromInputsx(&U_esa_rate,femmodel,SealevelUEsaRateEnum,VertexSIdEnum); … … 364 364 365 365 /*retrieve parameters:*/ 366 femmodel->parameters->FindParam(&horiz,S ealevelriseHorizEnum);366 femmodel->parameters->FindParam(&horiz,SolidearthSettingsHorizEnum); 367 367 femmodel->parameters->FindParam(&geometrydone,SealevelriseGeometryDoneEnum); 368 368 … … 470 470 471 471 /*Recover some parameters: */ 472 femmodel->parameters->FindParam(&max_nonlinear_iterations,S ealevelriseMaxiterEnum);473 femmodel->parameters->FindParam(&eps_rel,S ealevelriseReltolEnum);474 femmodel->parameters->FindParam(&eps_abs,S ealevelriseAbstolEnum);472 femmodel->parameters->FindParam(&max_nonlinear_iterations,SolidearthSettingsMaxiterEnum); 473 femmodel->parameters->FindParam(&eps_rel,SolidearthSettingsReltolEnum); 474 femmodel->parameters->FindParam(&eps_abs,SolidearthSettingsAbstolEnum); 475 475 476 476 /*computational flag: */ 477 femmodel->parameters->FindParam(&rotation,S ealevelriseRotationEnum);477 femmodel->parameters->FindParam(&rotation,SolidearthSettingsRotationEnum); 478 478 479 479 /*Figure out size of g-set deflection vector and allocate solution vector: */ … … 589 589 590 590 /*retrieve some parameters:*/ 591 femmodel->parameters->FindParam(&horiz,S ealevelriseHorizEnum);591 femmodel->parameters->FindParam(&horiz,SolidearthSettingsHorizEnum); 592 592 593 593 /*find size of vectors:*/ … … 637 637 638 638 /*retrieve some parameters:*/ 639 femmodel->parameters->FindParam(&frequency,S ealevelriseGeodeticRunFrequencyEnum);639 femmodel->parameters->FindParam(&frequency,SolidearthSettingsRunFrequencyEnum); 640 640 femmodel->parameters->FindParam(&dt,TimesteppingTimeStepEnum); 641 641 … … 990 990 if(VerboseSolution()) _printf0_(" computing earth mass transport\n"); 991 991 992 /*This module has to be recoded! We do not grab spc thicknesses from an earth model 993 * anymore! Thicknesses come from a mass transport module applied to each basin! 994 * Commeting out for now:*/ 995 _error_("EarthMassTransport error message: not supported anymore!"); 996 992 997 /*This mass transport module for the Earth is because we might have thickness variations as spcs 993 998 * specified in the md.slr class, outside of what we will get from the icecaps. That's why we get t … … 996 1001 /*No mass transport module was called, so we are just going to retrieve the geometry thickness 997 1002 * at this time step, at prior time step, and plug the difference as deltathickness: */ 998 femmodel->parameters->FindParam(&time,TimeEnum);1003 /*femmodel->parameters->FindParam(&time,TimeEnum); 999 1004 femmodel->parameters->FindParam(&dt,TimesteppingTimeStepEnum); 1000 1005 nv=femmodel->vertices->NumberOfVertices(); 1001 1006 1002 1007 GetVectorFromInputsx(&newthickness,femmodel,SealevelriseSpcthicknessEnum,VertexSIdEnum); 1003 GetVectorFromInputsx(&oldthickness,femmodel,SealevelriseSpcthicknessEnum,VertexSIdEnum,time-dt); 1008 GetVectorFromInputsx(&oldthickness,femmodel,SealevelriseSpcthicknessEnum,VertexSIdEnum,time-dt);*/ 1004 1009 1005 1010 /*compute deltathickness: */ 1006 deltathickness = new Vector<IssmDouble>(nv);1007 newthickness->Copy(deltathickness); deltathickness->AXPY(oldthickness,-1); 1011 /*deltathickness = new Vector<IssmDouble>(nv); 1012 newthickness->Copy(deltathickness); deltathickness->AXPY(oldthickness,-1); */ 1008 1013 1009 1014 /*plug into elements:*/ 1010 InputUpdateFromVectorx(femmodel,deltathickness,SealevelriseDeltathicknessEnum,VertexSIdEnum);1015 // InputUpdateFromVectorx(femmodel,deltathickness,SurfaceloadIceThicknessChangeEnum,VertexSIdEnum); 1011 1016 1012 1017 /*add to cumulated delta thickness: */ 1013 GetVectorFromInputsx(&cumdeltathickness,femmodel,SealevelriseCumDeltathicknessEnum,VertexSIdEnum);1018 /*GetVectorFromInputsx(&cumdeltathickness,femmodel,SealevelriseCumDeltathicknessEnum,VertexSIdEnum); 1014 1019 cumdeltathickness->AXPY(deltathickness,1); 1015 InputUpdateFromVectorx(femmodel,cumdeltathickness,SealevelriseCumDeltathicknessEnum,VertexSIdEnum); 1020 InputUpdateFromVectorx(femmodel,cumdeltathickness,SealevelriseCumDeltathicknessEnum,VertexSIdEnum);*/ 1016 1021 1017 1022 /*free ressources:*/ 1018 delete oldthickness;1023 /*delete oldthickness; 1019 1024 delete newthickness; 1020 1025 delete deltathickness; 1021 delete cumdeltathickness; 1026 delete cumdeltathickness;*/ 1022 1027 1023 1028 } /*}}}*/ -
issm/trunk-jpl/src/c/shared/Enum/Enum.vim
r25066 r25118 307 307 syn keyword cConstant RootPathEnum 308 308 syn keyword cConstant SaveResultsEnum 309 syn keyword cConstant S ealevelPlanetRadiusEnum310 syn keyword cConstant S ealevelPlanetAreaEnum309 syn keyword cConstant SolidearthPlanetRadiusEnum 310 syn keyword cConstant SolidearthPlanetAreaEnum 311 311 syn keyword cConstant SealevelEustaticEnum 312 syn keyword cConstant S ealevelriseAbstolEnum313 syn keyword cConstant SealevelriseAngularVelocityEnum314 syn keyword cConstant S ealevelriseElasticEnum312 syn keyword cConstant SolidearthSettingsAbstolEnum 313 syn keyword cConstant RotationalAngularVelocityEnum 314 syn keyword cConstant SolidearthSettingsElasticEnum 315 315 syn keyword cConstant SealevelriseGeometryDoneEnum 316 syn keyword cConstant SealevelriseEquatorialMoiEnum317 syn keyword cConstant SealevelriseFluidLoveEnum316 syn keyword cConstant RotationalEquatorialMoiEnum 317 syn keyword cConstant TidalLoveK2SecularEnum 318 318 syn keyword cConstant SealevelriseGRigidEnum 319 319 syn keyword cConstant SealevelriseGElasticEnum 320 syn keyword cConstant S ealevelriseGeodeticEnum321 syn keyword cConstant S ealevelriseGeodeticRunFrequencyEnum320 syn keyword cConstant SolidearthSettingsComputesealevelchangeEnum 321 syn keyword cConstant SolidearthSettingsRunFrequencyEnum 322 322 syn keyword cConstant SealevelriseHElasticEnum 323 syn keyword cConstant S ealevelriseHorizEnum324 syn keyword cConstant S ealevelriseMaxiterEnum325 syn keyword cConstant S ealevelriseOceanAreaScalingEnum326 syn keyword cConstant SealevelrisePolarMoiEnum327 syn keyword cConstant S ealevelriseReltolEnum323 syn keyword cConstant SolidearthSettingsHorizEnum 324 syn keyword cConstant SolidearthSettingsMaxiterEnum 325 syn keyword cConstant SolidearthSettingsOceanAreaScalingEnum 326 syn keyword cConstant RotationalPolarMoiEnum 327 syn keyword cConstant SolidearthSettingsReltolEnum 328 328 syn keyword cConstant SealevelriseRequestedOutputsEnum 329 syn keyword cConstant S ealevelriseRigidEnum330 syn keyword cConstant S ealevelriseRotationEnum329 syn keyword cConstant SolidearthSettingsRigidEnum 330 syn keyword cConstant SolidearthSettingsRotationEnum 331 331 syn keyword cConstant SealevelriseRunCountEnum 332 syn keyword cConstant SealevelriseTidalLoveHEnum333 syn keyword cConstant SealevelriseTidalLoveKEnum332 syn keyword cConstant TidalLoveHEnum 333 syn keyword cConstant TidalLoveKEnum 334 334 syn keyword cConstant SealevelriseTransitionsEnum 335 335 syn keyword cConstant SealevelriseUElasticEnum … … 687 687 syn keyword cConstant SealevelriseCumDeltathicknessEnum 688 688 syn keyword cConstant SealevelriseCumDeltathicknessOldEnum 689 syn keyword cConstant S ealevelriseDeltathicknessEnum690 syn keyword cConstant S ealevelriseSpcthicknessEnum691 syn keyword cConstant S ealevelriseHydroRateEnum689 syn keyword cConstant SurfaceloadOtherEnum 690 syn keyword cConstant SurfaceloadIceThicknessChangeEnum 691 syn keyword cConstant SurfaceloadWaterHeightChangeEnum 692 692 syn keyword cConstant SealevelriseIndicesEnum 693 693 syn keyword cConstant SealevelriseGEnum -
issm/trunk-jpl/src/c/shared/Enum/EnumDefinitions.h
r25066 r25118 301 301 RootPathEnum, 302 302 SaveResultsEnum, 303 S ealevelPlanetRadiusEnum,304 S ealevelPlanetAreaEnum,303 SolidearthPlanetRadiusEnum, 304 SolidearthPlanetAreaEnum, 305 305 SealevelEustaticEnum, 306 S ealevelriseAbstolEnum,307 SealevelriseAngularVelocityEnum,308 S ealevelriseElasticEnum,306 SolidearthSettingsAbstolEnum, 307 RotationalAngularVelocityEnum, 308 SolidearthSettingsElasticEnum, 309 309 SealevelriseGeometryDoneEnum, 310 SealevelriseEquatorialMoiEnum,311 SealevelriseFluidLoveEnum,310 RotationalEquatorialMoiEnum, 311 TidalLoveK2SecularEnum, 312 312 SealevelriseGRigidEnum, 313 313 SealevelriseGElasticEnum, 314 S ealevelriseGeodeticEnum,315 S ealevelriseGeodeticRunFrequencyEnum,314 SolidearthSettingsComputesealevelchangeEnum, 315 SolidearthSettingsRunFrequencyEnum, 316 316 SealevelriseHElasticEnum, 317 S ealevelriseHorizEnum,318 S ealevelriseMaxiterEnum,319 S ealevelriseOceanAreaScalingEnum,320 SealevelrisePolarMoiEnum,321 S ealevelriseReltolEnum,317 SolidearthSettingsHorizEnum, 318 SolidearthSettingsMaxiterEnum, 319 SolidearthSettingsOceanAreaScalingEnum, 320 RotationalPolarMoiEnum, 321 SolidearthSettingsReltolEnum, 322 322 SealevelriseRequestedOutputsEnum, 323 S ealevelriseRigidEnum,324 S ealevelriseRotationEnum,323 SolidearthSettingsRigidEnum, 324 SolidearthSettingsRotationEnum, 325 325 SealevelriseRunCountEnum, 326 SealevelriseTidalLoveHEnum,327 SealevelriseTidalLoveKEnum,326 TidalLoveHEnum, 327 TidalLoveKEnum, 328 328 SealevelriseTransitionsEnum, 329 329 SealevelriseUElasticEnum, … … 684 684 SealevelriseCumDeltathicknessEnum, 685 685 SealevelriseCumDeltathicknessOldEnum, 686 S ealevelriseDeltathicknessEnum,687 S ealevelriseSpcthicknessEnum,688 S ealevelriseHydroRateEnum,686 SurfaceloadOtherEnum, 687 SurfaceloadIceThicknessChangeEnum, 688 SurfaceloadWaterHeightChangeEnum, 689 689 SealevelriseIndicesEnum, 690 690 SealevelriseGEnum, -
issm/trunk-jpl/src/c/shared/Enum/EnumToStringx.cpp
r25066 r25118 309 309 case RootPathEnum : return "RootPath"; 310 310 case SaveResultsEnum : return "SaveResults"; 311 case S ealevelPlanetRadiusEnum : return "SealevelPlanetRadius";312 case S ealevelPlanetAreaEnum : return "SealevelPlanetArea";311 case SolidearthPlanetRadiusEnum : return "SolidearthPlanetRadius"; 312 case SolidearthPlanetAreaEnum : return "SolidearthPlanetArea"; 313 313 case SealevelEustaticEnum : return "SealevelEustatic"; 314 case S ealevelriseAbstolEnum : return "SealevelriseAbstol";315 case SealevelriseAngularVelocityEnum : return "SealevelriseAngularVelocity";316 case S ealevelriseElasticEnum : return "SealevelriseElastic";314 case SolidearthSettingsAbstolEnum : return "SolidearthSettingsAbstol"; 315 case RotationalAngularVelocityEnum : return "RotationalAngularVelocity"; 316 case SolidearthSettingsElasticEnum : return "SolidearthSettingsElastic"; 317 317 case SealevelriseGeometryDoneEnum : return "SealevelriseGeometryDone"; 318 case SealevelriseEquatorialMoiEnum : return "SealevelriseEquatorialMoi";319 case SealevelriseFluidLoveEnum : return "SealevelriseFluidLove";318 case RotationalEquatorialMoiEnum : return "RotationalEquatorialMoi"; 319 case TidalLoveK2SecularEnum : return "TidalLoveK2Secular"; 320 320 case SealevelriseGRigidEnum : return "SealevelriseGRigid"; 321 321 case SealevelriseGElasticEnum : return "SealevelriseGElastic"; 322 case S ealevelriseGeodeticEnum : return "SealevelriseGeodetic";323 case S ealevelriseGeodeticRunFrequencyEnum : return "SealevelriseGeodeticRunFrequency";322 case SolidearthSettingsComputesealevelchangeEnum : return "SolidearthSettingsComputesealevelchange"; 323 case SolidearthSettingsRunFrequencyEnum : return "SolidearthSettingsRunFrequency"; 324 324 case SealevelriseHElasticEnum : return "SealevelriseHElastic"; 325 case S ealevelriseHorizEnum : return "SealevelriseHoriz";326 case S ealevelriseMaxiterEnum : return "SealevelriseMaxiter";327 case S ealevelriseOceanAreaScalingEnum : return "SealevelriseOceanAreaScaling";328 case SealevelrisePolarMoiEnum : return "SealevelrisePolarMoi";329 case S ealevelriseReltolEnum : return "SealevelriseReltol";325 case SolidearthSettingsHorizEnum : return "SolidearthSettingsHoriz"; 326 case SolidearthSettingsMaxiterEnum : return "SolidearthSettingsMaxiter"; 327 case SolidearthSettingsOceanAreaScalingEnum : return "SolidearthSettingsOceanAreaScaling"; 328 case RotationalPolarMoiEnum : return "RotationalPolarMoi"; 329 case SolidearthSettingsReltolEnum : return "SolidearthSettingsReltol"; 330 330 case SealevelriseRequestedOutputsEnum : return "SealevelriseRequestedOutputs"; 331 case S ealevelriseRigidEnum : return "SealevelriseRigid";332 case S ealevelriseRotationEnum : return "SealevelriseRotation";331 case SolidearthSettingsRigidEnum : return "SolidearthSettingsRigid"; 332 case SolidearthSettingsRotationEnum : return "SolidearthSettingsRotation"; 333 333 case SealevelriseRunCountEnum : return "SealevelriseRunCount"; 334 case SealevelriseTidalLoveHEnum : return "SealevelriseTidalLoveH";335 case SealevelriseTidalLoveKEnum : return "SealevelriseTidalLoveK";334 case TidalLoveHEnum : return "TidalLoveH"; 335 case TidalLoveKEnum : return "TidalLoveK"; 336 336 case SealevelriseTransitionsEnum : return "SealevelriseTransitions"; 337 337 case SealevelriseUElasticEnum : return "SealevelriseUElastic"; … … 689 689 case SealevelriseCumDeltathicknessEnum : return "SealevelriseCumDeltathickness"; 690 690 case SealevelriseCumDeltathicknessOldEnum : return "SealevelriseCumDeltathicknessOld"; 691 case S ealevelriseDeltathicknessEnum : return "SealevelriseDeltathickness";692 case S ealevelriseSpcthicknessEnum : return "SealevelriseSpcthickness";693 case S ealevelriseHydroRateEnum : return "SealevelriseHydroRate";691 case SurfaceloadOtherEnum : return "SurfaceloadOther"; 692 case SurfaceloadIceThicknessChangeEnum : return "SurfaceloadIceThicknessChange"; 693 case SurfaceloadWaterHeightChangeEnum : return "SurfaceloadWaterHeightChange"; 694 694 case SealevelriseIndicesEnum : return "SealevelriseIndices"; 695 695 case SealevelriseGEnum : return "SealevelriseG"; -
issm/trunk-jpl/src/c/shared/Enum/StringToEnumx.cpp
r25066 r25118 315 315 else if (strcmp(name,"RootPath")==0) return RootPathEnum; 316 316 else if (strcmp(name,"SaveResults")==0) return SaveResultsEnum; 317 else if (strcmp(name,"S ealevelPlanetRadius")==0) return SealevelPlanetRadiusEnum;318 else if (strcmp(name,"S ealevelPlanetArea")==0) return SealevelPlanetAreaEnum;317 else if (strcmp(name,"SolidearthPlanetRadius")==0) return SolidearthPlanetRadiusEnum; 318 else if (strcmp(name,"SolidearthPlanetArea")==0) return SolidearthPlanetAreaEnum; 319 319 else if (strcmp(name,"SealevelEustatic")==0) return SealevelEustaticEnum; 320 else if (strcmp(name,"S ealevelriseAbstol")==0) return SealevelriseAbstolEnum;321 else if (strcmp(name," SealevelriseAngularVelocity")==0) return SealevelriseAngularVelocityEnum;322 else if (strcmp(name,"S ealevelriseElastic")==0) return SealevelriseElasticEnum;320 else if (strcmp(name,"SolidearthSettingsAbstol")==0) return SolidearthSettingsAbstolEnum; 321 else if (strcmp(name,"RotationalAngularVelocity")==0) return RotationalAngularVelocityEnum; 322 else if (strcmp(name,"SolidearthSettingsElastic")==0) return SolidearthSettingsElasticEnum; 323 323 else if (strcmp(name,"SealevelriseGeometryDone")==0) return SealevelriseGeometryDoneEnum; 324 else if (strcmp(name," SealevelriseEquatorialMoi")==0) return SealevelriseEquatorialMoiEnum;325 else if (strcmp(name," SealevelriseFluidLove")==0) return SealevelriseFluidLoveEnum;324 else if (strcmp(name,"RotationalEquatorialMoi")==0) return RotationalEquatorialMoiEnum; 325 else if (strcmp(name,"TidalLoveK2Secular")==0) return TidalLoveK2SecularEnum; 326 326 else if (strcmp(name,"SealevelriseGRigid")==0) return SealevelriseGRigidEnum; 327 327 else if (strcmp(name,"SealevelriseGElastic")==0) return SealevelriseGElasticEnum; 328 else if (strcmp(name,"S ealevelriseGeodetic")==0) return SealevelriseGeodeticEnum;329 else if (strcmp(name,"S ealevelriseGeodeticRunFrequency")==0) return SealevelriseGeodeticRunFrequencyEnum;328 else if (strcmp(name,"SolidearthSettingsComputesealevelchange")==0) return SolidearthSettingsComputesealevelchangeEnum; 329 else if (strcmp(name,"SolidearthSettingsRunFrequency")==0) return SolidearthSettingsRunFrequencyEnum; 330 330 else if (strcmp(name,"SealevelriseHElastic")==0) return SealevelriseHElasticEnum; 331 else if (strcmp(name,"S ealevelriseHoriz")==0) return SealevelriseHorizEnum;332 else if (strcmp(name,"S ealevelriseMaxiter")==0) return SealevelriseMaxiterEnum;333 else if (strcmp(name,"S ealevelriseOceanAreaScaling")==0) return SealevelriseOceanAreaScalingEnum;334 else if (strcmp(name," SealevelrisePolarMoi")==0) return SealevelrisePolarMoiEnum;335 else if (strcmp(name,"S ealevelriseReltol")==0) return SealevelriseReltolEnum;331 else if (strcmp(name,"SolidearthSettingsHoriz")==0) return SolidearthSettingsHorizEnum; 332 else if (strcmp(name,"SolidearthSettingsMaxiter")==0) return SolidearthSettingsMaxiterEnum; 333 else if (strcmp(name,"SolidearthSettingsOceanAreaScaling")==0) return SolidearthSettingsOceanAreaScalingEnum; 334 else if (strcmp(name,"RotationalPolarMoi")==0) return RotationalPolarMoiEnum; 335 else if (strcmp(name,"SolidearthSettingsReltol")==0) return SolidearthSettingsReltolEnum; 336 336 else if (strcmp(name,"SealevelriseRequestedOutputs")==0) return SealevelriseRequestedOutputsEnum; 337 else if (strcmp(name,"S ealevelriseRigid")==0) return SealevelriseRigidEnum;338 else if (strcmp(name,"S ealevelriseRotation")==0) return SealevelriseRotationEnum;337 else if (strcmp(name,"SolidearthSettingsRigid")==0) return SolidearthSettingsRigidEnum; 338 else if (strcmp(name,"SolidearthSettingsRotation")==0) return SolidearthSettingsRotationEnum; 339 339 else if (strcmp(name,"SealevelriseRunCount")==0) return SealevelriseRunCountEnum; 340 else if (strcmp(name," SealevelriseTidalLoveH")==0) return SealevelriseTidalLoveHEnum;341 else if (strcmp(name," SealevelriseTidalLoveK")==0) return SealevelriseTidalLoveKEnum;340 else if (strcmp(name,"TidalLoveH")==0) return TidalLoveHEnum; 341 else if (strcmp(name,"TidalLoveK")==0) return TidalLoveKEnum; 342 342 else if (strcmp(name,"SealevelriseTransitions")==0) return SealevelriseTransitionsEnum; 343 343 else if (strcmp(name,"SealevelriseUElastic")==0) return SealevelriseUElasticEnum; … … 704 704 else if (strcmp(name,"SealevelriseCumDeltathickness")==0) return SealevelriseCumDeltathicknessEnum; 705 705 else if (strcmp(name,"SealevelriseCumDeltathicknessOld")==0) return SealevelriseCumDeltathicknessOldEnum; 706 else if (strcmp(name,"S ealevelriseDeltathickness")==0) return SealevelriseDeltathicknessEnum;707 else if (strcmp(name,"S ealevelriseSpcthickness")==0) return SealevelriseSpcthicknessEnum;708 else if (strcmp(name,"S ealevelriseHydroRate")==0) return SealevelriseHydroRateEnum;706 else if (strcmp(name,"SurfaceloadOther")==0) return SurfaceloadOtherEnum; 707 else if (strcmp(name,"SurfaceloadIceThicknessChange")==0) return SurfaceloadIceThicknessChangeEnum; 708 else if (strcmp(name,"SurfaceloadWaterHeightChange")==0) return SurfaceloadWaterHeightChangeEnum; 709 709 else if (strcmp(name,"SealevelriseIndices")==0) return SealevelriseIndicesEnum; 710 710 else if (strcmp(name,"SealevelriseG")==0) return SealevelriseGEnum; -
issm/trunk-jpl/src/c/shared/io/Marshalling/IoCodeConversions.cpp
r25066 r25118 163 163 fieldname=xNew<char>((strlen(field)+1)); xMemCpy<char>(fieldname,field,(strlen(field)+1)); 164 164 } 165 else if(strcmp(string_in,"S ealevelriseDeltathickness")==0){166 const char* field = "md.s lr.deltathickness";167 input_enum = S ealevelriseDeltathicknessEnum;165 else if(strcmp(string_in,"SurfaceloadIceThicknessChange")==0){ 166 const char* field = "md.solidearth.surfaceload.iceheightchange"; 167 input_enum = SurfaceloadIceThicknessChangeEnum; 168 168 fieldname=xNew<char>((strlen(field)+1)); xMemCpy<char>(fieldname,field,(strlen(field)+1)); 169 169 } … … 178 178 fieldname=xNew<char>((strlen(field)+1)); xMemCpy<char>(fieldname,field,(strlen(field)+1)); 179 179 } 180 else if(strcmp(string_in,"S ealevelriseHydroRate")==0){181 const char* field = "md.s lr.hydro_rate";182 input_enum = S ealevelriseHydroRateEnum;180 else if(strcmp(string_in,"SurfaceloadWaterHeightChange")==0){ 181 const char* field = "md.solidearth.surfaceload.waterheightchange"; 182 input_enum = SurfaceloadWaterHeightChangeEnum; 183 183 fieldname=xNew<char>((strlen(field)+1)); xMemCpy<char>(fieldname,field,(strlen(field)+1)); 184 184 } -
issm/trunk-jpl/src/m/classes/model.m
r25074 r25118 23 23 rifts = 0; 24 24 dsl = 0; 25 s lr= 0;25 solidearth = 0; 26 26 27 27 debug = 0; … … 126 126 if isa(md.calving,'double'); md.calving=calving(); end 127 127 %2016 February 3 128 if isa(md.slr,'double'); md.s lr=slr(); end128 if isa(md.slr,'double'); md.solidearth=solidearth(); end 129 129 %2016 October 11 130 130 if isa(md.esa,'double'); md.esa=esa(); end … … 888 888 md.frontalforcings=extrude(md.frontalforcings,md); 889 889 md.hydrology = extrude(md.hydrology,md); 890 md.s lr = extrude(md.slr,md);890 md.solidearth = extrude(md.solidearth,md); 891 891 md.dsl = extrude(md.dsl,md); 892 892 … … 1209 1209 md.friction = friction(); 1210 1210 md.rifts = rifts(); 1211 md.s lr = slr();1211 md.solidearth = solidearth(); 1212 1212 md.dsl = dsl(); 1213 1213 md.timestepping = timestepping(); … … 1390 1390 disp(sprintf('%19s: %-22s -- %s','initialization' ,['[1x1 ' class(self.initialization) ']'],'initial guess/state')); 1391 1391 disp(sprintf('%19s: %-22s -- %s','rifts' ,['[1x1 ' class(self.rifts) ']'],'rifts properties')); 1392 disp(sprintf('%19s: %-22s -- %s','s lr' ,['[1x1 ' class(self.slr) ']'],'slr forcings'));1392 disp(sprintf('%19s: %-22s -- %s','solidearth' ,['[1x1 ' class(self.solidearth) ']'],'solidearth inputs and settings')); 1393 1393 disp(sprintf('%19s: %-22s -- %s','dsl' ,['[1x1 ' class(self.dsl) ']'],'dynamic sea-level ')); 1394 1394 disp(sprintf('%19s: %-22s -- %s','debug' ,['[1x1 ' class(self.debug) ']'],'debugging tools (valgrind, gprof)'));
Note:
See TracChangeset
for help on using the changeset viewer.