How to open GEOS-FP netcdf file through https url?

[zxdawn]

I'm trying to read the GEOS-FP wind data on this website. Here's one example:

import xarray as xr

url = 'https://portal.nccs.nasa.gov/datashare/gmao/geos-fp/das/Y2023/M10/D01/GEOS.fp.asm.tavg1_2d_slv_Nx.20231001_0030.V01.nc4'
xr.open_dataset(url)

I can open the url in browser and download it successfully. I also checked the downloading URL shown in the browser, that's the same url.

Any ideas are welcome :)

[jhamman]

Try using fsspec along with Xarray:

import fsspec
import xarray as xr

url = 'https://portal.nccs.nasa.gov/datashare/gmao/geos-fp/das/Y2023/M10/D01/GEOS.fp.asm.tavg1_2d_slv_Nx.20231001_0030.V01.nc4'
with fsspec.open(url).open() as f:
    ds = xr.open_dataset(f)
print(ds)

<xarray.Dataset>
Dimensions:   (lon: 1152, lat: 721, time: 1)
Coordinates:
  * lon       (lon) float64 -180.0 -179.7 -179.4 -179.1 ... 179.1 179.4 179.7
  * lat       (lat) float64 -90.0 -89.75 -89.5 -89.25 ... 89.25 89.5 89.75 90.0
  * time      (time) datetime64[ns] 2023-10-01T00:30:00
Data variables: (12/45)
    CLDPRS    (time, lat, lon) float32 ...
    CLDTMP    (time, lat, lon) float32 ...
    DISPH     (time, lat, lon) float32 ...
    H1000     (time, lat, lon) float32 ...
    H250      (time, lat, lon) float32 ...
    H500      (time, lat, lon) float32 ...
    ...        ...
    V250      (time, lat, lon) float32 ...
    V2M       (time, lat, lon) float32 ...
    V500      (time, lat, lon) float32 ...
    V50M      (time, lat, lon) float32 ...
    V850      (time, lat, lon) float32 ...
    TAITIME   (time) float64 ...
Attributes: (12/28)
    Contact:               http://gmao.gsfc.nasa.gov/
    Institution:           NASA Global Modeling and Assimilation Office
    History:               Original file generated: Sun Oct  1 09:57:03 2023 GMT
    Filename:              GEOS.fp.asm.tavg1_2d_slv_Nx.20231001_0030.V01.nc4
    Comment:               GMAO filename: f5295_fp.tavg1_2d_slv_Nx.20231001_0...
    Source:                GEOSadas-5_29_5_OPS experiment_id: f5295_fp
    ...                    ...
    NorthernmostLatitude:  90.0
    WesternmostLongitude:  -180.0
    EasternmostLongitude:  179.6875
    LatitudeResolution:    0.25
    LongitudeResolution:   0.3125
    DataResolution:        0.25 x 0.3125

[zxdawn]

Thanks, it works well!