GSFIT Data Format

From EOVSA Wiki
Revision as of 11:34, 11 May 2019 by Gnita (talk | contribs)
Jump to navigation Jump to search

GSFIT Log File Format

NOTE: Although such GSFIT/GSFITCP may appear to be traditional IDL *.sav, attempting to restore them via IDL command line using the native RESTORE procedure would only recover the first record logged in the file. To recover the full information stored in such files one may use one of the two support routines that are part of the GSFIT package

IDL> fitlist=gsfit_log2list('gsfitcp.log')
% Compiled module: MULTI_RESTORE.
% RESTORE: Portable (XDR) SAVE/RESTORE file.
% RESTORE: Save file written by gnita@baozi.hpcnet.campus.njit.edu, Fri Mar  1 18:15:09 2019.
% RESTORE: IDL version 8.5 (linux, x86_64).
% RESTORE: Restored variable: HEADER.
% RESTORE: Restored variable: TEMPLATE.
% RESTORE: Restored variable: HINT.
% RESTORE: Restored variable: KEY.
% RESTORE: Restored variable: XDR.
% GSFIT_LOG2LIST: Restoring 409599 additional records
IDL> help,fitlist
FITLIST         STRUCT    = -> <Anonymous> Array[409600]
IDL> help,fitlist,/str
** Structure <21fed180>, 24 tags, length=504, data length=504, refs=2:
   X               LONG                67
   Y               LONG                42
   T               LONG                 0
   FMIN            INT              0
   FMAX            INT             23
   DATA            FLOAT     Array[30]
   ERRDATA         FLOAT     Array[30]
   SPECFIT         FLOAT     Array[30]
   N_NTH           DOUBLE          0.00000000
   ERRN_NTH        DOUBLE          0.00000000
   B               DOUBLE          0.00000000
   ERRB            DOUBLE          0.00000000
   THETA           DOUBLE          0.00000000
   ERRTHETA        DOUBLE          0.00000000
   N_TH            DOUBLE          0.00000000
   ERRN_TH         DOUBLE          0.00000000
   DELTA           DOUBLE          0.00000000
   ERRDELTA        DOUBLE          0.00000000
   E_MAX           DOUBLE          0.00000000
   ERRE_MAX        DOUBLE          0.00000000
   T_E             DOUBLE          0.00000000
   ERRT_E          DOUBLE          0.00000000
   CHISQR          DOUBLE          0.00000000
   RESIDUAL        DOUBLE          0.00000000


GSFIT Parameter Map Structure

IDL> help, maps
** Structure <1e0e3430>, 28 tags, length=189883336, data length=189666533, refs=1:
   FITMAPS         STRUCT    -> <Anonymous> Array[30, 100]       ;[nfreq x ntimes] GSFIT solution maps
   DATAMAPS        STRUCT    -> <Anonymous> Array[30, 100]       ;[nfreq x ntimes] input data maps use to compute the GSFIT solutions
   ERRMAPS         STRUCT    -> <Anonymous> Array[30, 100]       ;[nfreq x ntimes] input error maps used to weight the input data 
   N_NTH           STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of estimated nonthermal electron densities
   ERRN_NTH        STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of 1-sigma fit errors of the estimated nonthermal electron densities
   B               STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of estimated absolute magnetic field strength
   ERRB            STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of 1-sigma fit errors of the estimated 
   THETA           STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of estimated LOS angle of the magnetic field vector 
   ERRTHETA        STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of 1-sigma fit errors of the estimated LOS angle of the magnetic field vector 
   N_TH            STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of estimated thermal electron densities
   ERRN_TH         STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of 1-sigma fit errors of the estimated thermal electron densities
   DELTA           STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of estimated distribution over energy power law slope
   ERRDELTA        STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of 1-sigma fit errors of the estimated distribution over energy power law slope
   E_MAX           STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of estimated maximum energy cutoff of nonthermal electrons
   ERRE_MAX        STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of 1-sigma fit errors of the estimated maximum energy cutoff of nonthermal electrons
   T_E             STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of estimated plasma temperature
   ERRT_E          STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of 1-sigma fit errors of the estimated plasma temperature
   CHISQR          STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of fit CHISQR
   RESIDUAL        STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of fit Residuals 
   PEAKFLUX        STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of estimated peak flux of the fit solution
   ERRPEAKFLUX     STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of 1-sigma fit errors of the estimated peak flux of the fit solution 
   PEAKFREQ        STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of estimated peak frequency of the fit solution
   ERRPEAKFREQ     STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of 1-sigma fit errors of the estimated peak frequency of the fit solution
   WB              STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps estimated magnetic field energy
   ERRWB           STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of 1-sigma fit errors of the estimated magnetic field energy
   WNTH            STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps estimated nonthermal electron energy
   ERRWNTH         STRUCT    -> <Anonymous> Array[100]           ;[ntimes] maps of 1-sigma fit errors of the estimated nonthermal electron energy 
   HEADER          STRUCT    -> <Anonymous> Array[1]             ;structure containing information about the GSFIT solutions from which these maps were created
