adjusting spitzer conductivity

#!/usr/bin/env python

from nuker_catalog import galaxy
import numpy as np
import dill

#gal=galaxy.NukerGalaxy.from_dir('NGC4551','/Users/aleksey/Second_Year_Project/hydro/batch_collected/NGC4551/vw_500.0')
gal=dill.load(open('/home/aleksey/Second_Year_Project/hydro/batch_collected/NGC4551/vw_500.0/grid.p','rb'))
gal.set_param('cond_scheme','spitzer')
# gal.set_param('eps',0.)
# gal.set_param('bdry','bp')
gal.set_param('bdry', 's_fixed')


# gal=dill.load(open('batch_collected/NGC4551/vw_500.0/grid.p','rb'))

gal.set_param('outdir', '/home/aleksey/Second_Year_Project/hydro/batch_collected/NGC4551/vw_500.0_spitzer_sfixed_2')
# gal.set_param('eps_cond',1.)
gal.set_param('tinterval',-1)
gal.set_param('sinterval',100)
# gal.set_temp(1.E7)

gal.set_param('eps_cond',0.)
gal.solve_adjust(3.*gal.tcross, 'eps_cond', 3.E-2, n=30)
# gal.solve_adjust(gal.tcross, 'eps_cond', 1.E-1)
# gal.solve_adjust(gal.tcross, 'eps_cond', 1.)