2017-02-23 39 views
1

我遇到了使用netcdf4函數寫入netcdf3文件的問題。我嘗試使用創建可變功能,但它給了我這個錯誤:創建NetCDF:試圖創建NetCDF-4操作上的NetCDF-3文件在python中使用模塊netcdf4寫入NetCDF3文件

nc = Dataset(root.fileName,'a',format="NETCDF4") 

Hycom_U = nc.createVariable('/variables/Hycom_U','float',('time','lat','lon',)) 

Hycom_V = nc.createVariable('/variables/Hycom_V','f4',('time','lat','lon',)) 

NC =

root group (NETCDF3_CLASSIC data model, file format NETCDF3): 
    netcdf_library_version: 4.1.3 
    format_version: HFRNet_1.0.0 
    product_version: HFRNet_1.1.05 
    Conventions: CF-1.0 
    title: Near-Real Time Surface Ocean Velocity, Hawaii, 
2 km Resolution 
    institution: Scripps Institution of Oceanography 
    source: Surface Ocean HF-Radar 
    history: 22-Feb-2017 00:55:46: NetCDF file created 
22-Feb-2017 00:55:46: Filtered U and V by GDOP < 1.25 ; 
FMRC Best Dataset 
    references: Terrill, E. et al., 2006. Data Management and Real-time 
Distribution in the HF-Radar National Network. Proceedings 
of the MTS/IEEE Oceans 2006 Conference, Boston MA, 
September 2006. 
    creator_name: Mark Otero 
    creator_email: [email protected] 
    creator_url: http://cordc.ucsd.edu/projects/mapping/ 
    summary: Surface ocean velocities estimated from HF-Radar are 
representative of the upper 0.3 - 2.5 meters of the 
ocean. The main objective of near-real time 
processing is to produce the best product from 
available data at the time of processing. Radial 
velocity measurements are obtained from individual 
radar sites through the U.S. HF-Radar Network. 
Hourly radial data are processed by unweighted 
least-squares on a 2 km resolution grid of Hawaii 
to produce near real-time surface current maps. 
    geospatial_lat_min: 20.487279892 
    geospatial_lat_max: 21.5720806122 
    geospatial_lon_min: -158.903594971 
    geospatial_lon_max: -157.490005493 
    grid_resolution: 2km 
    grid_projection: equidistant cylindrical 
    regional_description: Unites States, Hawaiian Islands 
    cdm_data_type: GRID 
    featureType: GRID 
    location: Proto fmrc:HFRADAR,_US_Hawaii,_2km_Resolution,_Hourly_RTV 
    History: Translated to CF-1.0 Conventions by Netcdf-Java CDM (NetcdfCFWriter) 
Original Dataset = fmrc:HFRADAR,_US_Hawaii,_2km_Resolution,_Hourly_RTV; Translation Date = Thu Feb 23 13:35:32 GMT 2017 
    dimensions(sizes): time(25), lat(61), lon(77) 
    variables(dimensions): float32 u(time,lat,lon), float64 time_run(time), float64 time(time), float32 lat(lat), float32 lon(lon), float32 v(time,lat,lon) 
    groups: 

什麼的NetCDF 3個行動,我可以用於將數據添加到文件中?我發現我可以手動添加數據,只需執行這個直接添加數據的nc.variables["Hycom_U"]=U2,但沒有別的。有一個更好的方法嗎?

回答

0

我想通了!我根本無法使用路徑作爲varname。

Hycom_U = nc.createVariable( 'Hycom_U', '浮動',( '時間', 'LAT', 'LON',))

它正確地創建的變量我。

1

我認爲問題是,你聲稱該文件是netCDF4格式:

nc = Dataset(root.fileName,'a',format="NETCDF4")` 

但你真的想表明它是netCDF3:

nc = Dataset(root.fileName,'a',format="NETCDF3_CLASSIC") 

的其他文檔可以發現here

+0

不幸的是,它不能解決問題。我嘗試了不同的格式,他們都給出了同樣的錯誤。 :( – EagleTamer

相關問題