;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; Side-Down relative Reflectivity Calibration ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; .run commonwcrset ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; openwcrf,'/sova4/cifex04/apr20a/W2004-04-20-13-07-55' readwcrdata,144905,145100,/timein openwcrf,'/sova4/cifex04/apr20a/W2004-04-20-15-53-37' readwcrdata,162850,163035,/timein ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; setxpswin,0,3 roll=44 rgh=5 rgv=rgh h=10.*alog10(hh[rgh,*]) v=10.*alog10(vv[rgv,*]) plot,h,yrange=[-65,0],/xst,psym=3 oplot,v & grid ; Exclude returns below -60dBm. This is for simplicity since the mean ; noise is not subtracted ind1=where(h gt -55.) junk=reform(v-h) ;plot,findgen(numrec),junk,psym=3 ,/xst ;oplot,ind1,junk[ind1],psym=5,symsize=0.5 & grid plot,time,junk,psym=3,xst=5 & timeaxis & grid oplot,time,(smooth(junk,51)),thick=3 coeffs=linfit(h,v) plot,h[ind1],v[ind1],psym=3,$ title='Slope='+num2str(coeffs[1],fr=2)+', Bias='+num2str(coeffs[0])+$ '; (V=slope*H+bias)', xtit='H-channel',ytit='V-channel' oplot,!x.crange,coeffs[0]+!x.crange*coeffs[1] & grid junk1=sqrt(radrng[rgh]^2+radrng[rgv]^2) ; dist between beams junk2=radrng[rgh]*sin(roll*!dtor)-radrng[rgv]*sin((90-roll)*!dtor); vertical beam separation ( h-v) print,';'+num2str(rgh,fr=0)+'/'+num2str(rgv,fr=0)+$ ' '+num2str(junk1,fr=0)+'/'+num2str(junk2,fr=0)+$ num2str(mean(junk[ind1]),fr=2,len=9)+ $ num2str(std(junk[ind1]),fr=2,len=10)+$ num2str(-coeffs[0]+(coeffs[1]-1)*50,fr=2,len=10)+num2str(coeffs[1],fr=2,len=8) ;------------------------- 25 March 2004 ------------------------------ ;W2004-03-25-19-24-52 19:25:30-19:27:05 side/down, 44deg right turn ;rgh/rgv HtoVdist std(Zv-Zh) Bias Slope ; hor/ver dB dB @-50dBm ; m ;2/2 190/-3 3.93 2.59 -10.76 0.68 ;3/3 233/-4 3.59 2.74 -9.77 0.69 ;4/4 275/-4 3.18 2.70 -8.53 0.71 ;5/5 318/-5 2.96 2.87 -9.26 0.64 ;7/7 403/-7 2.79 2.84 -8.00 0.66 ;10/10 530/-9 2.50 3.00 -7.39 0.61 ;------------------------- 20 April 2004 ------------------------------ ;W2004-04-20-13-07-55 14:49:05-14:51:00 side/down, 44deg right turn, TURBO OFF ;rgh/rgv HtoVdist std(Zv-Zh) Bias Slope ; hor/ver dB dB @-50dBm ; m ;2/2 190/-3 3.41 2.42 -4.09 0.93 ;3/3 233/-4 3.39 2.64 -4.16 0.90 ;4/4 275/-4 3.50 2.73 -4.27 0.87 ;5/5 318/-5 3.50 2.72 -4.05 0.88 ;7/7 403/-7 4.39 2.89 -4.46 0.96 ;9/9 487/-8 5.25 3.39 -5.43 1.10 ;W2004-04-20-15-53-37 16:28:50-16:30:35 side/down, 44deg right turn ;rgh/rgv HtoVdist std(Zv-Zh) Bias Slope ; hor/ver dB dB @-50dBm ; m ;2/2 190/-3 4.24 3.40 -5.81 0.76 ;3/3 233/-4 4.59 3.43 -5.12 0.88 ;4/4 275/-4 4.78 3.65 -4.64 1.05 ;5/5 318/-5 4.97 4.20 -4.73 1.19 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; Dual down relative Reflectivity Calibration ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; .run commonwcrset ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; openwcrf,'' readwcrdata ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; setxpswin,0,2 tas=100. th=30.*!dtor ; down-forward rgv=15 rgh=round(rgv/cos(th)) rgc=20.*alog10(radrng[rgh]/radrng[rgv]) ;h-range correction; not used for 500ns rgc=0. h=10.*alog10(hh[rgh,*])+rgc v=10.*alog10(vv[rgv,*]) ;plot,h,yrange=[-65,-20],/xst,psym=3 ;oplot,v & grid xcor=c_correlate(hh[rgh,*],vv[rgv,*],findgen(200)) plot,xcor,psym=1 & grid maxc=max(xcor,ind) plot,h[0:numrec-ind-1],yrange=[-65,-10],/xst,psym=3 oplot,v[ind:numrec-1] & grid ; Exclude returns below -60dBm. This is for simplicity since the mean ; noise is not subtracted ind1=where(h[0:numrec-ind-1] gt -60.) junk=h[0:numrec-ind-1]-v[ind:numrec-1] plot,findgen(numrec-ind),junk,psym=3 ,/xst oplot,ind1,junk[ind1],psym=5,symsize=0.5 & grid coeffs=linfit((h[0:numrec-ind-1])[ind1],(v[ind:numrec-1])[ind1]) plot,(h[0:numrec-ind-1])[ind1],(v[ind:numrec-1])[ind1],psym=3,$ title='Slope='+num2str(coeffs[1],fr=2)+', Bias='+num2str(coeffs[0])+$ '; (V=slope*H+bias)',xtit='H-channel',ytit='V-channel' oplot,!x.crange,coeffs[0]+!x.crange*coeffs[1] & grid junk1=tan(th)*radrng[rgv]/tas/timeinterval ; beam time separation assuming tas print,';'+num2str(rgh,fr=0)+'/'+num2str(rgv,fr=0)+num2str(maxc,fr=2,len=11)+$ num2str(ind,fr=0,len=5)+'/'+num2str(ind*timeinterval,fr=2)+$ num2str(junk1,fr=0,len=8)+num2str(mean(junk[ind1]),fr=2,len=9)+ $ num2str(std(junk[ind1]),fr=2,len=10)+$ num2str(-coeffs[0]+(coeffs[1]-1)*50,fr=2,len=10)+num2str(coeffs[1],fr=2,len=8) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; Dual side relative Reflectivity Calibration ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; .run commonwcrset ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; openwcrf,'' readwcrdata ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; setxpswin,0,2 tas=100. th=34.5*!dtor ; side-forward rgh=10 rgv=round(rgh/cos(th)) rgc=20.*alog10(radrng[rgv]/radrng[rgh]) ;v-range correction; not used for 500ns rgc=0. h=10.*alog10(hh[rgh,*]) v=10.*alog10(vv[rgv,*])+rgc xcor=c_correlate(vv[rgv,*],hh[rgh,*],findgen(200)) plot,xcor,psym=1 & grid maxc=max(xcor,ind) plot,h[ind:numrec-1],yrange=[-65,-10],/xst,psym=3 oplot,v[0:numrec-ind-1] & grid ; Exclude returns below -60dBm. This is for simplicity since the mean ; noise is not subtracted ind1=where(v[0:numrec-ind-1] gt -35.) junk=v[0:numrec-ind-1]-h[ind:numrec-1] plot,findgen(numrec-ind),junk,psym=3 ,/xst oplot,ind1,junk[ind1],psym=5,symsize=0.5 & grid coeffs=linfit((v[0:numrec-ind-1])[ind1],(h[ind:numrec-1])[ind1]) plot,(v[0:numrec-ind-1])[ind1],(h[ind:numrec-1])[ind1],psym=3,$ title='Slope='+num2str(coeffs[1],fr=2)+', Bias='+num2str(coeffs[0])+$ '; (H=slope*V+bias)',xtit='V-channel',ytit='H-channel' oplot,!x.crange,coeffs[0]+!x.crange*coeffs[1] & grid junk1=tan(th)*radrng[rgv]/tas/timeinterval print,';'+num2str(rgh,fr=0)+'/'+num2str(rgv,fr=0)+num2str(maxc,fr=2,len=11)+$ num2str(ind,fr=0,len=5)+'/'+num2str(ind*timeinterval,fr=2)+$ num2str(junk1,fr=0,len=8)+num2str(mean(junk[ind1]),fr=2,len=9)+ $ num2str(std(junk[ind1]),fr=2,len=10)+$ num2str(-coeffs[0]+(coeffs[1]-1)*55,fr=2,len=10)+num2str(coeffs[1],fr=2,len=8)