Problems expanding univariate liability threshold script to bivariate and multivariate script

Hi,
I'm fairly new to OpenMx and twin modelling. I have completed the twin modelling summer school at the SGDP this year, where we have learnt how to run a univariate liability threshold model. We have now tried to expand this script to a bivariate model, however we have only been able to make it work until the saturated model. Would anyone be able to let us know how to further expand the script so that it also works for the remaining ACE, AE etc models and what else we need to adjust so that it work for a multivariate (common pathway) liability threshold model.
Please find the script attached.
Thanks a lot for all your help. It is very much appreciated.
Best,
Isabel