And yes indeed, MATLAB is true, I have checked the slr class (slr.m), in the classdef slr, there is no such thing as steric_rate defined in here. What should I edit in the classdef slr? Should I include the 'steric_date' and attributing it the same as 'sealevel'? Here are the list of your classdef in slr.m:
fielddisplay(self,'deltathickness','thickness change: ice height equivalent [m]');
fielddisplay(self,'sealevel','current sea level (prior to computation) [m]');
fielddisplay(self,'spcthickness','thickness constraints (NaN means no constraint) [m]');
fielddisplay(self,'reltol','sea level rise relative convergence criterion, (default, NaN: not applied)');
fielddisplay(self,'abstol','sea level rise absolute convergence criterion, NaN: not applied');
fielddisplay(self,'maxiter','maximum number of nonlinear iterations');
fielddisplay(self,'love_h','load Love number for radial displacement');
fielddisplay(self,'love_k','load Love number for gravitational potential perturbation');
fielddisplay(self,'love_l','load Love number for horizontal displacements');
fielddisplay(self,'tide_love_k','tidal load Love number (deg 2)');
fielddisplay(self,'tide_love_h','tidal load Love number (deg 2)');
fielddisplay(self,'fluid_love','secular fluid Love number');
fielddisplay(self,'equatorial_moi','mean equatorial moment of inertia [kg m^2]');
fielddisplay(self,'polar_moi','polar moment of inertia [kg m^2]');
fielddisplay(self,'angular_velocity','mean rotational velocity of earth [per second]');
fielddisplay(self,'ocean_area_scaling','correction for model representation of ocean area [default: No correction]');
fielddisplay(self,'hydro_rate','rate of hydrological expansion (in mm/yr)');
fielddisplay(self,'Ngia','rate of viscous (GIA) geoid expansion (in mm/yr)');
fielddisplay(self,'Ugia','rate of viscous (GIA) bedrock uplift (in mm/yr)');
fielddisplay(self,'loop_increment','vector assembly (in the convolution) framentation');
fielddisplay(self,'geodetic','compute geodetic SLR? ( in addition to steric?) default 0');
fielddisplay(self,'geodetic_run_frequency','how many time steps we skip before we run SLR solver during transient (default: 1)');
fielddisplay(self,'rigid','rigid earth graviational potential perturbation');
fielddisplay(self,'elastic','elastic earth graviational potential perturbation');
fielddisplay(self,'rotation','earth rotational potential perturbation');
fielddisplay(self,'degacc','accuracy (default .01 deg) for numerical discretization of the Green''s functions');
fielddisplay(self,'transitions','indices into parts of the mesh that will be icecaps');
fielddisplay(self,'requested_outputs','additional outputs requested');
Let me know what should I do, I would happy to do it. (I'll consider this as my learning process). Thank you
Best,
Abdul