http://scv.bu.edu/computation/bluegene/IBMdocs/compiler/xlf-10.1/html/xlflr/where.htm
https://stackoverflow.com/questions/24882425/openmp-array-reductions-with-fortran