问一下如何把多年的ncep逐日资料整合成一个文件

puyiliang · 发表于 2014-6-10 21:51:27

longlivehj 发表于 2014-6-7 19:40
addfiles读取全部文件，然后用ncl官网examples下file i/o里output netcdf中提到的第一种方法写到文件就好了 ...

朋友，能不能帮我看看，这个转数据因为之前没弄过，自己用你介绍的方法试了一下，应该是可以的，但是我弄的这个脚本运行老是出错，能不能帮忙看看，哪里错了....
load "$NCARG_ROOt/lib/ncarg/nclscripts/csm/gsn_code.ncl"

load "$NCARG_ROOt/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

load "$NCARG_ROOt/lib/ncarg/nclscripts/csm/contributed.ncl"

load "$NCARG_ROOt/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl"

;******************** MAIN **********************************


begin
  ;========================
  ; get list of all files and open as "one big file"
  ;========================
   all_files = systemfunc ("ls /home/puyiliang/data/Uwnd_daily/uwnd.*.nc")
   fall    = addfiles (all_files, "r") ; note the "s" of addfile
  ;========================
  ; choose how files are combined and read in variable across files
  ;========================
   ListSetType (fall, "cat")          ; concatenate or "merge" (default)

   uwnd = fall[:]->uwnd          ; note syntax [:]

   time=fall[:]->time

printVarSummary(uwnd)

ncdf    = addfile("Uwnd_daily.nc" ,"c")  ; open output netCDF file

filedimdef(ncdf,"time",-1,True)
ncdf->uwnd  = uwnd
end
运行的时候老是说倒数第二，第三行有问题...
先谢了

longlivehj · 发表于 2014-6-11 07:40:53

puyiliang 发表于 2014-6-10 21:51
朋友，能不能帮我看看，这个转数据因为之前没弄过，自己用你介绍的方法试了一下，应该是可以的，但是我弄 ...

呵呵，看不出有什么问题。提示有什么错误啊？
另外，你那个NCARG_ROOt环境变量，最后一个t小写，挺奇怪的。环境变量一般都大写。

huayazi · 发表于 2014-7-30 15:13:34

{:5_235:}

wodeshiwo · 发表于 2014-8-1 16:19:20

用NCO可以

veronica · 发表于 2014-8-26 17:04:11

利用ncrcat命令　

aselly · 发表于 2014-10-21 20:29:47

请问楼主这个程序问题改好了吗，我也正在弄这个

xcy_NMC · 发表于 2014-10-22 08:27:45

可以用ＮＣＯ的，直接连接需要的变量

xigua · 发表于 2014-10-22 08:31:21

cdo 一个命令搞定

wyanghu · 发表于 2018-7-20 21:08:49

请问楼主最后解决了吗？可以分享下经验吗？谢谢

puyiliang · 发表于 2018-8-1 17:24:12

wyanghu 发表于 2018-7-20 21:08
请问楼主最后解决了吗？可以分享下经验吗？谢谢

选择自己要研究的区域，这样文件就会小一点，可以直接合并

		自动登录	找回密码
密码			立即注册

[经验总结] 问一下如何把多年的ncep逐日资料整合成一个文件