Diagnostic.cpp 9.43 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
/**
 *
 * Copyright (C) 2009-2014 INRA.
 * Copyright (C) 2009-2014 CETIOM.
 *
 *
 * TODO licence
 */

// @@tagdynamic@@
// @@tagdepends: vle.discrete-time @@endtagdepends


#include <math.h>
#include <vle/DiscreteTime.hpp>
#include <sunflo_utils.hpp>
#include <Phenologie.hpp>
#include <ParametresPlante.hpp>
#include <ParametresVariete.hpp>
20
#include <ParametresRendement.hpp>
21
22
23
24
25
26
27
28
29
30

namespace sunflo {

using namespace vle::discrete_time;

class Diagnostic : public DiscreteTimeDyn
{
public :
    ParametresVariete pv;
    ParametresPlante pp;
31
    ParametresRendement pr;
32
33
34
35
36
    bool first_compute;

    /*Sync*/ Var TT_A2;
    /*Sync*/ Var ETRETM;

37
38
39
40
41
42
43
////////////////////////////////////////////////////////////////////////
/////////////////////////////REMI ADDS//////////////////////////////////
////////////////////////////////////////////////////////////////////////

    /*Sync*/ Var PhasePhenoPlante;

////////////////////////////////////////////////////////////////////////
44
////////////////////////////////////////////////////////////////////////
45
46
47
48
49
50
51
52
53
54
55
56
57

    /// # Jours de stress jusqu'a la floraison (estimationTTentreSemisEtLevee_casPhaseSemisLeveeSimulee - F1)
    /// # (pas utilisé dans le modèle, sortie de Diagnostic)
    Var ISH1;

    /// # Jours de stress Autour de la floraison (F1 + 350°Cd)
    /// # (pas utilisé dans le modèle, sortie de Diagnostic)
    Var ISH2;

    /// # Jours de stress après la floraison (F1+350 °Cd – M3)
    /// # (pas utilisé dans le modèle, sortie de Diagnostic)
    Var ISH3;

58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82

////////////////////////////////////////////////////////////////////////
/////////////////////////////REMI ADDS//////////////////////////////////
////////////////////////////////////////////////////////////////////////

    /// # Jours entre maturation et récolte.
    Var D_MH;
    /// # Jours entre semis et germination.
    Var D_SE;

    /// # Somme ETPET entre maturation et récolte.
    Var SET_MH;

    /// # Moyenne ETPET entre maturation et récolte.
    Var MET_MH;

    /// # Somme ETPET entre germinatioan et initiation florale.
    Var SET_EF;

    /// # Moyenne ETPET entre germinatioan et initiation florale.
    Var MET_EF;

    // Thermal stress entre germinatioan et initiation florale.
    Var NHT28_EF;

RemiM's avatar
RemiM committed
83
    // Thermal stress during whole season
84
    Var NHT28;
RemiM's avatar
RemiM committed
85
86
87
88
89

    // Radiations
    Var RG;
    // Sum of radiations during flowering
    Var SGR_FM;
90
91
92
////////////////////////////////////////////////////////////////////////
////////////////////////////////////////////////////////////////////////

93
94
95
96
97
    Diagnostic(const vle::devs::DynamicsInit& model,
            const vle::devs::InitEventList& events) :
                DiscreteTimeDyn(model, events)
    {
        pv.initialiser(events);
98
        pr.initialiser(events);
99
100
101
102
103
104
105
        pp.initialiser(events);

        TT_A2.init(this, "TT_A2", events);
        ETRETM.init(this, "ETRETM", events);
        ISH1.init(this, "ISH1" , events);
        ISH2.init(this, "ISH2" , events);
        ISH3.init(this, "ISH3" , events);
106

107
108
109
110
111
112
113
114
115
116
117
118
////////////////////////////////////////////////////////////////////////	
/////////////////////////////REMI ADDS//////////////////////////////////
	D_MH.init(this, "D_MH" , events);
	D_SE.init(this, "D_SE", events);

        SET_MH.init(this, "SET_MH" , events);
        MET_MH.init(this, "MET_MH" , events);

        SET_EF.init(this, "SET_MH" , events);
        MET_EF.init(this, "MET_MH" , events);

	NHT28_EF.init(this, "NHT28_EF.init", events);
RemiM's avatar
RemiM committed
119
120
121
122
123


	SGR_FM.init(this, "SGR_FM", events);
	RG.init(this, "RG", events);

124
125
126
127
        PhasePhenoPlante.init(this, "PhasePhenoPlante", events);
////////////////////////////////////////////////////////////////////////
////////////////////////////////////////////////////////////////////////

128
        first_compute = true;
129
130
131
132
133
134
135
136
137
138
139
    }

    virtual ~Diagnostic() { }

