function jdmat=roms_time(fname,var);
% ROMS_TIME get time from Roms file
% Usage: jdmat=roms_time(fname,var);
%
%  fname = Roms netcdf file
%  var = variable name to find associated time variable for 
nc = netcdf(fname);
if nargin==1,
 time_variable='ocean_time';
else
if isempty(nc{var})
  error([ 'Variable ' var ' is not present in file ' file])
end

% get the time
time_variable = nc{var}.time(:);
if isempty(time_variable)
  time_variable = 'ocean_time'; 
end               
end
ot = nc{time_variable}(:);
time_units = nc{time_variable}.units(:);
if strcmp(lower(time_units(1:3)),'sec');
   ot = ot/86400;
end
jdmat=ot+datenum([1968 5 23 0 0 0]);
close(nc);