Will let you know in a few days. Hope to have a play this weekend.
I think I last had a play a year or so back and I did the numpy int16 thing (taking care that the values fit in 65536 bins) and it just wrote out a BITPIX=32. I have a feeling I have to force its hand with BSCALE and BZERO explicitly plus there is some option for in effect, “Don’t mess about with BZERO and BSCALE”. It may be you have to – yourself – make sure the image counts are in the range -32768 to 32767…
If you had code to hand I would be keen to see it.