This page was generated from nldas.ipynb. Interactive online version:

NLDAS2 Forcing

from pathlib import Path

import pynldas2 as nldas

from pygeohydro import WBD

The NLDAS2 database provides forcing data at 1/8th-degree grid spacing and range from 01 Jan 1979 to present. Let’s take a look at NLDAS2 grid mask that includes land, water, soil, and vegetation masks:

grid = nldas.get_grid_mask()
grid

<xarray.Dataset>
Dimensions:      (lon: 464, lat: 224, time: 1, bnds: 2)
Coordinates:
  * lon          (lon) float32 -124.9 -124.8 -124.7 ... -67.31 -67.19 -67.06
  * lat          (lat) float32 25.06 25.19 25.31 25.44 ... 52.69 52.81 52.94
  * time         (time) datetime64[ns] 2000-01-01
    spatial_ref  int64 0
Dimensions without coordinates: bnds
Data variables:
    time_bnds    (time, bnds) datetime64[ns] ...
    NLDAS_mask   (time, lat, lon) float32 ...
    CONUS_mask   (time, lat, lon) float32 ...
    NLDAS_veg    (time, lat, lon) float32 ...
    NLDAS_soil   (time, lat, lon) float32 ...
Attributes: (12/13)
    missing_value:          -9999.0
    time_definition:        constant
    title:                  NLDAS masks and predominant vegetation/soil
    institution:            NASA GSFC
    history:                created on date: Fri Mar  8 15:58:50 2019
    references:             Mitchell_etal_JGR_2004; Xia_etal_JGR_2012
    ...                     ...
    website:                https://ldas.gsfc.nasa.gov/nldas/
    MAP_PROJECTION:         EQUIDISTANT CYLINDRICAL
    SOUTH_WEST_CORNER_LAT:  25.0625
    SOUTH_WEST_CORNER_LON:  -124.9375
    DX:                     0.125
    DY:                     0.125

xarray.Dataset

• Dimensions:
  • lon: 464
  • lat: 224
  • time: 1
  • bnds: 2
• Coordinates: (4)
  • lon
    (lon)
    float32
    -124.9 -124.8 ... -67.19 -67.06

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    vmin :

    -124.9375

    vmax :

    -67.0625

    axis :

    X

    array([-124.9375, -124.8125, -124.6875, ...,  -67.3125,  -67.1875,  -67.0625],
          dtype=float32)

  • lat
    (lat)
    float32
    25.06 25.19 25.31 ... 52.81 52.94

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    vmin :

    25.0625

    vmax :

    52.9375

    axis :

    Y

    array([25.0625, 25.1875, 25.3125, ..., 52.6875, 52.8125, 52.9375],
          dtype=float32)

  • time
    (time)
    datetime64[ns]
    2000-01-01

    long_name :

    time

    time_increment :

    1

    bounds :

    time_bnds

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

  • spatial_ref
    ()
    int64
    0

    crs_wkt :

    GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4326"]]

    semi_major_axis :

    6378137.0

    semi_minor_axis :

    6356752.314245179

    inverse_flattening :

    298.257223563

    reference_ellipsoid_name :

    WGS 84

    longitude_of_prime_meridian :

    0.0

    prime_meridian_name :

    Greenwich

    geographic_crs_name :

    WGS 84

    grid_mapping_name :

    latitude_longitude

    spatial_ref :

    GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4326"]]

    GeoTransform :

    -125.0 0.125 0.0 25.0 0.0 0.125

    array(0)

