Versions Compared

Old Version 5

changes.mady.by.user Martin Janousek

Saved on

compared with

New Version 6

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 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)

...