import s3fs
import xarray as xr

#fs = s3fs.S3FileSystem(anon=True)
#aws_url = 'https://rgw.met.no/6c6a817bd3bb47588bcc9e1221cb52b5:escience2022/Ada/conc/1979/01/ice_conc_sh_ease2-250_cdr-v2p0_197901301200.nc'

import urllib
new_x = "temp.nc"
x = 'https://rgw.met.no/6c6a817bd3bb47588bcc9e1221cb52b5:escience2022/Ada/conc/1979/01/ice_conc_sh_ease2-250_cdr-v2p0_197901301200.nc'
urllib.request.urlretrieve(x, new_x)

('temp.nc', <http.client.HTTPMessage at 0x7f82524009a0>)

esacci = xr.open_dataset('temp.nc')
esacci

<xarray.Dataset>
Dimensions:                   (time: 1, nv: 2, xc: 432, yc: 432)
Coordinates:
  * time                      (time) datetime64[ns] 1979-01-30T12:00:00
  * xc                        (xc) float64 -5.388e+03 -5.362e+03 ... 5.388e+03
  * yc                        (yc) float64 5.388e+03 5.362e+03 ... -5.388e+03
    lat                       (yc, xc) float32 ...
    lon                       (yc, xc) float32 ...
Dimensions without coordinates: nv
Data variables:
    Lambert_Azimuthal_Grid    int32 -2147483647
    time_bnds                 (time, nv) datetime64[ns] 1979-01-30 1979-01-31
    ice_conc                  (time, yc, xc) float64 ...
    raw_ice_conc_values       (time, yc, xc) float64 ...
    total_standard_error      (time, yc, xc) float64 ...
    smearing_standard_error   (time, yc, xc) float64 ...
    algorithm_standard_error  (time, yc, xc) float64 ...
    status_flag               (time, yc, xc) float32 ...
Attributes: (12/40)
    title:                     Sea Ice Concentration Climate Data Record from...
    summary:                   This climate data record of sea ice concentrat...
    topiccategory:             Oceans ClimatologyMeteorologyAtmosphere
    keywords:                  Earth Science > Cryosphere > Sea Ice > Sea Ice...
    keywords_vocabulary:       GCMD Science Keywords
    geospatial_lat_min:        -90.0
    ...                        ...
    Conventions:               CF-1.6 ACDD-1.3
    standard_name_vocabulary:  CF Standard Name Table (v36, 21 September 2016)
    id:                        10.15770/EUM_SAF_OSI_0008
    product_version:           2.0
    history:                   2017-02-14T19:39:25Z creation\n2022-06-07: Cor...
    geospatial_bounds_crs:     EPSG:6932

xarray.Dataset

• Dimensions:
  • time: 1
  • nv: 2
  • xc: 432
  • yc: 432
• Coordinates: (5)
  • time
    (time)
    datetime64[ns]
    1979-01-30T12:00:00

    long_name :

    reference time of product

    standard_name :

    time

    axis :

    T

    bounds :

    time_bnds

    array(['1979-01-30T12:00:00.000000000'], dtype='datetime64[ns]')

  • xc
    (xc)
    float64
    -5.388e+03 -5.362e+03 ... 5.388e+03

    units :

    km

    long_name :

    x coordinate of projection (eastings)

    standard_name :

    projection_x_coordinate

    array([-5387.5, -5362.5, -5337.5, ...,  5337.5,  5362.5,  5387.5])

  • yc
    (yc)
    float64
    5.388e+03 5.362e+03 ... -5.388e+03

    units :

    km

    long_name :

    y coordinate of projection (northings)

    standard_name :

    projection_y_coordinate

    array([ 5387.5,  5362.5,  5337.5, ..., -5337.5, -5362.5, -5387.5])

  • lat
    (yc, xc)
    float32
    ...

    units :

    degrees_north

    long_name :

    latitude coordinate

    standard_name :

    latitude

    [186624 values with dtype=float32]

  • lon
    (yc, xc)
    float32
    ...

    units :

    degrees_east

    long_name :

    longitude coordinate

    standard_name :

    longitude

    [186624 values with dtype=float32]

