Versions Compared

Old Version 6

changes.mady.by.user Martin Janousek

Saved on

compared with

New Version 7

changes.mady.by.user Martin Janousek

Saved on

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

parameterstepsforecast values preprocessingscores

2-m temperature

0,3,6,...,69,72, 78,84,90,...,240

T2m extracted from model forecast field and then using the recent model orography reduced to the observation station level by

Code Block
T2m_reduced = T2m_model + (elev_model - elev_station)*0.0065

(assuming T2m are in K and elev_* in meters)

ME, MAE, RMSE

2-m dewpoint

0,3,6,...,69,72, 78,84,90,...,240extracted from model forecast fieldME, MAE, RMSE

2-m relative humidity

0,3,6,...,69,72, 78,84,90,...,240

computed from T2m and TD2m

Code Block
def magnusFormula(T):
    log10_e_w=-2937.4/T-4.9283*numpy.log10(T)+25.5518
    return 10.**log10_e_w

e=magnusFormula(Tdew)
e_s=magnusFormula(T)
RH=e/e_s


ME, MAE, RMSE

10-m wind direction

0,3,6,...,69,72, 78,84,90,...,240

computed from 10-m wind components

Code Block
dd10m = 180./pi * arctan2(u10m,v10m) - 180.

where the observed 10-m wind speed is less than 3m/s do not compute dd10m scores; however, that observation is still included into the total count of observations in the current month (which should reach at least 90% for the station score to be included in the report)

 

ME, MAE, RMSE

fc-ob differences computed as follows:

Code Block
diff = fc - ob
where (diff>180): diff-=360
where (diff<-180): diff+=360


10-m wind speed

0,3,6,...,69,72, 78,84,90,...,240

computed from 10-m wind components

Code Block
ff10m = sqrt(u10m**2+v10m**2)


ME, MAE, RMSE

CT (5,10,15)

total cloud cover

0,3,6,...,69,72, 78,84,90,...,240extracted from model forecast field (0.-1.) then converted to octas

ME, MAE, RMSE

CT (2,6) (warning) inconsistent 7)

(warning) please note incorrect thresholds (2,6) have been written in the earlier version of this document (corrected on 2 April 2019 with WMO guidelines (it will be corrected here soon)

prior CT computations the forecast values are rounded to the nearest octa

24-h precipitation24,48,72,...,240

extracted from model forecast as a difference of the precipitation accumulation fields, e.g.

tp24(step=120) = accum_precip(step=120) - accum_precip(step=96)

CT (1,10,50)
6-h precipitation6,12,18,...,240

extracted from model forecast as a difference of the precipitation accumulation fields, e.g.

tp6(step=120) = accum_precip(step=120) - accum_precip(step=112)

CT (1,5,25)

...