Changeset 6931f61 in flex_extract.git for For_developers/Sphinx/Source/Documentation/disagg.rst
- Timestamp:
- Jul 8, 2020, 10:00:37 PM (4 years ago)
- Branches:
- master, ctbto, dev
- Children:
- e0005c9
- Parents:
- 8a53342
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
For_developers/Sphinx/Source/Documentation/disagg.rst
r8a53342 r6931f61 13 13 14 14 .. csv-table:: Flux fields 15 :header: "Short Name", "Name", "Units", " Interpolation Type"15 :header: "Short Name", "Name", "Units", "Disaggregation" 16 16 :align: center 17 17 :widths: 5,15,5,10 18 18 19 LSP, "large-scale precipitation", "m", " modified linear interpolation"20 CP, "convective precipitation", "m", " modified linear interpolation"21 SSHF, "surface sensible heat flux", "J m :math:`^{-2}`", "bicubic interpolation"22 EWSS, "eastward turbulent surface stress", "N m :math:`^{-2}` s", "bicubic interpolation"23 NSSS, "northward turbulent surface stress", "N m :math:`^{-2}` s", "bicubic interpolation"24 SSR, "surface net solar radiation", "J m :math:`^{-2}`", "bicubic interpolation"19 LSP, "large-scale precipitation", "m", "precipitation " 20 CP, "convective precipitation", "m", "precipitation" 21 SSHF, "surface sensible heat flux", "J m :math:`^{-2}`", "flux" 22 EWSS, "eastward turbulent surface stress", "N m :math:`^{-2}` s", "flux " 23 NSSS, "northward turbulent surface stress", "N m :math:`^{-2}` s", "flux " 24 SSR, "surface net solar radiation", "J m :math:`^{-2}`", "flux " 25 25 26 26 … … 43 43 44 44 In the old versions (below 7.1), a relatively simple method processes the precipitation fields in a way that is consistent with the linear interpolation between times where input fields are available that is applied in ``FLEXPART`` for all variables. 45 This scheme (from Paul James) at first divides the accumulated values by the number of hours (i.e., 3 or 6). ???45 This scheme (from Paul James) at first divides the accumulated values by the number of hours (i.e., 3 or 6). 46 46 The best option for disaggregation, which was realised, is conservation within the interval under consideration plus the two adjacent ones. 47 47 Unfortunately, this leads to undesired temporal smoothing of the precipitation time series – maxima are damped and minima are raised. … … 161 161 162 162 The accumulated values for the other variables are first divided by the number of hours and 163 then interpolated to the exact times using a bicubic interpolation which conserves the integrals of the fluxes within each timespan. 164 Disaggregation is done for four adjacent time intervals (:math:`p_a, p_b, p_c, p_d`) which produces a new, disaggregated value that is the output at the central point of the four adjacent time intervals. 163 then disaggregates to the exact times and conserves the integrals of the fluxes within each timespan. 164 Disaggregation is done for four adjacent time intervals (:math:`F_0, F_1, F_2, F_3`) which produces a new, disaggregated value that is the output at the central point of the four adjacent time intervals. 165 This new point :math:`F` is used for linear interpolation of the complete timeseries afterwards. 165 166 166 167 .. math:: 168 167 169 168 p_a &= (a_3 - a_0 + 3. * (a_1 - a_2)) / 6. \\ 169 p_b &= (a_2 + a_0) / 2. - a_1 - 9. * p_a / 2. \\ 170 p_c &= a_1 - a_0 - 7. * p_a / 2. - 2. * p_b \\ 171 p_d &= a_0 - p_a / 4. - p_b / 3. - p_c / 2. 172 173 This new point :math:`p` is used for linear interpolation of the complete timeseries afterwards. 174 175 .. math:: 176 177 p = 8. * p_a + 4. * p_b + 2. * p_c + p_d 170 F = -\frac{1}{12}F_0 + \frac{7}{12}F_1 + \frac{7}{12}F_2 -\frac{1}{12}F_3 178 171 179 172
Note: See TracChangeset
for help on using the changeset viewer.