most of ic routines are now parallelized

however, there are not well checked and will probably fail
when giving directely an array of random numbers.

new scripts in examples/plots/star_formation_rate