setup.directorymaker(logger, targetdir)
setup.directorymaker(logger, '{0}{1}/'.format(targetdir, 'fullpostfiles'))
setup.directorymaker(logger, '{0}{1}/'.format(targetdir, 'truncpostfiles'))
#####First we'll calculate your isotherm and write a pressure.dat file in your target directory
satP = setup.pSat(logger, species[-1], T)
istm = setup.isothermcalculator(logger, satP, iso_length, minrelpress, maxrelpress)
setup.PressureFileWriter(logger, species[-1], T, satP, istm, targetdir) #goes in the directory above your individual one, gets symbolic linked later
##### Now to make some control files for you
setup.GcmcControlChanger(logger, species[-1], sorb_el_list, T, n, framework, '{0}'.format(targetdir), n_iterations, restart, ctrl_file_name, pressure)
for i, value in enumerate(istm): #now I make the extra .ctr files to give you final .xyz files for each simulaiton
setup.GcmcControlChanger(logger, species[-1], sorb_el_list, T, '1', framework, '{0}'.format(targetdir), '1', 'RESTARTFILE {0}.{1}.res.{2}'.format(framework, species[-1], i+1), '{0}kpa_restart.ctr'.format(value), value) #makes your control files for all of your pressure points
setup.PostControlChanger(logger, species[-1], T, n, framework, '{0}'.format(targetdir), '0')
setup.PostControlChanger(logger, species[-1], T, n, framework, '{0}'.format(targetdir), '60')
setup.GcmcRunWriter(logger, species[-1], T, framework, parentdir, '{0}'.format(targetdir), istm, '')#the last variable is the relative location of your interactions files
setup.AtmAtmMover(logger, sorb_el_list, forcefield, coultype, hicut, frameowrk, targetdir)
setup.SorbSorbWriter(logger, species[-1], sorb_el_list,framework,targetdir, pmap, emap)
setup.IntraWriter(logger, species[-1], sorb_el_list, framework, targetdir)
logger.info('''
##################################################################
So I've done the following:
Created the directory {0} which will run your simulation.
Calculated and written an isotherm between {1} and {2} kPa, then placed in it in {0}.
Created a gcmc.ctr file, 2 post control files, and a run.gcmc in {0}.
Finally I wrote an atom_atom_file, sorb_sorb_file, and intramolecular_file into {0}.
You're good to go!
##################################################################
'''.format(targetdir, satP*10**minrelpress, satP*10**maxrelpress))
setup.PressureFileWriter(logger, species[i], T, satP, istm, '{0}{1:02d}/'.format(targetdir, j)) #goes in the directory above your individual one, gets symbolic linked later
workingdir += n_per_species
##### Now to make some control files for you
workingdir = 1
for i in range(1, n_species+1):
for directory in range(workingdir, workingdir+n_per_species):
setup.GcmcControlChanger(logger, species[i], sorb_el_list, T, n, framework, '{0}{1:02d}/'.format(targetdir, directory), n_iterations, restart, ctrl_file_name, pressure)
for j, value in enumerate(istm): #now I make the extra .ctr files to give you final .xyz files for each simulaiton
setup.GcmcControlChanger(logger, species[i], sorb_el_list, T, '1', framework, '{0}{1:02d}/'.format(targetdir, directory), '1', 'RESTARTFILE {0}.{1}.res.{2}'.format(framework, species[-1], j+1), '{0}kpa_restart.ctr'.format(value), value) #makes your control files for all of your pressure points
setup.PostControlChanger(logger, species[i], T, n, framework, '{0}{1:02d}/'.format(targetdir, directory), '0')
setup.PostControlChanger(logger, species[i], T, n, framework, '{0}{1:02d}/'.format(targetdir, directory), '60')
setup.GcmcRunWriter(logger, species[i], T, framework, parentdir, '{0}{1:02d}/'.format(targetdir, directory), istm, '')#the last variable is the relative location of your interactions files
setup.AtmAtmMover(logger, sorb_el_list, forcefield, coultype, hicut, frameowrk, '{0}{1:02d}/'.format(targetdir, directory))
setup.SorbSorbWriter(logger, species[i], sorb_el_list,framework,'{0}{1:02d}/'.format(targetdir, directory), pmap, emap)
setup.IntraWriter(logger, species[i], sorb_el_list, framework, '{0}{1:02d}/'.format(targetdir, directory))
setup.TaskfarmRunWriter(logger, 'multiple', T, framework, parentdir, targetdir, Taskfarm_total)
logger.info('''
##################################################################
So I've done the following:
Created the directory {0} and {1} subdirectories for taskfarming.
Calculated and written isotherms for all {2} species you've got, {3} times over each.
Created a gcmc.ctr file, 2 post control files, and a run.gcmc file for each subdirectory.
Created a taskfarm document and a backup, and finally run.taskfarmer in {0}. It might complain about not having a python script, but don't worry about that.
Finally I wrote an atom_atom_file, sorb_sorb_file, and intramolecular_file into each subdirectory.
You're good to go!
##################################################################
'''.format(targetdir, Taskfarm_total, n_species, n_per_species))
import logging
import datetime
import argparse
import musicpy.setup as setup
import musicpy.Antoine as Antoine
import musicpy.Forcefield as Forcefield
now = datetime.datetime.now()
#################
####This section sorts out your messages from this script to the console and a log file
logger = logging.getLogger(__name__)
logger.setLevel(logging.DEBUG)
handler = logging.FileHandler('{0}/sorb_sorb_maker.log'.format('.'))
handler.setLevel(logging.DEBUG)
handler2 = logging.StreamHandler()
handler2.setLevel(logging.INFO)
logger.addHandler(handler)
logger.addHandler(handler2)
logger.debug('-----------------------------------')
logger.debug('Generating new files on {0}-{1}-{2} at {3}:{4}'.format(
now.year, now.month, now.day, now.hour, now.minute))
logger.debug('-----------------------------------')
##################
####This section sorts out your list of species' and frameworks to amke an intramolecular file for you
species = ['XXXX']
framework = 'YYYY'
sorb_el_list = ['''### your element list here###''']
targetdirectory = './experiments/{0}/298/'.format(species)
##################
setup.IntraWriter(logger, species, sorb_el_list, framework, targetdirectory)