Changeset 24750


Ignore:
Timestamp:
04/27/20 17:20:42 (5 years ago)
Author:
Eric.Larour
Message:

CHG: simplified test2002.m by skipping fields not stricly necessary for the slr solution.
In addition, created new 'hydro' material in the materials.m class. This triggered a series
of modifications to the model processor. The biggest was the fact that hmaterials hoooks in
the elements were pointing to materials that did not exist! So we did a post processing of the
hmaterials hooks to nullify them. This in trun triggered Tria.cpp modification of the Configure
routine to skip whenever hmaterial is NULL.

Location:
issm/trunk-jpl
Files:
12 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

    r24733 r24750  
    10841084        else this->hnodes = NULL;
    10851085        this->hvertices->configure(verticesin);
    1086         this->hmaterial->configure(materialsin);
     1086        if(this->hmaterial) this->hmaterial->configure(materialsin);
    10871087
    10881088        /*Now, go pick up the objects inside the hooks: */
    … …  
    10901090        else this->nodes=NULL;
    10911091        this->vertices = (Vertex**)this->hvertices->deliverp();
    1092         this->material = (Material*)this->hmaterial->delivers();
     1092        if(this->hmaterial)this->material = (Material*)this->hmaterial->delivers();
    10931093
    10941094        /*point parameters to real dataset: */

  • issm/trunk-jpl/src/c/classes/IoModel.cpp

    r24579 r24750  
    644644                                        if(strcmp(record_name,"md.hydrology.model")==0) integer = IoCodeToEnumHydrology(integer);
    645645                                        if(strcmp(record_name,"md.materials.type")==0) integer = IoCodeToEnumMaterials(integer);
     646                                        if(strcmp(record_name,"md.materials.nature")==0) integer = IoCodeToEnumNature(integer);
    646647                                        if(strcmp(record_name,"md.timestepping.type")==0) integer = IoCodeToEnumTimestepping(integer);
    647648                                        if(strcmp(record_name,"md.amr.type")==0) integer = IoCodeToEnumAmr(integer);

  • issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp

    r24335 r24750  
    171171                        //make sure materials that are not tied to elements come last:  for now, only Matlitho qualifies.
    172172                        for(int i=0;i<nnat;i++){
    173                                 if (IoCodeToEnumMaterials(nature[i])==MatlithoEnum){
     173                                if (IoCodeToEnumNature(nature[i])==MatlithoEnum){
    174174                                        int temp=nature[nnat-1];
    175175                                        nature[nnat-1]=nature[i];
    … …  
    180180                        //go through list of materials, and create them:
    181181                        for(int i=0;i<nnat;i++){
    182                                 switch(IoCodeToEnumMaterials(nature[i])){ //{{{
     182                                switch(IoCodeToEnumNature(nature[i])){ //{{{
    183183                                        case MaticeEnum:
    184184                                                iomodel->FetchDataToInput(inputs2,elements,"md.materials.rheology_B",MaterialsRheologyBEnum);
    … …  
    199199                                                materials->AddObject(new Matlitho(materials->Size()+1,iomodel));
    200200                                                iomodel->DeleteData(9,"md.materials.radius","md.materials.viscosity","md.materials.lame_lambda","md.materials.lame_mu","md.materials.burgers_viscosity","md.materials.burgers_mu","md.materials.isburgers","md.materials.issolid","md.materials.density");
     201                                                break;
     202
     203                                        case MathydroEnum:
     204                                                {
     205                                                        /*If we don't have any materials pointed to by elements (meaning, if we are running only litho or hydro),
     206                                                         * then we need to zero out the hmaterial pointers inside the elements dataset so that it won't error out
     207                                                         * during configuration: */
     208                                                        bool isice=false;
     209                                                        for (int j=0;j<nnat;j++){
     210                                                                if((IoCodeToEnumNature(nature[i])==MaticeEnum)||
     211                                                                                (IoCodeToEnumNature(nature[i])==MatenhancediceEnum)||
     212                                                                                (IoCodeToEnumNature(nature[i])==MatestarEnum)||
     213                                                                                (IoCodeToEnumNature(nature[i])==MatdamageiceEnum)){
     214                                                                        isice=true; break; }
     215                                                        }
     216                                                        if (!isice){
     217                                                                /*go through elements, and zero the hmaterials pointers: */
     218                                                                for(int j=0;j<elements->Size();j++){
     219                                                                        Element* element=xDynamicCast<Element*>(elements->GetObjectByOffset(j));
     220                                                                        switch(element->ObjectEnum()){
     221                                                                                case TriaEnum:
     222                                                                                        {
     223                                                                                                Tria* tria= xDynamicCast<Tria*>(element);
     224                                                                                                tria->hmaterial=NULL;
     225                                                                                        }
     226                                                                                        break;
     227                                                                                default:
     228                                                                                        _error_("Not implemented yet");
     229                                                                        }
     230                                                                }
     231                                                        }
     232                                                }
    201233                                                break;
    202234
    … …  
    252284
    253285                                        default:
    254                                                 _error_("Materials "<<EnumToStringx(IoCodeToEnumMaterials(nature[i]))<<" not supported");
     286                                                _error_("Materials nature type "<<EnumToStringx(IoCodeToEnumNature(nature[i]))<<" not supported");
    255287
    256288                                } //}}}

  • issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp

    r24531 r24750  
    341341                        iomodel->FetchData(&nature,&nnat,&dummy,"md.materials.nature");
    342342                        for(int i=0;i<nnat;i++){
    343                                 switch(IoCodeToEnumMaterials(nature[i])){
     343                                switch(IoCodeToEnumNature(nature[i])){
    344344                                        case MatlithoEnum:
    345345                                                break;
    … …  
    374374                                                parameters->AddObject(iomodel->CopyConstantObject("md.materials.earth_density",MaterialsEarthDensityEnum));
    375375                                                break;
     376                                        case MathydroEnum:
     377                                                parameters->AddObject(iomodel->CopyConstantObject("md.materials.rho_ice",MaterialsRhoIceEnum));
     378                                                parameters->AddObject(iomodel->CopyConstantObject("md.materials.rho_water",MaterialsRhoSeawaterEnum));
     379                                                parameters->AddObject(iomodel->CopyConstantObject("md.materials.earth_density",MaterialsEarthDensityEnum));
     380                                                break;
    376381                                }
    377382                        }

  • issm/trunk-jpl/src/c/shared/Enum/Enum.vim

    r24735 r24750  
    11501150syn keyword cConstant MaticeEnum
    11511151syn keyword cConstant MatlithoEnum
     1152syn keyword cConstant MathydroEnum
    11521153syn keyword cConstant MatrixParamEnum
    11531154syn keyword cConstant MaxAbsVxEnum

  • issm/trunk-jpl/src/c/shared/Enum/EnumDefinitions.h

    r24735 r24750  
    11501150        MaticeEnum,
    11511151        MatlithoEnum,
     1152        MathydroEnum,
    11521153        MatrixParamEnum,
    11531154        MaxAbsVxEnum,

  • issm/trunk-jpl/src/c/shared/Enum/EnumToStringx.cpp

    r24735 r24750  
    11521152                case MaticeEnum : return "Matice";
    11531153                case MatlithoEnum : return "Matlitho";
     1154                case MathydroEnum : return "Mathydro";
    11541155                case MatrixParamEnum : return "MatrixParam";
    11551156                case MaxAbsVxEnum : return "MaxAbsVx";

  • issm/trunk-jpl/src/c/shared/Enum/StringToEnumx.cpp

    r24735 r24750  
    11791179              else if (strcmp(name,"Matice")==0) return MaticeEnum;
    11801180              else if (strcmp(name,"Matlitho")==0) return MatlithoEnum;
     1181              else if (strcmp(name,"Mathydro")==0) return MathydroEnum;
    11811182              else if (strcmp(name,"MatrixParam")==0) return MatrixParamEnum;
    11821183              else if (strcmp(name,"MaxAbsVx")==0) return MaxAbsVxEnum;
    … …  
    12431244              else if (strcmp(name,"ProfilingCurrentFlops")==0) return ProfilingCurrentFlopsEnum;
    12441245              else if (strcmp(name,"ProfilingCurrentMem")==0) return ProfilingCurrentMemEnum;
    1245               else if (strcmp(name,"ProfilingSolutionTime")==0) return ProfilingSolutionTimeEnum;
    12461246         else stage=11;
    12471247   }
    12481248   if(stage==11){
    1249               if (strcmp(name,"Regionaloutput")==0) return RegionaloutputEnum;
     1249              if (strcmp(name,"ProfilingSolutionTime")==0) return ProfilingSolutionTimeEnum;
     1250              else if (strcmp(name,"Regionaloutput")==0) return RegionaloutputEnum;
    12501251              else if (strcmp(name,"Regular")==0) return RegularEnum;
    12511252              else if (strcmp(name,"RecoveryAnalysis")==0) return RecoveryAnalysisEnum;

  • issm/trunk-jpl/src/c/shared/io/Marshalling/IoCodeConversions.cpp

    r24469 r24750  
    262262                case 3: return MaticeEnum;
    263263                case 4: return MatenhancediceEnum;
    264                 case 5: return MatlithoEnum;
    265                 case 6: return MaterialsEnum;
     264                case 5: return MaterialsEnum;
    266265                default: _error_("Marshalled materials code \""<<enum_in<<"\" not supported yet");
     266        }
     267}/*}}}*/
     268int IoCodeToEnumNature(int enum_in){/*{{{*/
     269        switch(enum_in){
     270                case 1: return MatdamageiceEnum;
     271                case 2: return MatestarEnum;
     272                case 3: return MaticeEnum;
     273                case 4: return MatenhancediceEnum;
     274                case 5: return MaterialsEnum;
     275                case 6: return MatlithoEnum;
     276                case 7: return MathydroEnum;
     277                default: _error_("Marshalled materials nature code \""<<enum_in<<"\" not supported yet");
    267278        }
    268279}/*}}}*/

  • issm/trunk-jpl/src/c/shared/io/Marshalling/IoCodeConversions.h

    r23657 r24750  
    1010int IoCodeToEnumHydrology(int enum_in);
    1111int IoCodeToEnumMaterials(int enum_in);
     12int IoCodeToEnumNature(int enum_in);
    1213int IoCodeToEnumTimestepping(int enum_in);
    1314int IoCodeToEnumAmr(int enum_in);

  • issm/trunk-jpl/src/m/classes/materials.m

    r24060 r24750  
    1919                        %check this is acceptable:
    2020                        for i=1:length(self.nature),
    21                                 if ~(strcmpi(self.nature{i},'litho') | strcmpi(self.nature{i},'ice')),
    22                                         error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'')');
     21                                if ~(strcmpi(self.nature{i},'litho') | strcmpi(self.nature{i},'ice') | strcmpi(self.nature{i},'hydro')),
     22                                        error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'' or ''hydro'')');
    2323                                end
    2424                        end
    … …  
    5555                                        self.addprop('density');
    5656                                        self.addprop('issolid');
    57                                 otherwise
    58                                         error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'')');
     57                                case 'hydro'
     58                                        self.addprop('rho_ice');
     59                                        self.addprop('rho_water');
     60                                        self.addprop('earth_density');
     61                                otherwise
     62                                        error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'' or ''hydro'')');
    5963                                end
    6064                        end
    … …  
    126130                                        self.issolid=[1;1]; % is layer solid or liquid.
    127131
    128                                 otherwise
    129                                         error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'')');
     132                                case 'hydro'
     133                                        %ice density (kg/m^3)
     134                                        self.rho_ice=917.;
     135
     136                                        %ocean water density (kg/m^3)
     137                                        self.rho_water=1023.;
     138                                        %SLR
     139                                        self.earth_density= 5512;  % average density of the Earth, (kg/m^3)
     140
     141                                otherwise
     142                                        error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'' or ''hydro'')');
    130143                                end
    131144                        end
    … …  
    166179                                        fielddisplay(self,'density','array describing each layer''s density (numlayers) [kg/m^3]');
    167180                                        fielddisplay(self,'issolid','array describing whether the layer is solid or liquid (default 1) (numlayers)');
    168                                 otherwise
    169                                         error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'')');
     181                                case 'hydro'
     182                                        disp(sprintf('   \nHydro:'));
     183                                        fielddisplay(self,'rho_ice','ice density [kg/m^3]');
     184                                        fielddisplay(self,'rho_water','ocean water density [kg/m^3]');
     185                                        fielddisplay(self,'earth_density','Mantle density [kg/m^-3]');
     186
     187                                otherwise
     188                                        error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'' or ''hydro'')');
    170189                                end
    171190                        end
    … …  
    209228                                                        error(['Fluid layers detected at layers #', num2str(ind'), ', but having 2 or more adjacent fluid layers is not supported yet. Consider merging them.'])
    210229                                        end
    211                                 otherwise
    212                                         error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'')');
     230
     231                                case 'hydro'
     232                                        md = checkfield(md,'fieldname','materials.rho_ice','>',0);
     233                                        md = checkfield(md,'fieldname','materials.rho_water','>',0);
     234                                        md = checkfield(md,'fieldname','materials.earth_density','>',0,'numel',1);
     235
     236                                otherwise
     237                                        error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'' or ''hydro'')');
    213238                                end
    214239                        end
    … …  
    218243
    219244                        %1: MatdamageiceEnum 2: MatestarEnum 3: MaticeEnum 4: MatenhancediceEnum 5: MaterialsEnum
    220                         WriteData(fid,prefix,'name','md.materials.type','data',6,'format','Integer');
    221245                        WriteData(fid,prefix,'name','md.materials.nature','data',naturetointeger(self.nature),'format','IntMat','mattype',3);
     246                        WriteData(fid,prefix,'name','md.materials.type','data',5,'format','Integer'); %DANGER, this can evolve if you have classes.
    222247
    223248                        for i=1:length(self.nature),
    … …  
    251276                                        WriteData(fid,prefix,'object',self,'class','materials','fieldname','burgers_viscosity','format','DoubleMat','mattype',3);
    252277                                        WriteData(fid,prefix,'object',self,'class','materials','fieldname','burgers_mu','format','DoubleMat','mattype',3);
    253                                 otherwise
    254                                         error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'')');
     278                                case 'hydro'
     279                                        WriteData(fid,prefix,'object',self,'class','materials','fieldname','rho_ice','format','Double');
     280                                        WriteData(fid,prefix,'object',self,'class','materials','fieldname','rho_water','format','Double');
     281                                        WriteData(fid,prefix,'object',self,'class','materials','fieldname','earth_density','format','Double');
     282
     283                                otherwise
     284                                        error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'' or ''hydro'')');
    255285                                end
    256286                        end
    … …  
    299329                                        writejsdouble(fid,[modelname '.materials.burgers_viscosity'],self.burgers_viscosity);
    300330                                        writejsdouble(fid,[modelname '.materials.burgers_mu'],self.burgers_mu);
    301                                 otherwise
    302                                         error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'')');
     331
     332                                case 'hydro'
     333                                        writejsdouble(fid,[modelname '.materials.rho_ice'],self.rho_ice);
     334                                        writejsdouble(fid,[modelname '.materials.rho_water'],self.rho_water);
     335                                        writejsdouble(fid,[modelname '.materials.earth_density'],self.earth_density);
     336
     337                                otherwise
     338                                        error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'' or ''hydro'')');
    303339                                end
    304340                        end
    … …  
    313349                        for i=1:length(strnat),
    314350                                switch strnat{i},
    315                                 case 'damageice'
     351                                        case 'damageice'
    316352                                        intnat(i)=1;
    317353                                case 'estar'
    … …  
    321357                                case 'enhancedice'
    322358                                        intnat(i)=4;
    323                                 case 'litho'
    324                                         intnat(i)=5;
    325                                 otherwise
    326                                         error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'')');
    327                                 end
    328                         end
    329                 end % }}}
     359                                %case 'materials' %this case will never happen, kept to ensure equivalent of codes between IoCodeToMaterialsEnumm and IoCodeToNatureEnum
     360                                %       intnat(i)=5;
     361                                case 'litho'
     362                                        intnat(i)=6;
     363                                case 'hydro'
     364                                        intnat(i)=7;
     365                                otherwise
     366                                        error('materials constructor error message: nature of the material not supported yet! (''ice'' or ''litho'' or ''hydro'')');
     367                                end
     368                        end
     369                end % }}}

  • issm/trunk-jpl/test/NightlyRun/test2002.m

    r24482 r24750  
    5454md.slr.ocean_area_scaling=0;
    5555
    56 %geometry
    57 di=md.materials.rho_ice/md.materials.rho_water;
    58 md.geometry.thickness=ones(md.mesh.numberofvertices,1);
    59 md.geometry.surface=(1-di)*zeros(md.mesh.numberofvertices,1);
    60 md.geometry.base=md.geometry.surface-md.geometry.thickness;
    61 md.geometry.bed=md.geometry.base;
     56%Geometry for the bed, arbitrary:
     57md.geometry.bed=-ones(md.mesh.numberofvertices,1);
    6258
    63 %materials
    64 md.initialization.temperature=273.25*ones(md.mesh.numberofvertices,1);
    65 md.materials.rheology_B=paterson(md.initialization.temperature);
    66 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
     59%Materials:
     60md.materials=materials('hydro');
    6761
    6862%New stuff
Note: See TracChangeset for help on using the changeset viewer.