2017: Difference between revisions

From EOVSA Wiki
Jump to navigation Jump to search
 
(5 intermediate revisions by the same user not shown)
Line 1: Line 1:
== List of EOVSA Flares with Spectrogram Data ==
== List of EOVSA Flares with Spectrogram Data ==
=== July ===
=== September ===
{| class="wikitable"
{| class="wikitable"
!  Date    || Time (UT) || GOES Class || Position || Spectrogram || STIX Coverage || AIA Movie || EOVSA Images || Comment
!  Date    || Time (UT) || GOES Class || Position || Spectrogram || STIX Coverage || AIA Movie || EOVSA Images || Comment
|-
|-
| [http://ovsa.njit.edu/browser/?suntoday_date=2017-09-07 2017-09-07] || 14:20 || X1.3 || || [[File:eovsa.spec.flare_id_20190415193100.png|thumb|center|100px|]] [http://ovsa.njit.edu/events/2019/eovsa.spec.flare_id_20190415193100.fits plot data] || NA || [https://www.lmsal.com/hek/her?cmd=view-voevent&ivorn=ivo://helio-informatics.org/FL_SSWLatestEvents_20170907_082206_817 AIA] || [http://ovsa.njit.edu/SynopticImg/eovsamedia/eovsa-browser/2017/09/17/eovsa.lev1_mbd_12s.flare_id_20170917142000.mp4 Quicklook Movie] <br/> [http://ovsa.njit.edu/fits/flares/2017/09/07/20170907142000/ FITS Files]
| [http://ovsa.njit.edu/browser/?suntoday_date=2017-09-07 2017-09-07] || 14:20 || X1.3 || || [[File:eovsa.spec.flare_id_20170907142000.png|thumb|center|100px|]] [http://ovsa.njit.edu/events/2017/eovsa.spec.flare_id_20170907142000.fits plot data] || NA || [https://www.lmsal.com/hek/her?cmd=view-voevent&ivorn=ivo://helio-informatics.org/FL_SSWLatestEvents_20170907_082206_817 AIA] || [http://ovsa.njit.edu/SynopticImg/eovsamedia/eovsa-browser/2017/09/07/eovsa.lev1_mbd_12s.flare_id_20170907142000.mp4 Quicklook Movie] <br/> [http://ovsa.njit.edu/fits/flares/2017/09/07/20170907142000/ FITS Files]
|-
| [http://ovsa.njit.edu/browser/?suntoday_date=2017-09-08 2017-09-08] || 15:10 || M2.9 || || [[File:eovsa.spec.flare_id_20170908151000.png|thumb|center|100px|]] [http://ovsa.njit.edu/events/2019/eovsa.spec.flare_id_20170908151000.fits plot data] || NA || [https://www.lmsal.com/hek/her?cmd=view-voevent&ivorn=ivo://helio-informatics.org/FL_SSWLatestEvents_20170908_094755_252 AIA] || [http://ovsa.njit.edu/SynopticImg/eovsamedia/eovsa-browser/2017/09/08/eovsa.lev1_mbd_12s.flare_id_20170908151000.mp4 Quicklook Movie] <br/> [http://ovsa.njit.edu/fits/flares/2017/09/08/20170908151000/ FITS Files]
|-
|-
|}
|}

Latest revision as of 22:30, 13 August 2024

List of EOVSA Flares with Spectrogram Data

September

Date Time (UT) GOES Class Position Spectrogram STIX Coverage AIA Movie EOVSA Images Comment
2017-09-07 14:20 X1.3
Eovsa.spec.flare id 20170907142000.png
plot data
NA AIA Quicklook Movie
FITS Files
2017-09-08 15:10 M2.9
Eovsa.spec.flare id 20170908151000.png
plot data
NA AIA Quicklook Movie
FITS Files

Python code to read plotdata file

from __future__ import print_function
def rd_datfile(file):
    ''' Read EOVSA binary spectrogram file and return a dictionary with times 
        in Julian Date, frequencies in GHz, and cross-power data in sfu.
        
        Return Keys:
          'time'     Numpy array of nt times in JD format
          'fghz'     Numpy array of nf frequencies in GHz
          'data'     Numpy array of size [nf, nt] containing cross-power data
          
        Returns empty dictionary ({}) if file size is not compatible with inferred dimensions
    '''
    import struct
    import numpy as np
    def dims(file):
        # Determine time and frequency dimensions (assumes the file has fewer than 10000 times)
        f = open(file,'rb')
        tmp = f.read(83608)  # max 10000 times and 451 frequencies
        f.close()
        nbytes = len(tmp)
        tdat = np.array(struct.unpack(str(int(nbytes/8))+'d',tmp[:nbytes]))
        nt = np.where(tdat < 2400000.)[0]
        nf = np.where(np.logical_or(tdat[nt[0]:] > 18, tdat[nt[0]:] < 1))[0]
        return nt[0], nf[0]
    nt, nf = dims(file)
    f = open(file,'rb')
    tmp = f.read(nt*8)
    times = struct.unpack(str(nt)+'d',tmp)
    tmp = f.read(nf*8)
    fghz = struct.unpack(str(nf)+'d',tmp)
    tmp = f.read()
    f.close()
    if len(tmp) != nf*nt*4:
        print('File size is incorrect for nt=',nt,'and nf=',nf)
        return {}
    data = np.array(struct.unpack(str(nt*nf)+'f',tmp)).reshape(nf,nt)
    return {'time':times, 'fghz':fghz, 'data':data}

IDL code to read plotdata file

function rd_datfile,file
    ; Read EOVSA binary spectrogram file and return a structure with times 
    ; in Julian Date, frequencies in GHz, and cross-power data in sfu.
    ;    
    ; Return tags:
    ;      'time'     Array of nt times in JD format
    ;      'fghz'     Array of nf frequencies in GHz
    ;      'data'     Array of size [nf, nt] containing cross-power data
    ;
    ;    Returns empty dictionary ({}) if file size is not compatible with inferred dimensions
    openr,/get_lun,lun,file
    tmp = dblarr(10451)
    readu,lun,tmp
    free_lun,lun
    nt = (where(tmp lt 2400000.))[0]
    nf = (where(tmp[nt[0]:*] gt 18 or tmp[nt[0]:*] lt 1))[0]
    times = dblarr(nt)
    fghz = dblarr(nf)
    data = fltarr(nt, nf)
    openr,/get_lun,lun,file
    readu,lun,times
    readu,lun,fghz
    readu,lun,data
    free_lun,lun
    data = create_struct('time',times,'fghz',fghz,'data',transpose(data))
    return, data
    end