• Data variables: (5)
  • time_bnds
    (time, bnds)
    datetime64[ns]
    ...

    [2 values with dtype=datetime64[ns]]

  • NLDAS_mask
    (time, lat, lon)
    float32
    ...

    units :

    0 = water; 1 = land

    standard_name :

    land mask

    long_name :

    NLDAS UMD unified land mask

    cell_methods :

    time: mean

    vmin :

    0.0

    vmax :

    1.0

    [103936 values with dtype=float32]

  • CONUS_mask
    (time, lat, lon)
    float32
    ...

    units :

    0 = water; 1 = land

    standard_name :

    CONUS mask

    long_name :

    NLDAS CONUS mask

    cell_methods :

    time: mean

    vmin :

    0.0

    vmax :

    1.0

    [103936 values with dtype=float32]

  • NLDAS_veg
    (time, lat, lon)
    float32
    ...

    units :

    0 - 13

    standard_name :

    vegetation class

    long_name :

    NLDAS UMD predominant vegetation class

    cell_methods :

    time: mean

    vmin :

    0.0

    vmax :

    13.0

    [103936 values with dtype=float32]

  • NLDAS_soil
    (time, lat, lon)
    float32
    ...

    units :

    0 - 14

    standard_name :

    soil texture type

    long_name :

    NLDAS STATSGO predominant soil texture type

    cell_methods :

    time: mean

    vmin :

    0.0

    vmax :

    14.0

    [103936 values with dtype=float32]

• Indexes: (3)
  • lon
    PandasIndex

    PandasIndex(Float64Index([-124.9375, -124.8125, -124.6875, -124.5625, -124.4375, -124.3125,
                  -124.1875, -124.0625, -123.9375, -123.8125,
                  ...
                   -68.1875,  -68.0625,  -67.9375,  -67.8125,  -67.6875,  -67.5625,
                   -67.4375,  -67.3125,  -67.1875,  -67.0625],
                 dtype='float64', name='lon', length=464))

  • lat
    PandasIndex

    PandasIndex(Float64Index([25.0625, 25.1875, 25.3125, 25.4375, 25.5625, 25.6875, 25.8125,
                  25.9375, 26.0625, 26.1875,
                  ...
                  51.8125, 51.9375, 52.0625, 52.1875, 52.3125, 52.4375, 52.5625,
                  52.6875, 52.8125, 52.9375],
                 dtype='float64', name='lat', length=224))

  • time
    PandasIndex

    PandasIndex(DatetimeIndex(['2000-01-01'], dtype='datetime64[ns]', name='time', freq=None))

• Attributes: (13)

  missing_value :

  -9999.0

  time_definition :

  constant

  title :

  NLDAS masks and predominant vegetation/soil

  institution :

  NASA GSFC

  history :

  created on date: Fri Mar 8 15:58:50 2019

  references :

  Mitchell_etal_JGR_2004; Xia_etal_JGR_2012

  conventions :

  CF-1.6

  website :

  https://ldas.gsfc.nasa.gov/nldas/

  MAP_PROJECTION :

  EQUIDISTANT CYLINDRICAL

  SOUTH_WEST_CORNER_LAT :

  25.0625

  SOUTH_WEST_CORNER_LON :

  -124.9375

  DX :

  0.125

  DY :

  0.125

For example, let’s plot the vegetation mask.

ax = grid.NLDAS_veg.plot()
ax.figure.savefig(Path("_static", "nldas_grid.png"), facecolor="w", bbox_inches="tight")

Next, we use PyGeoHydro to get the geometry of a HUC8 with ID of 1306003:

huc8 = WBD("huc8")
geometry = huc8.byids("huc8", "13060003").geometry[0]

PyNLDAS2 allows us to get the data for a list of coordinates using pynldas2.get_bycoords or for a region as gridded data using pynldas2.get_bygeom. Here, we use the latter. Note that if we don’t pass any variables, all variables will be downloaded.

clm = nldas.get_bygeom(geometry, "2010-01-01", "2010-01-31", 4326)

clm

<xarray.Dataset>
Dimensions:      (time: 744, y: 11, x: 15)
Coordinates:
  * time         (time) datetime64[ns] 2010-01-01 ... 2010-01-31T23:00:00
  * y            (y) float64 33.56 33.69 33.81 33.94 ... 34.44 34.56 34.69 34.81
  * x            (x) float64 -105.6 -105.4 -105.3 ... -104.1 -103.9 -103.8
    spatial_ref  int64 0
