timespan,'7 7 27',2,/days
thm_load_sst
;
options,'*',spec=0,ylog=1 ; easier to see lineplot
;options,'*',spec=1,ylog=0   ; prettier to see specplot
tplot,'th?_psif_en th?_psir_en'
;tplot,'th?_psef_en th?_pser_en'