    virtual void compute(const vle::devs::Time& /*time*/ )
    {
        if (first_compute) {
            first_compute = false;
            ISH1 = 0.0; //codeMMpourMemo Initial Value = 0.0
            ISH2 = 0.0; //codeMMpourMemo Initial Value = 0.0
            ISH3 = 0.0; //codeMMpourMemo Initial Value = 0.0
140
141
142
143
144
145
146
147
148
149

/////////////////////////////REMI ADDS//////////////////////////////////
	    D_MH = 0.0; 
	    D_SE = 0.0; 
	    SET_MH = 0.0;
	    MET_MH = 0.0;
	    SET_EF = 0.0;
	    MET_EF = 0.0;
	    NHT_28 = 0.0;
            NHT28_MH = 0.0;
RemiM's avatar
RemiM committed
150
	    SGR_FM = 0.0;
151
////////////////////////////////////////////////////////////////////////
152
153
154
        } else {
            {
                double ddt = 0.0;
155
                if ( ETRETM() > pr.SeuilETRETM ){
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
                    ddt = 0.0;
                } else if ( TT_A2() > pv.date_TT_F1 ){
                    ddt = 0.0;
                } else if ( TT_A2() <= 0.0 ){
                    // traduit condition TT_A2 <= 0
                    ddt = 0.0;
                } else {
                    ddt = 1.0;
                }
                ISH1 = ISH1(-1) + ddt;
            }

            {
                double ddt = 0.0;

171
                if ( ETRETM() > pr.SeuilETRETM ){
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
                    ddt = 0.0;
                } else if ( TT_A2() > (pv.date_TT_F1 + pp.date_TT_F1M0) ){
                    ddt = 0.0;
                } else if ( TT_A2() < pv.date_TT_F1 ){
                    ddt = 0.0;
                } else {
                    ddt = 1.0;
                }
                ISH2 = ISH2(-1) + ddt;
            }


            {
                double ddt = 0.0;

187
                if ( ETRETM() > pr.SeuilETRETM ){
188
189
190
191
192
193
194
195
196
197
                    ddt = 0.0;
                } else if ( TT_A2() > pv.date_TT_M3 ){
                    ddt = 0.0;
                } else if ( TT_A2() < (pv.date_TT_F1 + pp.date_TT_F1M0) ){
                    ddt = 0.0;
                } else {
                    ddt = 1.0;
                }
                ISH3 = ISH3(-1) + ddt;
            }
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296

////////////////////////////////////////////////////////////////////////
/////////////////////////////REMI ADDS//////////////////////////////////
////////////////////////////////////////////////////////////////////////

            {
                double ddt = 0.0;
                if ( ( PhasePhenoPlante() < PHASEPHENOPLANTE_MATURATION )
                        // traduit condition TT_A2 < date_TT_M0
                        || ( PhasePhenoPlante() < PHASEPHENOPLANTE_JUVENILE )
                        || ( PhasePhenoPlante() >= PHASEPHENOPLANTE_RECOLTEE )){
                    // traduit condition TT_A2 = 0
                    ddt = 0.0;
                } else {
                    ddt = 1;
                }
                D_MH = D_MH(-1) + ddt;
            }

            {
                double ddt = 0.0;
                if ( PhasePhenoPlante() == PHASEPHENOPLANTE_JUVENILE){
                    ddt = 0.0;
                } else {
                    ddt = 1;
                }
                D_EF = D_EF(-1) + ddt;
            }

            {
                double ddt = 0.0;
                if ( PhasePhenoPlante() == PHASEPHENOPLANTE_GERMINATION ){
                    // traduit condition TT_A2 = 0
                    ddt = 0.0;
                } else {
                    ddt = 1;
                }
                D_SE = D_SE(-1) + ddt;
            }

            {
                double ddt = 0.0;
                if ( ( PhasePhenoPlante() < PHASEPHENOPLANTE_MATURATION )
                        // traduit condition TT_A2 < date_TT_M0
                        || ( PhasePhenoPlante() < PHASEPHENOPLANTE_JUVENILE )
                        || ( PhasePhenoPlante() >= PHASEPHENOPLANTE_RECOLTEE ) ){
                    // traduit condition TT_A2 = 0
                    ddt = 0.0;
                } else {
                    ddt = ETRETM();
                }
                SET_MH = SET_MH(-1) + ddt;
            }

            // calcul de MET_MH
            MET_MH = SET_MH() / D_MH();

            {
                double ddt = 0.0;
                if ( ( PhasePhenoPlante() == PHASEPHENOPLANTE_JUVENILE ) ){
                    // traduit condition TT_A2 = 0
                    ddt = 0.0;
                } else {
                    ddt = ETRETM();
                }
                SET_EF = SET_EF(-1) + ddt;
            }

            // calcul de MET_EF
            MET_EF = SET_EF() / D_EF();


            // calcul de NHT28_EF
            {
                double ddt = 0.0;
                if ( Tx() < 28.0 ){
                    ddt = 0.0;
                } else if (( PhasePhenoPlante() != PHASEPHENOPLANTE_JUVENILE )){
                    ddt = 0.0;
                } else {
                    ddt = 1.0;
                }

                NHT28_EF = NHT28_EF(-1) + ddt;
            }

            // calcul de NHT28
            {
                double ddt = 0.0;
                if ( Tx() < 28.0 ){
                    ddt = 0.0;
                } else if (( PhasePhenoPlante() >= PHASEPHENOPLANTE_RECOLTEE) || PhasePhenoPlante() == PHASEPHENOPLANTE_NONSEMEE){
                    ddt = 0.0;
                } else {
                    ddt = 1.0;
                }

                NHT28 = NHT28(-1) + ddt;
            }
RemiM's avatar
RemiM committed
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311

            // calcul de SGR_FM
            {
                double ddt = 0.0;
                if ( PhasePhenoPlante() != PHASEPHENOPLANTE_FLORAISON){
                    ddt = 0.0;
                } else {
                    ddt = RG();
                }

                SGR_FM = SGR_FM(-1) + ddt;
            }



312
313
////////////////////////////////////////////////////////////////////////
////////////////////////////////////////////////////////////////////////
314
315
316
317
318
319
320
321
        }
    }
};

} // namespace sunflo

DECLARE_DYNAMICS(sunflo::Diagnostic);