WRF-Python插值结果如何输入netcdf文件

AllenLu

如题，WRF-Python插值结果如何输入netcdf文件呢

from netCDF4 import Dataset
from wrf import getvar, interplevel, ALL_TIMES

wrfin = Dataset("wrfout_d03_2019-10-01_00:00:00")

RH = getvar(wrfin, "wspd_wdir",timeidx=ALL_TIMES)
height = getvar(wrfin, "height_agl",timeidx=ALL_TIMES)
pblh = getvar(wrfin, "PBLH", timeidx=ALL_TIMES)

RH_pblh = interplevel(RH , height , pblh)


要怎么把RH_pblh输入到.nc文件呢，萌新求助o(╥﹏╥)o

AllenLu

f = h5py.File("RH_pblh.nc",mode = 'w')
Dataset = f.create_dataset('ws',data = RH_pblh ,compression='gzip')

这样可以创建2D文件，但是想创建Geo2D文件要怎么做呢？
有没有大神知道呢？

AllenLu

1. # the Scientific Python netCDF 3 interface
   
2. # http://dirac.cnrs-orleans.fr/ScientificPython/
   
3. #from Scientific.IO.NetCDF import NetCDFFile as Dataset
   
4. # the 'classic' version of the netCDF4 python interface
   
5. # http://code.google.com/p/netcdf4-python/
   
6. from netCDF4 import Dataset
   
7. from numpy import arange, dtype # array module from http://numpy.scipy.org
   
8. """
   
9. This is an example program which writes some 4D pressure and
   
10. temperatures.
    
11. The companion program pres_temp_4D_rd.py shows how
    
12. to read the netCDF data file created by this program.
    
13. 
    
14. This example demonstrates the netCDF Python API.
    
15. It will work either with the Scientific Python NetCDF version 3 interface
    
16. (http://dirac.cnrs-orleans.fr/ScientificPython/)
    
17. of the 'classic' version of the netCDF4 interface.
    
18. (http://netcdf4-python.googlecode.com/svn/trunk/docs/netCDF4_classic-module.html)
    
19. To switch from one to another, just comment/uncomment the appropriate
    
20. import statements at the beginning of this file.
    
21. 
    
22. Jeff Whitaker <jeffrey.s.whitaker@noaa.gov> 20070202
    
23. """
    
24. # the netCDF variable will be nrecs x nlevs x nlats x nlons.
    
25. nrecs = 2; nlevs = 2; nlats = 6; nlons = 12
    
26. # open a new netCDF file for writing.
    
27. ncfile = Dataset('pres_temp_4D.nc','w')
    
28. # latitudes and longitudes of grid
    
29. lats_out = -25.0 + 5.0*arange(nlats,dtype='float32')
    
30. lons_out = -125.0 + 5.0*arange(nlons,dtype='float32')
    
31. # output data.
    
32. press_out = 900. + arange(nlevs*nlats*nlons,dtype='float32') # 1d array
    
33. press_out.shape = (nlevs,nlats,nlons) # reshape to 2d array
    
34. temp_out = 9. + arange(nlevs*nlats*nlons,dtype='float32') # 1d array
    
35. temp_out.shape = (nlevs,nlats,nlons) # reshape to 2d array
    
36. # create the lat and lon dimensions.
    
37. ncfile.createDimension('latitude',nlats)
    
38. ncfile.createDimension('longitude',nlons)
    
39. # create level dimension.
    
40. ncfile.createDimension('level',nlevs)
    
41. # create time dimension (record, or unlimited dimension)
    
42. ncfile.createDimension('time',None)
    
43. # Define the coordinate variables. They will hold the coordinate
    
44. # information, that is, the latitudes and longitudes.
    
45. # Coordinate variables only given for lat and lon.
    
46. lats = ncfile.createVariable('latitude',dtype('float32').char,('latitude',))
    
47. lons = ncfile.createVariable('longitude',dtype('float32').char,('longitude',))
    
48. # Assign units attributes to coordinate var data. This attaches a
    
49. # text attribute to each of the coordinate variables, containing the
    
50. # units.
    
51. lats.units = 'degrees_north'
    
52. lons.units = 'degrees_east'
    
53. # write data to coordinate vars.
    
54. lats[:] = lats_out
    
55. lons[:] = lons_out
    
56. # create the pressure and temperature variables
    
57. press = ncfile.createVariable('pressure',dtype('float32').char,('time','level','latitude','longitude'))
    
58. temp = ncfile.createVariable('temperature',dtype('float32').char,('time','level','latitude','longitude'))
    
