# make cubic membrane, visualize and write on file
import membrane.membrane as mem
m = mem.Membrane(unit_cell_size_=200)
# choose the form between "G", "P", "D", "micelle", "reveerse-micelle"
m.form = "G"
# start making the system with 100% glycoderamides
m.MakeMartini3(glyco=1.0)
m.Show()
m.WriteOnFileGMX("cubic_g_20_gly1.0.pdb")
#m.PrintInfo()
print "number of molecules: ", len(m.molecules)
print "sample_size=", m.sampling_size
# Making vescular systems
import membrane.membrane as mem
m=mem.Membrane(unit_cell_size_ = 150)
m.step_num = 16
m.micelle_radius = 50
m.MakeMartini8(1.0)
m.AddWater()
m.Show()
m.WriteOnFileGMX("vesicle_s15_5.pdb",make_TER=True)
import membrane.membrane as mem
import molar.pdb as pdb
m=mem.Membrane()
m.ReadFileGMX("vesicle_s15_5_no-water.pdb")
m.form = "vesicle"
m.side_length = 150
m.AddWater()
m.WriteOnFileGMX("vesicle_s15_5.pdb",make_TER=True)