Python main_qc Examples

Example #1

File: test_radar_qc_jetstream.py Project: gustfrontar/radar_qc

sys.path.append('./src/fortran')

import numpy as np
import matplotlib.pyplot as plt
import pyart
import netCDF4
import radar_qc_module as rqc
import os

import conf

#=======================================

conf.options['toporawdatapath']="/home/jruiz/Dropbox/DATA/radar_qc/data/terrain_data/raw/"
conf.options['toporadardatapath']="/home/jruiz/Dropbox/DATA/radar_qc/data/terrain_data/radar/"


# read in the file, create a RadarMapDisplay object
filename = './cfrad.20091117_174348.000_to_20091117_174737.000_PAR_SUR.nc'

#Performs QC operations based on options
[radar , qc_output] = rqc.main_qc( filename , conf.options )

print('End of QC')

plt.figure()
plt.pcolor(qc_output['x'][:,:,0],qc_output['y'][:,:,0],qc_output['smooth_rho'][:,:,0])
plt.colorbar()
plt.show()

Example #2

File: main_in.py Project: gustfrontar/radar_qc

#Obtenemos la lista de objetos radares.
radar_list = ot.read_multiple_files(  file_list , instrument_list )
my_updated_dirs =  []

my_updated_tars =  []

for radar in radar_list :

      print('=============================================================================')
      print(' QC')
      print('=============================================================================')
      print('')
 
      #Call QC routine

      [ radar , qc_output ] = rqc.main_qc( options , radar )

      #Save data in cfradial format

      if write_cfradial  :

         print('=============================================================================')
         print(' WRITING QC OUTPUT IN CFRADIAL FORMAT')
         print('=============================================================================')
         print('')


         ot.save_cfradial( datapath_out + '/cfradial/' , radar )

      print('')
      print('=============================================================================')

Example #3

File: test_radar_qc.py Project: gustfrontar/LETKF_WRF

#Attenuation parameters

options['attfiltertr']=20.0       #Attenuation threshold in dBZ
options['attcalerror']=1.0        #Calibration error
options['attfilter_save']=True    #Save filter fields

#Blocking parameters


#Detect missing parameters


#=======================================


# read in the file, create a RadarMapDisplay object
filename = '/home/jruiz/share/DATA/OBS/OBS_REAL_PARANA_20091117_CFRADIAL/cfrad.20091117_200345.000_to_20091117_200734.001_PAR_SUR.nc3'

#Performs QC operations based on options
[radar , qc_output] = rqc.main_qc( filename , options )

print('End of QC')

qc_output['echo_top'][ qc_output['echo_top'] == options['undef'] ] = 0.0

plt.figure()
plt.pcolor(qc_output['echo_top'][:,:,1])
plt.show()