nmc_met_io icon indicating copy to clipboard operation
nmc_met_io copied to clipboard

nmc_met_io/nmc_met_io/retrieve_cmadaas.py中cmadaas_model_3D_grids函数体中,example运行失败

Open stonemoonkey opened this issue 3 years ago • 1 comments

>>> valid_times = [6*i for i in range(13)]
>>> levels = [1000, 950, 925, 900, 850, 800, 700, 600, 500, 400, 300, 250, 200, 100]
>>> data = cmadaas_model_3D_grids("NAFP_FOR_FTM_HIGH_EC_ANEA", "2020021512", 24, 'TEM', levels, 100, units="C", scale_off=[1.0, -273.15], 
                                 levattrs={'long_name':'pressure_level', 'units':'hPa', '_CoordinateAxisType':'Pressure'})

参数24,似应改为valid_times,否则报错:int类型不能被迭代。

stonemoonkey avatar Apr 02 '22 08:04 stonemoonkey

例子能够正常运行, 24表示预报时效,不能为列表.

levels = [1000, 950, 925, 900, 850, 800, 700, 600, 500, 400, 300, 250, 200, 100]
data = cmadaas_model_3D_grid("NAFP_FOR_FTM_HIGH_EC_ANEA", "2020021512", 24, 'TEM', levels, 100, units="C", scale_off=[1.0, -273.15],  levattrs={'long_name':'pressure_level', 'units':'hPa', '_CoordinateAxisType':'Pressure'})
print(data)
<xarray.Dataset>
Dimensions:                  (time: 1, level: 14, lat: 281, lon: 361)
Coordinates:
  * time                     (time) datetime64[ns] 2020-02-16T12:00:00
  * level                    (level) int64 1000 950 925 900 ... 300 250 200 100
  * lat                      (lat) float64 60.0 59.75 59.5 ... -9.5 -9.75 -10.0
  * lon                      (lon) float64 60.0 60.25 60.5 ... 149.5 149.8 150.0
    forecast_reference_time  datetime64[ns] 2020-02-15T12:00:00
    forecast_period          (time) float64 24.0
Data variables:
    data                     (time, level, lat, lon) float32 -5.413 ... -79.57
Attributes:
    Conventions:  CF-1.6
    Origin:       CIMISS Server by MUSIC API

NMC-DAVE avatar Jun 01 '22 13:06 NMC-DAVE