1 / 4

Time Series Analysis of Particles and Fields data: SST corrections

Time Series Analysis of Particles and Fields data: SST corrections. In the tail or near Earth or during high speed flows, SST contribution to can be significant: Velocity, pressure corrections from SST Materials in: http://www.igpp.ucla.edu/public/vassilis/ESS265/20080602

dolph
Télécharger la présentation

Time Series Analysis of Particles and Fields data: SST corrections

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Time Series Analysis of Particles and Fields data: SST corrections • In the tail or near Earth or during high speed flows,SST contribution to can be significant: • Velocity, pressure corrections from SST Materials in: http://www.igpp.ucla.edu/public/vassilis/ESS265/20080602 class_notes_V_P_corr_SST.ppt joint_esa_sst_mom_v02.pro fix_dat.pro Time Series Analysis1

  2. Density Correction • Interpolate densities • Add • date='2008-03-01' • startdate = '2008-03-01/00:00' • timespan,startdate, 4.0, /hour • Trange=['08-03-01/00:00','08-03-01/04:00'] • Tzoom=['08-03-01/01:40','08-03-01/02:40'] • ;... select exact time interval to calculate join ESA/SST moments • tbeg = time_double(date+'/00:00') • tend = time_double(date+'/04:00') • ;select a probe • sc='b' • thm_load_state,probe=sc,coord='gsm',/get_support • thm_load_fit, level=1, probe=sc,datatype=['efs', 'fgs'],/verbose • thm_cotrans,strjoin('th'+sc+'_fgs'),out_suf='_gsm', in_c='dsl', out_c='gsm' • ; • ; SST now • thm_load_sst,probe=sc,lev=1 • thm_part_moments, probe = sc, instr= ['ps?f'], $ • moments = ['density', 'velocity', 't3'], $ • mag_suffix='_peir_magt3', $ • scpot_suffix='_peir_sc_pot';,/median • ; work in gsm • thm_cotrans,'th'+sc+'_ps?f_velocity',in_coord='dsl',out_coord='gsm',out_suffix='_gsm' • ; • ; ESA now • thm_load_esa,probe=sc • ; Interpolate densities • tinterpol_mxn,'th'+sc+'_peer_density','th'+sc+'_peir_density',/overwrite,/nan_extrapolate • tinterpol_mxn,'th'+sc+'_ps?f_density','th'+sc+'_peir_density',/overwrite,/nan_extrapolate • … • ; ...total ion density • totNi = sst_i_n.y + esa_i_n.y Ni Ne Time Series Analysis2

  3. Velocity Correction • Interpolate densities • Add flux • ; • ; • ; ...sst Flux • sstFi = sst_i_v.y*0. • sstFi[*,0] = sst_i_n.y*sst_i_v.y[*,0] • sstFi[*,1] = sst_i_n.y*sst_i_v.y[*,1] • sstFi[*,2] = sst_i_n.y*sst_i_v.y[*,2] • ; ...esa Flux • esaFi = esa_i_v.y*0. • esaFi[*,0] = esa_i_n.y*esa_i_v.y[*,0] • esaFi[*,1] = esa_i_n.y*esa_i_v.y[*,1] • esaFi[*,2] = esa_i_n.y*esa_i_v.y[*,2] • ; ...total ion density • totNi = sst_i_n.y + esa_i_n.y • store_data, 'th'+sc+'_Ni',$data={x:esa_i_n.x, y:totNi} • options, 'th'+sc+'_Ni', 'ytitle', $'Ni !C!C1/cm!U3' • ylim, 'th'+sc+'_Ni', 0.01, 1., 1 • ; ...total ion velocity (GSM) • totVi = esa_i_v.y*0. • totVi[*,0] = (sstFi[*,0]+esaFi[*,0])/totNi • totVi[*,1] = (sstFi[*,1]+esaFi[*,1])/totNi • totVi[*,2] = (sstFi[*,2]+esaFi[*,2])/totNi Time Series Analysis3

  4. Pressure Correction • Remove SST noise • Interpolate pressures • Then add • ; • ; SST now • ; SST now • thm_load_sst,probe=sc,lev=1 • thm_part_moments, probe = sc, instr= ['ps?f'], $ • moments = ['density', 'velocity', 't3'], $ • mag_suffix='_peir_magt3', $ • scpot_suffix='_peir_sc_pot';,/median • ; …interpolate • ; … add • ; ...pressure • ; ...SST: perpendicular temperature only • sst_Tperp = .5*(sst_i_t3.y[*,0]+sst_i_t3.y[*,1]) • sst_i_p_nPa = 0.16*.001*sst_i_n.y * sst_Tperp • ; perp. pressure in nPa • store_data, 'th'+sc+'_psif_p_perp_nPa', $ • data={x:sst_i_n.x, y:sst_i_p_nPa} • options, 'th'+sc+'_psif_p_perp_nPa', $ • 'ytitle', 'sst Pi !C!CnPa' • ; ...ESA: scalar temperature • esa_Ti = total(esa_i_T.y,2)/3. • store_data,'Ti_th'+sc+'_peir', $ • data={x:esa_i_n.x, y:esa_Ti} • ; ...ESA ion pressure: • esa_i_p_nPa = 0.16 *.001 * esa_i_n.y*esa_Ti • ; scalar pressure in nPa • store_data, 'th'+sc+'_peir_p_nPa', $ • data={x:esa_i_n.x, y:esa_i_p_nPa} • options, 'th'+sc+'_peir_p_nPa', $ • 'ytitle', 'esa Pi !C!CnPa' • ; ...Total ion pressure • totPi = sst_i_p_nPa + esa_i_p_nPa • store_data, 'th'+sc+'_i_p_nPa', $ • data={x:esa_i_n.x, y:totPi} • options, 'th'+sc+'_i_p_nPa', 'ytitle', 'Pi !C!CnPa' Time Series Analysis4

More Related