function [argos,argot,argop, argoo,argolat,argolon] = read_nc_write_mat(input_dir,datactr,time,callSign,profnum);
%[argos,argot,argop,argoo] = read_nc_write_mat(datactr,time,callSign,profnum);
%
%
if nargin==2
    filename = datactr;
else
    if profnum<10
        profstr = ['00',num2str(profnum)];
    elseif profnum<100
        profstr = ['0',num2str(profnum)];
    else
        profstr = num2str(profnum);
    end
    filename = [datactr,'_',time,num2str(callSign),'_',profstr];
end

if ~exist([input_dir,filename,'.nc'])
    disp([input_dir,filename,'.nc does not exist']);
end

% Temperature, pressure
argot = nc_varget([input_dir,filename,'.nc'],'temperature');
argop = nc_varget([input_dir,filename,'.nc'],'pressure')';

% If salinity exists, get it
issal = nc_isvar([input_dir,filename,'.nc'],'salinity');
if issal
    argos = nc_varget([input_dir,filename,'.nc'],'salinity');
else
    argos = NaN*argot;
end
% If oxygen exists, get it
issal = nc_isvar([input_dir,filename,'.nc'],'oxygen');
if issal
    argoo = nc_varget([input_dir,filename,'.nc'],'oxygen');
else
    argoo = NaN*argot;
end


%disp('saving')
eval(['save ',input_dir,filename,'.mat argos argot argop argoo'])