function	[d]=fmris_write_minc(d,Z,T);

d.dim=d.dim(length(d.dim):-1:1);

if nargin ==3
   if T(1)==1 & Z(1)==1
      newh=newimage(d.file_name,d.dim,d.parent_file,'float');
   else
      newh=openimage(d.file_name,'w');
   end
elseif nargin==1
   newh=newimage(d.file_name,d.dim,d.parent_file,'float');
   Z=1:d.dim(2);
   %T=1;
   T=1:(d.dim(1)+(d.dim(1)==0));
end

data=squeeze(reshape(d.data,d.dim(4)*d.dim(3),length(Z),length(T)));

%if (length(T)==1)&(T==1); T=[]; end
%putimages(newh,data,Z,T);

if d.dim(1)==0
   putimages(newh,data,Z);
elseif length(T)==1|length(Z)==1
   putimages(newh,data,Z,T);
else
   for t=T
      putimages(newh,squeeze(data(:,:,t)),Z,t);
   end
end

closeimage(newh);