function gifMovie(g,fname,clim,delay)
% function gifMovie(g,fname,clim,delay)
%
% g    data
% delay  time in seconds


if ndims(g) == 3
    % apparently GIF doesn't accept 3D images; 2D or 4D only
    dim = size(g);
    g = reshape(g,[dim(1:2) 1 dim(3)]);
end

if nargin < 3
    clim = [min(g(:)) max(g(:))];
end

if nargin < 4
    delay = 1/3;
end

% In this case, my clim was:
% clim = [-200 -50], g - offset / length
g = round((g - clim(1)) / diff(clim) *255);
% Values must be within range:
g(g<1) = 1; g(g>255) = 255;
% g 1:255, with 0 for nans
g = uint8(g);

cmap = jet(256);

imwrite(g,cmap,fname,'DelayTime',delay,'LoopCount',Inf,...
    'BackgroundColor',0,'TransparentColor',0,'DisposalMethod','restoreBG')