t42_grid.py 文件源码

def __init__(self, num_lons=128, num_lats=64, num_levels=1,
                 mask_file=None, description=''):

        levels = range(num_levels)

        self.type = 'Spectral'
        self.full_name = 'T42'

        if mask_file:
            with nc.Dataset(mask_file) as f:
                try:
                    mask = np.round(f.variables['WGOCN'][0, 0, :, :-1])
                except KeyError as e:
                    print("Error: var WGOCN not in {}.".format(mask_file),
                          file=sys.stderr)
                    raise e
        else:
            # Default is all unmasked.
            mask = np.zeros((num_lats, num_lons))

        assert mask.shape[0] == num_lats
        assert mask.shape[1] == num_lons

        super(T42Grid, self).__init__(num_lons, num_lats, mask_t=mask,
                                      levels=levels, description=description)