Changeset 22323


Ignore:
Timestamp:
01/03/18 21:02:33 (7 years ago)
Author:
seroussi
Message:

CHG: adding area correction to compute volumes and areas

Location:
issm/trunk-jpl/src
Files:
13 edited

Legend:

Unmodified
Added
Removed

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

    r22004 r22323  
    229229        iomodel->FetchData(6,"md.mesh.x","md.mesh.y","md.mesh.z","md.geometry.base","md.geometry.thickness","md.mask.ice_levelset");
    230230        if (iomodel->domaintype == Domain3DsurfaceEnum) iomodel->FetchData(3,"md.mesh.lat","md.mesh.long","md.mesh.r");
     231        else iomodel->FetchDataToInput(elements,"md.mesh.scale_factor",MeshScaleFactorEnum,1.);
    231232       
    232233        CreateNumberNodeToElementConnectivity(iomodel,solution_type);

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

    r22294 r22323  
    308308        MeshYEnum,
    309309        MeshZEnum,
     310        MeshScaleFactorEnum,
    310311        MeshElementsEnum,
    311312        MeshNumberofelementsEnum,
    … …  
    750751        SubelementMigration2Enum,
    751752        SubelementMigration3Enum,
     753        SubelementMigration4Enum,
    752754        ContactEnum,
    753755        GroundingOnlyEnum,

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

    r22294 r22323  
    314314                case MeshYEnum : return "MeshY";
    315315                case MeshZEnum : return "MeshZ";
     316                case MeshScaleFactorEnum : return "MeshScaleFactor";
    316317                case MeshElementsEnum : return "MeshElements";
    317318                case MeshNumberofelementsEnum : return "MeshNumberofelements";
    … …  
    740741                case SubelementMigration2Enum : return "SubelementMigration2";
    741742                case SubelementMigration3Enum : return "SubelementMigration3";
     743                case SubelementMigration4Enum : return "SubelementMigration4";
    742744                case ContactEnum : return "Contact";
    743745                case GroundingOnlyEnum : return "GroundingOnly";

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

    r22294 r22323  
    320320              else if (strcmp(name,"MeshY")==0) return MeshYEnum;
    321321              else if (strcmp(name,"MeshZ")==0) return MeshZEnum;
     322              else if (strcmp(name,"MeshScaleFactor")==0) return MeshScaleFactorEnum;
    322323              else if (strcmp(name,"MeshElements")==0) return MeshElementsEnum;
    323324              else if (strcmp(name,"MeshNumberofelements")==0) return MeshNumberofelementsEnum;
    … …  
    382383              else if (strcmp(name,"TimesteppingStartTime")==0) return TimesteppingStartTimeEnum;
    383384              else if (strcmp(name,"TimesteppingFinalTime")==0) return TimesteppingFinalTimeEnum;
    384               else if (strcmp(name,"TimesteppingCflCoefficient")==0) return TimesteppingCflCoefficientEnum;
    385385         else stage=4;
    386386   }
    387387   if(stage==4){
    388               if (strcmp(name,"TimesteppingTimeAdapt")==0) return TimesteppingTimeAdaptEnum;
     388              if (strcmp(name,"TimesteppingCflCoefficient")==0) return TimesteppingCflCoefficientEnum;
     389              else if (strcmp(name,"TimesteppingTimeAdapt")==0) return TimesteppingTimeAdaptEnum;
    389390              else if (strcmp(name,"TimesteppingTimeStep")==0) return TimesteppingTimeStepEnum;
    390391              else if (strcmp(name,"TimesteppingInterpForcings")==0) return TimesteppingInterpForcingsEnum;
    … …  
    505506              else if (strcmp(name,"SmbBNeg")==0) return SmbBNegEnum;
    506507              else if (strcmp(name,"SMBhenning")==0) return SMBhenningEnum;
    507               else if (strcmp(name,"SMBcomponents")==0) return SMBcomponentsEnum;
    508508         else stage=5;
    509509   }
    510510   if(stage==5){
    511               if (strcmp(name,"SmbAccumulation")==0) return SmbAccumulationEnum;
     511              if (strcmp(name,"SMBcomponents")==0) return SMBcomponentsEnum;
     512              else if (strcmp(name,"SmbAccumulation")==0) return SmbAccumulationEnum;
    512513              else if (strcmp(name,"SmbEvaporation")==0) return SmbEvaporationEnum;
    513514              else if (strcmp(name,"SmbRunoff")==0) return SmbRunoffEnum;
    … …  
    628629              else if (strcmp(name,"StringExternalResult")==0) return StringExternalResultEnum;
    629630              else if (strcmp(name,"Step")==0) return StepEnum;
    630               else if (strcmp(name,"Time")==0) return TimeEnum;
    631631         else stage=6;
    632632   }
    633633   if(stage==6){
    634               if (strcmp(name,"WaterColumnOld")==0) return WaterColumnOldEnum;
     634              if (strcmp(name,"Time")==0) return TimeEnum;
     635              else if (strcmp(name,"WaterColumnOld")==0) return WaterColumnOldEnum;
    635636              else if (strcmp(name,"SurfaceObservation")==0) return SurfaceObservationEnum;
    636637              else if (strcmp(name,"WeightsSurfaceObservation")==0) return WeightsSurfaceObservationEnum;
    … …  
    751752              else if (strcmp(name,"None")==0) return NoneEnum;
    752753              else if (strcmp(name,"AggressiveMigration")==0) return AggressiveMigrationEnum;
    753               else if (strcmp(name,"SoftMigration")==0) return SoftMigrationEnum;
    754754         else stage=7;
    755755   }
    756756   if(stage==7){
    757               if (strcmp(name,"SubelementMigration")==0) return SubelementMigrationEnum;
     757              if (strcmp(name,"SoftMigration")==0) return SoftMigrationEnum;
     758              else if (strcmp(name,"SubelementMigration")==0) return SubelementMigrationEnum;
    758759              else if (strcmp(name,"SubelementMigration2")==0) return SubelementMigration2Enum;
    759760              else if (strcmp(name,"SubelementMigration3")==0) return SubelementMigration3Enum;
     761              else if (strcmp(name,"SubelementMigration4")==0) return SubelementMigration4Enum;
    760762              else if (strcmp(name,"Contact")==0) return ContactEnum;
    761763              else if (strcmp(name,"GroundingOnly")==0) return GroundingOnlyEnum;
    … …  
    873875              else if (strcmp(name,"EsaRequestedOutputs")==0) return EsaRequestedOutputsEnum;
    874876              else if (strcmp(name,"EsaNumRequestedOutputs")==0) return EsaNumRequestedOutputsEnum;
    875               else if (strcmp(name,"TransientAmrFrequency")==0) return TransientAmrFrequencyEnum;
    876               else if (strcmp(name,"AmrType")==0) return AmrTypeEnum;
    877877         else stage=8;
    878878   }
    879879   if(stage==8){
    880               if (strcmp(name,"AmrRestart")==0) return AmrRestartEnum;
     880              if (strcmp(name,"TransientAmrFrequency")==0) return TransientAmrFrequencyEnum;
     881              else if (strcmp(name,"AmrType")==0) return AmrTypeEnum;
     882              else if (strcmp(name,"AmrRestart")==0) return AmrRestartEnum;
    881883              else if (strcmp(name,"AmrNeopz")==0) return AmrNeopzEnum;
    882884              else if (strcmp(name,"AmrLevelMax")==0) return AmrLevelMaxEnum;
    … …  
    996998              else if (strcmp(name,"HydrologyDCEfficientAnalysis")==0) return HydrologyDCEfficientAnalysisEnum;
    997999              else if (strcmp(name,"HydrologySommersAnalysis")==0) return HydrologySommersAnalysisEnum;
    998               else if (strcmp(name,"HydrologySolution")==0) return HydrologySolutionEnum;
    999               else if (strcmp(name,"MeltingAnalysis")==0) return MeltingAnalysisEnum;
    10001000         else stage=9;
    10011001   }
    10021002   if(stage==9){
    1003               if (strcmp(name,"MasstransportAnalysis")==0) return MasstransportAnalysisEnum;
     1003              if (strcmp(name,"HydrologySolution")==0) return HydrologySolutionEnum;
     1004              else if (strcmp(name,"MeltingAnalysis")==0) return MeltingAnalysisEnum;
     1005              else if (strcmp(name,"MasstransportAnalysis")==0) return MasstransportAnalysisEnum;
    10041006              else if (strcmp(name,"MasstransportSolution")==0) return MasstransportSolutionEnum;
    10051007              else if (strcmp(name,"FreeSurfaceBaseAnalysis")==0) return FreeSurfaceBaseAnalysisEnum;

  • issm/trunk-jpl/src/m/classes/mesh2d.js

    r21065 r22323  
    4444                        fielddisplay(this,"long","vertices longitude [degrees]");
    4545                        fielddisplay(this,"epsg","EPSG code (ex: 3413 for UPS Greenland, 3031 for UPS Antarctica)");
     46                        fielddisplay(self,"scale_factor","Projection correction for volume, area, etc. computation)");
    4647                } //}}}
    4748                this.classname = function () { //{{{
    … …  
    7071                        checkfield(md,'fieldname','mesh.average_vertex_connectivity','>=',9,'message',"'mesh.average_vertex_connectivity' should be at least 9 in 2d");
    7172                        checkfield(md,'fieldname','mesh.segments','NaN',1,'Inf',1,'>',0,'size',[NaN, 3]);
     73                        if(this.scale_factor.length>1) checkfield(md,'fieldname','mesh.scale_factor','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices, 1]);
    7274
    7375                        switch(solution){
    … …  
    9294                        WriteData(fid,prefix,'object',this,'class','mesh','fieldname','vertexonboundary','format','DoubleMat','mattype',1);
    9395                        WriteData(fid,prefix,'object',this,'class','mesh','fieldname','segments','format','DoubleMat','mattype',3);
     96                        WriteData(fid,prefix,'object',this,'class','mesh','fieldname','scale_factor','format','DoubleMat','mattype',1);
    9497                }//}}}
    9598                this.fix=function() { //{{{
    … …  
    104107                        this.lat=NullFix(this.lat,NaN);
    105108                        this.long=NullFix(this.long,NaN);
     109                        this.scale_factor=NullFix(this.scale_factor,NaN);
    106110                }//}}}
    107111
    … …  
    118122                this.long                        = NaN;
    119123                this.epsg                        = 0;
     124                this.scale_factor                = NaN;
    120125
    121126                this.vertexonboundary            = NaN;

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

    r22004 r22323  
    1616                long                        = NaN;
    1717                epsg                        = 0;
     18                scale_factor                = NaN;
    1819
    1920                vertexonboundary            = NaN;
    … …  
    9495                        md = checkfield(md,'fieldname','mesh.average_vertex_connectivity','>=',9,'message','''mesh.average_vertex_connectivity'' should be at least 9 in 2d');
    9596                        md = checkfield(md,'fieldname','mesh.segments','NaN',1,'Inf',1,'>',0,'size',[NaN 3]);
     97                        if numel(md.mesh.scale_factor)>1,
     98                                md = checkfield(md,'fieldname','mesh.scale_factor','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     99                        end
    96100
    97101                        if strcmp(solution,'ThermalSolution')
    … …  
    127131                        fielddisplay(self,'long','vertices longitude [degrees]');
    128132                        fielddisplay(self,'epsg','EPSG code (ex: 3413 for UPS Greenland, 3031 for UPS Antarctica)');
     133                        fielddisplay(self,'scale_factor','Projection correction for volume, area, etc. computation)');
    129134                end % }}}
    130135                function marshall(self,prefix,md,fid) % {{{
    … …  
    141146                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','vertexonboundary','format','DoubleMat','mattype',1);
    142147                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','segments','format','DoubleMat','mattype',3);
     148                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','scale_factor','format','DoubleMat','mattype',1);
    143149                end % }}}
    144150                function t = domaintype(self) % {{{
    … …  
    162168                        writejs1Darray(fid,[modelname '.mesh.long'],self.long);
    163169                        writejsdouble(fid,[modelname '.mesh.epsg'],self.epsg);
     170                        writejsdouble(fid,[modelname '.mesh.scale_factor'],self.scale_factor);
    164171                        writejs1Darray(fid,[modelname '.mesh.vertexonboundary'],self.vertexonboundary);
    165172                        writejs2Darray(fid,[modelname '.mesh.edges'],self.edges);

  • issm/trunk-jpl/src/m/classes/mesh2d.py

    r22150 r22323  
    2424                self.long                        = float('NaN');
    2525                self.epsg                        = 0;
     26                self.scale_factor                = float('NaN');
    2627
    2728                self.vertexonboundary            = float('NaN');
    … …  
    6869                string="%s\n%s"%(string,fielddisplay(self,"long","vertices longitude [degrees]"))
    6970                string="%s\n%s"%(string,fielddisplay(self,"epsg","EPSG code (ex: 3413 for UPS Greenland, 3031 for UPS Antarctica)"))
     71                string="%s\n%s"%(string,fielddisplay(self,"scale_factor","Projection correction for volume, area, etc. computation"))
    7072                return string
    7173                #}}}
    … …  
    9597                md = checkfield(md,'fieldname','mesh.average_vertex_connectivity','>=',9,'message',"'mesh.average_vertex_connectivity' should be at least 9 in 2d")
    9698                md = checkfield(md,'fieldname','mesh.segments','NaN',1,'Inf',1,'>',0,'size',[np.nan,3]);
     99                if(numel(self.scale_factor)>1):
     100                        md = checkfield(md,'fieldname','mesh.scale_factor','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     101               
    97102                if solution=='ThermalSolution':
    98103                        md.checkmessage("thermal not supported for 2d mesh")
    … …  
    122127                WriteData(fid,prefix,'object',self,'class','mesh','fieldname','vertexonboundary','format','DoubleMat','mattype',1)
    123128                WriteData(fid,prefix,'object',self,'class','mesh','fieldname','segments','format','DoubleMat','mattype',3)
     129                WriteData(fid,prefix,'object',self,'class','mesh','fieldname','scale_factor','format','DoubleMat','mattype',1)
    124130        # }}}

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

    r22004 r22323  
    1616                long                        = NaN;
    1717                epsg                        = NaN;
     18                scale_factor                = NaN;
    1819
    1920                vertexonboundary            = NaN;
    … …  
    9495                        md = checkfield(md,'fieldname','mesh.vertexonsurface','size',[md.mesh.numberofvertices 1],'values',[0 1]);
    9596                        md = checkfield(md,'fieldname','mesh.average_vertex_connectivity','>=',9,'message','''mesh.average_vertex_connectivity'' should be at least 9 in 2d');
     97                        if numel(md.mesh.scale_factor)>1,
     98                                md = checkfield(md,'fieldname','mesh.scale_factor','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     99                        end
    96100
    97101                        if strcmp(solution,'ThermalSolution')
    … …  
    125129                        fielddisplay(self,'long','vertices longitude [degrees]');
    126130                        fielddisplay(self,'epsg','EPSG code (ex: 3413 for UPS Greenland, 3031 for UPS Antarctica)');
     131                        fielddisplay(self,'scale_factor','Projection correction for volume, area, etc. computation)');
    127132                end % }}}
    128133                function marshall(self,prefix,md,fid) % {{{
    … …  
    139144                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','vertexonsurface','format','BooleanMat','mattype',1);
    140145                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','average_vertex_connectivity','format','Integer');
     146                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','scale_factor','format','DoubleMat','mattype',1);
    141147                end % }}}
    142148                function t = domaintype(self) % {{{

  • issm/trunk-jpl/src/m/classes/mesh2dvertical.py

    r22267 r22323  
    2424                self.long                        = float('NaN')
    2525                self.epsg                        = float('NaN')
     26                self.scale_factor                = float('NaN');
    2627
    2728                self.vertexonboundary            = float('NaN')
    … …  
    6667                string="%s\n%s"%(string,fielddisplay(self,"long","vertices longitude [degrees]"))
    6768                string="%s\n%s"%(string,fielddisplay(self,"epsg","EPSG code (ex: 3413 for UPS Greenland, 3031 for UPS Antarctica)"))
     69                string="%s\n%s"%(string,fielddisplay(self,"scale_factor","Projection correction for volume, area, etc. computation"))
    6870                return string
    6971                #}}}
    … …  
    9496                md = checkfield(md,'fieldname','mesh.vertexonsurface','size',[md.mesh.numberofvertices],'values',[0,1])
    9597                md = checkfield(md,'fieldname','mesh.average_vertex_connectivity','>=',9,'message',"'mesh.average_vertex_connectivity' should be at least 9 in 2d")
     98                if(numel(self.scale_factor)>1):
     99                        md = checkfield(md,'fieldname','mesh.scale_factor','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
    96100
    97101                if solution=='ThermalSolution':
    … …  
    128132                WriteData(fid,prefix,'object',self,'class','mesh','fieldname','vertexonsurface','format','BooleanMat','mattype',1)
    129133                WriteData(fid,prefix,'object',self,'class','mesh','fieldname','average_vertex_connectivity','format','Integer')
     134                WriteData(fid,prefix,'object',self,'class','mesh','fieldname','scale_factor','format','DoubleMat','mattype',1)
    130135        # }}}

  • issm/trunk-jpl/src/m/classes/mesh3dprisms.js

    r20970 r22323  
    3434                        checkfield(md,'fieldname','mesh.z','<=',md.geometry.surface+Math.pow(10, -10),'message','\'mesh.z\' higher than surface elevation');
    3535                        checkfield(md,'fieldname','mesh.average_vertex_connectivity','>=',24,'message','\'mesh.average_vertex_connectivity\' should be at least 24 in 3d');
     36                        if(this.scale_factor.length>1) checkfield(md,'fieldname','mesh.scale_factor','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices, 1]);
    3637                } // }}}
    3738                this.disp = function()  { // {{{
    … …  
    7576                        fielddisplay(this,'long','vertices longitude [degrees]');
    7677                        fielddisplay(this,'epsg','EPSG code (ex: 3413 for UPS Greenland, 3031 for UPS Antarctica)');
     78                        fielddisplay(self,"scale_factor","Projection correction for volume, area, etc. computation)");
    7779                } // }}}
    7880                this.marshall = function(md,prefix,fid) { // {{{
    … …  
    9597                        WriteData(fid,prefix,'object',this,'class','mesh','fieldname','numberofvertices2d','format','Integer');
    9698                        WriteData(fid,prefix,'object',this,'class','mesh','fieldname','numberofelements2d','format','Integer');
     99                        WriteData(fid,prefix,'object',this,'class','mesh','fieldname','scale_factor','format','DoubleMat','mattype',1);
    97100                } // }}}
    98101        this.fix=function() { //{{{
    … …  
    109112            this.lat=NaNFix(this.lat,NaN);
    110113            this.long=NaNFix(this.long,NaN);
     114                                this.scale_factor=NullFix(this.scale_factor,NaN);
    111115        }//}}}
    112116                this.domaintype = function() { // {{{)
    … …  
    135139        this.long                        = NaN;
    136140        this.epsg                        = 0;
     141                  this.scale_factor                = NaN;
    137142
    138143        this.vertexonbase                = NaN;

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

    r22004 r22323  
    1717                long                        = NaN;
    1818                epsg                        = 0;
     19                scale_factor                = NaN;
    1920
    2021                vertexonbase                = NaN;
    … …  
    120121                                md = checkmessage(md,'md.mesh.z is not consistent with md.geometry.thickness, you changed the geometry after extrusion !!');
    121122                        end
     123                        if numel(md.mesh.scale_factor)>1,
     124                                md = checkfield(md,'fieldname','mesh.scale_factor','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     125                        end
    122126                end % }}}
    123127                function disp(self) % {{{
    … …  
    161165                        fielddisplay(self,'long','vertices longitude [degrees]');
    162166                        fielddisplay(self,'epsg','EPSG code (ex: 3413 for UPS Greenland, 3031 for UPS Antarctica)');
     167                        fielddisplay(self,'scale_factor','Projection correction for volume, area, etc. computation)');
    163168                end % }}}
    164169                function marshall(self,prefix,md,fid) % {{{
    … …  
    181186                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','numberofvertices2d','format','Integer');
    182187                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','numberofelements2d','format','Integer');
     188                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','scale_factor','format','DoubleMat','mattype',1);
    183189                end % }}}
    184190                function type = domaintype(self) % {{{

  • issm/trunk-jpl/src/m/classes/mesh3dprisms.py

    r21410 r22323  
    2525                self.long                        = float('NaN');
    2626                self.epsg                        = 0;
     27                self.scale_factor                = float('NaN');
    2728
    2829                self.vertexonbase                = float('NaN');
    … …  
    9091                string="%s\n%s"%(string,fielddisplay(self,"long","vertices longitude [degrees]"))
    9192                string="%s\n%s"%(string,fielddisplay(self,"epsg","EPSG code (ex: 3413 for UPS Greenland, 3031 for UPS Antarctica)"))
     93                string="%s\n%s"%(string,fielddisplay(self,"scale_factor","Projection correction for volume, area, etc. computation"))
    9294                return string
    9395                #}}}
    … …  
    118120                md = checkfield(md,'fieldname','mesh.vertexonsurface','size',[md.mesh.numberofvertices],'values',[0,1])
    119121                md = checkfield(md,'fieldname','mesh.average_vertex_connectivity','>=',24,'message',"'mesh.average_vertex_connectivity' should be at least 24 in 3d")
     122                if(numel(self.scale_factor)>1):
     123                        md = checkfield(md,'fieldname','mesh.scale_factor','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
    120124
    121125                return md
    … …  
    149153                WriteData(fid,prefix,'object',self,'class','mesh','fieldname','numberofvertices2d','format','Integer')
    150154                WriteData(fid,prefix,'object',self,'class','mesh','fieldname','numberofelements2d','format','Integer')
     155                WriteData(fid,prefix,'object',self,'class','mesh','fieldname','scale_factor','format','DoubleMat','mattype',1)
    151156        # }}}

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

    r22004 r22323  
    1717                long                        = NaN;
    1818                epsg                        = 0;
     19                scale_factor                = NaN;
    1920
    2021                vertexonbase                = NaN;
    … …  
    109110                        md = checkfield(md,'fieldname','mesh.z','<=',md.geometry.surface+10^-10,'message','''mesh.z'' higher than surface elevation');
    110111                        md = checkfield(md,'fieldname','mesh.average_vertex_connectivity','>=',24,'message','''mesh.average_vertex_connectivity'' should be at least 24 in 3d');
     112                        if numel(md.mesh.scale_factor)>1,
     113                                md = checkfield(md,'fieldname','mesh.scale_factor','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     114                        end
     115
    111116                end % }}}
    112117                function disp(self) % {{{
    … …  
    150155                        fielddisplay(self,'long','vertices longitude [degrees]');
    151156                        fielddisplay(self,'epsg','EPSG code (ex: 3413 for UPS Greenland, 3031 for UPS Antarctica)');
     157                        fielddisplay(self,'scale_factor','Projection correction for volume, area, etc. computation)');
    152158                end % }}}
    153159                function marshall(self,prefix,md,fid) % {{{
    … …  
    166172                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','lowerelements','format','DoubleMat','mattype',2);
    167173                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','upperelements','format','DoubleMat','mattype',2);
     174                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','scale_factor','format','DoubleMat','mattype',1);
    168175                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','average_vertex_connectivity','format','Integer');
    169176                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','elements2d','format','DoubleMat','mattype',3);
Note: See TracChangeset for help on using the changeset viewer.