Data variables:
    prcp         (time, y, x) float64 nan nan nan nan nan ... 0.0 0.0 nan nan
    pet          (time, y, x) float64 nan nan nan nan ... 0.0143 0.0145 nan nan
    temp         (time, y, x) float64 nan nan nan nan ... 275.6 275.1 nan nan
    wind_u       (time, y, x) float64 nan nan nan nan nan ... 0.69 0.79 nan nan
    wind_v       (time, y, x) float64 nan nan nan nan nan ... 4.03 4.15 nan nan
    rlds         (time, y, x) float64 nan nan nan nan ... 241.8 239.1 nan nan
    rsds         (time, y, x) float64 nan nan nan nan ... 208.7 208.5 nan nan
    humidity     (time, y, x) float64 nan nan nan nan ... 0.004118 nan nan
Attributes:
    tz:       UTC

xarray.Dataset

• Dimensions:
  • time: 744
  • y: 11
  • x: 15
• Coordinates: (4)
  • time
    (time)
    datetime64[ns]
    2010-01-01 ... 2010-01-31T23:00:00

    array(['2010-01-01T00:00:00.000000000', '2010-01-01T01:00:00.000000000',
           '2010-01-01T02:00:00.000000000', ..., '2010-01-31T21:00:00.000000000',
           '2010-01-31T22:00:00.000000000', '2010-01-31T23:00:00.000000000'],
          dtype='datetime64[ns]')

  • y
    (y)
    float64
    33.56 33.69 33.81 ... 34.69 34.81

    axis :

    Y

    long_name :

    latitude

    standard_name :

    latitude

    units :

    degrees_north

    array([33.5625, 33.6875, 33.8125, 33.9375, 34.0625, 34.1875, 34.3125, 34.4375,
           34.5625, 34.6875, 34.8125])

  • x
    (x)
    float64
    -105.6 -105.4 ... -103.9 -103.8

    axis :

    X

    long_name :

    longitude

    standard_name :

    longitude

    units :

    degrees_east

    array([-105.5625, -105.4375, -105.3125, -105.1875, -105.0625, -104.9375,
           -104.8125, -104.6875, -104.5625, -104.4375, -104.3125, -104.1875,
           -104.0625, -103.9375, -103.8125])

  • spatial_ref
    ()
    int64
    0

    crs_wkt :

    GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4326"]]

    semi_major_axis :

    6378137.0

    semi_minor_axis :

    6356752.314245179

    inverse_flattening :

    298.257223563

    reference_ellipsoid_name :

    WGS 84

    longitude_of_prime_meridian :

    0.0

    prime_meridian_name :

    Greenwich

    geographic_crs_name :

    WGS 84

    grid_mapping_name :

    latitude_longitude

    spatial_ref :

    GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4326"]]

    GeoTransform :

    -105.625 0.125 0.0 33.5 0.0 0.125

    array(0)

