也说格点数据插值为站点数据

朝雾

df = "J:\\CMIP5-data\\RegCM-data1\\data\\rcp45.ctl"
myApp.OpenGrADSData(df)
#---- 获取地面格点数据
sTime = DateTime.Parse("2006-01-01 00:00")
eTime = DateTime.Parse("2006-12-30 00:00")
while sTime <= eTime:
  inFile = sTime.ToString("yyyyMMdd")
  print inFile
  gData = myApp.GetGridData("t2m")
  #---- 设置输入、输出站点数据文件
  inf = "J:\\CMIP5-data\\RegCM-data1\\data\\station.csv"
  outf = "J:\\CMIP5-data\\RegCM-data1\\data\\rcp45-t2m\\rcp45_"+sTime.ToString("yyyyMMdd")+".txt"
  print outf
  #---- 格点数据插值为站点数据，并输出为文件
  gData.ToStation(inf,outf)
  sTime = sTime.AddHours(24)

[b]未能实现循环，把第一个时次的数据循环处理了很多遍，不知道问题出在哪了？不好意思又来问了急！/b]

MeteoInfo

朝雾 发表于 2013-11-4 21:41

                               
登录/注册后可看大图

df = "J:\\CMIP5-data\\RegCM-data1\\data\\rcp45.ctl"
myApp.OpenGrADSData(df)
#---- 获取地面格点数据 ...

参考此贴中的脚本循环：MeteoInfo脚本示例：格点数据平均
http://bbs.06climate.com/forum.p ... 604&fromuid=106

hubert

能从nc数据中提取一个点的值吗，格点数据700*400*248

MeteoInfo

hubert 发表于 2013-11-7 22:57

                               
登录/注册后可看大图

能从nc数据中提取一个点的值吗，格点数据700*400*248

试试不就知道了

hubert

关键是不知道怎么试啊

MeteoInfo

hubert 发表于 2013-11-8 09:00

                               
登录/注册后可看大图

关键是不知道怎么试啊

帖子里不是都讲了吗

hubert

那个文件太大，不好上传啊

rainbownn

有用。谢谢分享

朝雾

使用MeteoInfo的过程中发现一个问题，grads使用一个描述文件统一描述多个时间序列原始数据无法正常读取数据文件。
如：ctl
dset  rcp45_%y4.dat
title RegCM daily output variables
options little_endian
options template
options  365_day_calendar
undef -1.e34
xdef  27 linear   102  0.5
ydef  31 linear   27 0.5
zdef 1 levels 1008.23
tdef   34310 linear 01jan2006 1dy
读取后只能每年1月1日的数据

MeteoInfo

朝雾 发表于 2013-11-13 14:19

                               
登录/注册后可看大图

使用MeteoInfo的过程中发现一个问题，grads使用一个描述文件统一描述多个时间序列原始数据无法正常读取数据 ...

对MeteoInfo库文件进行了更新尝试解决读取GrADS Template数据的bug，不过没有合适的数据没法做测试，你可以下载最新文件试试（见置顶帖子）。