Changeset 17946
- Timestamp:
- 05/06/14 08:30:09 (12 years ago)
- Location:
- issm/trunk-jpl/src
- Files:
-
- 1 added
- 1 deleted
- 16 edited
-
c/analyses/EnthalpyAnalysis.cpp (modified) (4 diffs)
-
c/analyses/HydrologyShreveAnalysis.cpp (modified) (1 diff)
-
c/analyses/MasstransportAnalysis.cpp (modified) (1 diff)
-
c/analyses/StressbalanceAnalysis.cpp (modified) (6 diffs)
-
c/analyses/ThermalAnalysis.cpp (modified) (3 diffs)
-
c/classes/Loads/Friction.cpp (modified) (2 diffs)
-
c/classes/Materials/Matpar.cpp (modified) (3 diffs)
-
c/shared/Enum/EnumDefinitions.h (modified) (1 diff)
-
c/shared/Enum/EnumToStringx.cpp (modified) (1 diff)
-
c/shared/Enum/StringToEnumx.cpp (modified) (1 diff)
-
m/classes/matdamageice.m (modified) (1 diff)
-
m/classes/matdamageice.py (modified) (1 diff)
-
m/classes/matice.m (modified) (2 diffs)
-
m/classes/matice.py (modified) (1 diff)
-
m/enum/EnumDefinitions.py (modified) (1 diff)
-
m/enum/MaterialsRhoSeawaterEnum.m (added)
-
m/enum/MaterialsRhoWaterEnum.m (deleted)
-
m/qmu/setupdesign/IsScaled.m (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
r17944 r17946 245 245 element->FindParam(&dt,TimesteppingTimeStepEnum); 246 246 element->FindParam(&stabilization,ThermalStabilizationEnum); 247 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRho WaterEnum);247 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 248 248 IssmDouble rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum); 249 249 IssmDouble gravity = element->GetMaterialParameter(ConstantsGEnum); … … 373 373 element->FindParam(&dt,TimesteppingTimeStepEnum); 374 374 IssmDouble gravity = element->GetMaterialParameter(ConstantsGEnum); 375 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRho WaterEnum);375 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 376 376 IssmDouble rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum); 377 377 IssmDouble heatcapacity = element->GetMaterialParameter(MaterialsHeatcapacityEnum); … … 648 648 Input* pressure_input=element->GetInput(PressureEnum); _assert_(pressure_input); 649 649 IssmDouble gravity = element->GetMaterialParameter(ConstantsGEnum); 650 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRho WaterEnum);650 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 651 651 IssmDouble rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum); 652 652 IssmDouble heatcapacity = element->GetMaterialParameter(MaterialsHeatcapacityEnum); … … 1103 1103 1104 1104 rho_ice=element->GetMaterialParameter(MaterialsRhoIceEnum); 1105 rho_water=element->GetMaterialParameter(MaterialsRho WaterEnum);1105 rho_water=element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 1106 1106 1107 1107 element->GetVerticesCoordinates(&xyz_list); -
issm/trunk-jpl/src/c/analyses/HydrologyShreveAnalysis.cpp
r17886 r17946 331 331 /*Retrieve all inputs and parameters*/ 332 332 IssmDouble rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum); 333 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRho WaterEnum);333 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 334 334 IssmDouble g = element->GetMaterialParameter(ConstantsGEnum); 335 335 IssmDouble mu_water = element->GetMaterialParameter(MaterialsMuWaterEnum); -
issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp
r17924 r17946 711 711 basalelement->FindParam(&hydroadjustment,MasstransportHydrostaticAdjustmentEnum); 712 712 rho_ice = basalelement->GetMaterialParameter(MaterialsRhoIceEnum); 713 rho_water = basalelement->GetMaterialParameter(MaterialsRho WaterEnum);713 rho_water = basalelement->GetMaterialParameter(MaterialsRhoSeawaterEnum); 714 714 715 715 for(i=0;i<numnodes;i++) { -
issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
r17943 r17946 1477 1477 Input* thickness_input = element->GetInput(ThicknessEnum); _assert_(thickness_input); 1478 1478 Input* base_input = element->GetInput(BaseEnum); _assert_(base_input); 1479 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRho WaterEnum);1479 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 1480 1480 IssmDouble rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum); 1481 1481 IssmDouble gravity = element->GetMaterialParameter(ConstantsGEnum); … … 1929 1929 Input* thickness_input = element->GetInput(ThicknessEnum); _assert_(thickness_input); 1930 1930 Input* base_input = element->GetInput(BaseEnum); _assert_(base_input); 1931 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRho WaterEnum);1931 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 1932 1932 IssmDouble rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum); 1933 1933 IssmDouble gravity = element->GetMaterialParameter(ConstantsGEnum); … … 2447 2447 /*Retrieve all inputs and parameters*/ 2448 2448 Input* surface_input = element->GetInput(SurfaceEnum); _assert_(surface_input); 2449 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRho WaterEnum);2449 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 2450 2450 IssmDouble rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum); 2451 2451 IssmDouble gravity = element->GetMaterialParameter(ConstantsGEnum); … … 3012 3012 element->FindParam(&dt,TimesteppingTimeStepEnum); 3013 3013 if(dt==0) dt=1.e+5; 3014 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRho WaterEnum);3014 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 3015 3015 IssmDouble gravity = element->GetMaterialParameter(ConstantsGEnum); 3016 3016 Input* surface_input = element->GetInput(SurfaceEnum); _assert_(surface_input); … … 3569 3569 element->GetVerticesCoordinatesBase(&xyz_list_base); 3570 3570 Input* base_input=element->GetInput(BaseEnum); _assert_(base_input); 3571 IssmDouble rho_water=element->GetMaterialParameter(MaterialsRho WaterEnum);3571 IssmDouble rho_water=element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 3572 3572 IssmDouble gravity =element->GetMaterialParameter(ConstantsGEnum); 3573 3573 … … 3643 3643 element->NormalSection(&normal[0],xyz_list_front); 3644 3644 Input* surface_input = element->GetInput(SurfaceEnum); _assert_(surface_input); 3645 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRho WaterEnum);3645 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 3646 3646 IssmDouble gravity = element->GetMaterialParameter(ConstantsGEnum); 3647 3647 -
issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp
r17943 r17946 178 178 element->FindParam(&dt,TimesteppingTimeStepEnum); 179 179 element->FindParam(&stabilization,ThermalStabilizationEnum); 180 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRho WaterEnum);180 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 181 181 IssmDouble rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum); 182 182 IssmDouble gravity = element->GetMaterialParameter(ConstantsGEnum); … … 307 307 element->FindParam(&dt,TimesteppingTimeStepEnum); 308 308 IssmDouble gravity = element->GetMaterialParameter(ConstantsGEnum); 309 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRho WaterEnum);309 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 310 310 IssmDouble rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum); 311 311 IssmDouble heatcapacity = element->GetMaterialParameter(MaterialsHeatcapacityEnum); … … 520 520 Input* pressure_input=element->GetInput(PressureEnum); _assert_(pressure_input); 521 521 IssmDouble gravity = element->GetMaterialParameter(ConstantsGEnum); 522 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRho WaterEnum);522 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 523 523 IssmDouble rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum); 524 524 IssmDouble heatcapacity = element->GetMaterialParameter(MaterialsHeatcapacityEnum); -
issm/trunk-jpl/src/c/classes/Loads/Friction.cpp
r17943 r17946 62 62 element->GetInputValue(&bed, gauss,BaseEnum); 63 63 element->GetInputValue(&drag_coefficient, gauss,FrictionCoefficientEnum); 64 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRho WaterEnum);64 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 65 65 IssmDouble rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum); 66 66 IssmDouble gravity = element->GetMaterialParameter(ConstantsGEnum); … … 124 124 element->GetInputValue(&bed, gauss,BaseEnum); 125 125 element->GetInputValue(&drag_coefficient, gauss,FrictionCoefficientEnum); 126 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRho WaterEnum);126 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum); 127 127 IssmDouble rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum); 128 128 IssmDouble gravity = element->GetMaterialParameter(ConstantsGEnum); -
issm/trunk-jpl/src/c/classes/Materials/Matpar.cpp
r17926 r17946 28 28 this->mid = matpar_mid; 29 29 iomodel->Constant(&this->rho_ice,MaterialsRhoIceEnum); 30 iomodel->Constant(&this->rho_water,MaterialsRho WaterEnum);30 iomodel->Constant(&this->rho_water,MaterialsRhoSeawaterEnum); 31 31 iomodel->Constant(&this->rho_freshwater,MaterialsRhoFreshwaterEnum); 32 32 iomodel->Constant(&this->mu_water,MaterialsMuWaterEnum); … … 164 164 this->rho_ice=constant; 165 165 break; 166 case MaterialsRho WaterEnum:166 case MaterialsRhoSeawaterEnum: 167 167 this->rho_water=constant; 168 168 break; … … 239 239 switch(enum_in){ 240 240 case MaterialsRhoIceEnum: return this->rho_ice; 241 case MaterialsRho WaterEnum: return this->rho_water;241 case MaterialsRhoSeawaterEnum: return this->rho_water; 242 242 case MaterialsRhoFreshwaterEnum: return this->rho_freshwater; 243 243 case MaterialsMuWaterEnum: return this->mu_water; -
issm/trunk-jpl/src/c/shared/Enum/EnumDefinitions.h
r17941 r17946 192 192 DamageEvolutionRequestedOutputsEnum, 193 193 MaterialsRhoIceEnum, 194 MaterialsRho WaterEnum,194 MaterialsRhoSeawaterEnum, 195 195 MaterialsRhoFreshwaterEnum, 196 196 MaterialsMuWaterEnum, -
issm/trunk-jpl/src/c/shared/Enum/EnumToStringx.cpp
r17941 r17946 200 200 case DamageEvolutionRequestedOutputsEnum : return "DamageEvolutionRequestedOutputs"; 201 201 case MaterialsRhoIceEnum : return "MaterialsRhoIce"; 202 case MaterialsRho WaterEnum : return "MaterialsRhoWater";202 case MaterialsRhoSeawaterEnum : return "MaterialsRhoSeawater"; 203 203 case MaterialsRhoFreshwaterEnum : return "MaterialsRhoFreshwater"; 204 204 case MaterialsMuWaterEnum : return "MaterialsMuWater"; -
issm/trunk-jpl/src/c/shared/Enum/StringToEnumx.cpp
r17941 r17946 203 203 else if (strcmp(name,"DamageEvolutionRequestedOutputs")==0) return DamageEvolutionRequestedOutputsEnum; 204 204 else if (strcmp(name,"MaterialsRhoIce")==0) return MaterialsRhoIceEnum; 205 else if (strcmp(name,"MaterialsRho Water")==0) return MaterialsRhoWaterEnum;205 else if (strcmp(name,"MaterialsRhoSeawater")==0) return MaterialsRhoSeawaterEnum; 206 206 else if (strcmp(name,"MaterialsRhoFreshwater")==0) return MaterialsRhoFreshwaterEnum; 207 207 else if (strcmp(name,"MaterialsMuWater")==0) return MaterialsMuWaterEnum; -
issm/trunk-jpl/src/m/classes/matdamageice.m
r17908 r17946 174 174 WriteData(fid,'enum',MaterialsEnum(),'data',MatdamageiceEnum(),'format','Integer'); 175 175 WriteData(fid,'object',obj,'class','materials','fieldname','rho_ice','format','Double'); 176 WriteData(fid,'object',obj,'class','materials','fieldname','rho_water',' format','Double');176 WriteData(fid,'object',obj,'class','materials','fieldname','rho_water','enum',MaterialsRhoSeawaterEnum(),'format','Double'); 177 177 WriteData(fid,'object',obj,'class','materials','fieldname','rho_freshwater','format','Double'); 178 178 WriteData(fid,'object',obj,'class','materials','fieldname','mu_water','format','Double'); -
issm/trunk-jpl/src/m/classes/matdamageice.py
r17765 r17946 129 129 WriteData(fid,'enum',MaterialsEnum(),'data',MatdamageiceEnum(),'format','Integer'); 130 130 WriteData(fid,'object',self,'class','materials','fieldname','rho_ice','format','Double') 131 WriteData(fid,'object',self,'class','materials','fieldname','rho_water',' format','Double')131 WriteData(fid,'object',self,'class','materials','fieldname','rho_water','enum',MaterialsRhoSeawaterEnum(),'format','Double') 132 132 WriteData(fid,'object',self,'class','materials','fieldname','rho_freshwater','format','Double') 133 133 WriteData(fid,'object',self,'class','materials','fieldname','mu_water','format','Double') -
issm/trunk-jpl/src/m/classes/matice.m
r17931 r17946 30 30 end 31 31 methods 32 function createxml(obj,fid) % {{{33 fprintf(fid, '\n\n');34 fprintf(fid, '<!-- materials -->\n');35 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rho_ice" type="',class(obj.rho_ice),'" default="',convert2str(obj.rho_ice),'">',' <section name="materials" />',' <help> ice density [kg/m^3] </help>','</parameter>');36 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rho_water" type="',class(obj.rho_water),'" default="',convert2str(obj.rho_water),'">',' <section name="materials" />',' <help> ocean water density [kg/m^3] </help>','</parameter>');37 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rho_freshwater" type="',class(obj.rho_freshwater),'" default="',convert2str(obj.rho_freshwater),'">',' <section name="materials" />',' <help> fresh water density [kg/m^3] </help>','</parameter>');38 39 40 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="mu_water" type="',class(obj.mu_water),'" default="',convert2str(obj.mu_water),'">',' <section name="materials" />',' <help> water viscosity [N s/m^2] </help>','</parameter>');41 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="heatcapacity" type="',class(obj.heatcapacity),'" default="',convert2str(obj.heatcapacity),'">',' <section name="materials" />',' <help> heat capacity [J/kg/K] </help>','</parameter>');42 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="latentheat" type="',class(obj.latentheat),'" default="',convert2str(obj.latentheat),'">',' <section name="materials" />',' <help> latent heat of fusion [J/kg] </help>','</parameter>');43 44 45 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="thermalconductivity" type="',class(obj.thermalconductivity),'" default="',convert2str(obj.thermalconductivity),'">',' <section name="materials" />',' <help> ice thermal conductivity [W/m/K] </help>','</parameter>');46 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="temperateiceconductivity" type="',class(obj.temperateiceconductivity),'" default="',convert2str(obj.temperateiceconductivity),'">',' <section name="materials" />',' <help> temperate ice thermal conductivity [W/m/K] </help>','</parameter>');47 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="meltingpoint" type="',class(obj.meltingpoint),'" default="',convert2str(obj.meltingpoint),'">',' <section name="materials" />',' <help> melting point of ice at 1atm in K </help>','</parameter>');48 49 50 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="beta" type="',class(obj.beta),'" default="',convert2str(obj.beta),'">',' <section name="materials" />',' <help> rate of change of melting point with pressure [K/Pa] </help>','</parameter>');51 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="mixed_layer_capacity" type="',class(obj.mixed_layer_capacity),'" default="',convert2str(obj.mixed_layer_capacity),'">',' <section name="materials" />',' <help> mixed layer capacity [W/kg/K] </help>','</parameter>');52 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="thermal_exchange_velocity" type="',class(obj.thermal_exchange_velocity),'" default="',convert2str(obj.thermal_exchange_velocity),'">',' <section name="materials" />',' <help> thermal exchange velocity [m/s] </help>','</parameter>');53 54 55 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rheology_B" type="',class(obj.rheology_B),'" default="',convert2str(obj.rheology_B),'">',' <section name="materials" />',' <help> flow law parameter [Pa/s^(1/n)] </help>','</parameter>');56 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rheology_n" type="',class(obj.rheology_n),'" default="',convert2str(obj.rheology_n),'">',' <section name="materials" />',' <help> Glens flow law exponent </help>','</parameter>');57 58 % rheology_law drop-down59 fprintf(fid,'%s\n%s\n%s\n%s\n','<parameter key ="rheology_law" type="alternative" optional="false">',' <section name="materials" />',' <help> law for the temperature dependance of the rheology: "None", "Paterson", "Arrhenius" or "LliboutryDuval" </help>');60 fprintf(fid,'%s\n',' <option value="None" type="string" default="true"> </option>');61 fprintf(fid,'%s\n',' <option value="Paterson" type="string" default="false"> </option>');62 fprintf(fid,'%s\n',' <option value="Arrhenius" type="string" default="false"> </option>');63 fprintf(fid,'%s\n%s\n',' <option value="LliboutryDuval" type="string" default="false"> </option>','</parameter>');64 65 66 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="lithosphere_shear_modulus" type="',class(obj.lithosphere_shear_modulus),'" default="',convert2str(obj.lithosphere_shear_modulus),'">',' <section name="materials" />',' <help> Lithosphere shear modulus [Pa] </help>','</parameter>');67 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="lithosphere_density" type="',class(obj.lithosphere_density),'" default="',convert2str(obj.lithosphere_density),'">',' <section name="materials" />',' <help> Lithosphere density [g/cm^-3] </help>','</parameter>');68 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="mantle_shear_modulus" type="',class(obj.mantle_shear_modulus),'" default="',convert2str(obj.mantle_shear_modulus),'">',' <section name="materials" />',' <help> Mantle shear modulus [Pa] </help>','</parameter>');69 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="mantle_density" type="',class(obj.mantle_density),'" default="',convert2str(obj.mantle_density),'">',' <section name="materials" />',' <help> Mantle density [g/cm^-3] </help>','</parameter>');70 71 72 end % }}}32 function createxml(obj,fid) % {{{ 33 fprintf(fid, '\n\n'); 34 fprintf(fid, '<!-- materials -->\n'); 35 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rho_ice" type="',class(obj.rho_ice),'" default="',convert2str(obj.rho_ice),'">',' <section name="materials" />',' <help> ice density [kg/m^3] </help>','</parameter>'); 36 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rho_water" type="',class(obj.rho_water),'" default="',convert2str(obj.rho_water),'">',' <section name="materials" />',' <help> ocean water density [kg/m^3] </help>','</parameter>'); 37 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rho_freshwater" type="',class(obj.rho_freshwater),'" default="',convert2str(obj.rho_freshwater),'">',' <section name="materials" />',' <help> fresh water density [kg/m^3] </help>','</parameter>'); 38 39 40 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="mu_water" type="',class(obj.mu_water),'" default="',convert2str(obj.mu_water),'">',' <section name="materials" />',' <help> water viscosity [N s/m^2] </help>','</parameter>'); 41 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="heatcapacity" type="',class(obj.heatcapacity),'" default="',convert2str(obj.heatcapacity),'">',' <section name="materials" />',' <help> heat capacity [J/kg/K] </help>','</parameter>'); 42 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="latentheat" type="',class(obj.latentheat),'" default="',convert2str(obj.latentheat),'">',' <section name="materials" />',' <help> latent heat of fusion [J/kg] </help>','</parameter>'); 43 44 45 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="thermalconductivity" type="',class(obj.thermalconductivity),'" default="',convert2str(obj.thermalconductivity),'">',' <section name="materials" />',' <help> ice thermal conductivity [W/m/K] </help>','</parameter>'); 46 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="temperateiceconductivity" type="',class(obj.temperateiceconductivity),'" default="',convert2str(obj.temperateiceconductivity),'">',' <section name="materials" />',' <help> temperate ice thermal conductivity [W/m/K] </help>','</parameter>'); 47 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="meltingpoint" type="',class(obj.meltingpoint),'" default="',convert2str(obj.meltingpoint),'">',' <section name="materials" />',' <help> melting point of ice at 1atm in K </help>','</parameter>'); 48 49 50 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="beta" type="',class(obj.beta),'" default="',convert2str(obj.beta),'">',' <section name="materials" />',' <help> rate of change of melting point with pressure [K/Pa] </help>','</parameter>'); 51 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="mixed_layer_capacity" type="',class(obj.mixed_layer_capacity),'" default="',convert2str(obj.mixed_layer_capacity),'">',' <section name="materials" />',' <help> mixed layer capacity [W/kg/K] </help>','</parameter>'); 52 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="thermal_exchange_velocity" type="',class(obj.thermal_exchange_velocity),'" default="',convert2str(obj.thermal_exchange_velocity),'">',' <section name="materials" />',' <help> thermal exchange velocity [m/s] </help>','</parameter>'); 53 54 55 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rheology_B" type="',class(obj.rheology_B),'" default="',convert2str(obj.rheology_B),'">',' <section name="materials" />',' <help> flow law parameter [Pa/s^(1/n)] </help>','</parameter>'); 56 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="rheology_n" type="',class(obj.rheology_n),'" default="',convert2str(obj.rheology_n),'">',' <section name="materials" />',' <help> Glens flow law exponent </help>','</parameter>'); 57 58 % rheology_law drop-down 59 fprintf(fid,'%s\n%s\n%s\n%s\n','<parameter key ="rheology_law" type="alternative" optional="false">',' <section name="materials" />',' <help> law for the temperature dependance of the rheology: "None", "Paterson", "Arrhenius" or "LliboutryDuval" </help>'); 60 fprintf(fid,'%s\n',' <option value="None" type="string" default="true"> </option>'); 61 fprintf(fid,'%s\n',' <option value="Paterson" type="string" default="false"> </option>'); 62 fprintf(fid,'%s\n',' <option value="Arrhenius" type="string" default="false"> </option>'); 63 fprintf(fid,'%s\n%s\n',' <option value="LliboutryDuval" type="string" default="false"> </option>','</parameter>'); 64 65 66 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="lithosphere_shear_modulus" type="',class(obj.lithosphere_shear_modulus),'" default="',convert2str(obj.lithosphere_shear_modulus),'">',' <section name="materials" />',' <help> Lithosphere shear modulus [Pa] </help>','</parameter>'); 67 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="lithosphere_density" type="',class(obj.lithosphere_density),'" default="',convert2str(obj.lithosphere_density),'">',' <section name="materials" />',' <help> Lithosphere density [g/cm^-3] </help>','</parameter>'); 68 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="mantle_shear_modulus" type="',class(obj.mantle_shear_modulus),'" default="',convert2str(obj.mantle_shear_modulus),'">',' <section name="materials" />',' <help> Mantle shear modulus [Pa] </help>','</parameter>'); 69 fprintf(fid,'%s%s%s%s%s\n%s\n%s\n%s\n','<parameter key ="mantle_density" type="',class(obj.mantle_density),'" default="',convert2str(obj.mantle_density),'">',' <section name="materials" />',' <help> Mantle density [g/cm^-3] </help>','</parameter>'); 70 71 72 end % }}} 73 73 function obj = matice(varargin) % {{{ 74 74 switch nargin … … 181 181 WriteData(fid,'enum',MaterialsEnum(),'data',MaticeEnum(),'format','Integer'); 182 182 WriteData(fid,'object',obj,'class','materials','fieldname','rho_ice','format','Double'); 183 WriteData(fid,'object',obj,'class','materials','fieldname','rho_water',' format','Double');183 WriteData(fid,'object',obj,'class','materials','fieldname','rho_water','enum',MaterialsRhoSeawaterEnum(),'format','Double'); 184 184 WriteData(fid,'object',obj,'class','materials','fieldname','rho_freshwater','format','Double'); 185 185 WriteData(fid,'object',obj,'class','materials','fieldname','mu_water','format','Double'); -
issm/trunk-jpl/src/m/classes/matice.py
r17497 r17946 129 129 WriteData(fid,'enum',MaterialsEnum(),'data',MaticeEnum(),'format','Integer'); 130 130 WriteData(fid,'object',self,'class','materials','fieldname','rho_ice','format','Double') 131 WriteData(fid,'object',self,'class','materials','fieldname','rho_water',' format','Double')131 WriteData(fid,'object',self,'class','materials','fieldname','rho_water','enum',MaterialsRhoSeawaterEnum(),'format','Double') 132 132 WriteData(fid,'object',self,'class','materials','fieldname','rho_freshwater','format','Double') 133 133 WriteData(fid,'object',self,'class','materials','fieldname','mu_water','format','Double') -
issm/trunk-jpl/src/m/enum/EnumDefinitions.py
r17941 r17946 192 192 def DamageEvolutionRequestedOutputsEnum(): return StringToEnum("DamageEvolutionRequestedOutputs")[0] 193 193 def MaterialsRhoIceEnum(): return StringToEnum("MaterialsRhoIce")[0] 194 def MaterialsRho WaterEnum(): return StringToEnum("MaterialsRhoWater")[0]194 def MaterialsRhoSeawaterEnum(): return StringToEnum("MaterialsRhoSeawater")[0] 195 195 def MaterialsRhoFreshwaterEnum(): return StringToEnum("MaterialsRhoFreshwater")[0] 196 196 def MaterialsMuWaterEnum(): return StringToEnum("MaterialsMuWater")[0] -
issm/trunk-jpl/src/m/qmu/setupdesign/IsScaled.m
r13646 r17946 3 3 4 4 switch variablename, 5 case {'MaterialsRhoIce','MaterialsRho Water','MaterialsHeatCapacity','MaterialsThermalConductivity','Gravity','MaxVel'},5 case {'MaterialsRhoIce','MaterialsRhoSeawater','MaterialsHeatCapacity','MaterialsThermalConductivity','Gravity','MaxVel'}, 6 6 7 7 status=0;
Note:
See TracChangeset
for help on using the changeset viewer.
![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)