issue #27: rhocs last dimension is allocated incorrectly in EMfields.cpp