ncl关于风场距平检验的新手问题

marklee

由于四月刚换题目，然后开始学ncl做论文，而今年居然提前交初稿，然后检验的问题基本存在于我每一张做回归系数并检验的图，求助学长学姐无果，作为一个菜鸟也实在不知道，只能求助于论坛了

关于风场距平的显著性检验：

我的想法是将通过显著性检验（概率<0.05即siglvl）的u或v的概率值 直接赋予一个小于0.05的定值数字（0.01），没通过的赋予confi_rcv缺测值，然后进行二次检验（必须，否则无阴影，为什么？），再将confi_rcv的坐标赋给新变量confi_rc，直接对新变量进行画图（其中confi_rcv与confi_rcu是student检验的结果 概率值，已转化为二维数组）

  confi_rc = where(confi_rcv.lt.siglvl.or.confi_rcu.lt.siglvl,0.01,confi_rcv@_FillValue)

  confi_rc= where(confi_rc.lt.siglvl,confi_rc,confi_rc@_FillValue)

  copy_VarCoords(confi_rcv,confi_rc)

结果是存在如下警告，但是能画图：

warning:Attempt to reference attribute (_FillValue) which is undefined

画出来的图见 图1，通过显著性检验的地区非常多

第二种语句是

  confi_rc = where(confi_rcv.lt.siglvl.or.confi_rcu.lt.siglvl,confi_rcv,confi_rcv@_FillValue)

  confi_rc= where(confi_rc.lt.siglvl,confi_rc,confi_rc@_FillValue)

  copy_VarCoords(confi_rcv,confi_rc_slp)

想法是：将通过显著性检验的u或v的概率值（<0.05）   赋予一个变量的值：confi_rcu或者confi_rcv，但是这个语句本身存在一定问题，就是除非confi_rcv或u本身不存在缺测，最终“.or.”才成立。但是和上面的语句相比，不都是赋的小于0.05的数值吗？画出来的图是这样的，且不存在警告：

图2	图3
赋值confi_rcv	赋值confi_rcu
通过显著性检验的区域明显不同，和前面赋值具体数值数值，差距很大

所用的画图语句都是：

  base = gsn_csm_contour_map(wks,confi_rc,res)

......

  base=gsn_contour_shade(base,siglvl,-999.,opt)   ;siglvl=0.05

其他无任何改动

如果说风场距平的检验不好说明这个问题，我对高度场回归系数进行相同的测试：

第一种语句：此时不存在风场距平的第一种语句警告

confi_rc = where(confi_rc_slp.lt.siglvl,0.01,confi_rc_slp@_FillValue)

  copy_VarCoords(confi_rc_slp,confi_rc)

其中confi_rc_slp仍为student检验的结果：概率值

第二种语句：

confi_rc = where(confi_rc_slp.lt.siglvl,confi_rc_slp,confi_rc_slp@_FillValue)

  copy_VarCoords(confi_rc_slp,confi_rc)

其他语句都相同，结果如下：

图4	图5
第一种语句	第二种语句
可以看到结果显著性差异十分显著

所以不知道哪种语句才是正确的？如果都不正确，这个检验要怎么写呢？

kathy-

首先warning:Attempt to reference attribute (_FillValue) which is undefined这行报错的意思是你没有定义你的变量的缺省值是什么。

marklee

kathy- 发表于 2019-4-20 19:17
首先warning:Attempt to reference attribute (_FillValue) which is undefined这行报错的意思是你没有定义 ...

嗯嗯谢谢，但是高度场基本是同样的语句就没有出现这个警告

vicky-zz

这两个图显著性好像刚好相反