If you add f_Tinz = 'mxxxx' and f_Sinz = 'mxxxx' to your namelist this will output the 4D temperature (nx,ny,ncat,nilyr) and salinity. Also, the temperature profile in the snow, f_Tsnz. However, one has to be careful when averaging across categories making sure to do volume weighted.