#29 psrsplit-like application for search mode data
Hi! I am not sure if this is a bug or a feature request. I am trying to use psrsplit to split search mode fits file by frequency. I tried using NenuFAR and Parkes (DFB4 data), both cases it complains about illegal HDU number (301) and did not write out any files. Using fold mode fits seem to be OK. The verbose screen output is attached at the end. Is there any way around this issue?
Cheers,
Cherry
(Sing_py3) cng@nancep7:$ psrsplit -V -c 256 PDFB4_Search_J0437-4715.sf
Pulsar::Archive::load with 12 Agents
Pulsar::Archive::load testing Timer
Timer::load convert endian
Timer::load not a timer header
Timer::is_timer invalid MJD=544109933
Pulsar::TimerArchive::Agent::advocate fail Timer::fload
Pulsar::Archive::load testing Baseband
Timer::load convert endian
Timer::load not a timer header
Timer::is_timer invalid MJD=544109933
Pulsar::Archive::load testing EPN
Pulsar::Archive::load testing ASP
IntegrationManager null constructor
Pulsar::Archive ctor this=0x7fff132403d0 instances=1
Pulsar::Archive::init
Pulsar::Archive::add_extension appending Receiver
Pulsar::Archive::add_extension appending Backend
Pulsar::BasicArchive dtor this=0x7fff132403d0
Pulsar::Archive dtor this=0x7fff132403d0 instances=0
Pulsar::IntegrationManager dtor this=0x7fff132403d0
ASP load failed due to:
Error::stack
Pulsar::ASPArchive::load_header
Error::InvalidState
Error::message
Unrecognized HDRVER=5.3 in file PDFB4_Search_J0437-4715.sf
Pulsar::Archive::load testing PSRFITS
Pulsar::FITSArchive::Agent::advocate test PDFB4_Search_J0437-4715.sf
FITSAgent::advocate test reading MJD
Pulsar::Archive::load using PSRFITS
IntegrationManager null constructor
Pulsar::Archive ctor this=0x561437e11d60 instances=1
Pulsar::Archive::init
FITSArchive default construct
Pulsar::Archive::add_extension appending FITSHdrExtension
FITSArchive::load_header fits_open_file (PDFB4_Search_J0437-4715.sf)
Pulsar::Archive::get<ext> name=FITSHdrExtension
Pulsar::Archive::get<ext> failed to find extension type=PN6Pulsar12ObsExtensionE
Pulsar::Archive::add_extension appending ObsExtension
Pulsar::Archive::get<ext> name=FITSHdrExtension
FITSArchive::load_header reading FITS header version
psrfits_read_key: HDRVER=5.3
Got: Version 5.3
FITSArchive::load_header reading file creation date
psrfits_read_key: DATE=2014-03-14T05:42:39
FITSArchive::load_header reading observer name
psrfits_read_key: OBSERVER=RNM
Got observer: RNM
FITSArchive::load_header reading project ID
psrfits_read_key: PROJID=P456
Got PID: P456
FITSArchive::load_header reading telescope name
psrfits_read_key: TELESCOP=PARKES
FITSArchive::load_header TELESCOP=PARKES
Pulsar::Archive::get<ext> name=FITSHdrExtension
Pulsar::Archive::get<ext> name=ObsExtension
Pulsar::Archive::get<ext> failed to find extension type=PN6Pulsar9TelescopeE
Pulsar::Archive::add_extension appending Telescope
FITSArchive::load_header calling Telescopes::set_telescope_info
psrfits_read_key: RA=04:37:15.800
psrfits_read_key: DEC=-47:15:08.600
FITSArchive::load_ITRFExtension
Pulsar::Archive::add_extension appending ITRFExtension
FITSArchive::load_Receiver
FITSArchive::load_Receiver reading receiver
psrfits_read_key: FRONTEND=MULTI
psrfits_read_key: NRCVR=1
FITSArchive::load_Receiver reading feed polarization
psrfits_read_key: FD_POLN=LIN
FITSArchive::load_Receiver reading XPOL_ANG
XPOL_ANG: keyword not found in header (202)
psrfits_read_key: using default=0
FITSArchive::load_Receiver reading YPOL_ANG
YPOL_ANG: keyword not found in header (202)
psrfits_read_key: using default=90
FITSArchive::load_Receiver reading CAL_ANG
CAL_ANG: keyword not found in header (202)
psrfits_read_key: using default=45
FITSArchive::load_Receiver reading FD_SANG
FITSArchive::load_Receiver reading FD_HAND
FITSArchive::load_Receiver reading FD_XYPH
FITSArchive::load_Receiver reading feed track mode
psrfits_read_key: FD_MODE=FA
FITSArchive::load_Receiver reading requested feed angle
FITSArchive::load_Receiver reading attenuator settings
FITSArchive::load_Receiver WARNING: fits_read_key ATTEN_A,B: keyword not found in header (202)
Pulsar::Archive::add_extension appending Receiver
FITSArchive::load_WidebandCorrelator
psrfits_read_key: BACKEND=PDFB4
FITSArchive::load_WidebandCorrelator BACKEND=PDFB4
psrfits_read_key: BECONFIG=srch_256_512
FITSArchive::load_WidebandCorrelator BECONFIG=srch_256_512
psrfits_read_key: BE_PHASE=-1
FITSArchive::load_WidebandCorrelator BE_PHASE=-1
psrfits_read_key: BE_DCC=0
FITSArchive::load_WidebandCorrelator BE_DCC=0
FITSArchive::load_WidebandCorrelator reading cycle time
psrfits_read_key: TCYCLE=10
FITSArchive::load_WidebandCorrelator reading NRCVR
psrfits_read_key: NRCVR=1
FITSArchive::load_WidebandCorrelator reading BE_DELAY
psrfits_read_key: BE_DELAY=0
Pulsar::Archive::add_extension appending WidebandCorrelator
Pulsar::Archive::get<ext> name=FITSHdrExtension
Pulsar::Archive::get<ext> name=ObsExtension
Pulsar::Archive::get<ext> name=Telescope
Pulsar::Archive::get<ext> name=ITRFExtension
Pulsar::Archive::get<ext> name=Receiver
Pulsar::Archive::get<ext> name=WidebandCorrelator
FITSArchive::load_header reading OBS_MODE
FITSArchive::load_header OBS_MODE='SEARCH'
FITSArchive::load_header search mode file
psrfits_read_key: OBSFREQ=1369
psrfits_read_key: OBSBW=256
psrfits_read_key: OBSNCHAN=512
psrfits_read_key: CHAN_DM=0
FITSArchive::load_header reading source name
FITSArchive::load_header reading coordinates
psrfits_read_key: COORD_MD=J2000
Got coordinate type: J2000
FITSArchive::load_header reading equinox
psrfits_read_key: EQUINOX=2000
psrfits_read_key: RA=04:37:15.800
psrfits_read_key: STT_CRD1=04:37:15.800
psrfits_read_key: DEC=-47:15:08.600
psrfits_read_key: STT_CRD2=-47:15:08.600
psrfits_read_key: STP_CRD1=04:37:15.800
psrfits_read_key: STP_CRD2=-47:15:08.600
FITSArchive::load_header reading track mode
psrfits_read_key: TRK_MODE=TRACK
FITSArchive::load_header reading observation date
psrfits_read_key: DATE-OBS=2014-03-14T05:42:39
FITSArchive::load_header DATE-0BS parsed into
date='2014-03-14'
time='05:42:39'
FITSArchive::load_header reading BPA
psrfits_read_key: BPA=0.
FITSArchive::load_header reading BMAJ
psrfits_read_key: BMAJ=0.248356464572681
FITSArchive::load_header reading BMIN
psrfits_read_key: BMIN=0.248356464572681
FITSArchive::load_header reading MJDs
psrfits_read_key: STT_IMJD=56730
psrfits_read_key: STT_SMJD=20569
psrfits_read_key: STT_OFFS=1
FITSArchive::load_header MJD=56730:054249.999999995343387
FITSArchive::load_header reading start LST
psrfits_read_key: STT_LST=10937.9
FITSArchive::load_header reading IBEAM
psrfits_read_key: IBEAM=
FITSArchive::load_header IBEAM=''
FITSArchive::load_header reading PNT_ID
psrfits_read_key: PNT_ID=
FITSArchive::load_header PNT_ID=''
FITSArchive::load_header finished with primary HDU
FITSArchive::load_ProcHistory entered
Pulsar::Archive::get<ext> name=FITSHdrExtension
Pulsar::Archive::get<ext> name=ObsExtension
Pulsar::Archive::get<ext> name=Telescope
Pulsar::Archive::get<ext> name=ITRFExtension
Pulsar::Archive::get<ext> name=Receiver
Pulsar::Archive::get<ext> name=FITSHdrExtension
Pulsar::Archive::get<ext> name=ObsExtension
Pulsar::Archive::get<ext> name=Telescope
Pulsar::Archive::get<ext> name=ITRFExtension
Pulsar::Archive::get<ext> name=Receiver
Pulsar::Archive::get<ext> name=WidebandCorrelator
load ProcHistory::row entered
psrfits_read_col: DATE_PRO=2014-03-14T05:42:39
load ProcHistory::row DATE_PRO=2014-03-14T05:42:39
psrfits_read_col: PROC_CMD=PDFB4:2014-03-03_15:56
load ProcHistory::row PROC_CMD=PDFB4:2014-03-03_15:56
psrfits_read_col: NSUB=17
psrfits_read_col: NBIN_PRD=0
psrfits_read_col: TBIN=9.6e-05
psrfits_read_col: DM=0
psrfits_read_col: RM=0
psrfits_read_col: PR_CORR=0
psrfits_read_col: FD_CORR=0
psrfits_read_col: BE_CORR=0
psrfits_read_col: RM_CORR=0
psrfits_read_col: DEDISP=0
psrfits_read_col: DDS_MTHD=NONE
psrfits_read_col: SC_MTHD=
psrfits_read_col: CAL_MTHD=NONE
psrfits_read_col: CAL_FILE=NONE
psrfits_read_col: RFI_MTHD=
psrfits_read_col: RM_MODEL=NONE
psrfits_read_col: AUX_RM_C=0
psrfits_read_col: DM_MODEL=
psrfits_read_col: AUX_DM_C=0
psrfits_read_col: SCALE=FluxDen
load ProcHistory::row exiting
FITSArchive::load_ProcHistory not calibrated
Pulsar::Archive::add_extension appending ProcHistory
Pulsar::Archive::get<ext> name=FITSHdrExtension
Pulsar::Archive::get<ext> name=ObsExtension
Pulsar::Archive::get<ext> name=Telescope
Pulsar::Archive::get<ext> name=ITRFExtension
Pulsar::Archive::get<ext> name=Receiver
Pulsar::Archive::get<ext> name=WidebandCorrelator
Pulsar::Archive::get<ext> name=ProcHistory
Pulsar::Archive::get<ext> failed to find extension type=PN6Pulsar13AuxColdPlasmaE
Pulsar::Archive::add_extension appending AuxColdPlasma
FITSArchive::load_ProcHistory exiting
load_ObsDescription entered
FITSArchive::load_ObsDescription fits_movnam_hdu OBSDESCR: illegal HDU number (301)
FITSArchive::load_digistat entered
Pulsar::FITSArchive::load_digistat no DIG_STAT HDU
FITSArchive::load_DigitiserCounts failed with error
Error::stack
psrfits_move_hdu
Error::FailedCall
Error::message
fits_movnam_hdu (DIG_CNTS): illegal HDU number (301)</ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext></ext>
FITSArchive::load_Passband entered
FITSArchive::load_Passband BANDPASS HDU found
FITSArchive::load_Passband BANDPASS NCH_ORIG=512 BP_NPOL=1
FITSArchive::load_Passband offsets read
FITSArchive::load_Passband scale factors read
FITSArchive::load_Passband data read
Pulsar::Archive::add_extension appending Passband
FITSArchive::load_Passband exiting
FITSArchive::load_CoherentDedispersion entered
FITSArchive::load_FluxCalibratorExtension entered
Pulsar::FITSArchive::load_FluxCalibratorExtension no FLUX_CAL HDU
FITSArchive::load_CalibratorStokes entered
Pulsar::FITSArchive::load_CalibratorStokes no CAL_POLN HDU
FITSArchive::load_PolnCalibratorExtension entered
Pulsar::FITSArchive::load_PolnCalibratorExtension no FEEDPAR HDU
FITSArchive::load_CalibrationInterpolatorExtension entered
FITSArchive::load_FITSSUBHdrExtension entered
psrfits_read_key: INT_UNIT=SEC
psrfits_read_key: INT_TYPE=TIME
psrfits_read_key: NBITS=2
NCH_STRT: keyword not found in header (202)
psrfits_read_key: using default=-1
psrfits_read_key: NSBLK=4096
psrfits_read_key: ZERO_OFF=1.5
psrfits_read_key: SIGNINT=0
Pulsar::Archive::add_extension appending FITSSUBHdrExtension
FITSArchive::load_CrossCovarianceMatrix entered
Pulsar::FITSArchive::load_CrossCovarianceMatrix : no COV_MAT HDU
Pulsar::FITSArchive::load_Predictor try POLYCO
load_polyco entered
load_polyco no POLYCO HDU
Pulsar::FITSArchive::load_Predictor try T2PREDICT
load_T2Predictor entered
load_T2Predictor error: fits_movnam_hdu T2PREDICT: illegal HDU number (301)
psrfits_read_key: INT_TYPE=TIME
FITSArchive::load_integration_state INT_TYPE=TIME
Pulsar::Archive::get<ext> name=FITSHdrExtension
Pulsar::Archive::get<ext> name=ObsExtension
Pulsar::Archive::get<ext> name=Telescope
Pulsar::Archive::get<ext> name=ITRFExtension
Pulsar::Archive::get<ext> name=Receiver
Pulsar::Archive::get<ext> name=WidebandCorrelator
Pulsar::Archive::get<ext> name=ProcHistory
FITSArchive::load_integration_state version=5.3 load SUBINT state
psrfits_read_key: SCALE=FluxDen
psrfits_read_key: POL_TYPE=AA+BB
psrfits_read_key: NCHAN=512
psrfits_read_key: CHAN_BW=-0.5
REFFREQ: keyword not found in header (202)
psrfits_read_key: using default=0
DM: bad string to double convert (409)
psrfits_read_key: using default=0
RM: bad string to double convert (409)
psrfits_read_key: using default=0
NAUX: keyword not found in header (202)
psrfits_read_key: using default=0
FITSArchive::interpret_scale SCALE=FluxDen
FITSArchive::interpret_pol_type POL_TYPE=AA+BB
FITSArchive::interpret_pol_type Signal::Intensity
NINFO: keyword not found in header (202)
psrfits_read_key: using default=0
FITSArchive::load_Pointing_columns no additional Pointing::Info
FITSArchive::load_header exit
Pulsar::Archive::correct 3 entries
Pulsar::Archive::correct testing CalSource
Pulsar::Archive::correct applying CalSource
Pulsar::CalSource::apply
Pulsar::Archive::correct testing Dedispersed
Pulsar::Archive::correct testing DeFaradayed
psrsplit::process PDFB4_Search_J0437-4715.sf
Pulsar::FITSArchive dtor this=0x561437e11d60
Pulsar::BasicArchive dtor this=0x561437e11d60
Pulsar::Archive dtor this=0x561437e11d60 instances=0
Pulsar::IntegrationManager dtor this=0x561437e11d60
Pulsar::Archive::Agent::~Agent
Pulsar::Archive::Agent::~Agent
Pulsar::Archive::Agent::~Agent
Pulsar::Archive::Agent::~Agent
Pulsar::Archive::Agent::~Agent
Pulsar::Archive::Agent::~Agent
Pulsar::Archive::Agent::~Agent
Pulsar::Archive::Agent::~Agent
Pulsar::Archive::Agent::~Agent
Pulsar::Archive::Agent::~Agent
Pulsar::Archive::Agent::~Agent
Pulsar::Archive::Agent::~Agent</ext></ext></ext></ext></ext></ext></ext>
Ticket moved from /p/psrchive/bugs/485/
hi Cherry,
Apart from psredit, which can read header parameters from PSRFITS search-mode data files, PSRCHIVE applications deal only with average pulse profile data. Splitting up into sub-bands sounds like a useful thing to do, so i've moved this to a DSPSR feature request.
Cheers,
Willem