from astrochem.wrapper import *
import numpy as np

p = phys()
initial_abundances = {}
c = cell(av=20, nh=1e3, tgas=10, tdust=10)
s = solver(c,  "network.chm", p , 1e-15, 1e-6, initial_abundances, 1e3, 0)