IDL> maps=gsfit_log2map(file)
% RESTORE: Portable (XDR) SAVE/RESTORE file.
% RESTORE: Save file written by gnita@baozi.hpcnet.campus.njit.edu, Fri Mar  1 18:15:09 2019.
% RESTORE: IDL version 8.5 (linux, x86_64).
% RESTORE: Restored variable: HEADER.
% RESTORE: Restored variable: TEMPLATE.
% RESTORE: Restored variable: HINT.
% RESTORE: Restored variable: KEY.
% RESTORE: Restored variable: XDR.
% GSFIT_LOG2LIST: Restoring 409599 additional record

IDL> help,map
** Structure <1fe1b690>, 28 tags, length=189883336, data length=189666533, refs=1:
   FITMAPS         STRUCT    -> <Anonymous> Array[30, 100]
   DATAMAPS        STRUCT    -> <Anonymous> Array[30, 100]
   ERRMAPS         STRUCT    -> <Anonymous> Array[30, 100]
   N_NTH           STRUCT    -> <Anonymous> Array[100]
   ERRN_NTH        STRUCT    -> <Anonymous> Array[100]
   B               STRUCT    -> <Anonymous> Array[100]
   ERRB            STRUCT    -> <Anonymous> Array[100]
   THETA           STRUCT    -> <Anonymous> Array[100]
   ERRTHETA        STRUCT    -> <Anonymous> Array[100]
   N_TH            STRUCT    -> <Anonymous> Array[100]
   ERRN_TH         STRUCT    -> <Anonymous> Array[100]
   DELTA           STRUCT    -> <Anonymous> Array[100]
   ERRDELTA        STRUCT    -> <Anonymous> Array[100]
   E_MAX           STRUCT    -> <Anonymous> Array[100]
   ERRE_MAX        STRUCT    -> <Anonymous> Array[100]
   T_E             STRUCT    -> <Anonymous> Array[100]
   ERRT_E          STRUCT    -> <Anonymous> Array[100]
   CHISQR          STRUCT    -> <Anonymous> Array[100]
   RESIDUAL        STRUCT    -> <Anonymous> Array[100]
   PEAKFLUX        STRUCT    -> <Anonymous> Array[100]
   ERRPEAKFLUX     STRUCT    -> <Anonymous> Array[100]
   PEAKFREQ        STRUCT    -> <Anonymous> Array[100]
   ERRPEAKFREQ     STRUCT    -> <Anonymous> Array[100]
   WB              STRUCT    -> <Anonymous> Array[100]
   ERRWB           STRUCT    -> <Anonymous> Array[100]
   WNTH            STRUCT    -> <Anonymous> Array[100]
   ERRWNTH         STRUCT    -> <Anonymous> Array[100]
   HEADER          STRUCT    -> <Anonymous> Array[1]