1 ! The following parameters are used as flags to indicate either a
2 ! MISSING data value (real number) in the "measurement" TYPE or an
3 ! "end of data" flag (two fields with end_data flag mean end of
4 ! measurements for one observation). The undefined values are when
5 ! a GTS observation is read but determined to be in error (not according
6 ! to the code). Geopotential height is referenced as a positive value
7 ! undefined, which explains the duplicate.
9 INTEGER , PARAMETER :: undefined1 = 888888
10 REAL , PARAMETER :: undefined1_r = 888888.
11 INTEGER , PARAMETER :: undefined2 = -888888
12 REAL , PARAMETER :: undefined2_r = -888888.
13 INTEGER , PARAMETER :: missing = -888888
14 REAL , PARAMETER :: missing_r = -888888.
15 INTEGER , PARAMETER :: end_data = -777777
16 REAL , PARAMETER :: end_data_r = -777777.
18 ! The cutoff height for Aircraft data:
19 ! h > aircraft_cut, to retrive the pressure from standard atmosphere;
20 ! h <= aircraft_cut, to keep Aircraft observed height.
22 REAL , PARAMETER :: aircraft_cut = 3000.
24 ! Define various ways for bad data to be flagged. Each of these is
25 ! 2^n, so that the errors may be summed.
27 INTEGER, PARAMETER :: &
29 outside_of_domain = 2 ** 18, & ! Data outside horizontal domain
30 ! or time window, data set to missing_r
31 wrong_direction = 2 ** 17, & ! Wind vector direction <0 or> 360
32 ! => direction set to missing_r
33 negative_spd = 2 ** 16, & ! Wind vector norm is negative
34 ! => norm set to missing_r
35 zero_spd = 2 ** 15, & ! Wind vector norm is zero
36 ! => norm set to missing_r
37 wrong_wind_data = 2 ** 14, & ! Spike in wind profile
38 ! =>direction and norm set to missing_r
39 zero_t_td = 2 ** 13, & ! t or td = 0 => t or td, rh and qv
40 ! are set to missing_r,
41 t_fail_supa_inver = 2 ** 12, & ! t, td, rh and qv fail vertical
42 ! consistency =>all are set to missing_r
43 wrong_t_sign = 2 ** 11, & ! superadiabatic and inversion
44 ! in temperature curve => no action
45 above_model_lid = 2 ** 10, & ! heigh above model lid
46 ! =>Datum can be removed
47 reference_atmosphere = 2 ** 7, & ! Missing h, p or t
48 ! =>Datum interpolated from standard atm
49 from_background = 2 ** 6, & ! Missing h, p or t
50 ! =>Datum interpolated from model ref st
51 convective_adjustment = 2 ** 2, & ! convective adjustement check
52 ! =>Apply correction on t, td, rh and qv
53 Hydrostatic_recover = 3 , & ! Height from hydrostaic assumption with
54 ! the OBS data calibration
55 surface_correction = 2 ** 1, & ! Surface datum
56 ! =>Apply correction on datum
57 Reference_OBS_scaled = 2 ** 0 ! Height from reference state with OBS
59 ! =>Apply correction on datum