masProjects idl routines

Last modified: Thu Mar 16 13:07:07 2023.


List of Routines


Routine Descriptions

MASRDRSAT - LOOK FOR RADAR SATURATION

[Next Routine] [List of Routines]
NAME:
masrdrsat - look for radar saturation 
SYNTAX: istat=masrdrsat(fnmiar,rfiFreq,aero=aero,faa=faa,$
                           remy=remy,bwsrch=bwsrch,rfiBaseLnFreq=rfiBaseLnFreq)
ARGS:
   FNMIAR[N]: {} array holding files to process. see masfilelist
   rfiFreq[m]: float cfr for each rfi that you are interested in.
                     see also rdr keywords.
KEYWORDS:
bwsrch: float Number Mhz around center of rdr to search for peak.
              default is larger of 1 Mhz or 1.5 channels
aero:  if set then include aerostat freq's at [1241.74,1256.5,1246.2,1261.25] 
faa  : if set then include faa radars at [1330, 1350]
remy : if set then include remy rdr at [1270,1290]
savFile: string Name for save file.. leave off the .sav 
savaz: if set then record az positions
savza: if set then record za positions
deffreq: if set then use default freq rather than the header values (for early datasets).
rfiBaseLnFreq: Mhz . If supplied then include the peak hold around this freq. You can
              use it for a baseline. Needed if the time sampling is slower than the
              rdr ipps.
verbose:      if set then print record numbers as we process them

DESCRIPTION:
   Assume the two 170 Mhz bands have lo's of 175 and 325 Mhz and that the
skycfr is 1375Mhz. The radar band will be in band1 and the clean band will
be in band0 (since first lo is high side).

(See /pkg/rsi/local/libao/phil/mas/Proj/masrdrsat.pro)


MASRDRSATDEC - DECIMATE, PEAK HOLD RADAR SATURATION DATASET

[Previous Routine] [Next Routine] [List of Routines]
NAME:
masrdrsatdec - decimate, peak hold radar saturation dataset
SYNTAX: masrdrsatdec,todec,tpcln,tprdr,tpclnDec,tpRdrDec,medClnDec,medRdrDec,nptsD,indAr
ARGS:
   todec: long  sample to peak hold/decimate
tpCln[npnt,nbeam]: float total power data clean band
tpRdr[npt,nbeam,nrdr]: float peak power from rdr band
RETURNS:
tpClnDec[npnt/todec,nbeam]: float smoothed, decimated clean band
tpRdrDec[npnt/todec,nbeam,nrdr]: float peak held rdr band
medClnDec[npnt/todec,nbeam]: float median value of data before smoothing
medRdrDec[npnt/todec,nbeam,nrdr]:  median value of data before smoothing
nptsD: long    number of points each beam after smooth,decimate.
indAr[npnt/todec]: long    indices into tpcln,tprdr (by beam) for the centers
               of each smooth,decimated point. Use this to extact the 
               za info from the za arr.
DESCRIPTION:
   boxcar smooth and then decimate the total power clean band. Do a 
peak hold over the same channels for the radar samples.

(See /pkg/rsi/local/libao/phil/mas/Proj/masrdrsatdec.pro)


MASRDRSATPLT - PLOT RADAR SATURATION DATA

[Previous Routine] [Next Routine] [List of Routines]
NAME:
masrdrsatplt - plot radar saturation data
SYNTAX: masrdrsatplt,tpCln,tpRdr,beamList,rdrFrqAr,hdr,tit=tit,smo=smo,$
           wait=wait,tm=tm,tpRdrN=tpRdrN,tpClnN=tpClnN,smo1inc=smo1inc,$
			indRdrBl=indRdrBl
ARGS:
     tpCln[n,nbeams]: float total power for clean band
tpRdr[n,nbeams,nrdr]: float peak power for radar frequencies
    beamList[nbeams]: int   beam numbers for beams used
      rdrFrqAr[nrdr]: float lower freq Mhz for each rdr freq used
              hdr{}:        header from save file

