compdyn / partmc

Particle-resolved stochastic atmospheric aerosol model
http://lagrange.mechse.illinois.edu/partmc/
GNU General Public License v2.0
28 stars 16 forks source link

input_n_files does not free unit #69

Open jcurtis2 opened 6 years ago

jcurtis2 commented 6 years ago

When doing processing of a large scenario library, we exceed the number of available IO units (which is hardcoded to be 200) when calling input_n_files for every scenario. In input_n_files, we call get_unit() but never actually free it. I don't think there is any reason to not free it here. This is solved by freeing the unit (call free_unit(unit)) before the end of the subroutine.

If we don't want to free this, people should hardcode their n_repeat and n_index for their large scenario library processing instead of calling input_n_files.