var _____WB$wombat$assign$function_____=function(name){return (self._wb_wombat && self._wb_wombat.local_init && self._wb_wombat.local_init(name))||self[name];};if(!self.__WB_pmw){self.__WB_pmw=function(obj){this.__WB_source=obj;return this;}}{
let window = _____WB$wombat$assign$function_____("window");
let self = _____WB$wombat$assign$function_____("self");
let document = _____WB$wombat$assign$function_____("document");
let location = _____WB$wombat$assign$function_____("location");
let top = _____WB$wombat$assign$function_____("top");
let parent = _____WB$wombat$assign$function_____("parent");
let frames = _____WB$wombat$assign$function_____("frames");
let opens = _____WB$wombat$assign$function_____("opens");
function get_latest_met(info)
% get_latest_met(info)
%
% Last modified:  Time-stamp: <2025-09-14 17:21:11 haines>
%
% Abstract: Function that loads, appends new data, and saves with all data 
%
% Usage: get_latest_met(info)
% Processing steps
% (1) load previous month and current month data from netCDF 
% (2) append current to previous
% (3) subset data for latest 48 hours
% (4) calc u and v, and air_pressure_at_sea_level, and other times
% (5) create new file of global attributes, dimensions, and variables
% (6) save (ncclose) file
% Author: Sara Haines, Sep 14, 2025
  have_this_month = 0;
  have_prev_month = 0;
  x = dir([info.ncdir '/*.nc']);
  for i=1:length(x)
    % x(i).name
    if strfind(x(i).name, [info.platform_name '_' info.month_str '.nc'])
      have_this_month = 1;
      ncfile_this_month = [info.ncdir '/' x(i).name];
    elseif strfind(x(i).name, [info.platform_name '_' info.prev_month_str '.nc'])
      have_prev_month = 1;
      ncfile_prev_month = [info.ncdir '/' x(i).name];
    end
  end
  % have_this_month
  % have_prev_month
  % (1) load previous and current months of data 
  if (have_this_month & have_prev_month)
    % if have both append this month data to previous month data
    nc_prev = netcdf(ncfile_prev_month, 'nowrite');
    nc_this = netcdf(ncfile_this_month, 'nowrite');
    % get variables
    % (2) append current to previous
    dn = [nc_prev{'time'}(:); nc_this{'time'}(:)];
    % size(dn)
    awspd = [nc_prev{'WSPDA'}(:); nc_this{'WSPDA'}(:)];
    awdir = [nc_prev{'WDIRA'}(:); nc_this{'WDIRA'}(:)];
    wgust = [nc_prev{'WGSTA'}(:); nc_this{'WGSTA'}(:)];
    atemp = [nc_prev{'ATEMP'}(:); nc_this{'ATEMP'}(:)];
    rh = [nc_prev{'RH'}(:); nc_this{'RH'}(:)];
    abaro = [nc_prev{'BARO'}(:); nc_this{'BARO'}(:)];
  elseif (have_this_month & ~have_prev_month)
    % if only this month just only use this month
    nc_this = netcdf(ncfile_this_month, 'nowrite');
    % get variables
    dn = nc_this{'time'}(:);
    % size(dn)
    awspd = [nc_prev{'WSPDA'}(:)];
    awdir = [nc_prev{'WDIRA'}(:)];
    wgust = [nc_prev{'WGSTA'}(:)];
    atemp = [nc_prev{'ATEMP'}(:)];
    rh = [nc_prev{'RH'}(:)];
    abaro = [nc_prev{'BARO'}(:)];
  else
    % if don't have, don't write out new netcdf file
    fprintf(2, ['... ... no current met data to read\n']);
  end
  % (3) subset data to within the last 48 hours (2 days)
  if have_this_month
    current_dn = info.current_time_dn;
    two_days_ago = current_dn-2; % datenum is in units of days
    which_latest = find(two_days_ago<=dn & dn<current_dn);
    % whos
    dn = dn(which_latest);
    % size(dn)
    awspd = awspd(which_latest);
    awdir = awdir(which_latest);
    wgust = wgust(which_latest);
    atemp = atemp(which_latest);
    rh    = rh(which_latest);
    abaro = abaro(which_latest);
    % (4) calculations and other time variables for output
    % (4a) dew point from relative humidity 
    % procedure from paroscientific
    % http://www.paroscientific.com/dewpoint.htm
    % limits
    % 0 < atemp < 60 deg C
    % 0.1< rh < 100 %
    % 0 < adewp < 50 deg C
    % The uncertainty in the calculated dew point temperature is +/- 0.4 deg C.   
    % constants
    a = 17.27;
    b = 237.7; % deg C
    alpha = ((a*atemp)./(b+atemp)) + log(rh);
    adewp = ((a*alpha)./(b-alpha));
    % (4a) uv wind components (northward and eastward)
    r=pi/180;
    WVE = awspd .*(sin((awdir-180)*r));
    WVN = awspd .*(cos((awdir-180)*r));
    % (4b) air pressure at sea level
    % this procedure from class notes from Grant Petty
    % http://rain.aos.wisc.edu/~gpetty/aos330/lab1_slp_adjust.pdf
    atempK = atemp+273.1; % units deg K
    Rd = 287.1; % gas constant for dry air (Rd) units of J/(kg K)
    g = 9.807; % gravitional force units of m/s^2
    z0 = info.baro_ht; % height of measurement of abaro
    H = (Rd.*atempK)/g; % scale height depending on air temperature
    % air pressure at sea level units of mbar quivalent to hPa (hecto-Pascal)
    abaro_sealevel = abaro.*(exp(z0./H));
    xxx=abaro_sealevel;
    % [abaro(1:10) abaro_sealevel(1:10)]
    % (4c) time variables for wind product exercise
    %     (i) seconds since  2025-1-1 00:00:00 UTC
    seacoos_time_1970 = (dn-datenum(1970, 1, 1, 0, 0, 0))*24*60*60;
    seacoos_time_1970 = round(seacoos_time_1970);
    %     (ii) matlab datenum (source_time)
    % don't need to do anything. 
    % dn
    %     (iii) date_time_string
    date_time_string = datestr(dn,'yyyy-mm-dd HH:MM:SS' );
  end
  ncfile_v20 = [info.out_nc_v20_dir '/nccoos-' info.platform_name '-' info.package_name '-latest.nc'];
  % delete old netcdf file of latest data
  if exist(ncfile_v20, 'file')
    delete(ncfile_v20);
  end
  % (5) create new file of global attributes, dimensions, and variables
  % if have data then create new netcdf file and populate
  if have_this_month
    data_type = 'met';
    nc = netcdf(ncfile_v20, 'clobber');     
    netcdf_create_v20;
    nc = close(nc);    
  end
  % 
  if (have_this_month)
    nc_this = close(nc_this);
  end
  if (have_prev_month)
    nc_prev = close(nc_prev);
  end
}