【已解决】ncl 关于整层水汽通量计算和画图的问题 附脚本

guoguohh

ZP盛夏光年 发表于 2023-2-7 19:29
请问水汽通量散度怎么计算呢

上面提到的官网脚本里有

kiki07

ps =short2flt(fps->sp(:,:,:));Pa=(kg/(m*s2))这个是什么数据啊？

guoguohh

kiki07 发表于 2023-4-15 16:01
ps =short2flt(fps->sp(:,:,:));Pa=(kg/(m*s2))这个是什么数据啊？

地表气压

喵喵是我

16度的可乐 发表于 2022-4-2 17:39
dp   = dpres_plevel_Wrap(p, ps, ptop, 0) 这里dp会出现很多nan，需要修改成fillvalue=-9999

你好！我同样利用官网给出的样例进行计算，在dp计算之后，使用print("Is dp all missing ? " + all(ismissing(dp)))进行排查发现结果为True（全为缺测）。我根据你的信息添加了dp@_FillValue=-9.99E+33。但是仍旧没有解决。请问你是怎么解决的呢？

这是dp的信息（未添加缺测值为-9999之前）：
Variable: dp
Type: float
Total Size: 23210496 bytes
            5802624 values
Number of Dimensions: 4
Dimensions and sizes:        [year | 69] x [plev | 8] x [lat | 73] x [lon | 144]
Coordinates:
            year: [1951..2019]
            plev: [100000..30000]
            lat: [-90..90]
            lon: [ 0..357.5]
Number Of Attributes: 3
  _FillValue :        -9.96921e+36
  long_name :        pressure thickness
  units :        Pa
(0)        pressure thickness (Pa) : min=-9.96921e+36   max=-9.96921e+36
(0)        -----

小米辣在月球

16度的可乐 发表于 2022-4-2 17:39
dp   = dpres_plevel_Wrap(p, ps, ptop, 0) 这里dp会出现很多nan，需要修改成fillvalue=-9999

用NCEP Monthly 的 pres.mon.mean.nc算dp会出很多NAN，修改成fillvalue=-9999也不行。但是用NCEP daily 的pres.sfc.2008.nc计算dp就不会出现NAN....