%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Load Delft3D gridfile in Matlab and export to PI XML format for FEWS
% input
% - Required availability of Delft3D-Quickplot
% - The number of rows and columns need to match with the grid size
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

d3d_qp('openfile','lp02.grd')

d3d_qp('selectfield','morphologic grid')

d3d_qp('loaddata')

fid = fopen('grids.xml','w');


fprintf(fid,'%s\n','<?xml version="1.0" encoding="UTF-8"?>');
fprintf(fid,'%s\n','<!-- edited with XMLSpy v2009 sp1 (http://www.altova.com) by ICT (Stichting Deltares) -->');
fprintf(fid,'%s\n','<grids xmlns="http://www.wldelft.nl/fews" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.wldelft.nl/fews grids.xsd">');
fprintf(fid,'%s\n','	<irregular locationId="Delft3DFLOW">');
fprintf(fid,'%s\n','		<description>Delft3D-FLOW model</description>');
fprintf(fid,'%s\n','		<rows>67</rows>');
fprintf(fid,'%s\n','		<columns>63</columns>');
fprintf(fid,'%s\n','		<geoDatum>SVY21</geoDatum>');


for yy=1:numel(data.X(1,:)),

for xx=1:numel(data.X(:,1)),

        fprintf(fid,'%s\n','		<cellCentre>');
        fprintf(fid,'%s\n',['			<x>',num2str(data.X(xx,yy)),'</x>']);
        fprintf(fid,'%s\n',['			<y>',num2str(data.Y(xx,yy)),'</y>']);
        fprintf(fid,'%s\n','		</cellCentre>');

    end
end

fprintf(fid,'%s\n','	</irregular>');
fprintf(fid,'%s\n','</grids>');

fclose all