Hello,
I think the unit specified in comments for rcp here is wrong. IMO it should be [J/(kg*K)], not [J/K]. The value is OK though.
This has no implications whatsoever on the actual model, but I think it has some for users that want to do some (potentially pen-and-paper) diagnostics. I typically use the code as a reference relying on units provided in the source comments.