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

minor

parent 5fd90aa9
No related branches found
No related tags found
No related merge requests found
import numpy as np import numpy as np
import h5py import h5py
from util.util import get_grid_values, get_grid_values_all, is_night from util.util import get_grid_values, get_grid_values_all, is_night, compute_lwc_iwc
import glob import glob
import os import os
from pathlib import Path from pathlib import Path
...@@ -240,6 +240,18 @@ def scan(directory): ...@@ -240,6 +240,18 @@ def scan(directory):
print(num_files, num_keep) print(num_files, num_keep)
def test_nlcomp(file):
h5f = h5py.File(file, 'r')
cld_phs = get_grid_values_all(h5f, 'cloud_phase', scale_factor_name=None, range_name=None)
reff = get_grid_values_all(h5f, 'cld_reff_nlcomp')
opd = get_grid_values_all(h5f, 'cld_opd_nlcomp')
cld_dz = get_grid_values_all(h5f, 'cld_geo_thick')
lwc_c, iwc_c = compute_lwc_iwc(cld_phs, reff, opd, cld_dz)
return lwc_c, iwc_c
# def run_mean_std(directory): # def run_mean_std(directory):
# #
# data_dct = {name: [] for name in mod_res_params} # data_dct = {name: [] for name in mod_res_params}
......
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