Lateral fill on model grid

This example shows how to apply lateral_fill() on a model dataset. We use a model grid dataset for illustration purposes.

Import packages

%matplotlib inline
import numpy as np
import xarray as xr

import pop_tools

ds = pop_tools.get_grid('POP_gx3v7')
ds

Downloading file 'inputdata/ocn/pop/gx3v7/grid/horiz_grid_20030806.ieeer8' from 'https://svn-ccsm-inputdata.cgd.ucar.edu/trunk/inputdata/ocn/pop/gx3v7/grid/horiz_grid_20030806.ieeer8' to '/home/docs/.pop_tools'.

/home/docs/checkouts/readthedocs.org/user_builds/pop-tools/conda/latest/lib/python3.9/site-packages/urllib3/connectionpool.py:1095: InsecureRequestWarning: Unverified HTTPS request is being made to host 'svn-ccsm-inputdata.cgd.ucar.edu'. Adding certificate verification is strongly advised. See: https://urllib3.readthedocs.io/en/latest/advanced-usage.html#tls-warnings
  warnings.warn(

  0%|                                               | 0.00/650k [00:00<?, ?B/s]

 10%|███▉                                   | 66.6k/650k [00:00<00:01, 547kB/s]

 88%|██████████████████████████████████▍    | 574k/650k [00:00<00:00, 2.66MB/s]

  0%|                                               | 0.00/650k [00:00<?, ?B/s]

100%|███████████████████████████████████████| 650k/650k [00:00<00:00, 1.01GB/s]

Downloading file 'inputdata/ocn/pop/gx3v7/grid/topography_20100105.ieeei4' from 'https://svn-ccsm-inputdata.cgd.ucar.edu/trunk/inputdata/ocn/pop/gx3v7/grid/topography_20100105.ieeei4' to '/home/docs/.pop_tools'.

/home/docs/checkouts/readthedocs.org/user_builds/pop-tools/conda/latest/lib/python3.9/site-packages/urllib3/connectionpool.py:1095: InsecureRequestWarning: Unverified HTTPS request is being made to host 'svn-ccsm-inputdata.cgd.ucar.edu'. Adding certificate verification is strongly advised. See: https://urllib3.readthedocs.io/en/latest/advanced-usage.html#tls-warnings
  warnings.warn(

  0%|                                              | 0.00/46.4k [00:00<?, ?B/s]

  0%|                                              | 0.00/46.4k [00:00<?, ?B/s]

100%|█████████████████████████████████████| 46.4k/46.4k [00:00<00:00, 54.8MB/s]

Downloading file 'inputdata/ocn/pop/gx3v7/grid/region_mask_20090831.ieeei4' from 'https://svn-ccsm-inputdata.cgd.ucar.edu/trunk/inputdata/ocn/pop/gx3v7/grid/region_mask_20090831.ieeei4' to '/home/docs/.pop_tools'.

/home/docs/checkouts/readthedocs.org/user_builds/pop-tools/conda/latest/lib/python3.9/site-packages/urllib3/connectionpool.py:1095: InsecureRequestWarning: Unverified HTTPS request is being made to host 'svn-ccsm-inputdata.cgd.ucar.edu'. Adding certificate verification is strongly advised. See: https://urllib3.readthedocs.io/en/latest/advanced-usage.html#tls-warnings
  warnings.warn(

  0%|                                              | 0.00/46.4k [00:00<?, ?B/s]

  0%|                                              | 0.00/46.4k [00:00<?, ?B/s]

100%|█████████████████████████████████████| 46.4k/46.4k [00:00<00:00, 39.0MB/s]

<xarray.Dataset>
Dimensions:      (nlat: 116, nlon: 100, z_t: 60, z_w: 60, z_w_bot: 60)
Coordinates:
  * z_t          (z_t) float64 500.0 1.5e+03 2.5e+03 ... 5.125e+05 5.375e+05
  * z_w          (z_w) float64 0.0 1e+03 2e+03 3e+03 ... 4.75e+05 5e+05 5.25e+05
  * z_w_bot      (z_w_bot) float64 1e+03 2e+03 3e+03 ... 5e+05 5.25e+05 5.5e+05
Dimensions without coordinates: nlat, nlon
Data variables: (12/13)
    TLAT         (nlat, nlon) float64 -79.5 -79.5 -79.5 ... 68.36 68.26 68.21
    TLONG        (nlat, nlon) float64 323.3 326.9 330.5 ... 317.8 319.3 320.8
    ULAT         (nlat, nlon) float64 -78.73 -78.73 -78.73 ... 68.49 68.42 68.4
    ULONG        (nlat, nlon) float64 325.1 328.7 332.3 ... 318.7 320.1 321.5
    DXT          (nlat, nlon) float64 6.842e+06 6.828e+06 ... 6.129e+06
    DYT          (nlat, nlon) float64 1.686e+07 1.686e+07 ... 4.524e+06
    ...           ...
    DYU          (nlat, nlon) float64 1.686e+07 1.686e+07 ... 1.069e+07
    TAREA        (nlat, nlon) float64 1.154e+14 1.151e+14 ... 2.773e+13
    UAREA        (nlat, nlon) float64 1.318e+14 1.318e+14 ... 6.256e+13 6.27e+13
    KMT          (nlat, nlon) int32 0 0 0 0 0 0 0 0 0 0 ... 0 0 0 0 0 0 0 0 0 0
    REGION_MASK  (nlat, nlon) int32 0 0 0 0 0 0 0 0 0 0 ... 0 0 0 0 0 0 0 0 0 0
    dz           (z_t) float64 1e+03 1e+03 1e+03 ... 2.499e+04 2.5e+04 2.5e+04
Attributes:
    lateral_dims:       [116, 100]
    vertical_dims:      60
    vert_grid_file:     gx3v7_vert_grid
    horiz_grid_fname:   inputdata/ocn/pop/gx3v7/grid/horiz_grid_20030806.ieeer8
    topography_fname:   inputdata/ocn/pop/gx3v7/grid/topography_20100105.ieeei4
    region_mask_fname:  inputdata/ocn/pop/gx3v7/grid/region_mask_20090831.ieeei4
    type:               dipole
    title:              POP_gx3v7 grid

xarray.Dataset

• Dimensions:
  • nlat: 116
  • nlon: 100
  • z_t: 60
  • z_w: 60
  • z_w_bot: 60
• Coordinates: (3)
  • z_t
    (z_t)
    float64
    500.0 1.5e+03 ... 5.375e+05

    units :

    cm

    long_name :

    depth from surface to midpoint of layer

    positive :

    down

    array([5.000000e+02, 1.500000e+03, 2.500000e+03, 3.500000e+03, 4.500000e+03,
           5.500000e+03, 6.500000e+03, 7.500000e+03, 8.500000e+03, 9.500000e+03,
           1.050000e+04, 1.150000e+04, 1.250000e+04, 1.350000e+04, 1.450000e+04,
           1.550000e+04, 1.650984e+04, 1.754790e+04, 1.862913e+04, 1.976603e+04,
           2.097114e+04, 2.225783e+04, 2.364088e+04, 2.513702e+04, 2.676542e+04,
           2.854837e+04, 3.051192e+04, 3.268680e+04, 3.510935e+04, 3.782276e+04,
           4.087847e+04, 4.433777e+04, 4.827367e+04, 5.277280e+04, 5.793729e+04,
           6.388626e+04, 7.075633e+04, 7.870025e+04, 8.788252e+04, 9.847059e+04,
           1.106204e+05, 1.244567e+05, 1.400497e+05, 1.573946e+05, 1.764003e+05,
           1.968944e+05, 2.186457e+05, 2.413972e+05, 2.649001e+05, 2.889385e+05,
           3.133405e+05, 3.379794e+05, 3.627670e+05, 3.876452e+05, 4.125768e+05,
           4.375393e+05, 4.625190e+05, 4.875083e+05, 5.125028e+05, 5.375000e+05])

  • z_w
    (z_w)
    float64
    0.0 1e+03 2e+03 ... 5e+05 5.25e+05

    units :

    cm

    positive :

    down

    long_name :

    depth from surface to top of layer

    array([     0.    ,   1000.    ,   2000.    ,   3000.    ,   4000.    ,
             5000.    ,   6000.    ,   7000.    ,   8000.    ,   9000.    ,
            10000.    ,  11000.    ,  12000.    ,  13000.    ,  14000.    ,
            15000.    ,  16000.    ,  17019.6808,  18076.1292,  19182.1243,
            20349.9313,  21592.3446,  22923.3124,  24358.4534,  25915.5793,
            27615.2589,  29481.4713,  31542.3736,  33831.2257,  36387.4728,
            39258.0478,  42498.885 ,  46176.6575,  50370.6883,  55174.9119,
            60699.6663,  67072.8582,  74439.803 ,  82960.6956,  92804.3538,
           104136.8196, 117104.0188, 131809.3626, 148290.0716, 166499.2064,
           186301.4408, 207487.3978, 229803.9076, 252990.4017, 276809.8509,
           301067.0677, 325613.847 , 350344.8607, 375189.1888, 400101.1634,
           425052.4544, 450026.0482, 475012.0091, 500004.6829, 525000.927 ])

  • z_w_bot
    (z_w_bot)
    float64
    1e+03 2e+03 ... 5.25e+05 5.5e+05

    units :

    cm

    positive :

    down

    long_name :

    depth from surface to bottom of layer

    array([  1000.    ,   2000.    ,   3000.    ,   4000.    ,   5000.    ,
             6000.    ,   7000.    ,   8000.    ,   9000.    ,  10000.    ,
            11000.    ,  12000.    ,  13000.    ,  14000.    ,  15000.    ,
            16000.    ,  17019.6808,  18076.1292,  19182.1243,  20349.9313,
            21592.3446,  22923.3124,  24358.4534,  25915.5793,  27615.2589,
            29481.4713,  31542.3736,  33831.2257,  36387.4728,  39258.0478,
            42498.885 ,  46176.6575,  50370.6883,  55174.9119,  60699.6663,
            67072.8582,  74439.803 ,  82960.6956,  92804.3538, 104136.8196,
           117104.0188, 131809.3626, 148290.0716, 166499.2064, 186301.4408,
           207487.3978, 229803.9076, 252990.4017, 276809.8509, 301067.0677,
           325613.847 , 350344.8607, 375189.1888, 400101.1634, 425052.4544,
           450026.0482, 475012.0091, 500004.6829, 525000.927 , 549999.0364])

• Data variables: (13)
  • TLAT
    (nlat, nlon)
    float64
    -79.5 -79.5 -79.5 ... 68.26 68.21

    units :

    degrees_north

    long_name :

    T-grid latitude

    array([[-79.49699172, -79.49699172, -79.49699172, ..., -79.49699172,
            -79.49699172, -79.49699172],
           [-77.98148831, -77.98148831, -77.98148831, ..., -77.98148831,
            -77.98148831, -77.98148831],
           [-76.4659849 , -76.4659849 , -76.4659849 , ..., -76.4659849 ,
            -76.4659849 , -76.4659849 ],
           ...,
           [ 67.40867121,  67.46350911,  67.57190569, ...,  67.57190569,
             67.46350911,  67.40867121],
           [ 67.80494094,  67.85778431,  67.96217313, ...,  67.96217313,
             67.85778431,  67.80494094],
           [ 68.20741553,  68.25795363,  68.35772038, ...,  68.35772038,
             68.25795363,  68.20741553]])

  • TLONG
    (nlat, nlon)
    float64
    323.3 326.9 330.5 ... 319.3 320.8

    units :

    degrees_east

    long_name :

    T-grid longitude

    array([[323.3       , 326.9       , 330.5       , ..., 312.5       ,
            316.1       , 319.7       ],
           [323.3       , 326.9       , 330.5       , ..., 312.5       ,
            316.1       , 319.7       ],
           [323.3       , 326.9       , 330.5       , ..., 312.5       ,
            316.1       , 319.7       ],
           ...,
           [322.33511948, 323.99745622, 325.6363243 , ..., 317.3636757 ,
            319.00254378, 320.66488052],
           [322.28951224, 323.86045218, 325.40739731, ..., 317.59260269,
            319.13954782, 320.71048776],
           [322.24063519, 323.71368525, 325.16235526, ..., 317.83764474,
            319.28631475, 320.75936481]])

  • ULAT
    (nlat, nlon)
    float64
    -78.73 -78.73 -78.73 ... 68.42 68.4

    units :

    degrees_north

    long_name :

    U-grid latitude

    array([[-78.73381714, -78.73381714, -78.73381714, ..., -78.73381714,
            -78.73381714, -78.73381714],
           [-77.21763809, -77.21763809, -77.21763809, ..., -77.21763809,
            -77.21763809, -77.21763809],
           [-75.70146057, -75.70146057, -75.70146057, ..., -75.70146057,
            -75.70146057, -75.70146057],
           ...,
           [ 67.61685002,  67.69760248,  67.83004184, ...,  67.69760248,
             67.61685002,  67.58971351],
           [ 68.01581668,  68.09335013,  68.22039887, ...,  68.09335013,
             68.01581668,  67.98975038],
           [ 68.42115055,  68.49497   ,  68.61581693, ...,  68.49497   ,
             68.42115055,  68.39632083]])

  • ULONG
    (nlat, nlon)
    float64
    325.1 328.7 332.3 ... 320.1 321.5

    units :

    degrees_east

    long_name :

    U-grid longitude

    array([[325.1       , 328.7       , 332.3       , ..., 314.3       ,
            317.9       , 321.5       ],
           [325.1       , 328.7       , 332.3       , ..., 314.3       ,
            317.9       , 321.5       ],
           [325.1       , 328.7       , 332.3       , ..., 314.3       ,
            317.9       , 321.5       ],
           ...,
           [323.12673427, 324.73742057, 326.31651111, ..., 318.26257943,
            319.87326573, 321.5       ],
           [323.03230674, 324.54824165, 326.03197642, ..., 318.45175835,
            319.96769326, 321.5       ],
           [322.93098401, 324.34538104, 325.72719339, ..., 318.65461896,
            320.06901599, 321.5       ]])

  • DXT
    (nlat, nlon)
    float64
    6.842e+06 6.828e+06 ... 6.129e+06

    units :

    cm

    long_name :

    x-spacing centered at T points

    coordinates :

    TLONG TLAT

    array([[6842420.99516022, 6827933.89860412, 6799440.21095021, ...,
            6799440.21092715, 6827933.89860412, 6842420.99513141],
           [8337574.41211318, 8337574.41210233, 8337574.41212115, ...,
            8337574.41207265, 8337574.41210488, 8337574.41205383],
           [9370383.38211568, 9370383.38210603, 9370383.38212533, ...,
            9370383.38207025, 9370383.38210858, 9370383.38204867],
           ...,
           [7149775.8730794 , 7118056.5288889 , 7055576.66663754, ...,
            7055576.66663754, 7118056.5288889 , 7149775.8730794 ],
           [6647264.25261398, 6616506.47992619, 6555955.11675041, ...,
            6555955.11675041, 6616506.47992619, 6647264.25261398],
           [6128684.48251535, 6099081.43768296, 6040838.58796237, ...,
            6040838.58796237, 6099081.43768296, 6128684.48251535]])

  • DYT
    (nlat, nlon)
    float64
    1.686e+07 1.686e+07 ... 4.524e+06

    units :

    cm

    long_name :

    y-spacing centered at T points

    coordinates :

    TLONG TLAT

    array([[16859734.13037167, 16859734.13085433, 16859734.12850002, ...,
            16859734.14009841, 16859734.13389394, 16859734.12882792],
           [16859723.94967546, 16859723.95149922, 16859723.95129658, ...,
            16859723.96102332, 16859723.95575467, 16859723.94967546],
           [16859706.97653924, 16859706.98059205, 16859706.98464485, ...,
            16859706.98930557, 16859706.98545541, 16859706.97957885],
           ...,
           [ 4384323.09541593,  4396093.3643232 ,  4419279.39434137, ...,
             4419279.39434137,  4396093.3643232 ,  4384323.09541593],
           [ 4451249.83211331,  4462715.34484552,  4485288.30535927, ...,
             4485288.30535927,  4462715.34484552,  4451249.83211331],
           [ 4523811.70176776,  4534913.64449504,  4556757.67311923, ...,
             4556757.67311923,  4534913.64449504,  4523811.70176776]])

  • DXU
    (nlat, nlon)
    float64
    7.82e+06 7.82e+06 ... 5.865e+06

    units :

    cm

    long_name :

    x-spacing centered at U points

    coordinates :

    ULONG ULAT

    array([[7819620.90819425, 7819620.90819713, 7819620.90819425, ...,
            7819620.90817408, 7819620.90816543, 7819620.9081712 ],
           [8855527.91602127, 8855527.91602635, 8855527.91602381, ...,
            8855527.91600345, 8855527.91599328, 8855527.91599582],
           [9885238.84820044, 9885238.848205  , 9885238.84820044, ...,
            9885238.84817537, 9885238.84816397, 9885238.84816853],
           ...,
           [6886738.79791196, 6840294.77238846, 6764456.77880701, ...,
            6840294.77238846, 6886738.79791196, 6902380.62231771],
           [6377031.93462821, 6332166.82428814, 6258963.18566492, ...,
            6332166.82428814, 6377031.93462821, 6392147.88291026],
           [5850733.98557009, 5807753.2013572 , 5737681.40530693, ...,
            5807753.2013572 , 5850733.98557009, 5865221.08212043]])

  • DYU
    (nlat, nlon)
    float64
    1.686e+07 1.686e+07 ... 1.069e+07

    units :

    cm

    long_name :

    y-spacing centered at U points

    coordinates :

    ULONG ULAT

    array([[16859729.04083413, 16859729.04151942, 16859729.03827718, ...,
            16859729.05035823, 16859729.03929038, 16859729.039213  ],
           [16859715.46249943, 16859715.46959184, 16859715.46634959, ...,
            16859715.47567105, 16859715.46553903, 16859715.46371527],
           [16859693.40165876, 16859693.4123987 , 16859693.41016965, ...,
            16859693.41786999, 16859693.40672477, 16859693.40429308],
           ...,
           [ 4420713.64365424,  4438095.06551448,  4466472.63418616, ...,
             4438095.06551448,  4420713.64365424,  4414859.283875  ],
           [ 4490374.59155667,  4507254.39778389,  4534791.58069461, ...,
             4507254.39778389,  4490374.59155667,  4484686.9423244 ],
           [10693172.13884091, 10701475.63650846, 10715016.1651108 , ...,
            10701475.64109181, 10693172.13729717, 10690373.69329851]])

  • TAREA
    (nlat, nlon)
    float64
    1.154e+14 1.151e+14 ... 2.773e+13

    units :

    cm^2

    long_name :

    area of T cells

    coordinates :

    TLONG TLAT

    array([[1.15361399e+14, 1.15117150e+14, 1.14636754e+14, ...,
            1.14636754e+14, 1.15117150e+14, 1.15361399e+14],
           [1.40569203e+14, 1.40569203e+14, 1.40569203e+14, ...,
            1.40569203e+14, 1.40569203e+14, 1.40569203e+14],
           [1.57981918e+14, 1.57981918e+14, 1.57981918e+14, ...,
            1.57981918e+14, 1.57981918e+14, 1.57981918e+14],
           ...,
           [3.13469275e+13, 3.12916411e+13, 3.11805646e+13, ...,
            3.11805646e+13, 3.12916411e+13, 3.13469275e+13],
           [2.95886339e+13, 2.95275850e+13, 2.94053488e+13, ...,
            2.94053488e+13, 2.95275850e+13, 2.95886339e+13],
           [2.77250146e+13, 2.76588076e+13, 2.75266376e+13, ...,
            2.75266376e+13, 2.76588076e+13, 2.77250146e+13]])

  • UAREA
    (nlat, nlon)
    float64
    1.318e+14 1.318e+14 ... 6.27e+13

    units :

    cm^2

    long_name :

    area of U cells

    coordinates :

    ULONG ULAT

    array([[1.31836690e+14, 1.31836690e+14, 1.31836690e+14, ...,
            1.31836690e+14, 1.31836690e+14, 1.31836690e+14],
           [1.49301681e+14, 1.49301681e+14, 1.49301681e+14, ...,
            1.49301681e+14, 1.49301681e+14, 1.49301681e+14],
           [1.66662096e+14, 1.66662096e+14, 1.66662096e+14, ...,
            1.66662096e+14, 1.66662096e+14, 1.66662096e+14],
           ...,
           [3.04443002e+13, 3.03578785e+13, 3.02132611e+13, ...,
            3.03578785e+13, 3.04443002e+13, 3.04730392e+13],
           [2.86352622e+13, 2.85406868e+13, 2.83830936e+13, ...,
            2.85406868e+13, 2.86352622e+13, 2.86667821e+13],
           [6.25629056e+13, 6.21515294e+13, 6.14793490e+13, ...,
            6.21515294e+13, 6.25629056e+13, 6.27014052e+13]])

  • KMT
    (nlat, nlon)
    int32
    0 0 0 0 0 0 0 0 ... 0 0 0 0 0 0 0 0

    long_name :

    k Index of Deepest Grid Cell on T Grid

    coordinates :

    TLONG TLAT

    array([[ 0,  0,  0, ...,  0,  0,  0],
           [34,  0,  0, ..., 36, 36, 36],
           [35, 36, 33, ..., 36, 36, 36],
           ...,
           [ 0,  0,  0, ...,  0,  0,  0],
           [ 0,  0,  0, ...,  0,  0,  0],
           [ 0,  0,  0, ...,  0,  0,  0]], dtype=int32)

  • REGION_MASK
    (nlat, nlon)
    int32
    0 0 0 0 0 0 0 0 ... 0 0 0 0 0 0 0 0

    long_name :

    basin index number (signed integers)

    coordinates :

    TLONG TLAT

    array([[0, 0, 0, ..., 0, 0, 0],
           [1, 0, 0, ..., 1, 1, 1],
           [1, 1, 1, ..., 1, 1, 1],
           ...,
           [0, 0, 0, ..., 0, 0, 0],
           [0, 0, 0, ..., 0, 0, 0],
           [0, 0, 0, ..., 0, 0, 0]], dtype=int32)

  • dz
    (z_t)
    float64
    1e+03 1e+03 ... 2.5e+04 2.5e+04

    units :

    cm

    long_name :

    thickness of layer k

    array([ 1000.    ,  1000.    ,  1000.    ,  1000.    ,  1000.    ,
            1000.    ,  1000.    ,  1000.    ,  1000.    ,  1000.    ,
            1000.    ,  1000.    ,  1000.    ,  1000.    ,  1000.    ,
            1000.    ,  1019.6808,  1056.4484,  1105.9951,  1167.807 ,
            1242.4133,  1330.9678,  1435.141 ,  1557.1259,  1699.6796,
            1866.2124,  2060.9023,  2288.8521,  2556.2471,  2870.575 ,
            3240.8372,  3677.7725,  4194.0308,  4804.2236,  5524.7544,
            6373.1919,  7366.9448,  8520.8926,  9843.6582, 11332.4658,
           12967.1992, 14705.3438, 16480.709 , 18209.1348, 19802.2344,
           21185.957 , 22316.5098, 23186.4941, 23819.4492, 24257.2168,
           24546.7793, 24731.0137, 24844.3281, 24911.9746, 24951.291 ,
           24973.5938, 24985.9609, 24992.6738, 24996.2441, 24998.1094])

• Indexes: (3)
  • z_t
    PandasIndex

    PandasIndex(Index([             500.0,             1500.0,             2500.0,
                       3500.0,             4500.0,             5500.0,
                       6500.0,             7500.0,             8500.0,
                       9500.0,            10500.0,            11500.0,
                      12500.0,            13500.0,            14500.0,
                      15500.0,         16509.8404,          17547.905,
                  18629.12675,         19766.0278,        20971.13795,
                   22257.8285,         23640.8829, 25137.016349999998,
                   26765.4191, 28548.365099999995,        30511.92245,
                  32686.79965,  35109.34924999999, 37822.760299999994,
            40878.46639999999,        44337.77125,  48273.67289999999,
            52772.80009999999, 57937.289099999995,  63886.26224999999,
            70756.33059999999,  78700.24929999998,  87882.52469999998,
            98470.58669999999, 110620.41919999999, 124456.69069999999,
           140049.71709999998,         157394.639,        176400.3236,
           196894.41929999998, 218645.65269999998, 241397.15464999998,
                  264900.1263,        288938.4593,       313340.45735,
           337979.35384999996,       362767.02475,  387645.1760999999,
           412576.80889999995,        437539.2513, 462519.02864999993,
           487508.34599999996, 512502.80494999996,  537499.9816999999],
          dtype='float64', name='z_t'))

  • z_w
    PandasIndex

    PandasIndex(Index([               0.0,             1000.0,             2000.0,
                       3000.0,             4000.0,             5000.0,
                       6000.0,             7000.0,             8000.0,
                       9000.0,            10000.0,            11000.0,
                      12000.0,            13000.0,            14000.0,
                      15000.0,            16000.0,         17019.6808,
                   18076.1292,         19182.1243,         20349.9313,
                   21592.3446,         22923.3124,         24358.4534,
           25915.579299999998, 27615.258899999997, 29481.471299999997,
                   31542.3736, 33831.225699999995, 36387.472799999996,
            39258.04779999999, 42498.884999999995, 46176.657499999994,
           50370.688299999994,  55174.91189999999,  60699.66629999999,
            67072.85819999999,  74439.80299999999,  82960.69559999998,
            92804.35379999998, 104136.81959999999, 117104.01879999999,
                  131809.3626,        148290.0716,        166499.2064,
           186301.44079999998, 207487.39779999998, 229803.90759999998,
                  252990.4017, 276809.85089999996, 301067.06769999996,
           325613.84699999995, 350344.86069999996, 375189.18879999995,
           400101.16339999996,        425052.4544, 450026.04819999996,
           475012.00909999997, 500004.68289999996,  525000.9269999999],
          dtype='float64', name='z_w'))

  • z_w_bot
    PandasIndex

    PandasIndex(Index([            1000.0,             2000.0,             3000.0,
                       4000.0,             5000.0,             6000.0,
                       7000.0,             8000.0,             9000.0,
                      10000.0,            11000.0,            12000.0,
                      13000.0,            14000.0,            15000.0,
                      16000.0,         17019.6808,         18076.1292,
                   19182.1243,         20349.9313,         21592.3446,
                   22923.3124,         24358.4534, 25915.579299999998,
           27615.258899999997, 29481.471299999997,         31542.3736,
           33831.225699999995, 36387.472799999996,  39258.04779999999,
           42498.884999999995, 46176.657499999994, 50370.688299999994,
            55174.91189999999,  60699.66629999999,  67072.85819999999,
            74439.80299999999,  82960.69559999998,  92804.35379999998,
           104136.81959999999, 117104.01879999999,        131809.3626,
                  148290.0716,        166499.2064, 186301.44079999998,
           207487.39779999998, 229803.90759999998,        252990.4017,
           276809.85089999996, 301067.06769999996, 325613.84699999995,
           350344.86069999996, 375189.18879999995, 400101.16339999996,
                  425052.4544, 450026.04819999996, 475012.00909999997,
           500004.68289999996,  525000.9269999999,  549999.0363999999],
          dtype='float64', name='z_w_bot'))

• Attributes: (8)

  lateral_dims :

  [116, 100]

  vertical_dims :

  60

  vert_grid_file :

  gx3v7_vert_grid

  horiz_grid_fname :

  inputdata/ocn/pop/gx3v7/grid/horiz_grid_20030806.ieeer8

  topography_fname :

  inputdata/ocn/pop/gx3v7/grid/topography_20100105.ieeei4

  region_mask_fname :

  inputdata/ocn/pop/gx3v7/grid/region_mask_20090831.ieeei4

  type :

  dipole

  title :

  POP_gx3v7 grid

Generate a field with some missing values

field = ds.KMT.copy() * 1.0
field = field.where(ds.KMT > 0)
field.values[20:40, :] = np.nan
field.plot()

<matplotlib.collections.QuadMesh at 0x7fa74d8880a0>

Fill the missing values

field_filled = pop_tools.lateral_fill(field, isvalid_mask=(ds.KMT > 0))
h = field_filled.plot()

%load_ext watermark
%watermark -d -iv -m -g -h

Compiler    : GCC 11.3.0
OS          : Linux
Release     : 5.15.0-1004-aws
Machine     : x86_64
Processor   : x86_64
CPU cores   : 2
Architecture: 64bit

Hostname: build-21213814-project-451810-pop-tools

Git hash: d3c80c0576ae4838c0e04a0157734eb0c977e613

xarray   : 2023.6.0
numpy    : 1.24.4
pop_tools: 2023.3.0.post2+dirty