function mod_ge_output(filename,output,pltparams, varargin)
% 
% % Please use the link below to view the documentation.

% % If you don't see the link which is automatically inserted by MATLAB, 
% % try removing the line break above, so that this help block is merged 
% % with the previous one.
% % Reference page in help browser: <a href="matlab:doc ge_output">ge_output</a>
%
% etm mods lines 36-44 to define a style to NOT get directions in the
% balloons
%  no time slider in this implementation

AuthorizedOptions = authoptions('ge_output' );

% Assign default values to parameters:
       name = filename;
msgToScreen = false;
    
parsepairs %script that parses Parameter/value pairs.

if msgToScreen
   disp(['Running ' mfilename '...']) 
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

fileoutput = filename( (end-2):end );

if strcmp(fileoutput,'kml')

    fid = fopen(filename, 'wt');
 
 % header includes style definitions for use by the icons & baloons,
 % and the USGS overlay, and a helpful note about using the time-slider
  header = ['<?xml version="1.0" encoding="UTF-8"?>',10,... 
             '<kml xmlns="http://www.opengis.net/kml/2.2"',10,...
             'xmlns:gx="http://www.google.com/kml/ext/2.2">',10,...
             '<Document>',10,...
             '<name>',10,name,10,'</name>',10,...
             '<Style id="displayInfo">',10,...
                '<BalloonStyle>',10,...
                 '<text>',10,...
                  '<![CDATA[',10,...
                  'summary of contents of $[name]:',10,... 
                  '$[description]',10,...
                  ']]>',10,...
                  '</text></BalloonStyle>',10,...
                  '</Style>',10,...
               '<ScreenOverlay>',10,...                
                  '<name>USGS Logo</name>',10,...
                  '<description>  Created by the USGS Coastal Marine Geology Progrm</description>',10,...
                  '<drawOrder> 1 </drawOrder>',10,...
                   '  <Icon>',10,...
                   '   <href>http://stellwagen.er.usgs.gov/images/header_graphic_usgsIdentifier_green.jpg</href>',10,...
                   '  </Icon>',10,...
                  '  <overlayXY x="0" y="1" xunits="fraction" yunits="fraction"/>',10,...
                  '  <screenXY x=".009" y=".982" xunits="fraction" yunits="fraction"/>',10,...
                  '  <size x="168" y="62" xunits="pixels" yunits="pixels"/>',10,...
                 '</ScreenOverlay>',10,...
                '<Folder id="inst" >',10,...
                 '<name>inst_id</name>',10,...
                 '  <open>1</open>' ];
         
    footer = [10, '</Folder>',10,...
        '</Document>',10,...
        '</kml>'];

    fprintf(fid,'%s',header);
    fprintf(fid,'%s',output);
    fprintf(fid,'%s',footer);
    fclose(fid);

elseif strcmp(fileoutput,'kmz')
  
    newfilename = filename(1:(end-4));
    newfilename = strcat(newfilename,'.kml');
    
    fid = fopen( newfilename, 'wt');

    header = ['<?xml version="1.0" encoding="UTF-8"?>',10,...
             '<kml xmlns="http://earth.google.com/kml/2.1">',10,...
             '<Document>',10,...
             '<name>',10,name,10,'</name>',10];
    footer = [10,'</Document>',10,...
        '</kml>',10];

    fprintf(fid,'%s',header);
    fprintf(fid,'%s',output);
    fprintf(fid,'%s',footer);
    fclose(fid);
    
    zip(filename,newfilename);
    
    if ispc
        eval(['!del ',newfilename] );
        eval(['!move ',filename, '.zip ', filename])
    else
        eval(['!rm -f ',newfilename]);
        eval(['!mv ',filename, '.zip ', filename])
    end
   
else
    error('unsupported file extension used in ge_output');
end


if msgToScreen
   disp(['Running ' mfilename '...Done']) 
end