• Data variables: (8)
  • Lambert_Azimuthal_Grid
    ()
    int32
    ...

    grid_mapping_name :

    lambert_azimuthal_equal_area

    longitude_of_projection_origin :

    0.0

    latitude_of_projection_origin :

    -90.0

    false_easting :

    0.0

    false_northing :

    0.0

    semi_major_axis :

    6378137.0

    inverse_flattening :

    298.257223563

    proj4_string :

    +proj=laea +lon_0=0 +datum=WGS84 +ellps=WGS84 +lat_0=-90.0

    array(-2147483647, dtype=int32)

  • time_bnds
    (time, nv)
    datetime64[ns]
    ...

    array([['1979-01-30T00:00:00.000000000', '1979-01-31T00:00:00.000000000']],
          dtype='datetime64[ns]')

  • ice_conc
    (time, yc, xc)
    float64
    ...

    long_name :

    fully filtered concentration of sea ice using atmospheric correction of brightness temperatures and open water filters

    standard_name :

    sea_ice_area_fraction

    units :

    %

    valid_min :

    0

    valid_max :

    10000

    grid_mapping :

    Lambert_Azimuthal_Grid

    ancillary_variables :

    total_standard_error status_flag

    comment :

    this field is the primary sea ice concentration estimate for this climate data record

    [186624 values with dtype=float64]

  • raw_ice_conc_values
    (time, yc, xc)
    float64
    ...

    long_name :

    sea ice concentration estimates as retrieved by the algorithm, and that were edited away by the various filters

    units :

    %

    grid_mapping :

    Lambert_Azimuthal_Grid

    comment :

    this field can be used in combination with 'ice_conc' to access un-bounded normal distribution of sea ice concentration estimates corresponding to the uncertainty estimate in 'total_standard_error'. It does also feature some amount of un-physical data (below 0% and above 100%).

    [186624 values with dtype=float64]

  • total_standard_error
    (time, yc, xc)
    float64
    ...

    long_name :

    total uncertainty (one standard deviation) of concentration of sea ice

    standard_name :

    sea_ice_area_fraction standard_error

    units :

    %

    grid_mapping :

    Lambert_Azimuthal_Grid

    [186624 values with dtype=float64]

  • smearing_standard_error
    (time, yc, xc)
    float64
    ...

    long_name :

    smearing uncertainty (one standard deviation) of concentration of sea ice

    units :

    %

    grid_mapping :

    Lambert_Azimuthal_Grid

    comment :

    this is one of the two components contributing to 'total_standard_error'

    [186624 values with dtype=float64]

  • algorithm_standard_error
    (time, yc, xc)
    float64
    ...

    long_name :

    algorithm uncertainty (one standard deviation) of concentration of sea ice

    units :

    %

    grid_mapping :

    Lambert_Azimuthal_Grid

    comment :

    this is one of the two components contributing to 'total_standard_error'

    [186624 values with dtype=float64]

  • status_flag
    (time, yc, xc)
    float32
    ...

    long_name :

    status flag bit array for sea ice concentration retrievals

    standard_name :

    sea_ice_area_fraction status_flag

    valid_min :

    0

    valid_max :

    32767

    grid_mapping :

    Lambert_Azimuthal_Grid

    flag_masks :

    [ 1 2 4 8 16 32 64 128]

    flag_meanings :

    land lake open_water_filtered land_spill_over high_t2m spatial_interp temporal_interp max_ice_climo

    flag_descriptions :

    all bits to 0 (flag 0): Nominal retrieval by the SIC algorithm bit 1 (flag 1): Position is over land bit 2 (flag 2): Position is lake bit 3 (flag 4): SIC is set to zero by the open water filter bit 4 (flag 8): SIC value is changed for correcting land spill-over effects bit 5 (flag 16): Handle with caution, the 2m air temperature is high at this position, and this might be false ice bit 6 (flag 32): Value is the result of spatial interpolation bit 7 (flag 64): Value is the result of temporal interpolation bit 8 (flag 128): SIC is set to zero since position is outside maximum sea ice climatology.

    comment :

    Flag values found in the map might be combinations of those listed above. For example flag value 80 (= 2^4 + 2^6 = 16 + 64) indicates high_t2m AND temporal_interp

    [186624 values with dtype=float32]

