;+ ;pro thm_sst_crib ; This is an example crib sheet that will load Solid State Telescope data. ; Open this file in a text editor and then use copy and paste to copy ; selected lines into an idl window. Or alternatively compile and run ; using the command: ; .RUN THM_SST_CRIB ;Author: Davin Larson ; ; $Id:$ ;- date = '2007-03-23' if 0 then setbp,/thisfile dprint,getdebug=db,setdebug=3,print_trace=3,print_dtime=1,/print_dlevel if keyword_set(date) then timespan,date,ndays else timespan,systime(1)-2*3600d*24 prb = ['a','b','c','d','e'] prb = 'c' thx = 'th'+prb loaded = 0 if not keyword_set(loaded) then begin thm_load_sst,probe=prb ;bpif keyword_set(dbg) mtypes = ['psif','psef'] ; SST thm_part_spec_calc,probe=prb,moments=['density','flux','mftens','velocity','T3'],instruments=mtypes thm_part_spec_calc,probe=prb,instruments=mtypes stop ;;thm_load_mom,probe=prb ;thm_load_hsk,probe=prb,varformat='th?_*sst* th?_imon_?5va th?_imon_?5va' ;thm_load_hsk,probe=prb,varformat='th?_*esa*' ;thm_load_fit, /get_support_data ;thm_cal_fit, /verbose ;thm_load_state2,probe=prb,/get_support,/polar ;tplot_options,var_label = 'th'+prb+'_state_pos_Re' ; mtypes = ['sif','sef'] ;thm_part_mom_calc,probe=prb,comps=['density','flux','mftens','velocity','T3'],types=mtypes ; ESA data: if 0 then begin ; ESA portion not working yet thm_load_esa_pkt,probe=prb mtypes = ['peif','peef'] ; ESA thm_part_spec_calc,probe=prb,moments=['density','flux','mftens','velocity','T3'],instruments=mtypes endif stop ; mom=moments_3d() ;store_data,'Tha_pexf_density',data=tnames('tha_pe?f_density') ;store_data,'Tha_pexf_velocity',data=tnames('tha_pe?f_velocity') ; ;store_data,'Tha_peix_density',data=tnames('tha_pei?_density') ;store_data,'Tha_peex_density',data=tnames('tha_pee?_density') ;store_data,'Tha_peix_flux',data=tnames('tha_pei?_flux') ;store_data,'Tha_peex_flux',data=tnames('tha_pee?_flux') ;store_data,'Tha_peix_mftens',data=tnames('tha_pei?_mftens') ;store_data,'Tha_peex_mftens',data=tnames('tha_pee?_mftens') ylim,'*_density',.01,100,1 ylim,'*_flux',-5e7,5e7,0 ylim,'*_mftens',-1000,5000 ylim,'*_eflux',-2e7,2e7 loaded = keyword_set(tnames()) endif ;tplot,'T*' tplot,'th?_ps?f_en' if 0 then begin ;tplot,'tha_pei[fm]*' ;bp ;tplot,/add,'T* ;tplot,'thb_pei[fm]*' ;bp ;tplot,'thc_pei[fm]*' ;bp ;tplot,'thd_pei[fm]*' ;bp ;tplot,'the_pei[fm]*' ;bp tplot,'tha_pee[fm]*' ;bp ;tplot,'thb_pee[fm]*' ;bp ;tplot,'thc_pee[fm]*' ;bp ;tplot,'thd_pee[fm]*' ;bp ;tplot,'the_pee[fm]*' ;bp ;mom= moments_3d() dtype = 'tha_peef' spec3d,thm_part_dist(dtype,gettime(/c)) plot3d_new,thm_part_dist(dtype,gettime(/c)) ;thm_load_mom_del,probe=prb ;thm_load_esa,probe=prb ;; for s=12,16 do for w=20,44,8 do print, (2 * 64 + ((s+16) mod 32) * 2)*256 + w,s,w, form='("1B",Z4," ; start=",i," length=",i)' endif dprint,setdebug=db dprint,'Done' end