编辑-地球科学江南体育网页版堆栈交换江南电子竞技平台

你不是登录。你的编辑将被放置在一个队列,直到同行评议。

我们欢迎编辑,使读者更容易理解和更有价值。因为社区成员审查编辑,请尽量使文章更好比你怎么发现,例如,通过修正语法或添加额外的资源和超链接。

牧师

标题

身体

我已经创建了一个x, y netcdf文件和现在想添加额外的时间能量变量。这为每个网格单元能量时间价值。我的脚本在Python中是这样的:value = loadtxt (“countries.txt”, skiprows = 0) %国家data = loadtxt (“energy.txt”, skiprows = 0) %能量的能量。txt文件如:[国家id]…51 52 54 55 56 59 60 62…[时间步1]…649年南南南1083 9979南13602…[时间步2]…5743 43444 6365 1283 9062 4026南南…(时间步n)…………………………等(直到87648步),感兴趣的细胞是第一行:“cellsOfInt =数据[0,:]”——对一些细胞我nan值以来没有值的时间步骤。 - Each row correspond to one time step so for cell of interest 51, in the first time step we have nan value, in the second time step we have 5743 and so on till 87648. - All cells of interest are in countries.txt file that is 950 X 1000 with cells of interest ids or nan when there is no value. My code: y = nc.createDimension("y", 950) #y x = nc.createDimension("x", 1000) #x time = nc.createDimension("time", 87648) #time latitude=nc.createVariable("y", "f8", ("y", )) longitude=nc.createVariable("x", "f8", ("x", )) time = nc.createVariable("time","f8",("time",)) countries = nc.createVariable("countries", "f8", ("y", "x")); energy = nc.createVariable("energy", "f8", ("time", "y", "x"), fill_value=-9999, chunksizes=(1, 950, 1000)); time.standard_name='time' time.units ='hours since 2006-01-01 00:00:00.0' time.calendar='proleptic_gregorian' lats = np.arange(5497500,747500,-5000) #y lons = np.arange(2502500,7502500,5000) #x # Fill variables latitude[:] = lats longitude[:] = lons countries [:]=value UNTIL HERE WORKS FINE I have my 2d country but then it does not give any result for energy: for i in xrange(0,950): #rows for j in xrange(0,1000): #columns for n in xrange(0,87648): #rows for m in xrange(0,1): #columns if data[1,m] == countries[i,j]: #correspondance nc.variables["energy"][n]=data[n,m] Any idea how to fix this so I can have in country.netcdf file time dimension with energy values at each time step? Thanks a lot!!

标签

编辑总结

取消

1

\ begingroup美元三个评论(我不熟悉python):(一)数组索引从0开始(1)在Python中,对吧?那么它应该如果数据[0,m] = =国家(i, j):因为你是感兴趣的第一行;(b)米从0来1(的m xrange (0,1):),但是它应该从0来number_of_countries这的所有列能源是迭代。(c)的顺序n和米for循环似乎是完全无效的。更好的是:米…:如果数据[1 m] = =…:n…:nc.variables…=…。或更换后一个for循环nc.variables[“能量”][:]=数据(:,m) \ endgroup美元
- - - - - -daniel.heydebreck
2017年10月22日在39
1

\ begingroup美元看一看xarray,我认为这可以让你设置更加容易。 \ endgroup美元
- - - - - -user2821
2017年10月22日,在十一21
\ begingroup美元我认为使用nco和/或其他工具将是一个更直接的方法创建和增加netcdf文件。例如netcdf厨房水槽linux.die.net/man/1/ncks。看到一个列表:unidata.ucar.edu/software/netcdf/software.html \ endgroup美元
- - - - - -f.thorpe ♦
2017年11月22日在17:26

添加一个评论 |

正确的小错误或错误
澄清的意思,而不更改它
添加相关资源或链接
总是尊重作者的意图
不使用编辑与作者回复