The fit:
Tsys(az,za,f=1400MHz) = 27.004869 + 0.081847676*za + 0.64943977*(za-15)2 - 0.12718878*(za-15)3 - 0.32397045*cos(az) - 0.19157567*sin(az) - 0.13125605*cos(2*az) + 0.11939164*sin(2*az) + 0.14413149*cos(3*az) - 0.22303953*sin(3*az)
Tsys(az,za,f=all) = 26.595679 + 0.094904966*za + 0.56760290*(za-15)2 - 0.10967648*(za-15)3 - 0.24777213*cos(az) - 0.11816161*sin(az) - 0.10067234*cos(2*az) + 0.085879611*sin(2*az) + 0.20042664*cos(3*az) - 0.23599053*sin(3*az)
Fit for f=all