【已解决】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....

一只秃头媛

楼主，冒昧打扰，想咨询一下你的代码中ps =short2flt(fps->sp(:,:,:));Pa=(kg/(m*s2))的含义，sp是每一层的气压么，还是surface pressure.

guoguohh

一只秃头媛 发表于 2024-12-5 16:33
楼主，冒昧打扰，想咨询一下你的代码中ps =short2flt(fps->sp(:,:,:));Pa=(kg/(m*s2))的含义，sp是每一层的 ...

是surface pressure

一只秃头媛

guoguohh 发表于 2024-12-6 16:09
是surface pressure

非常感谢您的回复。

FleurChim

楼主，在使用这个代码时候发现前面可以输出正常值，但在计算uq_dpg = uq*dpg
时就出现了输出值全为Nan的问题，不太清楚这种是什么原因导致的
相关输出：(0)     pressure thickness (Pa) : min=1   max=2149
(0)
(0)     Layer Mass Weighting (kg/m2) : min=0.101972   max=219.137
(0)
(0)     Zonal Moisture Flux [uq] ([m s**-1][g/kg]) : min=-244.222   max=241.68
(0)
(0)     min=nan   max=nan

FleurChim

FleurChim 发表于 2024-12-11 21:52
楼主，在使用这个代码时候发现前面可以输出正常值，但在计算uq_dpg = uq*dpg
时就出现了输出值全为Nan的问 ...

找到原因了，是数据的问题