def _check(self, append=False): iris_grib.save_grib2(self.cube, self.target, append=append) self.save_pairs_from_cube.called_once_with(self.cube) args, kwargs = self.save_messages.call_args self.assertEqual(len(args), 2) messages, target = args self.assertEqual(list(messages), self.messages) self.assertEqual(target, self.target) self.assertEqual(kwargs, dict(append=append))
def _check(self, append=False): iris_grib.save_grib2(self.cube, self.target, append=append) self.save_pairs_from_cube.called_once_with(self.cube) args, kwargs = self.save_messages.call_args self.assertEqual(len(args), 2) messages, target = args self.assertEqual(list(messages), self.messages) self.assertEqual(target, self.target) self.assertEqual(kwargs, dict(append=append))
def _grib_save(cube, target, append=False, **kwargs): # A simple wrapper for the grib save routine, which allows the saver to be # registered without having the grib implementation installed. try: from iris_grib import save_grib2 except ImportError: raise RuntimeError('Unable to save GRIB file - ' '"iris_grib" package is not installed.') save_grib2(cube, target, append, **kwargs)
def _grib_save(cube, target, append=False, **kwargs): # A simple wrapper for the grib save routine, which allows the saver to be # registered without having the grib implementation installed. try: import iris_grib as igrib except ImportError: try: import gribapi except ImportError: raise RuntimeError("Unable to save GRIB file - the ECMWF " "`gribapi` package is not installed.") from iris.fileformats import grib as igrib return igrib.save_grib2(cube, target, append, **kwargs)
def _grib_save(cube, target, append=False, **kwargs): # A simple wrapper for the grib save routine, which allows the saver to be # registered without having the grib implementation installed. try: import iris_grib as igrib except ImportError: try: import gribapi except ImportError: raise RuntimeError('Unable to save GRIB file - the ECMWF ' '`gribapi` package is not installed.') from iris.fileformats import grib as igrib return igrib.save_grib2(cube, target, append, **kwargs)
latitude = DimCoord(np.linspace(29.45, 20.70, 35), standard_name='latitude', units='degrees') longitude = DimCoord(np.linspace(97.33, 106.58, 37), standard_name='longitude', units='degrees') time = DimCoord(0, standard_name='time', units='hours since {}'.format( str(datetime.datetime(2017, 3, 2, 17)))) forecast_period = DimCoord(int(000), standard_name='forecast_period', units='hours') cube = Cube(data, dim_coords_and_dims=[(latitude, 0), (longitude, 1)], units='m s-1', standard_name='wind_speed') cube.add_aux_coord(time) cube.add_aux_coord(forecast_period) cube.coord(axis='X').coord_system = iris.coord_systems.GeogCS(654321) cube.coord(axis='Y').coord_system = iris.coord_systems.GeogCS(654321) # iris_grib.save_grib2(cube, "%s%s_EC-1000hpa-24h_temperature-CHINA_0.25deg.grib2" %(name, fh), append=False) iris_grib.save_grib2(cube, "2017030217%s%s_MICAPS_wind-speed_YUNNAN_27km.grib2" % (name, fh), append=False) # print cube.data
pd.to_datetime('1981-01-01')).days * 24 new_t_coord = iris.coords.DimCoord(date, standard_name='time', units='hours since 1981-01-01') data1.add_dim_coord(new_t_coord, time1_dims) lon_coord1 = data1.coords("longitude")[0] lat_coord1 = data1.coords("latitude")[0] lon_coord1.coord_system = iris.coord_systems.GeogCS( iris.fileformats.pp.EARTH_RADIUS) lat_coord1.coord_system = iris.coord_systems.GeogCS( iris.fileformats.pp.EARTH_RADIUS) data1.add_aux_coord( iris.coords.DimCoord(0, standard_name='forecast_period', units='hours')) data1.add_aux_coord(iris.coords.DimCoord(0, "height", units="m")) data1.var_name = "SST" data1.long_name = "Sea Surface Temperature" data1.standard_name = 'sea_surface_temperature' data1.units = 'K' data1.discipline = 0 data1.category = 0 data1.typeOfFirstFixedSurface = 1 time = file[0:4] + '-' + file[4:6] + '-' + file[6:8] time1 = file[0:4] + '-' + file[4:6] + '-' + file[6:8] print(data1) iris_grib.save_grib2(data1, filename)