##// END OF EJS Templates
schain-jc
RSS
26/03/2018
ebocanegra -
r1153:ea60e85af34e
parent child
Show More
Side by Side Unified
Context file:
r1153:ea60e85af34e Loading diff...:
Show file before | Show file after | Hide comments
@@ -274,6 +274,12 class SpectralFilters(Operation):
274 274 SPCcut[i,0:int(Breaker2R[0]),:] = dataOut.noise[i]
275 275 SPCcut[i,-int(Delta):,:] = dataOut.noise[i]
276 276
277 SPCcut[i]=SPCcut[i]- dataOut.noise[i]
278 SPCcut[ numpy.where( SPCcut<0 ) ] = 1e-20
279
280 SPCroll[i]=SPCroll[i]-dataOut.noise[i]
281 SPCroll[ numpy.where( SPCroll<0 ) ] = 1e-20
282
277 283 #self.spc[i, 0:int(Breaker1W[0]) ,:] = dataOut.noise[i]
278 284 #self.spc[i, int(Breaker2W[0]):self.Num_Bin ,:] = dataOut.noise[i]
279 285
@@ -283,6 +289,7 class SpectralFilters(Operation):
283 289
284 290
285 291
292
286 293 SPC_ch1 = SPCroll
287 294
288 295 SPC_ch2 = SPCcut
@@ -752,20 +759,20 class PrecipitationProc(Operation):
752 759
753 760 Numerator = ( (4*numpy.pi)**3 * aL**2 * 16 * numpy.log(2) )
754 761 Denominator = ( Pt * Gt * Gr * Lambda**2 * SPEED_OF_LIGHT * tauW * numpy.pi * ThetaT * ThetaR)
755 RadarConstant = 1/4.1396e+08# Numerator / Denominator #
762 RadarConstant = 1e-10 * Numerator / Denominator #
756 763 print '***'
757 764 print '*** RadarConstant' , RadarConstant, '****'
758 765 print '***'
759 766 ''' ============================= '''
760 767
761 self.spc[0] = self.spc[0]-dataOut.noise[0]
762 self.spc[1] = self.spc[1]-dataOut.noise[1]
763 self.spc[2] = self.spc[2]-dataOut.noise[2]
768 self.spc[0] = (self.spc[0]-dataOut.noise[0])
769 self.spc[1] = (self.spc[1]-dataOut.noise[1])
770 self.spc[2] = (self.spc[2]-dataOut.noise[2])
764 771
765 772 self.spc[ numpy.where(self.spc < 0)] = 0
766 773
767 SPCmean = numpy.mean(self.spc,0) - numpy.mean(dataOut.noise)
768 SPCmean[ numpy.where(SPCmean < 0)] = 1e-20
774 SPCmean = (numpy.mean(self.spc,0) - numpy.mean(dataOut.noise))
775 SPCmean[ numpy.where(SPCmean < 0)] = 0
769 776
770 777 ETAn = numpy.zeros([self.Num_Bin,self.Num_Hei])
771 778 ETAv = numpy.zeros([self.Num_Bin,self.Num_Hei])
@@ -823,7 +830,7 class PrecipitationProc(Operation):
823 830
824 831 Z[R] = numpy.nansum( N_dist[:,R] * (D_range[:,R])**6 )*1e-18
825 832
826 RR[R] = 6*10**-4.*numpy.pi * numpy.nansum( D_range[:,R]**3 * N_dist[:,R] * Velrange[0:self.Num_Bin] ) #Rainfall rate
833 RR[R] = 3.6e-6*1e-9*6*10**-4.*numpy.pi * numpy.nansum( D_range[:,R]**3 * N_dist[:,R] * Velrange[0:self.Num_Bin] ) #Rainfall rate
827 834
828 835 Ze[R] = (numpy.nansum( ETAn[:,R]) * Lambda**4) / ( numpy.pi**5 * Km)
829 836
@@ -841,25 +848,25 class PrecipitationProc(Operation):
841 848 dataOut.channelList = [0,1,2]
842 849
843 850 dataOut.data_param[0]=dBZ
844 dataOut.data_param[1]=dBZe
851 dataOut.data_param[1]=RR2
845 852 dataOut.data_param[2]=RR
846 853
847 854 #print 'VELRANGE', Velrange
848 855 print 'Range', len(Range)
849 856 print 'delv',del_V
850 857 #print 'DRANGE', D_range[:,50]
851 print 'NOISE', dataOut.noise[0]
858 print 'NOISE', dataOut.noise[0], 10*numpy.log10(dataOut.noise[0])
852 859 print 'radarconstant', RadarConstant
853 860 print 'Range', Range
854 print 'ETAn SHAPE', ETAn.shape
855 print 'ETAn ', numpy.nansum(ETAn, axis=0)
856 print 'ETAd ', numpy.nansum(ETAd, axis=0)
861 # print 'ETAn SHAPE', ETAn.shape
862 # print 'ETAn ', numpy.nansum(ETAn, axis=0)
863 # print 'ETAd ', numpy.nansum(ETAd, axis=0)
857 864 print 'Pr ', numpy.nansum(Pr, axis=0)
858 865 print 'dataOut.SPCparam[1]', numpy.nansum(dataOut.SPCparam[1][0], axis=0)
859 print 'Ze ', dBZe
860 print 'Z ', dBZ
861 #print 'RR2 ', RR2
862 #print 'RR ', RR
866 # print 'Ze ', dBZe
867 # print 'Z ', dBZ
868 print 'RR2 ', RR2
869 print 'RR ', RR
863 870 #print 'RR2 ', dBRR2
864 871 #print 'D_mean', D_mean
865 872 #print 'del_V', del_V
General Comments 0
You need to be logged in to leave comments. Login now