cd ../uncompressed
fn='tmpnc02.nc'
% for t=1:20,
% [Dnum(t), Dstr]=roms_get_date(fn,t,0,'ocean_time');
% end
% datevec(Dnum+datenum(1968,5,24))

% nc=netcdf(fn);
% time=nc{'ocean_time'}(:);
% close(nc)

time=nc_varget(fn,'ocean_time');

time(1)./(3600.*24)
jd=datenum(1968,5,24)+time./(3600.*24);

['beginning']
datevec(jd(1))

fn='tmpnc34.nc';
% for t=1:20,
% [Dnum(t), Dstr]=roms_get_date(fn,t,0,'ocean_time');
% end
% datevec(Dnum+datenum(1968,5,24))

% nc=netcdf(fn);
% time=nc{'ocean_time'}(:);
% close(nc)

time=nc_varget(fn,'ocean_time');


jd=datenum(1968,5,24)+time./(3600.*24);
['end']
datevec(jd(length(jd)))


cd ../analout