def _create_coord_list(self):
"""Read file coordinates into a CIS object"""
from cis.data_io.Coord import Coord, CoordList
from reame.utils import ncdf_read
try:
lon_data, lon_metadata = ncdf_read(self.filenames, "longitude")
lat_data, lat_metadata = ncdf_read(self.filenames, "latitude")
except IndexError:
lon_data, lon_metadata = ncdf_read(self.filenames, "lon")
lat_data, lat_metadata = ncdf_read(self.filenames, "lat")
lat = Coord(lat_data, lat_metadata, "Y")
lat.update_shape()
lat.update_range()
lon = Coord(lon_data, lon_metadata, "X")
lon.update_shape()
lat.update_range()
time_data, time_metadata = ncdf_read(self.filenames, "time")
# Ensure the standard name is set
time_metadata.standard_name = "time"
time = Coord(time_data, time_metadata, "T")
time.convert_TAI_time_to_std_time(ATSR_REFERENCE_TIME)
time.update_shape()
time.update_range()
return CoordList([lat, lon, time])
def AOD550(self):
from reame.utils import aod_without_prior, ncdf_read
aod_data, aod_metadata = ncdf_read(self.filenames, "AOD550")
unc_data, _ = ncdf_read(self.filenames, "AOD550_uncertainty")
typ_data, _ = ncdf_read(self.filenames, "aerosol_type")
new_aod = aod_without_prior(aod_data, unc_data, typ_data, [-1.0] * 10)
return new_aod, aod_metadata
def AOD870(self):
from reame.utils import aod_without_prior, ncdf_read
aod_data, aod_metadata = ncdf_read(self.filenames, "AOD870")
unc_data, _ = ncdf_read(self.filenames, "AOD870_uncertainty")
typ_data, _ = ncdf_read(self.filenames, "aerosol_type")
new_aod = aod_without_prior(aod_data, unc_data, typ_data, [
-0.951, -0.705, -0.711, -0.715, -0.717, -0.869, -0.960, -0.865,
-0.706, -0.407
])
return new_aod, aod_metadata
def _get_qcmask(self):
from reame.utils import ncdf_read
qcdata, _ = ncdf_read(self.filenames, "quality_flag")
return (qcdata != 1).filled(True)
def AOD1600_uncertainty(self):
from reame.utils import ncdf_read
return ncdf_read(self.filenames, "AOD1600_uncertainty")
def AOD1600(self):
from reame.utils import ncdf_read
return ncdf_read(self.filenames, "AOD1600")