ADDED   mttroot/mtt/bin/trans/ibg2abg_m
Index: mttroot/mtt/bin/trans/ibg2abg_m
==================================================================
--- mttroot/mtt/bin/trans/ibg2abg_m
+++ mttroot/mtt/bin/trans/ibg2abg_m
@@ -0,0 +1,122 @@
+#! /bin/sh
+
+     ###################################### 
+     ##### Model Transformation Tools #####
+     ######################################
+
+# Bourne shell script: ibg2abg_m
+# Based on rbg2abg_m
+
+while [ -n "`echo $1 | grep '^-'`" ]; do
+    case $1 in
+	-I )
+	    info=info ;;
+	-nounits )
+	    units=no ;;
+	*)
+	    echo "$1 is an invalid argument - ignoring" ;;
+    esac
+    shift
+done
+
+# Set up some vars
+sys=$1
+lbl_file=${sys}_lbl.txt
+abg_file=${sys}_abg.m
+err=mtt_error.txt
+log=ibg2abg_m.log
+
+# Remove the old log file
+rm -f ibg2abg_m.log
+rm -f ${abg_file}
+
+# Inform user
+echo Creating ${abg_file}
+
+# Use matrix manipulation to accomplish the transformation
+${MATRIX} > ${log} 2> ${err} <<EOF
+  name = '$1'
+  infofile = fopen('mtt_info.txt', 'w');
+  errorfile = fopen('mtt_error.txt', 'w');
+  
+  ## Interpret data from the ibg representation
+  ibonds = $1_ibg;
+  bonds = ibonds.bonds;
+
+  [bonds,components,n_vector_bonds] = ibg2abg(name,bonds, \
+					      infofile,errorfile)
+
+  ## Write the acausal bond graph m-file
+  write_abg(name,bonds,components,n_vector_bonds);
+EOF
+
+# Set the units for the abg file
+test_units()
+{
+      grep '^[\s]*[#|%]UNITS' < ${lbl_file} >/dev/null
+  if [ $? = "0" ]; then 
+    grep '^[\s]*[#|%]UNITS' < ${lbl_file} |\
+    gawk '{
+      printf("mtt_units.sh %s %s %s %s %s\n", sys, $2, $3, $4, $5)
+    }' sys=${sys}  | sh | grep ERROR
+    if [ $? = "0" ]; then
+      echo "    " *MTT_ERRROR: domains and units are not OK - exiting
+      exit 1  
+    else
+      echo "    " domains and units are OK 
+    fi
+  else
+    echo "  no domains or units declared"
+  fi
+} 
+
+check_ports_exist()
+{
+ declared_ports=`grep '^[\s]*[#|%]UNITS' < ${lbl_file} | gawk '{print $2}'`
+  for declared_port in $declared_ports; do
+    grep "${sys}\.ports\.${declared_port}\.type" ${abg_file} >/dev/null
+    if [ $? = "1" ]; then
+      echo "*MTT_ERROR: Units declared for non-existent port ${declared_port}"
+      exit 1
+    fi
+  done
+}
+
+set_units()
+{ 
+  grep '^[\s]*[#|%]UNITS' < ${lbl_file} |\
+  gawk '{
+    printf("  %s.ports.%s.domain = \"%s\";\n", sys, $2, $3);
+    printf("  %s.ports.%s.units.effort = \"%s\";\n", sys, $2, $4);
+    printf("  %s.ports.%s.units.flow = \"%s\";\n", sys, $2, $5);
+  }' sys=${sys}
+} 
+
+if [ -z "$units" ]; then
+  echo Checking port domains and units
+  check_ports_exist;
+  test_units;
+  echo "## Port domain and units" >> ${abg_file}
+  set_units >> ${abg_file}
+fi
+
+# Append any VAR declarations
+if [ -f "$1_lbl.txt" ]; then
+  echo "## Explicit variable declarations" >> ${abg_file}
+  grep '^[\s]*[%|#][V|P]AR' $1_lbl.txt | tr '%' '#' >> ${abg_file}    
+  grep '^[\s]*[%|#]NOT[V|P]AR' $1_lbl.txt | tr '%' '#' >> ${abg_file}    
+fi
+
+# Close off the function
+echo "endfunction" >> ${abg_file}
+
+# Errors and info
+if [ "$info" = "info" ]; then
+    cat mtt_info.txt
+fi
+
+if mtt_error mtt_error.txt; then
+    exit 0
+else
+    exit 1
+fi

