Commit 363f7306 authored by Espen Sollum's avatar Espen Sollum
Browse files

Added sza to netcdf files

parent cb211e16
......@@ -1374,6 +1374,12 @@ class Scene(object):
# atrack[:] = self.atrack
atrack[:] = tmp_atrack
sza = ncf.createVariable('sza',np.dtype('float').char,('npoints',))
sza.standard_name = 'sza'
sza.long_name = 'Satellite Zenith Angle'
sza.units = 'degrees'
sza[:] = self.szas
if hasattr(self, 'AshML'):
ML = ncf.createVariable('ML',np.dtype('float').char,('npoints',))
ML.standard_name = 'ash_mass_loading'
......
......@@ -28,7 +28,7 @@ importlib.reload(IASITools)
#IASI_PATH = "/sat_data_from_fred/eumetcast/EPS-11/"
# local files / 2010 Eyja
IASI_PATH = "/home/eso/repos/ash-iasi/data/"
IASI_PATH = "/home/eso/repos/ash-iasi/data/0416/"
### Paths to output
# locally for testing
......@@ -165,7 +165,7 @@ if __name__ == "__main__":
fns = sorted(glob(fs))
# Just test with 1 file
# fns = [fns[0]]
fns = [fns[0]]
# Invalid files (know to cause core dumps with message like
#GEOS_ERROR: TopologyException: found non-noded intersection between LINESTRING (..."
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment