python - netCDF のグループから変数を抽出する

Question

foll を含む netCDF ファイルがあります。ncdump:

netcdf test_nc {
dimensions:
    time = UNLIMITED ; // (20 currently)
    latitude = 360 ;
    longitude = 720 ;
    N = 3 ;
    strlen = 1 ;
variables:
    float data_array(time, latitude, longitude, N) ;
        data_array:_FillValue = -9999.f ;
        data_array:units = "1" ;
        data_array:long_name = "data_array" ;
    char N(N, strlen) ;
    double latitude(latitude) ;
        latitude:standard_name = "latitude" ;
        latitude:units = "degrees_north" ;
    double longitude(longitude) ;
        longitude:standard_name = "longitude" ;
        longitude:units = "degrees_east" ;
    double time(time) ;
        time:standard_name = "time" ;
        time:units = "days since 2000-01-01 00:00:00.0" ;
        time:calendar = "gregorian" ;
}

data_arrayグループ内の最後の変数を抽出するにはどうすればよいですか? を使用すると、次のようncksにすべてを抽出できます。data_array

ncks -v data_array test_nc.nc output_nc.nc

ただし、 で N=3 に対応する変数のみを抽出したいdata_array。python netCDF4 または nco ツールを使用するソリューションはすべて機能しますが (cdo は使用しません)、Windows でも機能する必要があります。