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)
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)
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)
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)