function [OKall,OKlegend] = nc_t_tide_test
%nc_t_tide_test  test for nc_t_tide
%
%See also: t_tide

rnd        = 1e-9; % dt variations are neglected if smaller than eps('single')
tstart     = datenum(1968,1,10); % [datenum]
nt         = 12;                 % [#] %% if 12 for 24hr series, canonical t_tide 1.3 crashes
tend       = tstart+1.333;       % [datenum] % test for short period not aligned with S2, to test centraltime
T          = 12/24;              % [day] % set to S2

   to  = linspace(tstart,tend,nt+(rem(nt,2)+1));
   zo  = sin(2*pi.*to./T);

   ncfile = 'nc_t_tide_scalar.nc';
   nc_t_tide(to,zo     ,'ncfile',ncfile,'units','m');
   nc_dump(ncfile,[],[filename(ncfile),'.cdl'])

   ncfile = 'nc_t_tide_vector.nc';
   nc_t_tide(to,zo+i*zo,'ncfile',ncfile,'units','m');
   nc_dump(ncfile,[],[filename(ncfile),'.cdl'])
   