In: Numerical Recipes in Fortran 90

ed of their indices. This is why the standard makes no statement about the order in which the individual operations in an array expression are executed; they might in fact be carried out simultaneously, on parallel hardware. By default, array expressions and assignments are performed for all elements of the same-shaped arrays referenced. This can be modified, however, by use of a where construction like this: where (harr > 0.0_sp) farr = 3.0_sp*garr + 1.0_sp