IDL有求涡度，散度的函数吗

15195775117

晓考拉 发表于 2018-9-20 20:05
不好意思，我上传完才发现附件下载要积分，然后我也不会弄成不要积分的了

没事，积分不值钱，交流的机会珍贵，我下载下来明天研究

15195775117

晓考拉 发表于 2018-9-20 20:05
不好意思，我上传完才发现附件下载要积分，然后我也不会弄成不要积分的了

NECP数据我没处理过，你的数据里有以下量，不知道用到的是哪些？分别是什么含义？
PRES_3_SFC      FLOAT     = Array[360, 181]
PRES_3_MWSL     FLOAT     = Array[360, 181]
PRES_3_TRO      FLOAT     = Array[360, 181]
PRES_3_PVL      FLOAT     = Array[360, 181, 2]
PRMSL_3_MSL     FLOAT     = Array[360, 181]
HGT_3_SFC       FLOAT     = Array[360, 181]
HGT_3_0DEG      FLOAT     = Array[360, 181]
HGT_3_MWSL      FLOAT     = Array[360, 181]
HGT_3_TRO       FLOAT     = Array[360, 181]
HGT_3_ISBL      FLOAT     = Array[360, 181, 26]
HGT_3_PVL       FLOAT     = Array[360, 181, 2]
HGT_3_HTFL      FLOAT     = Array[360, 181]
TOZONE_3_EATM   FLOAT     = Array[360, 181]
TMP_3_SFC       FLOAT     = Array[360, 181]
TMP_3_MWSL      FLOAT     = Array[360, 181]
TMP_3_TRO       FLOAT     = Array[360, 181]
TMP_3_ISBL      FLOAT     = Array[360, 181, 26]
TMP_3_GPML      FLOAT     = Array[360, 181, 3]
TMP_3_HTGL      FLOAT     = Array[360, 181]
TMP_3_SIGL      FLOAT     = Array[360, 181]
TMP_3_DBLY      FLOAT     = Array[360, 181, 4]
TMP_3_SPDY      FLOAT     = Array[360, 181]
TMP_3_PVL       FLOAT     = Array[360, 181, 2]
POT_3_SIGL      FLOAT     = Array[360, 181]
GPA_3_ISBL      FLOAT     = Array[360, 181, 2]
U_GRD_3_MWSL    FLOAT     = Array[360, 181]
U_GRD_3_TRO     FLOAT     = Array[360, 181]
U_GRD_3_ISBL    FLOAT     = Array[360, 181, 26]
U_GRD_3_GPML    FLOAT     = Array[360, 181, 3]
U_GRD_3_HTGL    FLOAT     = Array[360, 181]
U_GRD_3_SIGL    FLOAT     = Array[360, 181]
U_GRD_3_SPDY    FLOAT     = Array[360, 181]
U_GRD_3_PVL     FLOAT     = Array[360, 181, 2]
V_GRD_3_MWSL    FLOAT     = Array[360, 181]
V_GRD_3_TRO     FLOAT     = Array[360, 181]
V_GRD_3_ISBL    FLOAT     = Array[360, 181, 26]
V_GRD_3_GPML    FLOAT     = Array[360, 181, 3]
V_GRD_3_HTGL    FLOAT     = Array[360, 181]
V_GRD_3_SIGL    FLOAT     = Array[360, 181]
V_GRD_3_SPDY    FLOAT     = Array[360, 181]
V_GRD_3_PVL     FLOAT     = Array[360, 181, 2]
V_VEL_3_ISBL    FLOAT     = Array[360, 181, 21]
V_VEL_3_SIGL    FLOAT     = Array[360, 181]
ABS_V_3_ISBL    FLOAT     = Array[360, 181, 26]
SPF_H_3_HTGL    FLOAT     = Array[360, 181]
SPF_H_3_SPDY    FLOAT     = Array[360, 181]
R_H_3_0DEG      FLOAT     = Array[360, 181]
R_H_3_ISBL      FLOAT     = Array[360, 181, 21]
R_H_3_HTGL      FLOAT     = Array[360, 181]
R_H_3_SIGL      FLOAT     = Array[360, 181]
R_H_3_SIGY      FLOAT     = Array[360, 181, 4]
R_H_3_SPDY      FLOAT     = Array[360, 181]
R_H_3_EATM      FLOAT     = Array[360, 181]
R_H_3_HTFL      FLOAT     = Array[360, 181]
P_WAT_3_EATM    FLOAT     = Array[360, 181]
WEASD_3_SFC     FLOAT     = Array[360, 181]
T_CDC_3_CCY     FLOAT     = Array[360, 181]
C_WAT_3_EATM    FLOAT     = Array[360, 181]
LAND_3_SFC      FLOAT     = Array[360, 181]
ICE_C_3_SFC     FLOAT     = Array[360, 181]
LFT_X_3_SFC     FLOAT     = Array[360, 181]
VW_SH_3_TRO     FLOAT     = Array[360, 181]
VW_SH_3_PVL     FLOAT     = Array[360, 181, 2]
SOILW_3_DBLY    FLOAT     = Array[360, 181, 4]
CLWMR_3_ISBL    FLOAT     = Array[360, 181, 21]
O3MR_3_ISBL     FLOAT     = Array[360, 181, 6]
CIN_3_SFC       FLOAT     = Array[360, 181]
CIN_3_SPDY      FLOAT     = Array[360, 181]
CAPE_3_SFC      FLOAT     = Array[360, 181]
CAPE_3_SPDY     FLOAT     = Array[360, 181]
HPBL_3_SFC      FLOAT     = Array[360, 181]
LV_ISBL9        LONG      = Array[6]
LV_SIGY8_L1     LONG      = Array[4]
LV_SIGY8_L0     LONG      = Array[4]
LV_ISBL7        LONG      = Array[21]
LV_ISBL6        LONG      = Array[2]
LV_DBLY5_L1     LONG      = Array[4]
LV_DBLY5_L0     LONG      = Array[4]
LV_GPML4        LONG      = Array[3]
LV_ISBL3        LONG      = Array[26]
LV_PVL2         LONG      = Array[2]
LAT_3           FLOAT     = Array[181]
LON_3           FLOAT     = Array[360]

晓考拉

15195775117 发表于 2018-9-21 09:22
NECP数据我没处理过，你的数据里有以下量，不知道用到的是哪些？分别是什么含义？
PRES_3_SFC      FLOA ...

用这两个参量
U_GRD_3_ISBL    FLOAT     = Array[360, 181, 26]
V_GRD_3_ISBL    FLOAT     = Array[360, 181, 26]
应该就是用U、V风速画散度和涡度吧，26是表示垂直方向上划分了26层

15195775117

晓考拉 发表于 2018-9-22 17:39
用这两个参量
U_GRD_3_ISBL    FLOAT     = Array[360, 181, 26]
V_GRD_3_ISBL    FLOAT     = Array[3 ...

你看下对不对

晓考拉

15195775117 发表于 2018-9-26 17:28
你看下对不对

太厉害了，这个地图怎么加上去的呢，大神能把程序发给我么

15195775117

晓考拉 发表于 2018-10-17 09:54
太厉害了，这个地图怎么加上去的呢，大神能把程序发给我么

我已经把这个案例发布成帖子了：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=68750
当然我也不能保证一定对，帖子里有详细的解说，有问题我们再讨论

晓考拉

15195775117 发表于 2018-10-17 10:18
我已经把这个案例发布成帖子了：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=68750
当然 ...

非常感谢！