**Code and Data description**: Code was implemented in Wolfram Mathematica. Functions are written in Wolfram Language Package files (.wl) and they are executed in Notebook files (.nb). Folder **code**: **L05IIIS.wl** 1. **Filter**: function for expressing the variables (Eqs. (A10)-(A12)). 2. **L05IIIS**: L05-3 system (Eq. (A9)). **L05III0S.wl** 1. **Filter**: function for expressing the variables (Eqs. (A10)-(A12)). 2. **L05III0S**: L05-3 system (Eq. (A9)) with initial error. **L05IIS.wl** 1. **Filter**: function for expressing the variables (Eqs. (A6)-(A7)). 2. **L05IIS**: L05-2 system (Eq. (A8)). **L05II0S.wl** 1. **Filter**: function for expressing the variables (Eqs. (A6)-(A7)). 2. **L05II0S**: L05-2 system (Eq. (A8)) with initial error. **L05IS.wl** 1. **L05IS**: L05-1 system (Eq. (A1)). **AnsambleL05III.wl** 1. **Ansamble1IIIS**: Calculation of "reality" of L05-3 system (fourth-order Runge-Kutta method). 2. **Ansamble2IIIS**: Calculation of "prediction" of L05-3 system (fourth-order Runge-Kutta method). 3. **AproxRustCasIII**: Numerical solution of Eq. (4) (fourth-order Runge-Kutta method). 4. **AproxRustCasBezIII**: Numerical solution of Eq. (3) (fourth-order Runge-Kutta method). 5. **AproxRustCasLnIII**: Numerical solution of Eq. (6) (fourth-order Runge-Kutta method). 6. **AproxRustCasLnbezEIII**: Numerical solution of Eq. (5) (fourth-order Runge-Kutta method). 7. **AproxRustCasKvIII**: Numerical solution of Eq. (2) (fourth-order Runge-Kutta method). 8. **AproxRustCasExpIII**: Numerical solution of Eq. (1) (fourth-order Runge-Kutta method). **AnsambleL05II.wl** 1. **Ansamble1IIS**: Calculation of "reality" of L05-2 system (fourth-order Runge-Kutta method). 2. **Ansamble2IIS**: Calculation of "prediction" of L05-2 system (fourth-order Runge-Kutta method). **AnsambleL05I.wl** 1. **Ansamble1IS**: Calculation of "reality" of L05-1 system (fourth-order Runge-Kutta method). 2. **Ansamble2IS**: Calculation of "prediction" of L05-1 system (fourth-order Runge-Kutta method). 3. **EMEXP31**: Calculation of Fig. 12. 4. **OpakEMEXP31**: Calculation of Fig. 12. 5. **EMEXP31new**: Calculation of Fig. 12. 6. **RustCasI**: Calculation of E(t). 7. **RustDerI**: Calculation of dE/dt(E). 8. **RustDerCasI**: Calculation of dE/dt(t). 9. **AproxRustDerI**: Eq. (2). 10. **AproxRustDer2I**: Eq. (1). **Ansamble.wl** 1. **AproxRustDerModel33**: Eq. (6) with Lamda = 0.33 1/day. 2. **AproxRustDerModel2**: Eq. (6). 3. **AproxRustDer**: Eq. (4). **ECMWF.wl** 1. **RustDerECMWFPP**: Calculation of Fig. 16b. 2. **RustDerECMWFPPM**: Calculation of Fig. 16a. 3. **RustGraf**: Calculation of Fig. 15. 4. **AproxGraf**: Eq. (5). 5. **AproxDerECMWFPPModelE35**: Eq. (6). 6. **AproxDerECMWFPPModel**: Eq. (5). 7. **AproxDerECMWFPPModelEM35**: Eq. (6) with lamda = 0.35 1/day. 8. **AproxDerECMWFPPModelM25**: Eq. (5) with lamda = 0.25 1/day. Folder **Data/ECMWF**: 1. **ECMWFPPjg.wl**: data set (Fig. 16b). 2. **ECMWFPPjg.wl**: data set (Fig. 16a). 3. **AlphaQEI.wl**: data set (Fig. 15a, black curve). 4. **AlphaREI.wl**: data set (Fig. 15a, blue curve). 5. **AlphaREM.wl**: data set (Fig. 15a, red curve). 6. **BethaQEI.wl**: data set (Fig. 15b , red curve). 7. **BethaQEM**:data set (Fig. 15b, blue curve). Folder **Data/L05**: 1. **PPI360.wl**: data set (Fig. 1). 2. **PPII360.wl**: data set (Fig. 2). 3. **PPIII360.wl**: data set (Fig. 3). 4. **E0IIaI360.wl**: data set (Fig. 8, full black curve). 5. **E0IIIaI360.wl**: data set (Fig. 9, full black curve). 6. **EI360.wl**: data set (Fig. 5). 7. **EII360.wl**: data set (Fig. 6). 8. **EIIaI360.wl**: data set (Fig. 8, dashed black curve). 9. **EIIaI36002.wl**: data set (Fig. 8, dot-dashed black curve). 10. **EIII360.wl**: data set (Fig. 7). 11. **EIIIaI360.wl**: data set (Fig. 9, dashed black curve). 12. **EIIIaI36002.wl**: data set (Fig. 9, dot-dashed black curve). 13. **EMII360aI360.wl**: data set (Fig. 11a, 13a, 14a). 14. **EMIII360aI360.wl**: data set (Fig. 11b, 13b, 14b). Folder **Figs**: Procedures for displaying Fig. 1 - 16.