Skip to content
Snippets Groups Projects
Commit 7f91fe57 authored by tomrink's avatar tomrink
Browse files

minor...

parent 53b5ac7f
No related branches found
No related tags found
No related merge requests found
...@@ -9,6 +9,10 @@ half_width = 30 # search box centered on CALIOP profile (FGF coordinates) ...@@ -9,6 +9,10 @@ half_width = 30 # search box centered on CALIOP profile (FGF coordinates)
num_elems = 5424 num_elems = 5424
num_lines = 5424 num_lines = 5424
amv_fmwk_exmpl_file = '/ships19/cloud/scratch/4TH_AMV_INTERCOMPARISON/FMWK2_AMV/GOES16_ABI_2KM_FD_2019293_0020_34_WINDS_AMV_EN-14CT.nc'
amv_ops_exmpl_file = '/ships19/cloud/scratch/AMV_FAST_BIAS/OPS_AMVS/OR_ABI-L2-DMWF-M6C14_G16_s20201050000166_e20201050009474_c20201050023226.nc'
caliop_clvrx_exmpl_file = '/data/Personal/rink/clavrx_calipso/g16_s20201050200_06kmCLay.matchup.calipso.h5'
def get_search_box(cc, ll): def get_search_box(cc, ll):
...@@ -41,7 +45,6 @@ def match_calipso_clavrx_to_amvs(calipso_clavrx_ds, calipso_clavrx_file, amv_fil ...@@ -41,7 +45,6 @@ def match_calipso_clavrx_to_amvs(calipso_clavrx_ds, calipso_clavrx_file, amv_fil
nom_time = calipso_clavrx_ds.get_datetime(calipso_clavrx_file).timestamp() nom_time = calipso_clavrx_ds.get_datetime(calipso_clavrx_file).timestamp()
calipso_clavrx_nc4 = Dataset(calipso_clavrx_file) calipso_clavrx_nc4 = Dataset(calipso_clavrx_file)
# amv_fname, ftime, f_idx = amv_files.get_file_containing_time(nom_time)
amv_fname, ftime, f_idx = amv_files.get_file(nom_time) amv_fname, ftime, f_idx = amv_files.get_file(nom_time)
if f_idx is None: if f_idx is None:
return None return None
...@@ -155,10 +158,8 @@ def create_file(match_dct, output_path, target_filepath, caliop_clavrx_params, a ...@@ -155,10 +158,8 @@ def create_file(match_dct, output_path, target_filepath, caliop_clavrx_params, a
amvs_cnt = np.array(amvs_cnt) amvs_cnt = np.array(amvs_cnt)
# Sample file to retrieve and copy variable attributes # Sample file to retrieve and copy variable attributes
# rg_exmpl = Dataset('/ships19/cloud/scratch/4TH_AMV_INTERCOMPARISON/FMWK2_AMV/GOES16_ABI_2KM_FD_2019293_0020_34_WINDS_AMV_EN-14CT.nc', 'r') rg_exmpl = Dataset(amv_ops_exmpl_file)
rg_exmpl = Dataset('/ships19/cloud/scratch/AMV_FAST_BIAS/OPS_AMVS/OR_ABI-L2-DMWF-M6C14_G16_s20201050000166_e20201050009474_c20201050023226.nc') caliop_clavrx_exmpl = Dataset(caliop_clvrx_exmpl_file, 'r')
caliop_clavrx_exmpl = Dataset('/data/Personal/rink/clavrx_calipso/g16_s20201050200_06kmCLay.matchup.calipso.h5', 'r')
# the top level group for the output file # the top level group for the output file
rootgrp = Dataset(output_filepath, 'w', format='NETCDF4') rootgrp = Dataset(output_filepath, 'w', format='NETCDF4')
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment