ncl计算K指数问题

dh5002

计算秋季K指数,空间分布是空白图，时间序列也不对，请哪位大佬指点一下
begin
f1 = addfile("/mnt/d/text/reg/air.mon.range.nc","r")
f2 = addfile("/mnt/d/text/reg/rhum.mon.range.nc","r")

;;-----------------------------------Td850

T1 = f1->air(:,{850},:,:)
T850 = T1 - 273.16
copy_VarCoords(T1, T850)
rh1 = f2->rhum(:,:,:)
;;Td850=T850-((14.55+0.114*T850)*(1-0.01*rh1) + exp((2.5+0.007*T850)*(1-0.01*rh1))^exp(3) + (15.9+0.117*T850)*exp(1-0.01*rh1)^exp(14))
Td850 = dewtemp_trh(T850, rh1)
;;------------------------------------Td700
T2 = f1->air(:,{700},:,:)
T700 = T2 - 273.16
copy_VarCoords(T2, T700)
rh2 = f2->rhum(:,:,:)
;;Td700=T700-((14.55+0.114*T700)*(1-0.01*rh2) + exp((2.5+0.007*T700)*(1-0.01*rh2))^exp(3) + (15.9+0.117*T700)*exp(1-0.01*rh2)^exp(14))
Td700 = dewtemp_trh(T700, rh2)
;;------------------------------------T500
T3 = f1->air(:,{500},:,:)
T500 = T3 - 273.16
copy_VarCoords(T3, T500)
;;------------------------------------K指数
K=T850-T500+Td850-(T700-Td700)
;printVarSummary(K);[720] x [73] x [144]
copy_VarCoords(T1, K)
K_SON = month_to_season(K, "SON")

k_son_zhexian = dim_avg_n_Wrap(K_SON, (/1,2/))
print(k_son_zhexian)

k_son_avg = dim_avg_n_Wrap(K_SON, 0)
printMinMax(k_son_avg, 0)



这是出来的时间序列
(0)     -257.5296
(1)     -257.4906
(2)     -257.9269
(3)     -258.6134
(4)     -257.5385
(5)     -257.2753
(6)     -256.9373
(7)     -257.6821
(8)     -257.3702
(9)     -257.0979
(10)    -257.8559
(11)    -257.1031
(12)    -257.2643
(13)    -257.5393
(14)    -258.0568
(15)    -258.1897
(16)    -257.0003
(17)    -257.7199
(18)    -257.5693
(19)    -256.8394
(20)    -256.5392
(21)    -256.5257
(22)    -257.4327
(23)    -256.9323
(24)    -257.2335
(25)    -257.0083
(26)    -257.8497
(27)    -257.1943
(28)    -257.3214
(29)    -257.0381
(30)    -256.4925
(31)    -257.2253
(32)    -258.1451
(33)    -257.344
(34)    -257.1812
(35)    -256.4373
(36)    -256.2079
(37)    -256.6263
(38)    -256.5399
(39)    -256.9391
(40)    -256.4896
(41)    -255.9097
(42)    -255.8889
(43)    -255.9853
(44)    -255.7345
(45)    -255.9702
(46)    -256.0774
(47)    -256.1321
(48)    -255.7826
(49)    -256.2544
(50)    -256.1766
(51)    -255.6105
(52)    -256.1145
(53)    -255.9849
(54)    -256.0902
(55)    -254.8394
(56)    -255.5782
(57)    -256.0936
(58)    -255.5096
(59)    -255.7727

爱地理不爱物理

你这个相对湿度是哪一层的啊，为啥rh1和rh2是一样的

爱地理不爱物理

并且露点求出来还要减去273.16

dh5002

爱地理不爱物理 发表于 2021-12-30 10:11
你这个相对湿度是哪一层的啊，为啥rh1和rh2是一样的

最开始下的rh数据只有一层，后面重新下载，计算出来啦

dh5002

重新写的结果，最后的结果也不知道对不对，然后era5里面有K指数数据
begin
f1 = addfile("/mnt/d/text/reg/air.mon.range.nc","r")
f2 = addfile("/mnt/d/text/reg/rhum.mon.range.nc","r")

;;-----------------------------------Td850

T1 = f1->air(:,{850},:,:)
T850 = T1 + 273.15
copy_VarCoords(T1, T850)
rh1 = f2->rhum(:,{850},:,:)
;；Td850=T850-((14.55+0.114*T850)*(1-0.01*rh1) + exp((2.5+0.007*T850)*(1-0.01*rh1))^exp(3) + (15.9+0.117*T850)*exp(1-0.01*rh1)^exp(14))
Td850 = dewtemp_trh(T850, rh1) - 273.15
copy_VarCoords(T1, Td850)
;;------------------------------------Td700
T2 = f1->air(:,{700},:,:)
T700 = T2 + 273.15
copy_VarCoords(T2, T700)
rh2 = f2->rhum(:,{700},:,:)
;；Td700=T700-((14.55+0.114*T700)*(1-0.01*rh2) + exp((2.5+0.007*T700)*(1-0.01*rh2))^exp(3) + (15.9+0.117*T700)*exp(1-0.01*rh2)^exp(14))
Td700 = dewtemp_trh(T700, rh2) - 273.15
copy_VarCoords(T2, Td700)
;;------------------------------------T500
T3 = f1->air(:,{500},:,:)

;;------------------------------------K指数
K=T1-T3+Td850-(T2-Td700)