pointer to pyart

MeteoSwiss · Jan 22, 2024 · 7088012 · 7088012
1 parent 7b05689
commit 7088012
Show file tree

Hide file tree

Showing 5 changed files with 17 additions and 4 deletions.
diff --git a/src/pyart b/src/pyart
diff --git a/src/pyrad_proc/pyrad/__init__.py b/src/pyrad_proc/pyrad/__init__.py
@@ -3,7 +3,6 @@
 =====================================
 
 """
-
 # Detect if we're being called as part of Pyrad's setup procedure
 try:
     __PYRAD_SETUP__
@@ -19,7 +18,7 @@
     # Make sure that deprecation warnings get printed by default
     import warnings as _warnings
     _warnings.simplefilter("always", DeprecationWarning)
-
+    
     # print out helpful message if build fails or importing from source tree
     # fvj built not checked for the moment
     # from . import __check_build

diff --git a/src/pyrad_proc/pyrad/graph/plots_aux.py b/src/pyrad_proc/pyrad/graph/plots_aux.py
@@ -449,7 +449,7 @@ def get_norm(field_name, field_dict={}, isxarray=False):
     ticklabs = None
 
     ref_dict = pyart.config.get_metadata(field_name)
-    cmap = mpl.cm.get_cmap(pyart.config.get_field_colormap(field_name))
+    cmap = mpl.colormaps.get_cmap(pyart.config.get_field_colormap(field_name))
 
     if field_dict is not None and 'boundaries' in field_dict:
         if isxarray:

diff --git a/src/pyrad_proc/pyrad/proc/process_monitoring.py b/src/pyrad_proc/pyrad/proc/process_monitoring.py
@@ -138,6 +138,7 @@ def process_selfconsistency_kdp_phidp(procstatus, dscfg, radar_list=None):
         return None, None
 
     fzl = None
+    temp_ref ='fixed_fzl'
     # determine which freezing level reference
     if temp is not None:
         if temp in radar.fields:
@@ -399,6 +400,7 @@ def process_selfconsistency_bias(procstatus, dscfg, radar_list=None):
         return None, None
 
     fzl = None
+    temp_ref ='fixed_fzl'
     # determine which freezing level reference
     if temp is not None:
         if temp in radar.fields:
@@ -719,6 +721,7 @@ def process_selfconsistency_bias2(procstatus, dscfg, radar_list=None):
         return None, None
 
     fzl = None
+    temp_ref ='fixed_fzl'
     # determine which freezing level reference
     if temp is not None:
         if temp in radar.fields:
@@ -1068,6 +1071,7 @@ def process_rhohv_rain(procstatus, dscfg, radar_list=None):
         return None, None
 
     fzl = None
+    temp_ref ='fixed_fzl'
     # determine which freezing level reference
     if temp_field is not None:
         if temp_field in radar.fields:
@@ -1238,6 +1242,7 @@ def process_zdr_precip(procstatus, dscfg, radar_list=None):
     # if data in and above the melting layer has to be filtered determine the
     # field to use
     fzl = None
+    temp_ref = 'fixed_fzl'
     ml_filter = True
     if 'ml_filter' in dscfg:
         ml_filter = dscfg['ml_filter']

diff --git a/src/pyrad_proc/pyrad/proc/process_phase.py b/src/pyrad_proc/pyrad/proc/process_phase.py
@@ -21,6 +21,7 @@
     process_attenuation
 
 """
+import warnings
 
 from copy import deepcopy
 from warnings import warn
@@ -32,6 +33,11 @@
 from ..io.io_aux import get_datatype_fields
 from ..io.read_data_sensor import read_fzl_igra
 
+# Ignore warning in process_phidp_kdp_Maesaka
+warnings.filterwarnings("ignore", 
+    message=".*'partition' will ignore the 'mask' of the MaskedArray.*", 
+    category=UserWarning)
+
 def process_correct_phidp0(procstatus, dscfg, radar_list=None):
     """
     corrects phidp of the system phase
@@ -438,6 +444,7 @@ def process_phidp_kdp_Maesaka(procstatus, dscfg, radar_list=None):
     thickness = dscfg.get('ml_thickness', 700.)
 
     fzl = None
+    temp_ref = 'fixed_fzl'
     # determine which freezing level reference
     if temp_field is not None:
         if temp_field in radar.fields:
@@ -624,6 +631,7 @@ def process_phidp_kdp_lp(procstatus, dscfg, radar_list=None):
         return None, None
 
     fzl = None
+    temp_ref = 'fixed_fzl'
     # determine which freezing level reference
     if temp_field is not None:
         if temp_field in radar.fields:
@@ -1169,6 +1177,7 @@ def process_attenuation(procstatus, dscfg, radar_list=None):
 
     # determine which freezing level reference
     fzl = None
+    temp_ref = 'fixed_fzl'
     if temp is not None:
         if temp in radar.fields:
             temp_ref = 'temperature'
+3 −1		doc/environment.yml
+5 −4		doc/source/conf.py
+1 −0		doc/source/notebook-gallery.rst
+546 −0		doc/source/notebooks/read_mch_metranet_data.ipynb
+106 −506		pyart/__check_build/_check_build.c
+3 −2		pyart/config.py
+711 −1,662		pyart/correct/_fast_edge_finder.c
+452 −1,340		pyart/correct/_unwrap_1d.c
+454 −1,348		pyart/correct/_unwrap_2d.c
+454 −1,348		pyart/correct/_unwrap_3d.c
+13 −4		pyart/correct/attenuation.py
+32 −12		pyart/correct/bias_and_noise.py
+1,594 −2,973		pyart/io/_sigmetfile.c
+469 −1,415		pyart/io/nexrad_interpolate.c
+801 −2,349		pyart/map/_gate_to_grid_map.c
+456 −1,380		pyart/map/_load_nn_field_data.c
+863 −1,912		pyart/map/ckdtree.c
+395 −906		pyart/retrieve/_gecsx_functions_cython.c
+463 −1,394		pyart/retrieve/_kdp_proc.c
+1 −1		pyart/retrieve/wind.py