parametersout.update_parameters(parametersspawn)
# How much time slots do we need
tm = parametersout.get_timemanager()
slots = tm.compute_number_saves()
# Second IOM for output data of the spawning simulation
iomout = IOManager()
iomout.create_file(parametersout, filename="simulation_results_spawn.hdf5")
# Some data in the global data block
iomout.add_grid(parametersout, blockid="global")
iomout.save_grid(iomin.load_grid(blockid="global"), blockid="global")
# Allocate all the data blocks
for i in xrange(len(parametersout["spawn_components"])):
gid = iomout.create_group()
bid1 = iomout.create_block(groupid=gid)
bid2 = iomout.create_block(groupid=gid)
# Block for remainder / mother after spawning
iomout.add_wavepacket(parametersin, blockid=bid1)
# Block for spawned packet
iomout.add_wavepacket(parametersout, blockid=bid2)
# Really do the aposteriori spawning simulation
aposteriori_spawning(iomin, iomout, parametersin, parametersout)
# Close the inpout/output files
iomin.finalize()
iomout.finalize()