Add zonal mean output capability to simple models.

### What is the feature/what would you like to discuss?

With the addition of the `phys_grid_ctem` module it is now possible to output zonal mean diagnostics when using non-FV dycores. However the functionality has only been added for output from the `cam[4-6]` and `cam_dev` physics packages, and not for the simple physics models. It would be useful to have this capability added for simple physics. It just requires the addition of a couple of `phys_grid_ctem` subroutine calls to `src/physics/simple/physpkg.F90`, analogous to the current calls in `src/physics/cam/physpkg.F90`.

### Is there anyone in particular you want to be part of this conversation?

_No response_

### Will this change (regression test) answers?

No

### Will you be implementing this enhancement yourself?

Any CAM SE can do this