From 27308d25dc814cf63c104b3b96a0161d85837034 Mon Sep 17 00:00:00 2001
From: rink <rink@ssec.wisc.edu>
Date: Wed, 14 Oct 2020 11:56:34 -0500
Subject: [PATCH] remove dead code
---
modules/aeolus/aeolus_amv.py | 91 +-----------------------------------
1 file changed, 1 insertion(+), 90 deletions(-)
diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py
index 2f02e0f2..81c863f6 100644
--- a/modules/aeolus/aeolus_amv.py
+++ b/modules/aeolus/aeolus_amv.py
@@ -14,11 +14,6 @@ num_elems = 5424
num_lines = 5424
-first_time = True
-ftimes = []
-flist = None
-
-
class MyGenericException(Exception):
def __init__(self, message):
self.message = message
@@ -141,70 +136,6 @@ class CarrStereo(AMVFiles):
return self.meta_dict
-def get_datetime(pathname):
- fname = os.path.split(pathname)[1]
- toks = fname.split('_')
- dstr = toks[4]
- tstr = toks[5]
- dtstr = dstr+tstr
- dto = datetime.datetime.strptime(dtstr, '%Y%j%H%M').replace(tzinfo=timezone.utc)
-
- return dto
-
-
-def get_datetime_ops(pathname):
- fname = os.path.split(pathname)[1]
- toks = fname.split('_')
- dtstr = toks[3]
- dtstr = dtstr[:-3]
- dto = datetime.datetime.strptime(dtstr, 's%Y%j%H%M').replace(tzinfo=timezone.utc)
-
- return dto
-
-
-def get_datetime_carr(pathname):
- fname = os.path.split(pathname)[1]
- toks = fname.split('_')
- dtstr = toks[3]
- dto = datetime.datetime.strptime(dtstr, '%Y%j.%H%M.ch').replace(tzinfo=timezone.utc)
-
- return dto
-
-
-def get_file_containing_time(timestamp, files_path, file_time_span, amv_source='OPS', band='14'):
- global first_time, ftimes, flist
-
- if first_time is True:
-
- if amv_source == 'OPS':
- flist = glob.glob(files_path + 'OR_ABI-L2-DMWF*'+'C'+band+'*.nc')
- dto_func = get_datetime_ops
- elif amv_source == 'SSEC':
- flist = glob.glob(files_path + '*WINDS_AMV_EN-'+band+'*.nc')
- dto_func = get_datetime
- elif amv_source == 'CARR':
- flist = glob.glob(files_path + 'tdw_qc_GOES*'+'ch_'+band+'.nc')
- dto_func = get_datetime_carr
-
- for pname in flist: # TODO: make better with regular expressions (someday)
- dto = dto_func(pname)
- dto_start = dto
- dto_end = dto + datetime.timedelta(minutes=file_time_span)
- ftimes.append((dto_start.timestamp(), dto_end.timestamp()))
-
- first_time = False
-
- k = -1
- for i in range(len(ftimes)):
- if (timestamp >= ftimes[i][0]) and (timestamp < ftimes[i][1]):
- k = i
- break
- if k < 0:
- return None, None, None
-
- return flist[k], ftimes[k], k
-
-
# imports the S4 NOAA output
# filename: full path as a string, '/home/user/filename'
# returns a dict: time -> list of profiles (a profile is a list of levels)
@@ -380,26 +311,6 @@ def get_search_box(nav, lon, lat):
return c_rng, l_rng
-# Framework
-amv_dir_name = 'Wind_Dir'
-amv_spd_name = 'Wind_Speed'
-amv_lon_name = 'Longitude'
-amv_lat_name = 'Latitude'
-amv_elem_name = 'Element'
-amv_line_name = 'Line'
-amv_press_name = 'MedianPress'
-# -------------------------------
-# Carr stereo
-amv_lon_name = 'Lon'
-amv_lat_name = 'Lat'
-amv_press_name = 'pres'
-amv_h3d_name = 'H_3D'
-amv_alt_name = 'Altitude'
-# -------------------------------
-sub_lon = -137.0 # GOES-17
-# sub_lon = -75.0 # GOES-16
-
-
# aeolus_dict: time -> profiles
# amv_files_path: directory containing AMVs, '/home/user/amvdir/'
# return dict: aeolus time -> tuple (amv_lon, amv_lat, amv_pres, amv_spd, amv_dir)
@@ -592,7 +503,7 @@ def create_amv_to_aeolus_match_file(aeolus_files_dir, amv_files_dir, outfile=Non
# match_file: pathname for the product file
# dt_str_0: start time (YYYY-MM-DD_HH:MM)
# dt_str_1: end time (YYYY-MM-DD_HH:MM)
-# amv_var_names: list of amv parameters (see match_file) to subset
+# amv_var_names: list of amv parameters (see match_file ncdump) to subset
# returns: Xarray.DataArray
# amvs[nprofs, max_num_amvs_per_prof, num_of_params], profs[nprofs, max_num_levs_per_prof, num_of_params],
# prof_locs[nprofs, (lon, lat)
--
GitLab