【已解决】NCL对FNL风场平均

channam

---------------2020/04/10----------已解决----------------
目的：对2020年1月的FNL数据批量读取，并对风场UV进行平均（每个经纬度点上 时间维度的平均）
下面是从网上看到的批量读取的方法：
        dir   = "./fnl/"
        files = systemfunc("ls " + dir + "fnl_202001*")     ;批量读取数据
        f     = addfiles(files,"r")
        u = f[:]->UGRD_P0_L103_GLL0({10},：,：)     （其实我也不太懂f[:]->和f->的区别是什么，猜想是[:]是指所有fnl数据，f是单个fnl数据）
        v = f[:]->VGRD_P0_L103_GLL0({10},：,：)
是应该增加一个时间维度，再dim_avg_n对时间维度进行平均吗？有没有做过类似工作的坛友可以提点意见呢？

一大碗年糕

[:]是读取列表中多个文件的同一变量
NCL有两种方式读取同一变量，分别是cat和join，通过ListSetType设定，默认的处理方式是cat，只增加变量的长度而不增加其维数，如果你要用dim_avg_n对时间维平均的话最好用join读取方式,比如
ListSetType(f,"join")
v = f[:]->VGRD_P0_L103_GLL0({10},：,：)
这样V的维数就比原来的VGRD_P0_L103_GLL0多一维了变成v（读取的文件数量，{10}，：，：}

channam

一大碗年糕 发表于 2020-4-10 12:46
[:]是读取列表中多个文件的同一变量
NCL有两种方式读取同一变量，分别是cat和join，通过ListSetType设定， ...

哇~  又学了一招，谢谢你啦，试了一下解决了困扰我的问题。

channam

出现新问题，自己顶一下。