--- a/fourier/plotfft.m
+++ b/fourier/plotfft.m
@@ -41,6 +41,10 @@
 %
 %     'posfreq'    Display only the positive frequencies.
 %
+%     'dim',dim  If `coef` is multidimensional, `dim` indicates the 
+%                dimension along which are the individual channels oriented.
+%                Value 1 indicates columns, value 2 rows.
+%
 %
 %   In addition to these parameters, `plotfft` accepts any of the flags
 %   from |normalize|. The coefficients will be normalized as specified
@@ -52,9 +56,9 @@
   error('%s: Too few input parameters.',upper(mfilename));
 end;
 
-if ~isvector(coef)>1
-  error('Input is multidimensional.');
-end;
+% if ~isvector(coef)
+%   error('%s: Can only plot vectors.',upper(mfilename));
+% end;
 
 definput.import={'ltfattranslate','normalize'};
 definput.importdefaults={'null'};
@@ -66,15 +70,12 @@
 definput.keyvals.dynrange=[];
 
 definput.keyvals.opts={};
-
+definput.keyvals.dim=[];
 [flags,kv,fs]=ltfatarghelper({'fs','dynrange'},definput,varargin);
 
+[coef,~,N]=assert_sigreshape_pre(coef,[],kv.dim,upper(mfilename));
+
 coef=normalize(coef,flags.norm);
-
-if ~isvector(coef)
-  error('%s: Can only plot vectors.',upper(mfilename));
-end;
-N=length(coef);
 
 % Apply transformation to coefficients.
 if flags.do_db
@@ -119,7 +120,7 @@
 else
 
   N2=floor(N/2)+1;
-  coef=coef(1:N2);
+  coef=coef(1:N2,:);
   xr=(0:N2-1)*2/N;  
 end;