59. # set the units attribute.
    
60. press.units =  'hPa'
    
61. temp.units = 'celsius'
    
62. # write data to variables along record (unlimited) dimension.
    
63. # same data is written for each record.
    
64. for nrec in range(nrecs):
    
65.    press[nrec,:,::] = press_out
    
66.    temp[nrec,:,::] = temp_out
    
67. # close the file.
    
68. ncfile.close()
    
69. print '*** SUCCESS writing example file pres_temp_4D.nc'

# the Scientific Python netCDF 3 interface<br /> # http://dirac.cnrs-orleans.fr/ScientificPython/<br /> #from Scientific.IO.NetCDF import NetCDFFile as Dataset<br /> # the 'classic' version of the netCDF4 python interface<br /> # http://code.google.com/p/netcdf4-python/<br /> from netCDF4 import Dataset<br /> from numpy import arange, dtype # array module from http://numpy.scipy.org<br /> """<br /> This is an example program which writes some 4D pressure and<br /> temperatures. <br /> The companion program pres_temp_4D_rd.py shows how<br /> to read the netCDF data file created by this program.<br /> <br /> This example demonstrates the netCDF Python API.<br /> It will work either with the Scientific Python NetCDF version 3 interface<br /> (http://dirac.cnrs-orleans.fr/ScientificPython/)<br /> of the 'classic' version of the netCDF4 interface. <br /> (http://netcdf4-python.googlecode.com/svn/trunk/docs/netCDF4_classic-module.html)<br /> To switch from one to another, just comment/uncomment the appropriate<br /> import statements at the beginning of this file.<br /> <br /> Jeff Whitaker <<a href=\"mailto:jeffrey.s.whitaker@noaa.gov\">jeffrey.s.whitaker@noaa.gov</a>> 20070202<br /> """<br /> # the netCDF variable will be nrecs x nlevs x nlats x nlons.<br /> nrecs = 2; nlevs = 2; nlats = 6; nlons = 12<br /> # open a new netCDF file for writing.<br /> ncfile = Dataset('pres_temp_4D.nc','w') <br /> # latitudes and longitudes of grid<br /> lats_out = -25.0 + 5.0*arange(nlats,dtype='float32')<br /> lons_out = -125.0 + 5.0*arange(nlons,dtype='float32')<br /> # output data.<br /> press_out = 900. + arange(nlevs*nlats*nlons,dtype='float32') # 1d array<br /> press_out.shape = (nlevs,nlats,nlons) # reshape to 2d array<br /> temp_out = 9. + arange(nlevs*nlats*nlons,dtype='float32') # 1d array<br /> temp_out.shape = (nlevs,nlats,nlons) # reshape to 2d array<br /> # create the lat and lon dimensions.<br /> ncfile.createDimension('latitude',nlats)<br /> ncfile.createDimension('longitude',nlons)<br /> # create level dimension.<br /> ncfile.createDimension('level',nlevs)<br /> # create time dimension (record, or unlimited dimension)<br /> ncfile.createDimension('time',None)<br /> # Define the coordinate variables. They will hold the coordinate<br /> # information, that is, the latitudes and longitudes.<br /> # Coordinate variables only given for lat and lon.<br /> lats = ncfile.createVariable('latitude',dtype('float32').char,('latitude',))<br /> lons = ncfile.createVariable('longitude',dtype('float32').char,('longitude',))<br /> # Assign units attributes to coordinate var data. This attaches a<br /> # text attribute to each of the coordinate variables, containing the<br /> # units.<br /> lats.units = 'degrees_north'<br /> lons.units = 'degrees_east'<br /> # write data to coordinate vars.<br /> lats[:] = lats_out<br /> lons[:] = lons_out<br /> # create the pressure and temperature variables <br /> press = ncfile.createVariable('pressure',dtype('float32').char,('time','level','latitude','longitude'))<br /> temp = ncfile.createVariable('temperature',dtype('float32').char,('time','level','latitude','longitude'))<br /> # set the units attribute.<br /> press.units =  'hPa'<br /> temp.units = 'celsius'<br /> # write data to variables along record (unlimited) dimension.<br /> # same data is written for each record.<br /> for nrec in range(nrecs):<br />    press[nrec,:,::] = press_out<br />    temp[nrec,:,::] = temp_out<br /> # close the file.<br /> ncfile.close()<br /> print '*** SUCCESS writing example file pres_temp_4D.nc'

这个好像可以，但是他是自己做的数据，我wrf-python插值得到的数据要怎么才能输出到nc文件呢