% do_agg2ncml: Script to generate NcML aggregation of WHSC Time Series Data

server='http://stellwagen.er.usgs.gov/cgi-bin/nph-dods/'; % URL part 1
prefix='DATAFILES/MBAY_LT/';  % URL part 2
url_base=[server prefix];

ncml_file='nrsfc_ts.ncml';  %output NcML file 

% grab all the files in this directory...
% files=get_files_from_url(url_base);

%... or use specific files
%files={'3392-a1h.nc','3412-a1h.nc','3482-a1h.nc'};

% all needed for MBAY_LT 5m Salinity
 files={'3391ats-a1h.nc','3411sc-a1h.nc','3481sc-a1h.nc','3591sc-a1h.nc',...
     '3751sc-a1h.nc','3841sc-a1h.nc','3971sc-a1h.nc','4011sc-a1h.nc',...
     '4081sc-a1h.nc','4141sc-a1h.nc','4211sc-a1h.nc','4291sc-a1h.nc',...
     '4311sc-a1h.nc','4522sc-a1h.nc','4591sc-a1h.nc','4631sc-a1h.nc',...
     '4661sc-a1h.nc','4681sc-a1h.nc','4771sc-a1h_d1.nc','5291mc-a1h_d1.nc',...
     '5411mc-a1h.nc','5511mc-a1h.nc','5701mc-a1h.nc','5901mc-a1h.nc',...
     '6121mc-a1h.nc','6241mc-a1h.nc','6311mc-a1h.nc','6371mc-a1h.nc',...
     '6441mc-a1h.nc','6641mc-a1h.nc',...
     '6821mc-a1h.nc','6891mc-a1h.nc','6961mc-a1h.nc','7071mc-a1h.nc',...
     '7161mc-a1h.nc','7551mc-a1h.nc','7661mc-a1h.nc','7741mc-a1h.nc',...
     '7761mc-a1h.nc','7871mc-a1h.nc','7921mc-a1h.nc'};

% these were all thought to be non uniform
% files={'5701mc-a1h.nc','5901mc-a1h.nc','6121mc-a1h.nc','7071mc-a1h.nc'};

fid=fopen(ncml_file,'wt');
% Could write the dataset info here, but for now, we find it easier to just
% write the NetCDF dataset aggregation so it can be loaded as NcML without
% THREDDS, and add the <dataset> part later.

fprintf(fid,'<netcdf xmlns="http://www.unidata.ucar.edu/namespaces/netcdf/ncml-2.2">\n');
fprintf(fid,'<aggregation dimName="time" type="joinExisting">\n');

for i=1:length(files)
    stafile=char(files(i))
    %%
    %if ~isempty(strfind(lower(stafile),'a1h'));  %only filenames with
    % 'a1h' in them (hourly averaged)
    url=[url_base stafile];
    % [mvals]=get_whfc_meta(url);
    [mvals]=get_meta_ncml(url);     % uses different uniformity assessment
    if mvals.timebase_uniform
     
        agg_sta2ncml(fid,server,prefix,stafile,mvals);
    else
        disp(['non_uniform timebase for ',stafile])
    end
    %end
    %%
end
fprintf(fid,'</aggregation>\n');
fprintf(fid,'<variable name="depth">\n');
fprintf(fid,'  <attribute name="positive" value="down"/>\n');
fprintf(fid,'</variable>\n');
fprintf(fid,'<variable name="lon">\n');
fprintf(fid,'  <attribute name="units" value="degrees_east"/>\n');
fprintf(fid,'</variable>\n');
fprintf(fid,'<variable name="lat">\n');
fprintf(fid,'  <attribute name="units" value="degrees_north"/>\n');
fprintf(fid,'</variable>\n');
fprintf(fid,'<attribute name="Conventions" value="CF-1.0" />\n');
fprintf(fid,'</netcdf>\n');

fclose(fid);