Changeset 12301


Ignore:
Timestamp:
05/29/12 14:28:15 (13 years ago)
Author:
Mathieu Morlighem
Message:

Added md.materials.rho_freshwater

Location:
issm/trunk/src
Files:
2 added
8 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/EnumDefinitions/EnumDefinitions.h

    r12297 r12301  
    103103        MaterialsRhoIceEnum,
    104104        MaterialsRhoWaterEnum,
     105        MaterialsRhoFreshwaterEnum,
    105106        MaterialsMuWaterEnum,
    106107        MaterialsThermalExchangeVelocityEnum,
     
    160161        SurfaceforcingsMassBalanceEnum,
    161162        SurfaceforcingsIspddEnum,
    162         SurfaceforcingsMonthlytemperatures,
     163        SurfaceforcingsMonthlytemperaturesEnum,
    163164        ThermalMaxiterEnum,
    164165        ThermalPenaltyFactorEnum,
  • issm/trunk/src/c/modules/EnumToStringx/EnumToStringx.cpp

    r12297 r12301  
    107107                case MaterialsRhoIceEnum : return "MaterialsRhoIce";
    108108                case MaterialsRhoWaterEnum : return "MaterialsRhoWater";
     109                case MaterialsRhoFreshwaterEnum : return "MaterialsRhoFreshwater";
    109110                case MaterialsMuWaterEnum : return "MaterialsMuWater";
    110111                case MaterialsThermalExchangeVelocityEnum : return "MaterialsThermalExchangeVelocity";
     
    164165                case SurfaceforcingsMassBalanceEnum : return "SurfaceforcingsMassBalance";
    165166                case SurfaceforcingsIspddEnum : return "SurfaceforcingsIspdd";
     167                case SurfaceforcingsMonthlytemperaturesEnum : return "SurfaceforcingsMonthlytemperatures";
    166168                case ThermalMaxiterEnum : return "ThermalMaxiter";
    167169                case ThermalPenaltyFactorEnum : return "ThermalPenaltyFactor";
  • issm/trunk/src/c/modules/ModelProcessorx/Prognostic/UpdateElementsPrognostic.cpp

    r12297 r12301  
    7171          iomodel->FetchDataToInput(elements,ThermalSpctemperatureEnum);
    7272          iomodel->FetchDataToInput(elements,SurfaceforcingsPrecipitationEnum);
    73           iomodel->FetchDataToInput(elements,SurfaceforcingsMonthlytemperatures);
     73          iomodel->FetchDataToInput(elements,SurfaceforcingsMonthlytemperaturesEnum);
    7474        }
    7575        else{
  • issm/trunk/src/c/modules/StringToEnumx/StringToEnumx.cpp

    r12297 r12301  
    108108              else if (strcmp(name,"MaterialsRhoIce")==0) return MaterialsRhoIceEnum;
    109109              else if (strcmp(name,"MaterialsRhoWater")==0) return MaterialsRhoWaterEnum;
     110              else if (strcmp(name,"MaterialsRhoFreshwater")==0) return MaterialsRhoFreshwaterEnum;
    110111              else if (strcmp(name,"MaterialsMuWater")==0) return MaterialsMuWaterEnum;
    111112              else if (strcmp(name,"MaterialsThermalExchangeVelocity")==0) return MaterialsThermalExchangeVelocityEnum;
     
    136137              else if (strcmp(name,"PrognosticMinThickness")==0) return PrognosticMinThicknessEnum;
    137138              else if (strcmp(name,"PrognosticPenaltyFactor")==0) return PrognosticPenaltyFactorEnum;
    138               else if (strcmp(name,"PrognosticSpcthickness")==0) return PrognosticSpcthicknessEnum;
    139139         else stage=2;
    140140   }
    141141   if(stage==2){
    142               if (strcmp(name,"PrognosticStabilization")==0) return PrognosticStabilizationEnum;
     142              if (strcmp(name,"PrognosticSpcthickness")==0) return PrognosticSpcthicknessEnum;
     143              else if (strcmp(name,"PrognosticStabilization")==0) return PrognosticStabilizationEnum;
    143144              else if (strcmp(name,"PrognosticVertexPairing")==0) return PrognosticVertexPairingEnum;
    144145              else if (strcmp(name,"QmuIsdakota")==0) return QmuIsdakotaEnum;
     
    168169              else if (strcmp(name,"SurfaceforcingsMassBalance")==0) return SurfaceforcingsMassBalanceEnum;
    169170              else if (strcmp(name,"SurfaceforcingsIspdd")==0) return SurfaceforcingsIspddEnum;
     171              else if (strcmp(name,"SurfaceforcingsMonthlytemperatures")==0) return SurfaceforcingsMonthlytemperaturesEnum;
    170172              else if (strcmp(name,"ThermalMaxiter")==0) return ThermalMaxiterEnum;
    171173              else if (strcmp(name,"ThermalPenaltyFactor")==0) return ThermalPenaltyFactorEnum;
     
    258260              else if (strcmp(name,"IntInput")==0) return IntInputEnum;
    259261              else if (strcmp(name,"IntParam")==0) return IntParamEnum;
    260               else if (strcmp(name,"IntVecParam")==0) return IntVecParamEnum;
    261               else if (strcmp(name,"MacAyeal2dIceFront")==0) return MacAyeal2dIceFrontEnum;
    262262         else stage=3;
    263263   }
    264264   if(stage==3){
    265               if (strcmp(name,"MacAyeal3dIceFront")==0) return MacAyeal3dIceFrontEnum;
     265              if (strcmp(name,"IntVecParam")==0) return IntVecParamEnum;
     266              else if (strcmp(name,"MacAyeal2dIceFront")==0) return MacAyeal2dIceFrontEnum;
     267              else if (strcmp(name,"MacAyeal3dIceFront")==0) return MacAyeal3dIceFrontEnum;
    266268              else if (strcmp(name,"Matice")==0) return MaticeEnum;
    267269              else if (strcmp(name,"Matpar")==0) return MatparEnum;
     
    381383              else if (strcmp(name,"StressTensorzz")==0) return StressTensorzzEnum;
    382384              else if (strcmp(name,"IceVolume")==0) return IceVolumeEnum;
    383               else if (strcmp(name,"P0")==0) return P0Enum;
    384               else if (strcmp(name,"P1")==0) return P1Enum;
    385385         else stage=4;
    386386   }
    387387   if(stage==4){
    388               if (strcmp(name,"P1DG")==0) return P1DGEnum;
     388              if (strcmp(name,"P0")==0) return P0Enum;
     389              else if (strcmp(name,"P1")==0) return P1Enum;
     390              else if (strcmp(name,"P1DG")==0) return P1DGEnum;
    389391              else if (strcmp(name,"SaveResults")==0) return SaveResultsEnum;
    390392              else if (strcmp(name,"BoolElementResult")==0) return BoolElementResultEnum;
  • issm/trunk/src/c/objects/Elements/Tria.cpp

    r12300 r12301  
    23162316   /*Get material parameters :*/
    23172317   rho_ice=matpar->GetRhoIce();
    2318    //rho_freshwater=matpar->GetRhoWater();
     2318   rho_water=matpar->GetRhoFreshwater();
    23192319   
    2320    sconv=(1000/rho_ice)/12.; //rhow_rain/rhoi / 12 months
     2320   sconv=(rho_water/rho_ice)/12.; //rhow_rain/rhoi / 12 months
    23212321     
    23222322     /*PDD constant*/
  • issm/trunk/src/c/objects/Materials/Matpar.cpp

    r11527 r12301  
    2525Matpar::Matpar(int matpar_mid, IoModel* iomodel){
    2626
    27         this->mid                       = matpar_mid;
     27        this->mid = matpar_mid;
    2828        iomodel->Constant(&this->rho_ice,MaterialsRhoIceEnum);
    2929        iomodel->Constant(&this->rho_water,MaterialsRhoWaterEnum);
     30        iomodel->Constant(&this->rho_freshwater,MaterialsRhoFreshwaterEnum);
    3031        iomodel->Constant(&this->mu_water,MaterialsMuWaterEnum);
    3132        iomodel->Constant(&this->heatcapacity,MaterialsHeatcapacityEnum);
     
    6061        printf("   rho_ice: %g\n",rho_ice);
    6162        printf("   rho_water: %g\n",rho_water);
     63        printf("   rho_freshwater: %g\n",rho_freshwater);
    6264        printf("   mu_water: %g\n",mu_water);
    6365        printf("   heatcapacity: %g\n",heatcapacity);
     
    7678void Matpar::DeepEcho(void){
    7779
    78         printf("Matpar:\n");
    79         printf("   mid: %i\n",mid);
    80         printf("   rho_ice: %g\n",rho_ice);
    81         printf("   rho_water: %g\n",rho_water);
    82         printf("   mu_water: %g\n",mu_water);
    83         printf("   heatcapacity: %g\n",heatcapacity);
    84         printf("   thermalconductivity: %g\n",thermalconductivity);
    85         printf("   latentheat: %g\n",latentheat);
    86         printf("   beta: %g\n",beta);
    87         printf("   meltingpoint: %g\n",meltingpoint);
    88         printf("   referencetemperature: %g\n",referencetemperature);
    89         printf("   mixed_layer_capacity: %g\n",mixed_layer_capacity);
    90         printf("   thermal_exchange_velocity: %g\n",thermal_exchange_velocity);
    91         printf("   g: %g\n",g);
    92         return;
     80        this->Echo();
    9381}               
    9482/*}}}1*/
     
    122110        memcpy(marshalled_dataset,&rho_ice,sizeof(rho_ice));marshalled_dataset+=sizeof(rho_ice);
    123111        memcpy(marshalled_dataset,&rho_water,sizeof(rho_water));marshalled_dataset+=sizeof(rho_water);
     112        memcpy(marshalled_dataset,&rho_freshwater,sizeof(rho_freshwater));marshalled_dataset+=sizeof(rho_freshwater);
    124113        memcpy(marshalled_dataset,&mu_water,sizeof(mu_water));marshalled_dataset+=sizeof(mu_water);
    125114        memcpy(marshalled_dataset,&heatcapacity,sizeof(heatcapacity));marshalled_dataset+=sizeof(heatcapacity);
     
    143132                sizeof(rho_ice)+
    144133                sizeof(rho_water)+
     134                sizeof(rho_freshwater)+
    145135                sizeof(mu_water)+
    146136                sizeof(heatcapacity)+
     
    170160        memcpy(&rho_ice,marshalled_dataset,sizeof(rho_ice));marshalled_dataset+=sizeof(rho_ice);
    171161        memcpy(&rho_water,marshalled_dataset,sizeof(rho_water));marshalled_dataset+=sizeof(rho_water);
     162        memcpy(&rho_freshwater,marshalled_dataset,sizeof(rho_freshwater));marshalled_dataset+=sizeof(rho_freshwater);
    172163        memcpy(&mu_water,marshalled_dataset,sizeof(mu_water));marshalled_dataset+=sizeof(mu_water);
    173164        memcpy(&heatcapacity,marshalled_dataset,sizeof(heatcapacity));marshalled_dataset+=sizeof(heatcapacity);
     
    243234                        this->rho_ice=constant;
    244235                        break;
    245 
    246236                case MaterialsRhoWaterEnum:
    247237                        this->rho_water=constant;
    248238                        break;
    249 
     239                case MaterialsRhoFreshwaterEnum:
     240                        this->rho_freshwater=constant;
     241                        break;
    250242                case MaterialsMuWaterEnum:
    251243                        this->mu_water=constant;
    252244                        break;
    253 
    254245                case MaterialsHeatcapacityEnum:
    255246                        this->heatcapacity=constant;
    256247                        break;
    257 
    258248                case MaterialsThermalconductivityEnum:
    259249                        this->thermalconductivity=constant;
    260250                        break;
    261 
    262251                case  MaterialsLatentheatEnum:
    263252                        this->latentheat=constant;
    264253                        break;
    265 
    266254                case  MaterialsBetaEnum:
    267255                        this->beta=constant;
    268256                        break;
    269 
    270257                case  MaterialsMeltingpointEnum:
    271258                        this->meltingpoint=constant;
    272259                        break;
    273 
    274260                case  ConstantsReferencetemperatureEnum:
    275261                        this->referencetemperature=constant;
    276262                        break;
    277 
    278263                case  MaterialsMixedLayerCapacityEnum:
    279264                        this->mixed_layer_capacity=constant;
    280265                        break;
    281 
    282266                case  MaterialsThermalExchangeVelocityEnum:
    283267                        this->thermalconductivity=constant;
    284268                        break;
    285 
    286269                case  ConstantsGEnum:
    287270                        this->g=constant;
    288271                        break;
    289 
    290272                default:
    291273                        break;
     
    362344double Matpar::GetRhoWater(){
    363345        return rho_water;
     346}
     347/*}}}1*/
     348/*FUNCTION Matpar::GetRhoFreshwater {{{1*/
     349double Matpar::GetRhoFreshwater(){
     350        return rho_freshwater;
    364351}
    365352/*}}}1*/
  • issm/trunk/src/c/objects/Materials/Matpar.h

    r11995 r12301  
    1515
    1616        private:
    17                 int     mid;
    18                 double  rho_ice;
    19                 double  rho_water;
     17                int       mid;
     18                double  rho_ice;
     19                double  rho_water;
     20                double  rho_freshwater;
    2021                double  mu_water;
    2122                double  heatcapacity;
     
    7778                double GetRhoIce();
    7879                double GetRhoWater();
     80                double GetRhoFreshwater();
    7981                double GetMuWater();
    8082                double GetMixedLayerCapacity();
  • issm/trunk/src/m/classes/materials.m

    r11995 r12301  
    88                rho_ice                    = 0;
    99                rho_water                  = 0;
     10                rho_freshwater             = 0;
    1011                mu_water                   = 0;
    1112                heatcapacity               = 0;
     
    3435                        obj.rho_ice=917;
    3536
    36                         %water density (kg/m^3)
     37                        %ocean water density (kg/m^3)
    3738                        obj.rho_water=1023;
     39
     40                        %fresh water density (kg/m^3)
     41                        obj.rho_freshwater=1000;
    3842
    3943                        %water viscosity (N.s/m^2)
     
    6872                        checkfield(md,'materials.rho_ice','>',0);
    6973                        checkfield(md,'materials.rho_water','>',0);
     74                        checkfield(md,'materials.rho_freshwater','>',0);
    7075                        checkfield(md,'materials.mu_water','>',0);
    7176                        checkfield(md,'materials.rheology_B','>',0,'size',[md.mesh.numberofvertices 1]);
     
    7782
    7883                        fielddisplay(obj,'rho_ice','ice density [kg/m^3]');
    79                         fielddisplay(obj,'rho_water','water density [kg/m^3]');
     84                        fielddisplay(obj,'rho_water','ocean water density [kg/m^3]');
     85                        fielddisplay(obj,'freshrho_water','fresh water density [kg/m^3]');
    8086                        fielddisplay(obj,'mu_water','water viscosity [N s/m^2]');
    8187                        fielddisplay(obj,'heatcapacity','heat capacity [J/kg/K]');
     
    9399                        WriteData(fid,'object',obj,'fieldname','rho_ice','format','Double');
    94100                        WriteData(fid,'object',obj,'fieldname','rho_water','format','Double');
     101                        WriteData(fid,'object',obj,'fieldname','rho_freshwater','format','Double');
    95102                        WriteData(fid,'object',obj,'fieldname','mu_water','format','Double');
    96103                        WriteData(fid,'object',obj,'fieldname','heatcapacity','format','Double');
Note: See TracChangeset for help on using the changeset viewer.