• Data variables: (8)
  • prcp
    (time, y, x)
    float64
    nan nan nan nan ... 0.0 0.0 nan nan

    nldas_name :

    APCPsfc

    long_name :

    Precipitation hourly total

    units :

    mm

    array([[[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 0.    , 0.    ,    nan],
            ...,
            [   nan, 0.    , 0.    , ..., 0.    ,    nan,    nan],
            [   nan,    nan, 0.    , ..., 0.    ,    nan,    nan],
            [   nan,    nan,    nan, ..., 0.    ,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 0.    , 0.    ,    nan],
            ...,
            [   nan, 0.    , 0.    , ..., 0.    ,    nan,    nan],
            [   nan,    nan, 0.    , ..., 0.    ,    nan,    nan],
            [   nan,    nan,    nan, ..., 0.    ,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 0.    , 0.    ,    nan],
            ...,
    ...
            ...,
            [   nan, 0.    , 0.    , ..., 0.    ,    nan,    nan],
            [   nan,    nan, 0.    , ..., 0.    ,    nan,    nan],
            [   nan,    nan,    nan, ..., 0.    ,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 0.    , 0.    ,    nan],
            ...,
            [   nan, 0.    , 0.    , ..., 0.    ,    nan,    nan],
            [   nan,    nan, 0.    , ..., 0.    ,    nan,    nan],
            [   nan,    nan,    nan, ..., 0.    ,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 0.    , 0.    ,    nan],
            ...,
            [   nan, 0.    , 0.    , ..., 0.    ,    nan,    nan],
            [   nan,    nan, 0.    , ..., 0.    ,    nan,    nan],
            [   nan,    nan,    nan, ..., 0.    ,    nan,    nan]]])

  • pet
    (time, y, x)
    float64
    nan nan nan nan ... 0.0145 nan nan

    nldas_name :

    PEVAPsfc

    long_name :

    Potential evaporation

    units :

    mm

    array([[[       nan,        nan,        nan, ...,        nan,
                    nan,        nan],
            [       nan,        nan,        nan, ...,        nan,
                    nan,        nan],
            [       nan,        nan,        nan, ...,  1.546e-01,
              1.481e-01,        nan],
            ...,
            [       nan,  9.240e-02,  8.700e-02, ...,  1.086e-01,
                    nan,        nan],
            [       nan,        nan,  7.610e-02, ...,  1.044e-01,
                    nan,        nan],
            [       nan,        nan,        nan, ...,  9.760e-02,
                    nan,        nan]],
    
           [[       nan,        nan,        nan, ...,        nan,
                    nan,        nan],
            [       nan,        nan,        nan, ...,        nan,
                    nan,        nan],
            [       nan,        nan,        nan, ...,  3.700e-03,
              3.700e-03,        nan],
    ...
            [       nan,  1.582e-01,  1.492e-01, ...,  1.320e-02,
                    nan,        nan],
            [       nan,        nan,  1.711e-01, ...,  1.340e-02,
                    nan,        nan],
            [       nan,        nan,        nan, ...,  1.450e-02,
                    nan,        nan]],
    
           [[       nan,        nan,        nan, ...,        nan,
                    nan,        nan],
            [       nan,        nan,        nan, ...,        nan,
                    nan,        nan],
            [       nan,        nan,        nan, ...,  1.997e-01,
              1.799e-01,        nan],
            ...,
            [       nan,  1.582e-01,  1.492e-01, ...,  1.320e-02,
                    nan,        nan],
            [       nan,        nan,  1.711e-01, ...,  1.340e-02,
                    nan,        nan],
            [       nan,        nan,        nan, ...,  1.450e-02,
                    nan,        nan]]])

  • temp
    (time, y, x)
    float64
    nan nan nan nan ... 275.1 nan nan

    nldas_name :

    TMP2m

    long_name :

    2-m above ground temperature

    units :

    K

    array([[[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 276.53, 275.95,    nan],
            ...,
            [   nan, 270.57, 271.01, ..., 274.9 ,    nan,    nan],
            [   nan,    nan, 270.52, ..., 274.47,    nan,    nan],
            [   nan,    nan,    nan, ..., 273.71,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 275.98, 275.37,    nan],
            ...,
            [   nan, 269.83, 270.3 , ..., 274.01,    nan,    nan],
            [   nan,    nan, 269.77, ..., 273.55,    nan,    nan],
            [   nan,    nan,    nan, ..., 272.77,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 275.43, 274.8 ,    nan],
            ...,
    ...
            ...,
            [   nan, 278.83, 279.11, ..., 278.01,    nan,    nan],
            [   nan,    nan, 279.21, ..., 277.09,    nan,    nan],
            [   nan,    nan,    nan, ..., 275.85,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 283.03, 282.06,    nan],
            ...,
            [   nan, 278.27, 278.59, ..., 277.47,    nan,    nan],
            [   nan,    nan, 278.64, ..., 276.62,    nan,    nan],
            [   nan,    nan,    nan, ..., 275.48,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 283.02, 282.01,    nan],
            ...,
            [   nan, 277.72, 278.06, ..., 276.92,    nan,    nan],
            [   nan,    nan, 278.06, ..., 276.14,    nan,    nan],
            [   nan,    nan,    nan, ..., 275.11,    nan,    nan]]])

  • wind_u
    (time, y, x)
    float64
    nan nan nan nan ... 0.79 nan nan

    nldas_name :

    UGRD10m

    long_name :

    10-m above ground zonal wind

    units :

    m/s

    array([[[  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ..., -0.82, -0.82,   nan],
            ...,
            [  nan,  3.12,  3.05, ...,  0.27,   nan,   nan],
            [  nan,   nan,  3.3 , ...,  0.62,   nan,   nan],
            [  nan,   nan,   nan, ...,  1.01,   nan,   nan]],
    
           [[  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ..., -0.57, -0.56,   nan],
            ...,
            [  nan,  3.31,  3.29, ...,  0.47,   nan,   nan],
            [  nan,   nan,  3.5 , ...,  0.81,   nan,   nan],
            [  nan,   nan,   nan, ...,  1.4 ,   nan,   nan]],
    
           [[  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ..., -0.33, -0.31,   nan],
            ...,
    ...
            ...,
            [  nan,  6.83,  6.35, ...,  0.35,   nan,   nan],
            [  nan,   nan,  6.9 , ...,  0.46,   nan,   nan],
            [  nan,   nan,   nan, ...,  0.99,   nan,   nan]],
    
           [[  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,  0.44,  0.42,   nan],
            ...,
            [  nan,  5.8 ,  5.61, ...,  0.15,   nan,   nan],
            [  nan,   nan,  6.16, ...,  0.31,   nan,   nan],
            [  nan,   nan,   nan, ...,  0.89,   nan,   nan]],
    
           [[  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,  0.18,  0.16,   nan],
            ...,
            [  nan,  4.78,  4.87, ..., -0.05,   nan,   nan],
            [  nan,   nan,  5.42, ...,  0.16,   nan,   nan],
            [  nan,   nan,   nan, ...,  0.79,   nan,   nan]]])

  • wind_v
    (time, y, x)
    float64
    nan nan nan nan ... 4.15 nan nan

    nldas_name :

    VGRD10m

    long_name :

    10-m above ground meridional wind

    units :

    m/s

    array([[[  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ..., -2.93, -2.94,   nan],
            ...,
            [  nan, -1.71, -1.6 , ..., -1.82,   nan,   nan],
            [  nan,   nan, -1.33, ..., -1.61,   nan,   nan],
            [  nan,   nan,   nan, ..., -1.29,   nan,   nan]],
    
           [[  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ..., -1.43, -1.45,   nan],
            ...,
            [  nan, -1.43, -1.15, ..., -0.49,   nan,   nan],
            [  nan,   nan, -1.01, ..., -0.33,   nan,   nan],
            [  nan,   nan,   nan, ..., -0.08,   nan,   nan]],
    
           [[  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,  0.07,  0.04,   nan],
            ...,
    ...
            ...,
            [  nan,  1.08,  1.61, ...,  3.44,   nan,   nan],
            [  nan,   nan,  1.17, ...,  3.7 ,   nan,   nan],
            [  nan,   nan,   nan, ...,  4.31,   nan,   nan]],
    
           [[  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,  2.86,  2.9 ,   nan],
            ...,
            [  nan,  0.8 ,  1.36, ...,  3.45,   nan,   nan],
            [  nan,   nan,  0.99, ...,  3.69,   nan,   nan],
            [  nan,   nan,   nan, ...,  4.23,   nan,   nan]],
    
           [[  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,   nan,   nan,   nan],
            [  nan,   nan,   nan, ...,  2.54,  2.58,   nan],
            ...,
            [  nan,  0.52,  1.12, ...,  3.45,   nan,   nan],
            [  nan,   nan,  0.81, ...,  3.67,   nan,   nan],
            [  nan,   nan,   nan, ...,  4.15,   nan,   nan]]])

  • rlds
    (time, y, x)
    float64
    nan nan nan nan ... 239.1 nan nan

    nldas_name :

    DLWRFsfc

    long_name :

    Surface DW longwave radiation flux

    units :

    W/m^2

    array([[[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 208.65, 206.03,    nan],
            ...,
            [   nan, 176.63, 178.32, ..., 198.21,    nan,    nan],
            [   nan,    nan, 175.97, ..., 195.56,    nan,    nan],
            [   nan,    nan,    nan, ..., 191.74,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 208.65, 206.03,    nan],
            ...,
            [   nan, 176.62, 178.32, ..., 198.2 ,    nan,    nan],
            [   nan,    nan, 175.96, ..., 195.55,    nan,    nan],
            [   nan,    nan,    nan, ..., 191.73,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 208.64, 206.03,    nan],
            ...,
    ...
            ...,
            [   nan, 240.4 , 242.99, ..., 247.73,    nan,    nan],
            [   nan,    nan, 240.92, ..., 243.96,    nan,    nan],
            [   nan,    nan,    nan, ..., 239.08,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 277.78, 270.8 ,    nan],
            ...,
            [   nan, 240.4 , 242.99, ..., 247.73,    nan,    nan],
            [   nan,    nan, 240.92, ..., 243.96,    nan,    nan],
            [   nan,    nan,    nan, ..., 239.08,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 277.78, 270.8 ,    nan],
            ...,
            [   nan, 240.39, 242.99, ..., 247.72,    nan,    nan],
            [   nan,    nan, 240.92, ..., 243.95,    nan,    nan],
            [   nan,    nan,    nan, ..., 239.07,    nan,    nan]]])

  • rsds
    (time, y, x)
    float64
    nan nan nan nan ... 208.5 nan nan

    nldas_name :

    DSWRFsfc

    long_name :

    Surface DW shortwave radiation flux

    units :

    W/m^2

    array([[[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,   0.  ,   0.  ,    nan],
            ...,
            [   nan,   0.  ,   0.  , ...,   0.  ,    nan,    nan],
            [   nan,    nan,   0.  , ...,   0.  ,    nan,    nan],
            [   nan,    nan,    nan, ...,   0.  ,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,   0.  ,   0.  ,    nan],
            ...,
            [   nan,   0.  ,   0.  , ...,   0.  ,    nan,    nan],
            [   nan,    nan,   0.  , ...,   0.  ,    nan,    nan],
            [   nan,    nan,    nan, ...,   0.  ,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,   0.  ,   0.  ,    nan],
            ...,
    ...
            ...,
            [   nan, 495.9 , 488.08, ..., 474.7 ,    nan,    nan],
            [   nan,    nan, 487.34, ..., 477.7 ,    nan,    nan],
            [   nan,    nan,    nan, ..., 477.09,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 348.46, 349.09,    nan],
            ...,
            [   nan, 362.16, 357.74, ..., 341.92,    nan,    nan],
            [   nan,    nan, 362.  , ..., 337.31,    nan,    nan],
            [   nan,    nan,    nan, ..., 341.71,    nan,    nan]],
    
           [[   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ...,    nan,    nan,    nan],
            [   nan,    nan,    nan, ..., 214.72, 214.48,    nan],
            ...,
            [   nan, 228.54, 225.13, ..., 209.14,    nan,    nan],
            [   nan,    nan, 227.55, ..., 206.06,    nan,    nan],
            [   nan,    nan,    nan, ..., 208.49,    nan,    nan]]])

  • humidity
    (time, y, x)
    float64
    nan nan nan ... 0.004118 nan nan

    nldas_name :

    SPFH2m

    long_name :

    2-m above ground specific humidity

    units :

    kg/kg

    array([[[      nan,       nan,       nan, ...,       nan,       nan,
                   nan],
            [      nan,       nan,       nan, ...,       nan,       nan,
                   nan],
            [      nan,       nan,       nan, ..., 0.003091 , 0.003105 ,
                   nan],
            ...,
            [      nan, 0.002713 , 0.0027498, ..., 0.0028554,       nan,
                   nan],
            [      nan,       nan, 0.0027574, ..., 0.0027586,       nan,
                   nan],
            [      nan,       nan,       nan, ..., 0.0026542,       nan,
                   nan]],
    
           [[      nan,       nan,       nan, ...,       nan,       nan,
                   nan],
            [      nan,       nan,       nan, ...,       nan,       nan,
                   nan],
            [      nan,       nan,       nan, ..., 0.0030743, 0.0030971,
                   nan],
    ...
            [      nan, 0.0044181, 0.0043933, ..., 0.0040693,       nan,
                   nan],
            [      nan,       nan, 0.0043753, ..., 0.0040481,       nan,
                   nan],
            [      nan,       nan,       nan, ..., 0.0040677,       nan,
                   nan]],
    
           [[      nan,       nan,       nan, ...,       nan,       nan,
                   nan],
            [      nan,       nan,       nan, ...,       nan,       nan,
                   nan],
            [      nan,       nan,       nan, ..., 0.004287 , 0.0042138,
                   nan],
            ...,
            [      nan, 0.004371 , 0.0043666, ..., 0.004127 ,       nan,
                   nan],
            [      nan,       nan, 0.0043278, ..., 0.0041002,       nan,
                   nan],
            [      nan,       nan,       nan, ..., 0.0041182,       nan,
                   nan]]])

