get_NDBC_air.m

function [aa payload] = get_NDBC_air(fname,pdir)
%
%   read NODC netCDF4 file to get met information
%   created 07/12 by TJ. Hesser
% 
fprintf(1,'Analyizing Station %5s\n',fname(end-21:end-17));
%Get attributes distributed
fileatt = ncinfo(fname);
fatt = {fileatt.Attributes.Name};

% Get time out of the file and convert to year mon day hour min sec
times = double(ncread(fname,'/time'));
[year,mon,day,hour,min,sec] = datevec( ...
        datenum(1970,1,1)+times./86400);
    
% Get lat and lon out of file and convert to deg min sec  
latm = strcmp(fatt,'nominal_latittude');
lat = ncreadatt(fname,'/','nominal_latitude');
[lad,lam,las] = ddeg2dms(abs(lat));
lam = abs(lam);las = round(las);las = abs(las);
lad = repmat(lad,[length(year) 1]);lam = repmat(lam,[length(year) 1]);
las = repmat(las,[length(year) 1]);
lon = ncreadatt(fname,'/','nominal_longitude');
[lod,lom,los] = ddeg2dms(abs(lon));
lom = abs(lom);los = abs(los);los = round(los);
lod = repmat(lod,[length(year) 1]);lom = repmat(lom,[length(year) 1]);
los = repmat(los,[length(year) 1]);
if lon < 0
    lod = -lod;
end
if lat < 0
    lat = -lat;
end


% Get depth out of file
try
    dep = ncreadatt(fname,'/','sea_floor_depth_below_sea_level');
catch
    dep = -999.00;
end
dep = repmat(dep,[length(year) 1]);

% Get windspeed, windgust, wind direction
try
    wspd = ncread(fname,['/payload_1/anemometer_1/wind_speed']);
    payload = 'payload_1';
catch
    try
        wspd = ncread(fname,['/payload_2/anemometer_1/wind_speed']);
        payload = 'payload_2';
    catch
        payload = 'none';aa = 0;
        return
    end
        
end
try
    payatt = ncreadatt(fname,['/',payload],'description');
    np = length(payatt);sp = ['%',num2str(np),'s\n'];
    fprintf(1,'Using %9s \n',payload);
    fprintf(1,sp,payatt);
catch
    fprintf(1,'No Description of %9s \n',payload);
end

wspd_q = ncread(fname,['/',payload,'/anemometer_1/wind_speed_qc']);
ii = wspd_q ~= 0;
wspd(ii) = -99.99;
ig = wspd_q == 0;
if isempty(wspd(ig))
  wspd2 = ncread(fname,['/',payload,'/anemometer_2/wind_speed']);
  wspd_q2 = ncread(fname,['/',payload,'/anemometer_2/wind_speed_qc']);
  jj = wspd_q2(ii) == 0;
  wspd(ii(jj)) = wspd2(ii(jj));
end
wspd_info = ncinfo(fname,['/',payload,'/anemometer_1/wind_speed']);
ii = wspd == wspd_info.FillValue;
wspd(ii) = -99.99;

% Get anemometer height
try
    anemh = ncreadatt(fname,['/',payload,'/anemometer_1','height_of_instrument']);
catch
    anemh = 'NaN';
    re = 0;
end
if ~isnan(anemh)
    dp = [pdir,'alt_height.dat'];
    if exist(dp,'file')
        fid2 = fopen(dp);
        data = textscan(fid2,'%s%f');
        ii = strcmp(data{1},fname(end-21:end-17));
        anemh = data{2}(ii);
        fclose(fid2);
        if isempty(anemh)
            %ques = ['What is the anemometer height for buoy ',fname(end-21:end-17),': '];
            %anemh = input(ques);
            anemh = 5.0;
            fid2 = fopen(dp,'a');
            fprintf(fid2,'%5s %7.2f\n',fname(end-21:end-17),anemh);
            fclose(fid2);
        end
    else
        %ques = ['What is the anemometer height for buoy ',fname(end-21:end-17),': '];
        %anemh = input(ques);
        anemh = 5.0;
        fid2 = fopen(dp,'w');
        fprintf(fid2,'%5s %7.2f\n',fname(end-21:end-17),anemh);
        fclose(fid2);
    end
end
anemh = repmat(anemh,[length(year) 1]);

