pro read_h2 ;+ ; ; wkp 3/28/08 read in the names of a h2 cdf files in the ; current directory and plot out all the variables ; see the assoicated cdf skeleton table for information about ; the variables. ; ;- ; cmnd='ls -l po_h2*.cdf >! cdf_filelist' spawn, cmnd openr,lun,'cdf_filelist',/get_lun file_name='' while not eof(lun) do begin readf,lun,file_name ;-rw-rw-r-- 1 plasmas 2606101 Apr 2 13:46 /timas/flowdb/1996/po_h1_tim_19960317_v01.cdf ;1234567890123456789012345678901234567890123456789012345678901234 infile=strmid(file_name,54,41) print, infile cdfid=cdf_open(infile) print,'CDFID ',cdfid info=cdf_inquire(cdfid) help,/structure, info print, 'Dimensions: ',info.dim print, '' print,'Z Variable names and Types' for i=0, info.nzvars-1 do begin vinfo= cdf_varinq(cdfid,i,/zvariable) print, vinfo.name, ' Type:', vinfo.datatype, ' DIMS:', vinfo.DIM endfor cdf_control,cdfid,variable=0,/zvariable,get_var_info=einfo max_records=einfo.maxrec print, 'NOTE: The NO_PADVALUE_SPECIFIED: error is a feature of idl. no problem/wkp.' ;% CDF_CONTROL: Function completed but: NO_PADVALUE_SPECIFIED: A pad value has ; not been specified. print,'MAXRECORDS= ',max_records if max_records eq 0 then begin print, 'No valid Summary Data for this UT day' print_null,cdfid,einfo ; call to print_null.pro ???? print,cdfname stop endif if max_records lt 0 then begin print, 'ERRORR NO FILL RECORD IN THIS CDF' print, 'KEY PARAMETER CDF NOT CLOSED when generated' print, ' ..... or .......' print, 'you forgot to specify bin when you used ftp ' print, 'to transfer the file here! ' print, cdfname stop endif i=0 ; make the record star===0 (insruance statement status='' ; epoch ; ;H_Epoch Type:CDF_EPOCH DIMS: 0 cdf_varget,cdfid,'H_Epoch',H_Epoch, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for H_Epoch call/',Result,'/', status ;O_Epoch Type:CDF_EPOCH DIMS: 0 cdf_varget,cdfid,'O_Epoch',O_Epoch, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for O_Epoch call/',Result,'/', status ;He1_Epoch Type:CDF_EPOCH DIMS: 0 cdf_varget,cdfid,'He1_Epoch',He1_Epoch, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for He1_Epoch call/',Result,'/', status ;He2_Epoch Type:CDF_EPOCH DIMS: 0 cdf_varget,cdfid,'He2_Epoch',He2_Epoch, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for He2_Epoch call/',Result,'/', status ;H_Fluence Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'H_Fluence',H_Fluence, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for H_Fluence call/',Result,'/', status ;O_Fluence Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'O_Fluence',O_Fluence, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for O_Fluence call/',Result,'/', status ;He1_Fluence Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'He1_Fluence',He1_Fluence, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for He1_Fluence call/',Result,'/', status ;He2_Fluence Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'He2_Fluence',He2_Fluence, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for He2_Fluence call/',Result,'/', status ;H_sigma Type:CDF_BYTE DIMS: 0 cdf_varget,cdfid,'H_sigma',H_sigma, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for H_sigma call/',Result,'/', status ;O_sigma Type:CDF_BYTE DIMS: 0 cdf_varget,cdfid,'O_sigma',O_sigma, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for O_sigma call/',Result,'/', status ;He1_sigma Type:CDF_BYTE DIMS: 0 cdf_varget,cdfid,'He1_sigma',He1_sigma, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for He1_sigma call/',Result,'/', status ;He2_sigma Type:CDF_BYTE DIMS: 0 cdf_varget,cdfid,'He2_sigma',He2_sigma, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for He2_sigma call/',Result,'/', status ;H_Eflux Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'H_Eflux',H_Eflux, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for H_Eflux call/',Result,'/', status ;O_Eflux Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'O_Eflux',O_Eflux, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for O_Eflux call/',Result,'/', status ;He1_Eflux Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'He1_Eflux',He1_Eflux, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for He1_Eflux call/',Result,'/', status ;He2_Eflux Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'He2_Eflux',He2_Eflux, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for He2_Eflux call/',Result,'/', status ;H_efluxsigma Type:CDF_BYTE DIMS: 0 cdf_varget,cdfid,'H_efluxsigma',H_efluxsigma, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for H_efluxsigma call/',Result,'/', status ;O_efluxsigma Type:CDF_BYTE DIMS: 0 cdf_varget,cdfid,'O_efluxsigma',O_efluxsigma, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for H_efluxsigma call/',Result,'/', status ;He1_efluxsigma Type:CDF_BYTE DIMS: 0 cdf_varget,cdfid,'He1_efluxsigma',He1_efluxsigma, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for He1_efluxsigma call/',Result,'/', status ;He2_efluxsigma Type:CDF_BYTE DIMS: 0 cdf_varget,cdfid,'He2_efluxsigma',He2_efluxsigma, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for H_efluxsigma call/',Result,'/', status ;H_spins Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'H_spins',H_spins, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for H_spins call/',Result,'/', status ;O_spins Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'O_spins',O_spins, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for O_spins call/',Result,'/', status ;He1_spins Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'He1_spins',He1_spins, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for He1_spins call/',Result,'/', status ;He2_spins Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'He2_spins',He2_spins, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for He2_spins call/',Result,'/', status ;H_quality Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'H_quality',H_quality, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for H_quality call/',Result,'/', status ;O_quality Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'O_quality',O_quality, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for O_quality call/',Result,'/', status ;He1_quality Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'He1_quality',He1_quality, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for He1_quality call/',Result,'/', status ;He2_quality Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'He2_quality',He2_quality, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for He2_quality call/',Result,'/', status ; plot one H_flux e_flux quality spinss flux sigma eflux sigma ; also for o, he1, and he2 !p.multi=[0,2,3] start_year=0 start_month=0 start_day=0 stop_year=0 stop_month=0 stop_day=0 cdf_epoch, H_epoch(0), start_year, start_month, start_day,/breakdown cdf_epoch, H_epoch(N_elements(H_fluence)-1),stop_year,stop_month, stop_day,hour,min,sec,/breakdown cdf_epoch, begin_of_day,start_year,start_month,start_day,/compute print, 'FILE: ' ,infile print, 'START: ',start_year, start_month, start_day print, 'STOP: ',stop_year, stop_month, stop_day, hour, min, sec plot,(h_epoch-begin_of_day)/3600000.,H_fluence,/ylog, $ title='Upward H+ Fluence', xrange=[0,24], xstyle=1, yrange=[1e2, 1e10],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) plot,(h_epoch-begin_of_day)/3600000.,H_Eflux,/ylog, $ title='Upward H+ E_Flux', xrange=[0,24], xstyle=1, yrange=[1e2, 1e10],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;H_sigma Type:CDF_BYTE DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,H_sigma, $ title='H+ Fluence Sigma (%)', xrange=[0,24], xstyle=1, yrange=[0,150],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;H_efluxsigma Type:CDF_BYTE DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,H_efluxsigma, $ title='H+ E_Flux Sigma (%)', xrange=[0,24], xstyle=1, yrange=[0,150],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) plot,(h_epoch-begin_of_day)/3600000.,H_quality, $ title='H+ quality', xrange=[0,24], xstyle=1, yrange=[-5,30],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) plot,(h_epoch-begin_of_day)/3600000.,H_spins, $ title='H+ Spins', xrange=[0,24], xstyle=1, yrange=[0,8],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) dum='' print,' Hit return to continue' read,dum ; O+ plot,(o_epoch-begin_of_day)/3600000.,o_fluence,/ylog, $ title='Upward O+ Fluence', xrange=[0,24], xstyle=1, yrange=[1e2, 1e10],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) plot,(o_epoch-begin_of_day)/3600000.,o_Eflux,/ylog, $ title='Upward O+ E_Flux', xrange=[0,24], xstyle=1, yrange=[1e2, 1e10],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;H_sigma Type:CDF_BYTE DIMS: 0 plot,(o_epoch-begin_of_day)/3600000.,o_sigma, $ title='O+ Fluence Sigma (%)', xrange=[0,24], xstyle=1, yrange=[0,150],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;H_efluxsigma Type:CDF_BYTE DIMS: 0 plot,(o_epoch-begin_of_day)/3600000.,o_efluxsigma, $ title='O+ E_Flux Sigma (%)', xrange=[0,24], xstyle=1, yrange=[0,150],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) plot,(o_epoch-begin_of_day)/3600000.,o_quality, $ title='O+ quality', xrange=[0,24], xstyle=1, yrange=[-5,30],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) plot,(o_epoch-begin_of_day)/3600000.,o_spins, $ title='O+ Spins', xrange=[0,24], xstyle=1, yrange=[0,8],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) dum='' print,' Hit return to continue' read,dum ; He+ plot,(he1_epoch-begin_of_day)/3600000.,He1_fluence,/ylog, $ title='Upward He+ Fluence', xrange=[0,24], xstyle=1, yrange=[1e2, 1e10],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) plot,(he1_epoch-begin_of_day)/3600000.,He1_Eflux,/ylog, $ title='Upward He+ E_Flux', xrange=[0,24], xstyle=1, yrange=[1e2, 1e10],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;H_sigma Type:CDF_BYTE DIMS: 0 plot,(he1_epoch-begin_of_day)/3600000.,He1_sigma, $ title='He+ Fluence Sigma (%)', xrange=[0,24], xstyle=1, yrange=[0,150],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;H_efluxsigma Type:CDF_BYTE DIMS: 0 plot,(he1_epoch-begin_of_day)/3600000.,He1_efluxsigma, $ title='He+ E_Flux Sigma (%)', xrange=[0,24], xstyle=1, yrange=[0,150],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) plot,(he1_epoch-begin_of_day)/3600000.,He1_quality, $ title='He+ quality', xrange=[0,24], xstyle=1, yrange=[-5,30],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) plot,(he1_epoch-begin_of_day)/3600000.,He1_spins, $ title='He+ Spins', xrange=[0,24], xstyle=1, yrange=[0,8],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) dum='' print,' Hit return to continue' read,dum ; He++ plot,(he2_epoch-begin_of_day)/3600000.,He2_fluence,/ylog, $ title='Upward He++ Fluence', xrange=[0,24], xstyle=1, yrange=[1e2, 1e10],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) plot,(he2_epoch-begin_of_day)/3600000.,He2_Eflux,/ylog, $ title='Upward He++ E_Flux', xrange=[0,24], xstyle=1, yrange=[1e2, 1e10],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;H_sigma Type:CDF_BYTE DIMS: 0 plot,(he2_epoch-begin_of_day)/3600000.,He2_sigma, $ title='He++ Fluence Sigma (%)', xrange=[0,24], xstyle=1, yrange=[0,150],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;H_efluxsigma Type:CDF_BYTE DIMS: 0 plot,(he2_epoch-begin_of_day)/3600000.,He2_efluxsigma, $ title='He++ E_Flux Sigma (%)', xrange=[0,24], xstyle=1, yrange=[0,150],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) plot,(he2_epoch-begin_of_day)/3600000.,He2_quality, $ title='He++ quality', xrange=[0,24], xstyle=1, yrange=[-5,30],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) plot,(he2_epoch-begin_of_day)/3600000.,He2_spins, $ title='He++ Spins', xrange=[0,24], xstyle=1, yrange=[0,8],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) dum='' print,' Hit return to continue' read,dum ;for i=2000,3000 do print,i, (he2_epoch(i)-begin_of_day)/3600000.,he2_fluence(i),he2_quality(i),he2_spins(i),HE2_sigma(i) ;SZA Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'SZA',SZA, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for SZA call/',Result,'/', status ; empty for 3/17/96 ;Conjugate_SZA Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'Conjugate_SZA',Conjugate_SZA, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for Conjugate_SZA call/',Result,'/', status ;SZA_status Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'SZA_status',SZA_status, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for SZA_status call/',Result,'/', status ;Conjugate_SZA_status Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'Conjugate_SZA_status',Conjugate_SZA_status, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for Conjugate_SZA_status call/',Result,'/', status ;R Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'R',R, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for R call/',Result,'/', status ;L Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'L',L, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for L call/',Result,'/', status ;ILA Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'ILA',ILA, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for ILA call/',Result,'/', status ;MLT Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'MLT',MLT, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for MLT call/',Result,'/', status ;Orbit_status Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'Orbit_status',Orbit_status, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for Orbit_status call/',Result,'/', status ; plot 2 sza consza status status r l orbit_status invl !p.multi=[0,2,4] ;SZA plot,(h_epoch-begin_of_day)/3600000.,SZA, $ title='SZA', xrange=[0,24], xstyle=1, yrange=[0,180],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;Conjugate_SZA plot,(h_epoch-begin_of_day)/3600000.,Conjugate_SZA, $ title='Conjugate_SZA', xrange=[0,24], xstyle=1, yrange=[0,180],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;SZA_status plot,(h_epoch-begin_of_day)/3600000.,SZA_status, $ title='SZA_status', xrange=[0,24], xstyle=1, yrange=[0,110],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;Conjugate_SZA_status plot,(h_epoch-begin_of_day)/3600000.,Conjugate_SZA_status, $ title='Conjugate_SZA_status', xrange=[0,24], xstyle=1, yrange=[0,110],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;R plot,(h_epoch-begin_of_day)/3600000.,R, $ title='R (Re)', xrange=[0,24], xstyle=1, yrange=[1,10],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;L plot,(h_epoch-begin_of_day)/3600000.,L, $ title='L (Re)', xrange=[0,24], xstyle=1, yrange=[2,200],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;Orbit_status Type:CDF_UINT1 DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,Orbit_status, $ title='Orbit_status', xrange=[0,24], xstyle=1, yrange=[-10,110],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;ILA plot,(h_epoch-begin_of_day)/3600000.,ILA, $ title='ILA (degrees)', xrange=[0,24], xstyle=1, yrange=[30,90],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) dum='' print,'Hit return to continue' read,dum ; read plot these: ;Fec Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'Fec',Fec, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for Fec call/',Result,'/', status ;Bcr Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'Bcr',Bcr, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for Bcr call/',Result,'/', status ;Energy_range Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'Energy_range',Energy_range, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for Energy_range call/',Result,'/', status ;Intervals Type:CDF_UINT1 DIMS: 0;MLT cdf_varget,cdfid,'Intervals',Intervals, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for Intervals call/',Result,'/', status ;LocT Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'LocT',LocT, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for LocT call/',Result,'/', status ;Mlat Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'Mlat',Mlat, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for Mlat call/',Result,'/', status ;MagField Type:CDF_REAL4 DIMS: 3 cdf_varget,cdfid,'MagField',MagField, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for MagField call/',Result,'/', status ;Mag_status Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'Mag_status',Mag_status, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for Mag_status call/',Result,'/', status !p.multi=[0,2,5] ; plot these: ;Fec Type:CDF_REAL4 DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,Fec,/ylog, $ title='Fec (counts/spin)', xrange=[0,24], xstyle=1, yrange=[10,1e6],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;Bcr Type:CDF_REAL4 DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,Bcr,/ylog, $ title='Bcr (counts/spin)', xrange=[0,24], xstyle=1, yrange=[1,1e5],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;Energy_range Type:CDF_UINT1 DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,Energy_range, $ title='Energy_range', xrange=[0,24], xstyle=1, yrange=[-1,4],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;Intervals Type:CDF_UINT1 DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,Intervals, $ title='Intervals', xrange=[0,24], xstyle=1, yrange=[-1,9],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;MLT plot,(h_epoch-begin_of_day)/3600000.,MLT, $ title='MLT (hrs)', xrange=[0,24], xstyle=1, yrange=[0,24],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;LocT Type:CDF_REAL4 DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,LocT, $ title='LocT (hrs)', xrange=[0,24], xstyle=1, yrange=[0,24],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;Mlat Type:CDF_REAL4 DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,Mlat, $ title='Mlat (degrees)', xrange=[0,24], xstyle=1, yrange=[-90.,90.],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;MagField Type:CDF_REAL4 DIMS: 3 plot,(h_epoch-begin_of_day)/3600000.,MagField(0,*), $ title='MagField (GSM x,y,z)', xrange=[0,24], xstyle=1, yrange=[-10000,10000.],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) oplot,(h_epoch-begin_of_day)/3600000.,MagField(1,*) oplot,(h_epoch-begin_of_day)/3600000.,MagField(2,*) ;Mag_status Type:CDF_UINT1 DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,Mag_status, $ title='MagField (GSM x,y,z)', xrange=[0,24], xstyle=1, yrange=[-1 ,10.],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) dum='' print,' Hit return to continue' read,dum ;read these ;SC_Potential Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'Mag_status',SC_Potential, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for SC_Potential call/',Result,'/', status ;SCPot_status Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'Mag_status',SCPot_status, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for SCPot_status call/',Result,'/', status ;Imf_Gsm Type:CDF_REAL4 DIMS: 3 cdf_varget,cdfid,'Imf_Gsm',Imf_Gsm, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for Imf_Gsm call/',Result,'/', status ;IMF_status Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'IMF_status',IMF_status, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for IMF_status call/',Result,'/', status ;SW_pressure Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'SW_pressure',SW_pressure, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for SW_pressure call/',Result,'/', status ;SW_density Type:CDF_REAL4 DIMS: 0 cdf_varget,cdfid,'SW_density',SW_density, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for SW_density call/',Result,'/', status ;SW_velocity Type:CDF_REAL4 DIMS: 3 cdf_varget,cdfid,'SW_velocity',SW_velocity, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for SW_velocity call/',Result,'/', status ;SW_status Type:CDF_UINT1 DIMS: 0 cdf_varget,cdfid,'SW_status',SW_status, rec_count=max_records,/zvariable Result=CDF_ERROR(status) if result ne 'CDF_OK: Function completed successfully.' then $ print, 'Status after cdf_varget for SW_status call/',Result,'/', status ; plot these !p.multi=[0,2,4] ;SC_Potential Type:CDF_REAL4 DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,SC_Potential, $ title='SC_Potential from EFI (eV)', xrange=[0,24], xstyle=1, yrange=[-50,50.],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;SCPot_status Type:CDF_UINT1 DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,SCPot_status, $ title='SCPot_status ', xrange=[0,24], xstyle=1, yrange=[-10,130.],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;Imf_Gsm Type:CDF_REAL4 DIMS: 3 plot,(h_epoch-begin_of_day)/3600000.,Imf_Gsm(0,*), $ title='IMF Field (GSM x,y,z)', xrange=[0,24], xstyle=1, yrange=[-10000,10000.],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) oplot,(h_epoch-begin_of_day)/3600000.,Imf_Gsm(1,*) oplot,(h_epoch-begin_of_day)/3600000.,Imf_Gsm(2,*) ;IMF_status Type:CDF_UINT1 DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,IMF_status, $ title='IMF_status', xrange=[0,24], xstyle=1, yrange=[-10,130.],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;SW_pressure Type:CDF_REAL4 DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,SW_pressure, $ title='SW_pressure (nP)', xrange=[0,24], xstyle=1, yrange=[0,100.],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;SW_density Type:CDF_REAL4 DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,SW_density, $ title='SW_density (cm^-3)', xrange=[0,24], xstyle=1, yrange=[0,100.],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;SW_status Type:CDF_UINT1 DIMS: 0 plot,(h_epoch-begin_of_day)/3600000.,SW_status, $ title='SW_status', xrange=[0,24], xstyle=1, yrange=[-10,130.],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) ;SW_velocity Type:CDF_REAL4 DIMS: 3 plot,(h_epoch-begin_of_day)/3600000.,SW_velocity(0,*), $ title='SW_velocity (GSM x,y,z)', xrange=[0,24], xstyle=1, yrange=[-1000,1000.],ystyle=1, $ Xtitle='Year/Month/Day ='+string(start_year)+string(start_month)+string(Start_day) oplot,(h_epoch-begin_of_day)/3600000.,SW_velocity(1,*) oplot,(h_epoch-begin_of_day)/3600000.,SW_velocity(2,*) dum='' print,' Hit return to continue' read,dum CDF_CLOSE,cdfid print,'CDF FILE CLOSED' !p.multi=[0,1,1] endwhile ; end of the big loop stop end