Could you add a comment, saying that you are assuming the DEFAULT_GMF_REALIZATIONS is equal to 1 (or equivalently that ```gmf.ground_motion_fields``` is returning a list with a single 2d numpy array where the first dimension is equal to 1)?