heating=4.74e4
varfile='var.dat'
pc_read_var,obj=obj,varfile=varfile,varia=['tt','pp','rho'],/magic,/add,/trimall
;plot,obj.TT,obj.pp,ps=1,yst=3
plot_oi,obj.TT,obj.cooling-heating,ps=1,yst=3,xr=[100,1e4]
oplot,obj.TT,obj.TT*0.
END