For some reason my secs1d_dd_gummel_map code is not responding to changes in temperature enacted on the secs1d_physical_constants file (both are in the same directory). In other words when I change the temperature, the current densities remain unchanged. I was wondering if you could take a quick look at my slightly modified secs1d_dd_gummel_map.m file and see why this is. Any help is appreciated, thank you.
secs1d_dd_gummel_map_example.m