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

get clavrx parameters from a saved pickle file

parent 14c4f572
No related branches found
No related tags found
No related merge requests found
...@@ -3,28 +3,32 @@ from datetime import timezone ...@@ -3,28 +3,32 @@ from datetime import timezone
import glob import glob
import numpy as np import numpy as np
import re import re
import pickle
from netCDF4 import Dataset from netCDF4 import Dataset
from pathlib import Path from pathlib import Path
from util.util import GenericException from util.util import GenericException
from util.lon_lat_grid import LonLatGrid from util.lon_lat_grid import LonLatGrid
from util.geos_nav import GEOSNavigation from util.geos_nav import GEOSNavigation
from util.util import homedir from util.setup import ancillary_path
def get_parameters_clavrx(filename=homedir+'data/clavrx/clavrx_OR_ABI-L1b-RadF-M6C01_G16_s20192930000343.level2.nc'): # def get_parameters_clavrx(filename=homedir+'data/clavrx/clavrx_OR_ABI-L1b-RadF-M6C01_G16_s20192930000343.level2.nc'):
rg = Dataset(filename, 'r') def get_parameters_clavrx(filename=ancillary_path+'clavrx_parameters.pkl'):
var_s = rg.variables # rg = Dataset(filename, 'r')
var_names = list(var_s.keys()) # var_s = rg.variables
var_names_2d = [] # var_names = list(var_s.keys())
# var_names_2d = []
for str in var_names: #
v = var_s[str] # for str in var_names:
if len(v.shape) == 2: # v = var_s[str]
if not (str.find('latitude') != -1 or str.find('longitude') != -1): # if len(v.shape) == 2:
var_names_2d.append(str) # if not (str.find('latitude') != -1 or str.find('longitude') != -1):
# var_names_2d.append(str)
rg.close() #
# rg.close()
var_names_2d = pickle.load(open(filename, 'rb'))
return var_names_2d return var_names_2d
......
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