气象ts评分_TS平分
盟
8 G; e$
z. a3 U( Q
如果给定的数据是各次预报和观测的值，和
ets
评分的阈值，则用如下代码，其中多出一段判断命中、空报、漏
报和正确否定次数的代码，对
TS
评分也一样。
气象
数值
模式
气象数
值模式
数值模
式
气象
模式
海洋模
式
气候
模式
,WRF
MM5
,Grapes,HYCOM,LI
COM,
同化
集
合预报
模式联
盟
模
式论坛
) p' g4 [3 K7
B! E$
M+ K) x
subroutine QETS(ets,yuezhi,N,forc,obse)
气象
数值
模
式
气象
数值模
式
数值
模式
气象模
式
海洋
模式
气
候模式
,W
RF
,MM5,Grapes,HYCO
M,LICO
M,
同化
集合预
报
模
式联盟
模式论
坛
. ]' d& }6 t+ W9 {
implicit none
integer ::H,F,M,C !
命中，空报，漏报，正确否定
www.mn
muc.org; U( A; Z G+ U# a" Q7 l9 G# Q
integer ::Fo,Ob
，
N !
发布警报次数，事件发生次数，总预报次数
Meteorological Numerical Model Union of China (MNMUC)
中
国气象
数值模
式联盟
l6 e+ M6 b8
integer ::i
www.mnmuc.org$ K8 B, _2 E5 t8 ]
real ::forc(N),obse(N)!
预报和观测的降水序列
Meteorological Numerical Model Union of China (MN
MUC)
中
国气象
数值模
式联盟
4 s0 L! ^' d. ]2 n1
J9 w
real ::yuezhi,ets !
降水事件的评定阈值，
ets
得分
Meteorological Numerical Model Un
ion of
China (MN
MUC)
中
国气象
数值模
式联
盟
7 Z1
R3 H4 G!
~. V' G3
) Q
H=0
气象
数值
模式
气象数
值模式
数值
模式
气
象模式
海洋模
式
气
候模式
,W
RF
,MM5,Grapes,HYCO
M,LICOM,
同化
集合预
报
模式
联盟
模
式论
坛
7 g; A+ Y8 ?+ Y6 H
F=0
Meteorological Numerical Model Union of China (MNMU
C)
中
国气象
数值模
式联盟
5 v7 Y5 J& a- P. V, _% @
M=0
Meteorological Numerical Model Union
of China (MNMU
C)
中国
气象数
值模式
联盟
. j1 Z! I. `*
O) n. r, H6 ~. U
C=0
do i=1,N
if(forc(i)>=yuezhi.and.obse(i)>=yuezhi)H=H+1
if(forc(i)>=yuezhi.and.obse(i)
if(forc(i)=yuezhi)M=M+1
气象
数值
模式
气
象数值
模式
数
值模式
气象模
式
海洋
模式
气候模
式
,WRF
,MM5,Grapes,HYCO
M,LICO
M,
同
化
集合
预报
模
式联盟
模式论
坛
' I. S9
k2 C8
_+ O. a
if(forc(i)
气象
数值
模式
气
象数值
模式
数
值模式
气象
模式
海
洋模式
气候模
式
,WRF
MM5,Grapes,HY
COM,LI
COM,
同
化
集
合预报
模式联
盟
模式
论坛
% s" J0 x
/ s. x
5 n/
M1 ~' j
enddo
/ D! m4 K/ P& ? ~/ L
5 b( Q
Fo=H+F
j3 A% m" [0 p1 `. r2 V4
Ob=H+M
N=H+F+M+C
ets=(H-Fo*Ob/N)/(H+F+M-Fo*Ob/N)
return
气象
数值
模式
气
象数值
模式
数值模
式
气象
模式
海
洋模式
气候模
式
,WRF
MM5,Grapes,HY
COM,LI
COM,
同
化
集
合预报
模式联
盟
模式
论坛
, j8 B!
z2
C8
end subroutine
Meteoro
logical Numerical Model Un
ion of
China (MN
MUC)
中国气
象数值
模式联
盟
8 G; e$
z. a3 U( Q