Fix number of work array components for spherical harmonics transform
Radial grid output is only called on rank 0
Reduce file access for memory conservation version of gen_sph_grid