Make stats.F90 handle single-repeat analysis

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

Make stats.F90 handle single-repeat analysis #27

Closed mwest1066 closed 7 years ago

mwest1066 commented 9 years ago

Running 1_urban_plume with n_repeat 1 currently results in:

ERROR (PartMC-359779741): assertion failed

This is due to stats.F90 calculating student_t_95_coeff() with just one data point.

The statistics generation should be changed to do something sensible in this case, so that the same post-processing code can handle both multi- and single-repeat runs.

mwest1066 commented 7 years ago

Fixed in ffd6550b9126dbf7402290324757eab6bdd0efbb