Changeset 24778


Ignore:
Timestamp:
05/04/20 10:00:58 (5 years ago)
Author:
Eric.Larour
Message:

CHG: added proj field towards removing epsg field in the future. Left epsg
field for now.

Location:
issm/trunk-jpl/src/m/classes
Files:
5 edited

Legend:

Unmodified
Added
Removed

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

    r22879 r24778  
    1616                long                        = NaN;
    1717                epsg                        = 0;
     18                proj                        = 0;
    1819                scale_factor                = NaN;
    1920
    … …  
    4445                                        if strcmpi(oldself.hemisphere,'n'),
    4546                                                self.epsg=3413;
     47                                                self.proj=epsg2proj(3413);
    4648                                        else
    4749                                                self.epsg=3031;
     50                                                self.proj=epsg2proj(3031);
    4851                                        end
    4952                                end
    … …  
    131134                        fielddisplay(self,'long','vertices longitude [degrees]');
    132135                        fielddisplay(self,'epsg','EPSG code (ex: 3413 for UPS Greenland, 3031 for UPS Antarctica)');
     136                        fielddisplay(self,'proj','PROJ.4 compatible projection string');
    133137                        fielddisplay(self,'scale_factor','Projection correction for volume, area, etc. computation)');
    134138                end % }}}
    … …  
    172176                        writejs1Darray(fid,[modelname '.mesh.long'],self.long);
    173177                        writejsdouble(fid,[modelname '.mesh.epsg'],self.epsg);
     178                        writejsdouble(fid,[modelname '.mesh.proj'],self.proj);
    174179                        writejsdouble(fid,[modelname '.mesh.scale_factor'],self.scale_factor);
    175180                        writejs1Darray(fid,[modelname '.mesh.vertexonboundary'],self.vertexonboundary);

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

    r22879 r24778  
    1616                long                        = NaN;
    1717                epsg                        = NaN;
     18                proj                        = '';
    1819                scale_factor                = NaN;
    1920
    … …  
    4344                                        if strcmpi(oldself.hemisphere,'n'),
    4445                                                self.epsg=3413;
     46                                                self.proj=epsg2proj(3413);
    4547                                        else
    4648                                                self.epsg=3031;
     49                                                self.proj=epsg2proj(3031);
    4750                                        end
    4851                                end
    … …  
    129132                        fielddisplay(self,'long','vertices longitude [degrees]');
    130133                        fielddisplay(self,'epsg','EPSG code (ex: 3413 for UPS Greenland, 3031 for UPS Antarctica)');
     134                        fielddisplay(self,'proj','PROJ.4 compatible projection string');
    131135                        fielddisplay(self,'scale_factor','Projection correction for volume, area, etc. computation)');
    132136                end % }}}

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

    r23981 r24778  
    1717                long                        = NaN;
    1818                epsg                        = 0;
     19                proj                        = '';
    1920                scale_factor                = NaN;
    2021
    … …  
    5455                                        if strcmpi(oldself.hemisphere,'n'),
    5556                                                self.epsg=3413;
     57                                                self.proj=epsg2proj(3413);
    5658                                        else
    5759                                                self.epsg=3031;
     60                                                self.proj=epsg2proj(3031);
    5861                                        end
    5962                                end
    … …  
    166169                        fielddisplay(self,'long','vertices longitude [degrees]');
    167170                        fielddisplay(self,'epsg','EPSG code (ex: 3413 for UPS Greenland, 3031 for UPS Antarctica)');
     171                        fielddisplay(self,'proj','PROJ.4 compatible projection string');
    168172                        fielddisplay(self,'scale_factor','Projection correction for volume, area, etc. computation)');
    169173                end % }}}
    … …  
    216220                        writejs1Darray(fid,[modelname '.mesh.long'],self.long);
    217221                        writejs1Darray(fid,[modelname '.mesh.epsg'],self.epsg);
     222                        writejs1Darray(fid,[modelname '.mesh.proj'],self.proj);
    218223                        writejs1Darray(fid,[modelname '.mesh.vertexonbase'],self.vertexonbase);
    219224                        writejs1Darray(fid,[modelname '.mesh.vertexonsurface'],self.vertexonsurface);

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

    r24156 r24778  
    4040                                %Assign property values from struct
    4141                                self=structtoobj(mesh3dsurface(),oldself);
    42                                 if isfield(oldself,'hemisphere'),
    43                                         disp('md.mesh.hemisphere has been automatically converted to EPSG code');
    44                                         if strcmpi(oldself.hemisphere,'n'),
    45                                                 self.epsg=3413;
    46                                         else
    47                                                 self.epsg=3031;
    48                                         end
    49                                 end
     42                                %if isfield(oldself,'hemisphere'), %NOTE: not ever needed, epsg is not a field anymore.
     43                                %       disp('md.mesh.hemisphere has been automatically converted to EPSG code');
     44                                %       if strcmpi(oldself.hemisphere,'n'),
     45                                %               self.epsg=3413;
     46                                %       else
     47                                %               self.epsg=3031;
     48                                %       end
     49                                %end
    5050                        end
    5151

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

    r22879 r24778  
    1717                long                        = NaN;
    1818                epsg                        = 0;
     19                proj                        = '';
    1920                scale_factor                = NaN;
    2021
    … …  
    5455                                        if strcmpi(oldself.hemisphere,'n'),
    5556                                                self.epsg=3413;
     57                                                self.proj=epsg2proj(3413);
    5658                                        else
    5759                                                self.epsg=3031;
     60                                                self.proj=epsg2proj(3031);
    5861                                        end
    5962                                end
    … …  
    155158                        fielddisplay(self,'long','vertices longitude [degrees]');
    156159                        fielddisplay(self,'epsg','EPSG code (ex: 3413 for UPS Greenland, 3031 for UPS Antarctica)');
     160                        fielddisplay(self,'proj','PROJ.4 compatible projection string');
    157161                        fielddisplay(self,'scale_factor','Projection correction for volume, area, etc. computation)');
    158162                end % }}}
Note: See TracChangeset for help on using the changeset viewer.