cl()
r(e='rootname',c='rootname',s='root.dat')
sv('root.new')


#e = constants file
#c = c file
#s = solutions files which will be save
#sv = append the current run to fvdpmap4.new

# Second run if the first run was NOT adequate
cl()
r(e='rootname',c='rootname',s='root.new')
ap(root.new)
# ap saves the current run and appends the results to root.new
# To run this, you need to change the IRS command in the constants file
# See the HowtoRun file