def mini_topsim(config_file=None):
"""
Loads parameters from config_file, starts the sim, plots and writes to file
:param config_file: config_file with simulation parameters
Loads parameters from config_file.
If no config_file is passed passed, None is returned.
Creates a Surface Object and starts the simulation.
The correct timestep is calculated with the timestep function
from the advance module.
If a *.srf_save file with the same filename exists, the plot function with
both surfaces is called.
"""
print('Running miniTopSim ...')
if config_file is None:
if len(sys.argv) > 1:
config_filename = sys.argv[1]
else:
sys.exit('No Config file passed')
#config_filename = 'cosine.cfg'
config_file = config_filename
if not config_file.endswith('.cfg'):
print('Error: Incorrect config.')
sys.exit()
filename = os.path.splitext(config_file)[0] + '.srf'
if os.path.exists(filename):
os.remove(filename)
par.load_parameters(config_file)
dir_path = os.path.dirname(os.path.realpath(config_file))
par.INITIAL_SURFACE_FILE = os.path.join(dir_path, par.INITIAL_SURFACE_FILE)
tend = par.TOTAL_TIME
dt = par.TIME_STEP
surface = Surface()
sputter.init_sputtering()
time = 0
start_simulation_time = currenttime()
while time < tend:
surface.write(time, filename)
dtime = timestep(dt, time, tend)
advance(surface, dtime)
surface.eliminate_overhangs()
time += dtime
stop_simulation_time = currenttime()
simulation_time = stop_simulation_time - start_simulation_time
print('The Simulation took: {}s'.format(float(simulation_time)))
surface.write(time, filename)
filename_save = filename + '_save'
if par.PLOT_SURFACE:
if os.path.exists(filename_save):
print('*.srf_save file exists... plotting both!')
plot.plot(filename, filename_save)
else:
plot.plot(filename)
config_filename="test_overhangs.cfg"
config_file = os.path.join(os.path.dirname(__file__), config_filename)
if not config_file.endswith('.cfg'):
print('Error: Incorrect config.')
sys.exit()
filename = os.path.splitext(config_file)[0] + '.srf'
if os.path.exists(filename):
os.remove(filename)
par.load_parameters(config_file)
#Example Surface
surface = Surface()
surface.xvals=[0.,1.,2.,1.,1.5,3.,4.,5.,4.,5.,6.,7.,8.,7.,7.5,8.,9.]
surface.yvals=[10.,10.,10.,15.,17.,18.,18.,18.,12.,13.,13.,12.,11.,10.,9.,8.,8.]
plt.title('Eliminate Overhangs Test')
plt.ylabel('Surface Y values')
plt.xlabel('Surface X values')
plt.plot(surface.xvals, surface.yvals, 'b-', label='With Overhangs')
surface.eliminate_overhangs()
plt.plot(surface.xvals, surface.yvals, 'r-', label='Without Overhangs')
plt.legend()
plt.show()
