Ignore:
Timestamp:
04/12/16 21:32:01 (9 years ago)
Author:
Mathieu Morlighem
Message:

merged trunk-jpl and trunk for revision 20497

Location:
issm/trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk

  • issm/trunk/src

  • issm/trunk/src/c/classes/Materials/Matpar.cpp

    r19105 r20500  
    1919Matpar::Matpar(int matpar_mid, IoModel* iomodel){/*{{{*/
    2020
     21        rho_ice=0;
     22        rho_water=0;
     23        rho_freshwater=0;
     24        mu_water=0;
     25        heatcapacity=0;
     26        thermalconductivity=0;
     27        temperateiceconductivity=0;
     28        latentheat=0;
     29        beta=0;
     30        meltingpoint=0;
     31        referencetemperature=0;
     32        mixed_layer_capacity=0;
     33        thermal_exchange_velocity=0;
     34        g=0;
     35        omega=0;
     36        desfac=0;
     37        rlaps=0;
     38        rlapslgm=0;
     39        dpermil=0;
     40
     41        albedo_snow=0;
     42        albedo_ice=0;
     43
     44        sediment_compressibility=0;
     45        sediment_porosity=0;
     46        sediment_thickness=0;
     47        water_compressibility=0;
     48
     49        epl_compressibility=0;
     50        epl_porosity=0;
     51        epl_init_thickness=0;
     52        epl_colapse_thickness=0;
     53        epl_max_thickness=0;
     54        epl_conductivity=0;
     55
     56        lithosphere_shear_modulus=0;
     57        lithosphere_density=0;
     58        mantle_shear_modulus=0;
     59        mantle_density=0;
     60       
     61        earth_density=0;
     62
     63        poisson=0;
     64        young_modulus=0;
     65        ridging_exponent=0;
     66        cohesion=0;
     67        internal_friction_coef=0;
     68        compression_coef=0;
     69        traction_coef=0;
     70        time_relaxation_stress=0;
     71        time_relaxation_damage=0;
     72
    2173        bool isefficientlayer;
    2274        int  hydrology_model,smb_model,materials_type;
    2375        iomodel->Constant(&hydrology_model,HydrologyModelEnum);
    24         iomodel->Constant(&smb_model,SurfaceforcingsEnum);
     76        iomodel->Constant(&smb_model,SmbEnum);
    2577        iomodel->Constant(&materials_type,MaterialsEnum);
    2678
    … …  
    4698
    4799                        switch(smb_model){
    48                                 case SMBEnum:
     100                                case SMBforcingEnum:
    49101                                        /*Nothing to add*/
    50102                                        break;
     103                                case SMBgembEnum:
     104                                        iomodel->Constant(&this->albedo_ice,SmbAIceEnum);
     105                                        iomodel->Constant(&this->albedo_snow,SmbASnowEnum);
     106                                        break;
    51107                                case SMBpddEnum:
    52                                         iomodel->Constant(&this->desfac,SurfaceforcingsDesfacEnum);
    53                                         iomodel->Constant(&this->s0p,SurfaceforcingsS0pEnum);
    54                                         iomodel->Constant(&this->s0t,SurfaceforcingsS0tEnum);
    55                                         iomodel->Constant(&this->rlaps,SurfaceforcingsRlapsEnum);
    56                                         iomodel->Constant(&this->rlapslgm,SurfaceforcingsRlapslgmEnum);
    57                                         break;
     108                                        iomodel->Constant(&this->desfac,SmbDesfacEnum);
     109                                        iomodel->Constant(&this->rlaps,SmbRlapsEnum);
     110                                        iomodel->Constant(&this->rlapslgm,SmbRlapslgmEnum);
     111                                        break;
     112                                case SMBd18opddEnum:
     113                                        iomodel->Constant(&this->desfac,SmbDesfacEnum);
     114                                        iomodel->Constant(&this->rlaps,SmbRlapsEnum);
     115                                        iomodel->Constant(&this->rlapslgm,SmbRlapslgmEnum);
     116                                        iomodel->Constant(&this->dpermil,SmbDpermilEnum);                                       
    58117                                case SMBgradientsEnum:
     118                                        break;
    59119                                        /*Nothing to add*/
    60120                                        break;
    … …  
    90150                                /*Nothing to add*/
    91151                        }
     152                        else if(hydrology_model==HydrologysommersEnum){
     153                                /*Nothing to add*/
     154                        }
    92155                        else{
    93156                                _error_("Hydrology model "<<EnumToStringx(hydrology_model)<<" not supported yet");
    … …  
    99162                        iomodel->Constant(&this->mantle_shear_modulus,MaterialsMantleShearModulusEnum);
    100163                        iomodel->Constant(&this->mantle_density,MaterialsMantleDensityEnum);
     164
     165                        /*slr:*/
     166                        iomodel->Constant(&this->earth_density,MaterialsEarthDensityEnum);
     167
    101168                        break;
    102169                default:
    … …  
    130197        _printf_("   g: " << g << "\n");
    131198        _printf_("   desfac: " << desfac << "\n");
    132         _printf_("   s0p: " << s0p << "\n");
    133         _printf_("   s0t: " << s0t << "\n");
    134199        _printf_("   rlaps: " << rlaps << "\n");
    135200        _printf_("   rlapslgm: " << rlapslgm << "\n");
     201        _printf_("   dpermil: " << dpermil << "\n");
    136202        return;
    137203}
    … …  
    142208}               
    143209/*}}}*/
    144 int    Matpar::Id(void){ return mid; }/*{{{*/
    145 /*}}}*/
    146 int Matpar::ObjectEnum(void){/*{{{*/
     210int  Matpar::Id(void){ return mid; }/*{{{*/
     211/*}}}*/
     212int  Matpar::ObjectEnum(void){/*{{{*/
    147213
    148214        return MatparEnum;
    … …  
    175241        matpar->g=this->g;
    176242        matpar->desfac=this->desfac;
    177         matpar->s0p=this->s0p;
    178         matpar->s0t=this->s0t;
    179243        matpar->rlaps=this->rlaps;
    180244        matpar->rlapslgm=this->rlapslgm;
     245        matpar->dpermil=this->dpermil;
    181246
    182247        matpar->sediment_compressibility=this->sediment_compressibility;
    … …  
    196261        matpar->mantle_shear_modulus=this->mantle_shear_modulus;
    197262        matpar->mantle_density=this->mantle_density;
     263       
     264        matpar->earth_density=this->earth_density;
    198265
    199266        return matpar;
     267}
     268/*}}}*/
     269void Matpar::Marshall(char** pmarshalled_data,int* pmarshalled_data_size, int marshall_direction){ /*{{{*/
     270
     271        MARSHALLING_ENUM(MatparEnum);
     272
     273        MARSHALLING(mid);
     274        MARSHALLING(rho_ice);
     275        MARSHALLING(rho_water);
     276        MARSHALLING(rho_freshwater);
     277        MARSHALLING(mu_water);
     278        MARSHALLING(heatcapacity);
     279        MARSHALLING(thermalconductivity);
     280        MARSHALLING(temperateiceconductivity);
     281        MARSHALLING(latentheat);
     282        MARSHALLING(beta);
     283        MARSHALLING(meltingpoint);
     284        MARSHALLING(referencetemperature);
     285        MARSHALLING(mixed_layer_capacity);
     286        MARSHALLING(thermal_exchange_velocity);
     287        MARSHALLING(g);
     288        MARSHALLING(omega);
     289        MARSHALLING(desfac);
     290        MARSHALLING(rlaps);
     291        MARSHALLING(rlapslgm);
     292        MARSHALLING(dpermil);
     293
     294        //hydrology Dual Porous Continuum:
     295        MARSHALLING(sediment_compressibility);
     296        MARSHALLING(sediment_porosity);
     297        MARSHALLING(sediment_thickness);
     298        MARSHALLING(water_compressibility);
     299
     300        MARSHALLING(epl_compressibility);
     301        MARSHALLING(epl_porosity);
     302        MARSHALLING(epl_init_thickness);
     303        MARSHALLING(epl_colapse_thickness);
     304        MARSHALLING(epl_max_thickness);
     305        MARSHALLING(epl_conductivity);
     306
     307        //gia:
     308        MARSHALLING(lithosphere_shear_modulus);
     309        MARSHALLING(lithosphere_density);
     310        MARSHALLING(mantle_shear_modulus);
     311        MARSHALLING(mantle_density);
     312       
     313        //slr:
     314        MARSHALLING(earth_density);
     315
     316        //Sea ice:
     317        MARSHALLING(poisson);
     318        MARSHALLING(young_modulus);
     319        MARSHALLING(ridging_exponent);
     320        MARSHALLING(cohesion);
     321        MARSHALLING(internal_friction_coef);
     322        MARSHALLING(compression_coef);
     323        MARSHALLING(traction_coef);
     324        MARSHALLING(time_relaxation_stress);
     325        MARSHALLING(time_relaxation_damage);
     326
    200327}
    201328/*}}}*/
    … …  
    259386                        this->g=constant;
    260387                        break;
    261                 case  SurfaceforcingsDesfacEnum:
     388                case  SmbDesfacEnum:
    262389                        this->desfac=constant;
    263390                        break;
    264                 case SurfaceforcingsS0pEnum:
    265                         this->s0p=constant;
    266                         break;
    267                 case SurfaceforcingsS0tEnum:
    268                         this->s0t=constant;
    269                         break;
    270                 case SurfaceforcingsRlapsEnum:
     391                case SmbRlapsEnum:
    271392                        this->rlaps=constant;
    272393                        break;
    273                 case SurfaceforcingsRlapslgmEnum:
     394                case SmbRlapslgmEnum:
    274395                        this->rlapslgm=constant;
     396                        break;
     397                case  SmbDpermilEnum:
     398                        this->dpermil=constant;
    275399                        break;
    276400                default:
    … …  
    395519                case HydrologydcWaterCompressibilityEnum:    return this->water_compressibility;
    396520                case ConstantsGEnum:                         return this->g;
    397                 case SurfaceforcingsDesfacEnum:              return this->desfac;
    398                 case SurfaceforcingsS0pEnum:                 return this->s0p;
    399                 case SurfaceforcingsS0tEnum:                 return this->s0t;
    400                 case SurfaceforcingsRlapsEnum:               return this->rlaps;
    401                 case SurfaceforcingsRlapslgmEnum:            return this->rlapslgm;
     521                case SmbDesfacEnum:              return this->desfac;
     522                case SmbRlapsEnum:               return this->rlaps;
     523                case SmbRlapslgmEnum:            return this->rlapslgm;
     524                case SmbDpermilEnum:             return this->dpermil;
    402525                case MaterialsLithosphereShearModulusEnum:   return this->lithosphere_shear_modulus;
    403526                case MaterialsLithosphereDensityEnum:        return this->lithosphere_density;
    404527                case MaterialsMantleDensityEnum:             return this->mantle_density;
    405528                case MaterialsMantleShearModulusEnum:        return this->mantle_shear_modulus;
     529                case MaterialsEarthDensityEnum:              return this->earth_density;
    406530                default: _error_("Enum "<<EnumToStringx(enum_in)<<" not supported yet");
    407531        }
Note: See TracChangeset for help on using the changeset viewer.