function [jd]=ecomtime(cdf,ind)
% ECOMTIME Julian time data from POM.CDF or ECOM.CDF type file
% function [jd]=ecomtime(cdf,[ind])
%    cdf= ecom.cdf file name 
%
mexcdf('setopts',0);
ncid=mexcdf('open',cdf,'nowrite');
[nam,nt]=mexcdf('diminq',ncid,'time');
base_date=[1968 5 23 0 0 0];
if(nargin==1),
   t=mexcdf('varget',ncid,'ocean_time',0,nt);
else
   t=mexcdf('varget1',ncid,'ocean_time',ind-1);
end
jd0=julian(base_date);
jd=jd0+t/(3600*24);
mexcdf('close',ncid);