function dbn = mk_orig_water_dbn
% Converted by Frank Hutter:
%   Provided in Phrog format by Xavier Boyen.
%   Manually converted into BIF.
%   Converted to BNT from BIF by Web-based bif2bnt (2004-01-30T05:28:10)
%   Manually converted into function creating DBN.
%   Manually changed the node numbering s.t. A-B-C-D-E-F-G-H from the BK paper correspond to 1-2-3-4-5-6-7-8

node = struct('C_NI_12_ANT', 1, ...
              'CKNI_12_ANT', 2, ...
              'CBODD_12_ANT', 3, ...
              'CNOD_12_ANT', 4, ...
              'CBODN_12_ANT', 5, ...
              'CNON_12_ANT', 6, ...
              'CKND_12_ANT', 7, ...
              'CKNN_12_ANT', 8, ...
              'C_NI_12_OBS', 9, ...
              'CKNI_12_OBS', 10, ...
              'CBODD_12_OBS', 11, ...
              'CNOD_12_OBS', 12, ...
              'CBODN_12_OBS', 13, ...
              'CNON_12_OBS', 14, ...
              'CKND_12_OBS', 15, ...
              'CKNN_12_OBS', 16, ...
              'C_NI_12_ULT', 17, ...
              'CKNI_12_ULT', 18, ...
              'CBODD_12_ULT', 19, ...
              'CNOD_12_ULT', 20, ...
              'CBODN_12_ULT', 21, ...
              'CNON_12_ULT', 22, ...
              'CKND_12_ULT', 23, ...
              'CKNN_12_ULT', 24);

adjacency = zeros(24);
adjacency([node.C_NI_12_ANT], node.C_NI_12_OBS) = 1;
adjacency([node.CKNI_12_ANT], node.CKNI_12_OBS) = 1;
adjacency([node.CBODD_12_ANT], node.CBODD_12_OBS) = 1;
adjacency([node.CKND_12_ANT], node.CKND_12_OBS) = 1;
adjacency([node.CNOD_12_ANT], node.CNOD_12_OBS) = 1;
adjacency([node.CBODN_12_ANT], node.CBODN_12_OBS) = 1;
adjacency([node.CKNN_12_ANT], node.CKNN_12_OBS) = 1;
adjacency([node.CNON_12_ANT], node.CNON_12_OBS) = 1;
adjacency([node.C_NI_12_ANT], node.C_NI_12_ULT) = 1;
adjacency([node.CKNI_12_ANT], node.CKNI_12_ULT) = 1;
adjacency([node.CBODN_12_ANT node.CNOD_12_ANT node.CBODD_12_ANT node.CKNI_12_ANT node.C_NI_12_ANT], node.CBODD_12_ULT) = 1;
adjacency([node.CKNN_12_ANT node.CKND_12_ANT node.CKNI_12_ANT], node.CKND_12_ULT) = 1;
adjacency([node.CNON_12_ANT node.CNOD_12_ANT node.CBODD_12_ANT], node.CNOD_12_ULT) = 1;
adjacency([node.CNON_12_ANT node.CBODN_12_ANT node.CBODD_12_ANT], node.CBODN_12_ULT) = 1;
adjacency([node.CKNN_12_ANT node.CKND_12_ANT], node.CKNN_12_ULT) = 1;
adjacency([node.CNON_12_ANT node.CKNN_12_ANT node.CBODN_12_ANT node.CNOD_12_ANT], node.CNON_12_ULT) = 1;

ss = 16;
dnodes = 1:ss;
ant = 1:8;
onodes = 9:16;
ult = 17:24;
intra = adjacency(1:ss, 1:ss);
inter_real = adjacency(ant, ult);
inter = zeros(ss);
inter(ant,1:length(ult)) = inter_real;