Changeset 17946
- Timestamp:
- 05/06/14 08:30:09 (11 years ago)
- Location:
- issm/trunk-jpl/src
- Files:
-
- 1 added
- 1 deleted
- 16 edited
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 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 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.