关于NCL添加多文件的一些方法

刘哲

以下是在linux下 运行的load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
diri ="文件的绝对路径(上一级路径)"
fils = systemfunc("ls " + diri + "文件名（模糊的文件名）*. grib2") // 这个是把括号里的命令传递给linux， 然后把结果以字符串的形式传递给fils
f = addfiles(fils, "r")
ListSetType(f , "cat") // 如何排列添加的多文件，这里需要说明后面的参数，有两个参数（join ，cat） ，其中join 是把多个文件并列排布，就把每个文   件当成一个数组， join就是按文件添加的顺序（源文件的排列顺序）， 并排或者顺着排列，这样就f 相比元数据多出一维，多出的维就是文件的序号，例如f[0,:,:,:] ,这个就是读取第一个文件的所有数据。 cat 参数是把所有文件的元数据都串联起来，生成一个和元数据一样维数的数组，数据排列也是按照添加的顺序，这样的话，数据数量，就是文件数*元数据第一维的数量。
v = f[:]-> UN // UN代表你需要添加元数据
// 下面求多文件的累加和平均
v1 = dim_avg_n(f , n) // 后面一个n代表你要对f的哪个维求平均，这里要说明的是，就是 对那个维求平均的话，例如对time变量平均，就是把time上的数据点，累加在平均，放到它对应的经纬度的栅格点上（图在后面）。
   v2 = dim_sum _n(f ,n ) // 后面一个n代表你要对f的哪个维求和，和avg 很相似只是把time轴上的所有值相加，放到它对应的经纬度的栅格点上（图在后面）

三维栅格点图

残月

这个好，以前都是用变量加循环完成的，赞。

longlivehj

元数据有专门的定义，楼主可以查一下。感觉帖子中使用元数据的说法不妥。

奔奔

这个好诶。。。多谢多谢

刘哲

longlivehj 发表于 2014-8-9 17:36
元数据有专门的定义，楼主可以查一下。感觉帖子中使用元数据的说法不妥。

恩，好的！主要我也是新手，把自己的遇到的问题，及解决方法，与大家分享！师傅，您多指导

crystalone

给你支持下~

刘哲

crystalone 发表于 2014-8-18 17:16
给你支持下~

谢谢，我是一个新手，希望大家能互相帮忙

levisk

ncl 6> fils = systemfunc("ls " + "/home/data/" + "auxhist3_d02_2017-04-*-00-00")
fatal:Number of dimensions on right hand side do not match number of dimension in left hand side
请问加载文件报这个错怎么解决呢