I'm an undergrad and as a part of my final project I need to simulate Cooke and Deserno's model for self assembling bilayer membranes (2005). I've been trying to get the example script for simple bilayer membrane in the mbtools package to run and received the following error.
:: > loaded version 1.0.0 of analysis
:: > loaded version 1.3.3 of cmdline
:: > loaded version 1.0.0 of system_generation
:: > using paramsfile: simplebilayer.tcl
:: > loading parameter file .. done
::mbtools::utils > Warning: no checkpoint named ./simplebilayer//checkpoint.latest.gz
::mbtools::utils > setting up ./simplebilayer/ for simplebilayer.tcl
::mbtools::utils > set interaction: 0 FENE 30 1.5
::mbtools::utils > set interaction: 1 harmonic 10.0 4.0
::mbtools::utils > set interaction: 0 0 lennard-jones 1.0 0.95 1.066375 0.25 0.0
::mbtools::utils > set interaction: 0 1 lennard-jones 1.0 0.95 1.066375 0.25 0.0
::mbtools::utils > set interaction: 0 2 lennard-jones 1.0 0.95 1.066375 0.25 0.0
::mbtools::utils > set interaction: 1 1 lj-cos2 1.0 1.0 0.0 1.6
::mbtools::utils > set interaction: 1 2 lj-cos2 1.0 1.0 0.0 1.6
::mbtools::utils > set interaction: 2 2 lj-cos2 1.0 1.0 0.0 1.6
::mbtools::system_generation > setting up system
::mbtools::system_generation > creating a topology for moltype 0
::mbtools::system_generation::flat > placing lipids in a flat bilayer
::mbtools::system_generation > Warning: no trappedmols defined
::mbtools::utils > Setting the random seed to clock in seconds: 356729693
::mbtools::utils > initializing vmd to : offline
:: > warming up at 1.1
::mbtools::utils > warming up 20 times 100 timesteps
::mbtools::utils > wrote file ./simplebilayer//warm.0000.gz
::mbtools::utils > run 0 of 20 at time=0.20000000000000015 (cap=54.75)
::mbtools::utils > run 1 of 20 at time=0.4000000000000003 (cap=104.5)
::mbtools::utils > run 2 of 20 at time=0.6000000000000004 (cap=154.25)
::mbtools::utils > run 3 of 20 at time=0.8000000000000006 (cap=204.0)
::mbtools::utils > run 4 of 20 at time=1.0000000000000007 (cap=253.75)
::mbtools::utils > run 5 of 20 at time=1.2000000000000008 (cap=303.5)
::mbtools::utils > run 6 of 20 at time=1.400000000000001 (cap=353.25)
::mbtools::utils > run 7 of 20 at time=1.6000000000000012 (cap=403.0)
::mbtools::utils > run 8 of 20 at time=1.8000000000000014 (cap=452.75)
::mbtools::utils > run 9 of 20 at time=2.0000000000000013 (cap=502.5)
::mbtools::utils > wrote file ./simplebilayer//warm.0001.gz
::mbtools::utils > run 10 of 20 at time=2.1999999999999793 (cap=552.25)
::mbtools::utils > run 11 of 20 at time=2.3999999999999573 (cap=602.0)
::mbtools::utils > run 12 of 20 at time=2.5999999999999353 (cap=651.75)
::mbtools::utils > run 13 of 20 at time=2.7999999999999132 (cap=701.5)
::mbtools::utils > run 14 of 20 at time=2.999999999999891 (cap=751.25)
::mbtools::utils > run 15 of 20 at time=3.199999999999869 (cap=801.0)
::mbtools::utils > run 16 of 20 at time=3.399999999999847 (cap=850.75)
::mbtools::utils > run 17 of 20 at time=3.599999999999825 (cap=900.5)
::mbtools::utils > run 18 of 20 at time=3.799999999999803 (cap=950.25)
::mbtools::utils > run 19 of 20 at time=3.999999999999781 (cap=1000.0)
::mbtools::utils > uncapping forces
::mbtools::system_generation > Warning: no user fixed particles defined
:: > warming up again at 1.1
::mbtools::utils > warming up 1 times 0 timesteps
::mbtools::utils > Warning: warmup steps are zero
:: > starting integration: run 100 times 100 steps
::mbtools::analysis > setting up analysis
The operation "get_lipid_orients" you requested is not implemented.
while executing
"analyze get_lipid_orients "
(procedure "setup_flipflop" line 30)
invoked from within
"setup_flipflop"
("eval" body line 1)
invoked from within
"eval $setupcommand"
(procedure "::mbtools::analysis::setup_analysis" line 61)
invoked from within
"::mbtools::analysis::setup_analysis $analysis_flags -outputdir $outputdir -g $mgrid -str $stray_cut_off"
(file "./scripts/main.tcl" line 236)