• Indexes: (3)
  • time
    PandasIndex

    PandasIndex(DatetimeIndex(['2010-01-01 00:00:00', '2010-01-01 01:00:00',
                   '2010-01-01 02:00:00', '2010-01-01 03:00:00',
                   '2010-01-01 04:00:00', '2010-01-01 05:00:00',
                   '2010-01-01 06:00:00', '2010-01-01 07:00:00',
                   '2010-01-01 08:00:00', '2010-01-01 09:00:00',
                   ...
                   '2010-01-31 14:00:00', '2010-01-31 15:00:00',
                   '2010-01-31 16:00:00', '2010-01-31 17:00:00',
                   '2010-01-31 18:00:00', '2010-01-31 19:00:00',
                   '2010-01-31 20:00:00', '2010-01-31 21:00:00',
                   '2010-01-31 22:00:00', '2010-01-31 23:00:00'],
                  dtype='datetime64[ns]', name='time', length=744, freq=None))

  • y
    PandasIndex

    PandasIndex(Float64Index([33.5625, 33.6875, 33.8125, 33.9375, 34.0625, 34.1875, 34.3125,
                  34.4375, 34.5625, 34.6875, 34.8125],
                 dtype='float64', name='y'))

  • x
    PandasIndex

    PandasIndex(Float64Index([-105.5625, -105.4375, -105.3125, -105.1875, -105.0625, -104.9375,
                  -104.8125, -104.6875, -104.5625, -104.4375, -104.3125, -104.1875,
                  -104.0625, -103.9375, -103.8125],
                 dtype='float64', name='x'))

• Attributes: (1)

  tz :

  UTC

ax = clm.humidity.isel(time=slice(10, 16)).plot(col="time", col_wrap=3)
ax.fig.savefig(Path("_static", "nldas_humidity.png"), facecolor="w", bbox_inches="tight")