NCL函数的小心得

贫道敬孔

NCL中的函数要求最左维是什么参数，那这一维必须是什么参数。这时候如果需要用循环，即把第一维变为循环变量，第二维是原最左维的参数的话，也不行。
这时候不能用循环，只能用数组一个个算。

如使用gsn_csm_pres_hgt函数，要求变量最左边是lev，而如果有多图需要画，那首先：


1、多图的变量不能是第一维，第一维只能是层数

2、变量只能在第一维以外，如第二维

3、变量如果在第二/三维，这时候就不能用循环，会报错

这时候要画多图只能老老实实一个变量一个变量来画。
比如
        plot                      = new(7,graphic)  

        plot(0) = gsn_csm_pres_hgt(wks ,vd_Td(:,0,:)  , res)
        plot(1) = gsn_csm_pres_hgt(wks ,vd_Td(:,1,:)  , res)
        plot(2) = gsn_csm_pres_hgt(wks ,vd_Td(:,2,:)  , res)
        plot(3) = gsn_csm_pres_hgt(wks ,vd_Td(:,3,:)  , res)
        plot(4) = gsn_csm_pres_hgt(wks ,vd_Td(:,4,:)  , res)
        plot(5) = gsn_csm_pres_hgt(wks ,vd_Td(:,5,:)  , res)
        plot(6) = gsn_csm_pres_hgt(wks ,vd_Td_mean    , res)


目前心得是这样，大家还有更好的办法的话欢迎交流。

yangzi88552

你已经给出了分开写的，例如你给的例子就可以写成
do i = 0,6
     plot(i) = gsn_csm_pres_hgt(wks ,vd_Td(:,i,:)  , res)
end do
你要多种变量的话，甚至你可以把变量名存为一个新的数组，然后做循环一样的。
ncl的变量名不能用来f->var提取数据，但是用来画图是可以的

贫道敬孔

yangzi88552 发表于 2022-8-3 10:27
你已经给出了分开写的，例如你给的例子就可以写成
do i = 0,6
     plot(i) = gsn_csm_pres_hgt(wks ,vd_ ...

您说的这个方法我试过，但是plot()这个函数貌似循环只能在第一维呢，放在第二维就出错了

yangzi88552

贫道敬孔 发表于 2022-8-3 14:29
您说的这个方法我试过，但是plot()这个函数貌似循环只能在第一维呢，放在第二维就出错了

这个你需要查查你两种方式提取变量后的，变量具体的格式，这个和循环放在第几个维度是没有关系的，函数中看的只是你输入的变量是否符合规则

LDW

这种要求最左边或者最后边维度的，可以试试把维度调整一下，比如原来是x1（time,pressure,lat）,x = x1(pressure|:,time|:,lat|:)

贫道敬孔

LDW 发表于 2022-8-4 09:41
这种要求最左边或者最后边维度的，可以试试把维度调整一下，比如原来是x1（time,pressure,lat）,x = x1(pre ...

嗯呢，本来数组的最左维是时间，因为要求最左维是层数，所以已经调整维度了呢

贫道敬孔

yangzi88552 发表于 2022-8-3 17:40
这个你需要查查你两种方式提取变量后的，变量具体的格式，这个和循环放在第几个维度是没有关系的，函数中 ...

之前一层循环不在最外围的都没有成功过，我以为是NCL做单层循环必须放到最外围。我再试试，谢谢回复

wrry22

谢谢分享！！~~

任性小王子

非常感谢楼主分享，学习了。

贫道敬孔

wrry22 发表于 2022-8-6 16:58
谢谢分享！！~~

不客气哒，我也在大家的分享中学习了很多