• Attributes: (40)

  title :

  Sea Ice Concentration Climate Data Record from EUMETSAT OSI SAF

  summary :

  This climate data record of sea ice concentration is obtained from coarse resolution passive microwave satellite data over the polar regions (SMMR, SSM/I, and SSMIS). The processing chain features: 1) dynamic tuning of tie-points and algorithms, 2) correction of atmospheric noise using a Radiative Transfer Model, 3) computation of per-pixel uncertainties, and 4) an optimal hybrid sea ice concentration algorithm. This dataset was generated by the EUMETSAT Ocean and Sea Ice Satellite Application Facility (OSI SAF). The ESA CCI Programme contributed with Research and Development on the algorithms.

  topiccategory :

  Oceans ClimatologyMeteorologyAtmosphere

  keywords :

  Earth Science > Cryosphere > Sea Ice > Sea Ice Concentration Earth Science > Oceans > Sea Ice > Sea Ice Concentration Earth Science > Climate Indicators > Cryospheric Indicators > Sea Ice Concentration,Geographic Region > Southern Hemisphere,Vertical Location > Sea Surface,NO/MET > Norwegian Meteorological Institute

  keywords_vocabulary :

  GCMD Science Keywords

  geospatial_lat_min :

  -90.0

  geospatial_lat_max :

  -16.62393

  geospatial_lon_min :

  -180.0

  geospatial_lon_max :

  180.0

  geospatial_vertical_min :

  0.0

  geospatial_vertical_max :

  0.0

  sensor :

  SMMR,SSM/I,SSMIS

  platform :

  Nimbus-7,DMSP-F<08,10,11,13,14,15>,DMSP-F<16,17,18>

  source :

  FCDR of SMMR / SSMI / SSMIS Brightness Temperatures R3 (doi: 10.5676/EUM_SAF_CM/FCDR_MWI/V003),ERA-Interim fields from ECMWF (doi: 10.1002/qj.828)

  time_coverage_start :

  1979-01-30T00:00:00Z

  time_coverage_end :

  1979-01-31T00:00:00Z

  time_coverage_duration :

  P1D

  time_coverage_resolution :

  P1D

  project :

  OSI SAF - EUMETSAT

  institution :

  OSISAF - EUMETSAT

  creator_name :

  OSISAF - EUMETSAT

  creator_type :

  institution

  creator_url :

  http://www.osi-saf.org

  creator_email :

  osi-saf.helpdesk@meteo.fr

  license :

  All intellectual property rights of the Ocean and Sea Ice SAF products belong to EUMETSAT. The use of these products is granted to every user, free of charge. If users wish to use these products, EUMETSAT's copyright credit must be shown by displaying the words 'Copyright EUMETSAT' under each of the products shown. EUMETSAT offers no warranty and accepts no liability in respect of the Ocean and Sea Ice SAF products. EUMETSAT neither commits to nor guarantees the continuity, availability, or quality or suitability for any purpose of, the Ocean and Sea Ice SAF products.

  references :

  Product User Manual v1p0 (Dec 2016),Algorithm Theoretical Basis Document v1.2 (Dec 2016),http://osisaf.met.no,http://www.osi-saf.org

  date_created :

  2017-02-14

  cdm_data_type :

  Grid

  spatial_resolution :

  25.0 km grid spacing

  algorithm :

  SICCI2LF (19V, 37V, 37H)

  contributor_name :

  Leif Toudal Pedersen, Gorm Dybkjaer, Steinar Eastwood, Mari Anne Killie, John Lavelle, Thomas Lavergne, Helge Pfeiffer, Atle Soerensen, Rasmus Tonboe

  contributor_role :

  PrincipalInvestigator,author,author,author,author,author,author,author,author

  tracking_id :

  e4a4ec71-7c69-4521-8937-f430bb38469e

  naming_authority :

  int.eumetsat

  Conventions :

  CF-1.6 ACDD-1.3

  standard_name_vocabulary :

  CF Standard Name Table (v36, 21 September 2016)

  id :

  10.15770/EUM_SAF_OSI_0008

  product_version :

  2.0

  history :

  2017-02-14T19:39:25Z creation 2022-06-07: Corrected EPSG

  geospatial_bounds_crs :

  EPSG:6932