Axisymmetric models for neutron star merger remnants with realistic thermal and rotational profiles