在Python中读取netCDF文件时用户输入文件名

Question

我有一套从1953年到2014年的土壤湿度数据文件。它们都是土的形式_YYYY.nc版本（其中YYYY是其中一个年份）。有没有一种方法可以让我要求用户输入用户希望查看的年份，并让我的程序打开相应的功能？目前，我在gedit中手动更改年份，并编写函数来获取每个变量（土壤湿度随时间、纬度、经度的变化）：

    import netCDF4 as nc
    import numpy as np
    import numpy.ma as ma
    import csv as csv

    fid=nc.MFDataset('/data/reu_data/soil_moisture/cpc_soil_1957.nc','r')
    fid.close()
    ncf='/data/reu_data/soil_moisture/cpc_soil_1957.nc'

    def read_var(ncfile, varname):
        fid=nc.Dataset(ncfile, 'r')
        out=fid.variables[varname][:]   
        fid.close()
        return out

    time=read_var(ncf, 'time')
    lat=read_var(ncf, 'lat')
    lon=read_var(ncf, 'lon')
    soil=read_var(ncf, 'soilw')

Answer 1

您可以使用input()要求用户输入年份。然后可以使用它来生成文件路径。你知道吗

...
year = input("Enter year: "))
filename = '/data/reu_data/soil_moisture/cpc_soil_%s.nc' % (year,)
fid=nc.MFDataset(filename,'r')
fid.close()
...

您应该进行错误检查，以确保用户输入的值实际上是一年，并且在您的数据范围内。你知道吗

您可以阅读更多关于Pythonhere中输入/输出的内容。你知道吗