KEYWORDS:
   tit: string  to add to standard title
   smo: long    amount to smooth tpCln on first page. default:101
  wait:         if set then wait for keypress between plots
 expTm: float   time for expanded plot
durExp: float   duration in secs for expanded plot
smo1inc: float  increment for smoothed clean band. def:.008
indRdrBl: int   if present then this is the entry in the rdrFrqAr that should be used
                to baseline the rdr data. Normally the sampling is fast enough that Tsys
                can be found in each tpRdr time series. For slow sampling you should
                call masrdrsat with rfiBaseLnFreq=rfiBaseLnFreq and then use that index
                here. eq. rfiBaseLnFreq=1300. and then find index in rdrFrqar that=1300

RETURNS:
                tm[n]: float time (secs) from start of scan
tpRdrN[n,nbeams,nrdr]: float normalized to Tsys
     tpClnN[n,nbeams]: float normalized to Tsys
       frqRdr[lenfft]: float freq array for radar band
       frqCln[lenfft]: float freq array for clean band

(See /pkg/rsi/local/libao/phil/mas/Proj/masrdrsatplt.pro)


MASRDRSATZAPLT - PLOT ZA STRIPS FROM RADAR SATURATION DATA.

[Previous Routine] [List of Routines]
NAME:
masrdrsatzaplt - plot za strips from radar saturation data.
SYNTAX: masrdrsatzaplt,hdr,tpCln,tpRdr,za,rdrFrqAr,beamList,$
        wait=wait,deg=deg,zaStart=zaStart,zaStep=zaStep,tit=tit,$
        rdrFrqOrd=rdrFrqOrd,vmaxRdr=vmaxRdr,medRdr=medRdr,medCln=medCln
ARGS:
hdr{}:     fits binary header for this file
tpCln[nsmp,nbeam]: float the power data from the clean band
tpRdr[nsmp,nbeam,nrdr]: float the peak data from the radar band.
za[nsmp]:  the za for each sample
rdrFrqAr[nrdr]: float start freq each rdr band (see masrdrsat()).
beamList[nbeam]: int beam number for each array index
KEYWORDS:

wait   :        if set then wait for return between each plot.
deg:   : long   order of polynomial to fit when removing za dependence
                 over zastep points. Default is 1.
zaStart: float  first za to  use (deg) def=min value in dataset
zaStep:  float  za length plot each page. default is 2.5 deg
tit:  ''        string to add to title
rdrFrqOrd[nrdr]: int order to plot the radar data. This is relative to 
                 how the dat is stored.
vmaxRdr: float   max value for the rdr plots in db's above Tsys.
medCln[nsmp,nbeam]:float  if supplied then the median value at each sample.
                 use this to normalize the cln band data.
medRdr[nsmp,nbeam,nrdr]:float  if supplied then the median value at each sample.
                 use this to normalize the rdr data.

DESCRIPTION:
   make plots of radar saturation za strips. Data comes from masrdrsat().
Plot zaStep degrees per page (default is 2.5 degrees). On each page
plot the avg total power for the clean band followed by the
peak values in the radar band.

EXAMPLE:
   savFile='rdrsat.sav'
   nfound=masrdrsat(fnmiar,rfifreq,/aero,/faa,/remy,savFile=savFile,$
                    /verb,/savza)

   restore,savFile,/verb
   todec=1000
   masrdrsatdec(todec,,tpCln,tpRdr,tpClnD,tpRdrD,nptDec,indAr
   zaStart=2.
   zaStep=2.5
   masrdrsatzaplt,hdrD,tpClnD,tpRdrD,zaAr[indAr],rdrFrqAr,beamList,wait=wait,$
       zaStart=zaStart,zaStep=zaStep,tit=tit,deg=2

(See /pkg/rsi/local/libao/phil/mas/Proj/masrdrsatzaplt.pro)