wgus = ncread(fname,['/',payload,'/anemometer_1/wind_gust']);
wgus_q = ncread(fname,['/',payload,'/anemometer_1/wind_gust_qc']);
ii = wgus_q ~= 0;
wgus(ii) = -99.99;
ig = wgus_q == 0;
if isempty(wgus(ig))
  wgus2 = ncread(fname,['/',payload,'/anemometer_2/wind_gust']);
  wgus_q2 = ncread(fname,['/',payload,'/anemometer_2/wind_gust_qc']);
  jj = wgus_q2(ii) == 0;
  wgus(ii(jj)) = wgus2(ii(jj));
end
wgus_info = ncinfo(fname,['/',payload,'/anemometer_1/wind_gust']);
ii = wgus == wgus_info.FillValue;
wgus(ii) = -99.99;

wdir = ncread(fname,['/',payload,'/anemometer_1/wind_direction']);
wdir_q = ncread(fname,['/',payload,'/anemometer_1/wind_direction_qc']);
ii = wdir_q ~= 0;
wdir(ii) = -999.0;
ig = wdir_q == 0;
if isempty(wdir(ig))
  wdir2 = ncread(fname,['/',payload,'/anemometer_2/wind_direction']);
  wdir_q2 = ncread(fname,['/',payload,'/anemometer_2/wind_direction_qc']);
  jj = wdir_q2(ii) == 0;
  wdir(ii(jj)) = wdir2(ii(jj));
end
wdir_info = ncinfo(fname,['/',payload,'/anemometer_1/wind_direction']);
ii = wdir == wdir_info.FillValue;
wdir(ii) = -999.0;
wdir = single(wdir);

% Get Air temperature and Sea tempurature
atemp = ncread(fname,['/',payload,'/air_temperature_sensor_1/air_temperature']);
atemp_q = ncread(fname,['/',payload,'/air_temperature_sensor_1/air_temperature_qc']);
atemp_info = ncinfo(fname,['/',payload,'/air_temperature_sensor_1/air_temperature']);
ii = atemp == atemp_info.FillValue;
atemp = atemp - 272.15;
atemp(ii) = -99.99;
ii = atemp_q ~= 0;
atemp(ii) = -99.99;
ig = atemp_q == 0;
try
    if isempty(atemp(ig))
    atemp2 = ncread(fname,['/',payload,'/air_temperature_sensor_2/air_temperature']);
    atemp_q2 = ncread(fname,['/',payload,'/air_temperature_sensor_2/air_temperature']);
    jj = atemp_q2(ii) == 0;
    atemp(ii(jj)) = atemp2(ii(jj));
    end
catch
end

try
    stemp = ncread(fname,['/',payload,'/ocean_temperature_sensor_1/sea_surface_temperature']);
    stemp_q = ncread(fname,['/',payload,'/ocean_temperature_sensor_1/sea_surface_temperature_qc']);
    stemp_info = ncinfo(fname,['/',payload,'/ocean_temperature_sensor_1/sea_surface_temperature']);
    ii = stemp == stemp_info.FillValue;
    stemp = stemp - 272.15;
    stemp(ii) = -99.99;
    ii = stemp_q ~= 0;
    stemp(ii) = -99.99;
catch
    stemp = repmat(-99.99,[length(times) 1]);
end
% Get barometric pressure from file
barp = ncread(fname,['/',payload,'/barometer_1/air_pressure_at_sea_level']);
barp_q = ncread(fname,['/',payload,'/barometer_1/air_pressure_at_sea_level_qc']);
barp_info = ncinfo(fname,['/',payload,'/barometer_1/air_pressure_at_sea_level']);
ii = barp == barp_info.FillValue;
barp = double(barp)./100;
barp(ii) = -999.;
ii = barp_q ~= 0;
barp(ii) = -999.;
ig = barp_q == 0;
if isempty(barp(ig))
    try
        barp2 = ncread(fname,['/',payload,'/barometer_2/air_pressure_at_sea_level']);
        barp_q2 = ncread(fname,['/',payload,'/barometer_2/air_pressure_at_sea_level_qc']);
        jj = barp_q2(ii) == 0;
        barp(ii(jj)) = barp2(ii(jj));
        ii = barp == barp_info.FillValue;
        barp = barp./100;
        barp(ii) = -999.;
    catch
    end
end

aa = [year,mon,day,hour,min,lad,lam,las,lod,lom,los,dep,anemh, ...
    double(wspd),double(wgus),double(wdir),double(atemp),double(stemp), ...
    double(barp)];