ADDED   mttroot/mtt/bin/trans/m/ibg2abg.m
Index: mttroot/mtt/bin/trans/m/ibg2abg.m
==================================================================
--- mttroot/mtt/bin/trans/m/ibg2abg.m
+++ mttroot/mtt/bin/trans/m/ibg2abg.m
@@ -0,0 +1,171 @@
+## -*-octave-*-
+
+function [bonds,components,n_vector_bonds] = ibg2abg(name,bonds, \
+						     infofile, \
+						     errorfile)
+
+  ## Provide useful data for log
+  struct_levels_to_print = 5;
+
+  ## Find number of bonds
+  [n_bonds, junk] = size(struct_elements(bonds));
+
+  ## Get components
+  max_bonds_on_component = 0;
+
+  for i = 1:n_bonds
+    eval(sprintf("bond = bonds.b%i", i));
+
+    ## heads
+    head_type_name = bond.head.component;
+    head_type_name = deblank(split(head_type_name, ":"));
+    head_type = head_type_name(1,:);
+    head_name = head_type_name(2,:);
+    if (exist("component_struct"))
+      if (struct_contains(component_struct, head_name))
+	## create a copy to work on - clearer than eval(sprintf(...))
+	## (pointers would be better still!)
+	eval(sprintf("head = component_struct.%s;", head_name));
+      endif
+    endif
+    head.type = head_type;
+    if (!struct_contains(head,"bonds"))
+      head.n_bonds = 1;
+      head.bonds = [+i];
+      head.named_ports = !strcmp(bond.head.ports, "[]");
+    else
+      head.n_bonds = head.n_bonds + 1;
+      head.bonds = [head.bonds, +i];
+      head.named_ports = head.named_ports + !strcmp(bond.head.ports, \
+						    "[]");
+    endif
+    eval(sprintf("head.ports.bond%i.name = '%s'", i, bond.head.ports));
+    eval(sprintf("head.ports.bond%i.sign = '[in]'", i));
+    max_bonds_on_component = max(max_bonds_on_component, \
+				 head.n_bonds)
+    eval(sprintf("component_struct.%s = head", head_name));
+    clear head;
+
+    ## tails
+    tail_type_name = bond.tail.component;
+    tail_type_name = deblank(split(tail_type_name, ":"));
+    tail_type = tail_type_name(1,:);
+    tail_name = tail_type_name(2,:);
+    if (exist("component_struct"))
+      if (struct_contains(component_struct, tail_name))
+	## create a copy to work on - clearer than eval(sprintf(...))
+	tail = eval(sprintf("tail = component_struct.%s", tail_name));
+      endif
+    endif
+    tail.type = tail_type;
+    if (!struct_contains(tail,"bonds"))
+      tail.n_bonds = 1;
+      tail.bonds = [-i];
+      tail.named_ports = !strcmp(bond.tail.ports, "[]");
+    else
+      tail.n_bonds = tail.n_bonds + 1
+      tail.bonds = [tail.bonds, -i];
+      tail.named_ports = tail.named_ports + !strcmp(bond.tail.ports, \
+						    "[]");
+    endif
+    eval(sprintf("tail.ports.bond%i.name = '%s'", i, bond.tail.ports));
+    eval(sprintf("tail.ports.bond%i.sign = '[out]'", i));
+    max_bonds_on_component = max(max_bonds_on_component, \
+				 tail.n_bonds);
+    eval(sprintf("component_struct.%s = tail", tail_name));
+    clear tail;
+
+  endfor
+  
+  ## Get component indices from sys_cmp.m
+  [n_components, junk] = size(struct_elements(component_struct));
+  for i = 1:n_components
+    [comp_type, comp_name] = eval(sprintf("%s_cmp(%i)", name, i));
+    eval(sprintf("component_struct.%s.index = %i", comp_name, i));
+  endfor
+
+  ## Create the connections matrix (components)
+  components = zeros(n_components, max_bonds_on_component);
+  for [comp, comp_name] = component_struct
+    components(comp.index,1:comp.n_bonds) = comp.bonds;
+  endfor
+  
+  for [comp, comp_name] = component_struct
+    ## Unalias and/or default all ports on this component
+    if ((strcmp(deblank(comp.type), "0")) || (strcmp(deblank(comp.type), "1")))
+      disp("---- default junctions ----");
+      if (comp.named_ports == 1) # one named port
+	for [port, bond_number] = comp.ports
+	  if (!strcmp(port.name,"[]"))
+	    junction_port_name = port.name;
+	  endif
+	endfor
+	mtt_info(sprintf("Defaulting all ports on junction %s to %s", \
+			 comp_name, junction_port_name));
+	for [port, bond_number] = comp.ports
+	  port.name = junction_port_name;
+	endfor	
+      elseif ((comp.named_ports != 0) && (comp.named_ports != \
+					  comp.n_bonds)) # not allowed
+	mtt_error(sprintf("Junction %s must have 0, 1 or %i port labels", \
+			  comp_name, comp.n_bonds));
+      endif
+    else			# Not a junction
+      unlabelled_ports = comp.n_bonds - comp.named_ports;
+      if (unlabelled_ports == 1)
+	## find unlabelled ports
+	for [port, bond_number] = comp.ports
+	  if (strcmp(deblank(port.name), "[]"))
+	    ## found it - default to "in" or "out"
+	    eval(sprintf("comp.ports.%s.name = port.sign", bond_number));
+	    mtt_info(["Defaulting port name [" port.sign "]\t on component " \
+		      comp_name " (" deblank(comp.type) ")"], infofile);
+	  endif
+	endfor
+      elseif (unlabelled_ports == 2)
+	## find unlabelled ports
+	for [port, bond_number] = comp.ports
+	  if (strcmp(deblank(port.name), "[]"))
+	    ## got one
+	    if (exist("unlabelled_port1"))
+	      unlabelled_port2 = port;
+	      bond_number2 = bond_number;
+	    else
+	      unlabelled_port1 = port;
+	      bond_number1 = bond_number;
+	    endif
+	  endif
+	endfor
+	if (strcmp(unlabelled_port1.sign, "[in]") && \
+	    strcmp(unlabelled_port2.sign, "[in]"))
+	  mtt_error(["More than one unlabelled INport on component " \
+		     comp_name " (" deblank(comp.type) ")"], errorfile);
+	elseif (strcmp(unlabelled_port1.sign, "[out]") && \
+		strcmp(unlabelled_port2.sign, "[out]"))
+	  mtt_error(["More than one unlabelled OUTport on component " \
+		     comp_name " (" deblank(comp.type) ")"], errorfile);
+	else
+	  eval(sprintf("comp.ports.%s.name = \
+	  unlabelled_port1.sign", bond_number1));
+	  mtt_info(["Defaulting port name [" unlabelled_port1.sign
+		    "]\t on component "  comp_name " (" \
+		    deblank(comp.type) ")"], infofile);
+	  eval(sprintf("comp.ports.%s.name = \
+	  unlabelled_port2.sign", bond_number2));
+	  mtt_info(["Defaulting port name [" unlabelled_port2.sign
+		    "]\t on component "  comp_name " (" \
+		    deblank(comp.type) ")"], infofile);
+	endif
+	## ???	
+      endif
+    ## ???	
+    endif
+    ## ???
+    eval(sprintf("component_struct.%s = comp", comp_name));    
+  endfor
+component_struct;
+  ## ???
+  bonds;
+  components;
+  n_vector_bonds = 1;
+endfunction