Most recent update to network.f90 broke the indexing of the vdiff array in the thermal desorption branch.
Vdiff has length nspecies but we only need it for surface species.
Currently, if a grain species is at position i in the grain array and j in the species array, it's diffusion rate is placed at position i in the vdiff array.
Then, in rates.f90, it is accessed at position j so we get the wrong value.
Most recent update to network.f90 broke the indexing of the vdiff array in the thermal desorption branch.
Vdiff has length nspecies but we only need it for surface species.
Currently, if a grain species is at position i in the grain array and j in the species array, it's diffusion rate is placed at position i in the vdiff array.
Then, in rates.f90, it is accessed at position j so we get the wrong value.