wolffd@0: function ss = set(s,varargin)
wolffd@0: % SET Set properties for the MIRspectrum object
wolffd@0: % and return the updated object
wolffd@0: 
wolffd@0: propertyArgIn = varargin;
wolffd@0: ph = s.phase;
wolffd@0: log = s.log;
wolffd@0: pow = s.pow;
wolffd@0: xs = s.xscale;
wolffd@0: d = mirdata(s);
wolffd@0: d = set(d,'Title',get(s,'Title'),'Abs',get(s,'Abs'),'Ord',get(s,'Ord'));
wolffd@0: while length(propertyArgIn) >= 2,
wolffd@0:    prop = propertyArgIn{1};
wolffd@0:    val = propertyArgIn{2};
wolffd@0:    propertyArgIn = propertyArgIn(3:end);
wolffd@0:    switch prop
wolffd@0:        case 'Magnitude'
wolffd@0:            d = set(d,'Data',val);
wolffd@0:        case 'Frequency'
wolffd@0:            d = set(d,'Pos',val);
wolffd@0:        case 'Phase'
wolffd@0:            ph = val;
wolffd@0:        case 'log'
wolffd@0:            log = val;
wolffd@0:        case 'XScale'
wolffd@0:            xs = val;
wolffd@0:        case 'Power'
wolffd@0:            pow = val;
wolffd@0:        otherwise
wolffd@0:            d = set(d,prop,val);
wolffd@0:    end
wolffd@0: end
wolffd@0: ss.phase = ph;
wolffd@0: ss.log = log;
wolffd@0: ss.xscale = xs;
wolffd@0: ss.pow = pow;
wolffd@0: ss = class(ss,'mirspectrum',d);