function sta2ncml_tdsmeta_grid(fid,server,prefix,stafile,mvals,site);
% sta2ncml_tdsmeta: Write NcML to make time series file CF-compliant
% Usage: sta2ncml(fid,server,prefix,stafile,start,dt);
%  Inputs:
%   fid     = file id of text file *already opened for writing*
%   server  = first part of URL. For example:
%    'http://stellwagen.er.usgs.gov/thredds/dodsC/TSdata/'  % thredds OPeNDAP style
%    'http://stellwagen.er.usgs.gov/opendap/' % Hyrax OPeNDAP style
%    'http://stellwagen.er.usgs.gov/'                  % HTTP style
%   prefix  = The time series directory (next part of URL).
%             'MBAY_LT/';
%             'MYRTLEBEACH/';
%   stafile = name of NetCDF file:   '7201adc-a.nc';
%   mvals = a structure containing start stop times, lat, lon and timestep
%     mvals may be obtained by calling get_whfc_time()
%   site = a string containing the site ID (made up of mooring ID)
%
% Rich Signell (rsignell@usgs.gov)
% Notes on making our Time Series Data CF-compliant:
%   Time: Our time series data use non CF-compliant time convention,
%   with "time" variable containing Julian days and "time2" containing
%   milliseconds since midnight.  To make these CF-compliant, we use
%   NcML to construct a double precision time vector in Modified Julian Day,
%   one of the few internationally recognized time conventions that starts
%   at midnight. (See http://tycho.usno.navy.mil/mjd.html for more info.)
%
%   We specify %.16f precision, which will give us accuracy of:
%   - better than 3.0e-5 milliseconds until August 31, 2132 and
%   - better than 3.0e-4 milliseconds until October 12, 4596!
%
%   Lon,Lat and Depth: We give "lon" and "lat" CF-Compatible units of
%   "degrees_east" and "degrees_north" and give "depth" the attribute of
%   "positive=down" so that it recognized as a vertical coordinate.
start=(mvals.mjd_start-3600*24*1000*datenum(1858,11,17,0,0,0));  % convert datenum to MJD
location=[server prefix stafile];
stub=[prefix stafile];
fprintf(fid,'<dataset name="%s" ID="%s" urlPath="%s">\n',stafile,stub,stub);
% this line isn't needed if the experiment has this inherited
fprintf(fid,'<serviceName>manyservices</serviceName>\n');
% this inserts the THREDDS metadata
fprintf(fid,'<metadata inherited="false">\n');
fprintf(fid,'<geospatialCoverage>\n');
fprintf(fid,'<northsouth>\n');
fprintf(fid,'<start> %f </start>\n',mvals.lat);
fprintf(fid,'<size>0</size>\n');
fprintf(fid,'<units>degrees_north</units>\n');
fprintf(fid,'</northsouth>\n');
fprintf(fid,'<eastwest>\n');
fprintf(fid,'<start> %f </start>\n',mvals.lon);
fprintf(fid,'<size>0</size>\n');
fprintf(fid,'<units>degrees_west</units>\n');
fprintf(fid,'</eastwest>\n');
fprintf(fid,'</geospatialCoverage>\n');
fprintf(fid,'<timeCoverage>\n');
% in the TDS metadata, these must NOT be inside quotes
fprintf(fid,'<start>%s </start>\n',datestr(mvals.jd_start,31));
fprintf(fid,'<end> %s </end>\n',datestr(mvals.jd_stop,31));
fprintf(fid,'</timeCoverage>\n');
fprintf(fid,'<variables vocabulary="CF-1.4">\n');
for ik=1:length(mvals.varnames)
    % this inserts the cf_name obtained in lookup_cf if available
    %  if units are missing, it's a ' ' (space), so this patches most cases
    % with current GI-cat usage both should probably be cf_name
    % old way looks like:
    %  <variable name="S_40" vocabulary_name="SALINITY (PPT)" units="PPT">sea_water_salinity</variable>
    %    for it to find "sea_water_salinity", that string has to be the last element
    % based on the THREDDS metadata documentation, "sea_water_salinity" should
    %  really be in the vocabulary name, so probably it should go in both fields
    % 
    % this is the "broken" way that used to work in gi-cat
%     try
%     if isspace(mvals.units{ik})
%         fprintf(fid,'<variable name="%s" vocabulary_name="%s" units="%s">%s</variable>\n',...
%             mvals.varnames{ik},deblank(mvals.long_name{ik}),'unknown',deblank(mvals.cf_name{ik}));
%         
%     else
%         fprintf(fid,'<variable name="%s" vocabulary_name="%s" units="%s">%s</variable>\n',...
%             mvals.varnames{ik},deblank(mvals.long_name{ik}),mvals.units{ik},deblank(mvals.cf_name{ik}));
%     end
    %
    %  this is the swapped way
    try
    if isspace(mvals.units{ik})
        fprintf(fid,'<variable name="%s" vocabulary_name="%s" units="%s">%s</variable>\n',...
            mvals.varnames{ik},deblank(mvals.cf_name{ik}),'unknown',deblank(mvals.long_name{ik}));       
    else
        fprintf(fid,'<variable name="%s" vocabulary_name="%s" units="%s">%s</variable>\n',...
            mvals.varnames{ik},deblank(mvals.cf_name{ik}),mvals.units{ik},deblank(mvals.long_name{ik}));
    end
    
    catch ME
        disp ('we are here because it did not like the long_name')
        keyboard
    end
end
fprintf(fid,'</variables>');
fprintf(fid,'</metadata>\n');

% the netcdf line *MUST* follow the metadata
fprintf(fid,'<netcdf xmlns="http://www.unidata.ucar.edu/namespaces/netcdf/ncml-2.2" location="%s">\n',location);

% this part adds the cf-time
fprintf(fid,'<variable name="time" shape="time" type="double">\n');
fprintf(fid,'  <attribute name="units" value="milliseconds since 1858-11-17 00:00:00 UTC"/>\n');
fprintf(fid,'  <attribute name="long_name" value="Modified Julian Day"/>\n');
fprintf(fid,'  <values start="%.0f" increment="%.0f"/>\n',start,mvals.dtms);
fprintf(fid,'</variable>\n');
% get rid of old time2
fprintf(fid,'<remove type="variable" name="time2" />\n');
fprintf(fid,'<variable name="depth">\n');
fprintf(fid,'  <attribute name="positive" value="down"/>\n');
fprintf(fid,'  <attribute name="standard_name" value="depth"/>\n');
% standard_name of station_altitude is needed for point obs- don't use yet
fprintf(fid,'</variable>\n');
% lon lat (and depth for single depth instruments) should be shape ""
fprintf(fid,'<variable name="lon">\n');
fprintf(fid,'  <attribute name="units" value="degree_east"/>\n');
fprintf(fid,'  <attribute name="standard_name" value="longitude"/>\n');
fprintf(fid,'</variable>\n');
fprintf(fid,'<variable name="lat">\n');
fprintf(fid,'  <attribute name="units" value="degree_north"/>\n');
fprintf(fid,'  <attribute name="standard_name" value="latitude"/>\n');
fprintf(fid,'</variable>\n');
% now add the site id to allow it to by "station"
fprintf(fid,'<variable name="site" shape="" type="String">\n');
fprintf(fid,'<attribute name="standard_name" value="station_id"/>\n');
fprintf(fid,' <values>"%s"</values>\n', site);
fprintf(fid,' </variable>\n');
% for each of the non-coordinate variables, add the standard_name if
% available cf_sub is 1 when there was a CF name equivalant- do nothing
% when cf_sub is 0
for ik=1:length(mvals.varnames)
    if mvals.cf_sub{ik}==1
        eval(['fprintf(fid,''<variable name="' mvals.varnames{ik} '">\n'');'])
        eval(['fprintf(fid,''<attribute name="standard_name" value="' mvals.cf_name{ik} '"/>\n'');'])
        fprintf(fid,'</variable>\n');
    end
end
% add these necessary items for identification
fprintf(fid,'<attribute name="Conventions" value="CF-1.4" />\n');
% put the next two back in if needed as appropriate for "grid" 
%fprintf(fid,'<attribute name="CF:\\FeatureType" value="stationTimeSeries" />\n');
%fprintf(fid,'<attribute name="cdm_data_type" value="Station" />\n');
%
%here's more from the unidata dataset discovery attributes that are already
%in the thredds dataset, but that aren't already inherited
eval(['fprintf(fid,''<attribute name="id" value="' stub '" />\n'');'])
fprintf(fid,'<attribute name="naming_authority" value="gov.usgs.whsc" />\n');
fprintf(fid,'<attribute name="institution" value="U.S.Geological Survey, Woods Hole Coastal and Marine Science Center" />\n');
fprintf(fid,'<attribute name="creator_name" value="Ellyn Montgomery" />\n');
fprintf(fid,'<attribute name="publisher_name" value="USGS Coastal Marine Geology Program Sediment Transport group" />\n');
fprintf(fid,'<attribute name="publisher_url" value="http://stellwagen.er.usgs.gov" />\n');
fprintf(fid,'<attribute name="processing_level" value="Best Basic Version" />\n');
fprintf(fid,'<attribute name="project" value="USGS Coastal Marine Geology Program Sediment Transport group" />\n');
eval(['fprintf(fid,''<attribute name="contributor_name" value="' mvals.scipi '" />\n'');'])
fprintf(fid,'<attribute name="contributor_role" value="Project Investigator" />\n');
eval(['fprintf(fid,''<attribute name="title" value="' mvals.title '" />\n'');'])
eval(['fprintf(fid,''<attribute name="summary" value="' mvals.summary '" />\n'');'])
fprintf(fid,'<attribute name="standard_name_vocabulary" value="CF-1.4" />\n');
fprintf(fid,'<attribute name="Metadata_Conventions" value="Unidata Dataset Discovery v1.0" />\n');
eval(['fprintf(fid,''<attribute name="date_created" value="' datestr(now) '" />\n'');'])
% re-write start_time and stop time attributes in ISO 8601 format
fprintf(fid,'<attribute name="start_time" value="%sZ" />\n',datestr(mvals.jd_start,31));
fprintf(fid,'<attribute name="stop_time" value="%sZ" />\n',datestr(mvals.jd_stop,31));
fprintf(fid,'</netcdf>\n');
fprintf(fid,'</dataset>\n');