source: XIOS/trunk/src/calendar_util.cpp @ 424

Last change on this file since 424 was 424, checked in by ymipsl, 11 years ago

Add new timestep unit : ts, for duration

YM

File size: 7.2 KB
Line 
1#include "calendar_util.hpp"
2
3namespace xios
4{
5      /// ////////////////////// Définitions ////////////////////// ///
6
7      CDuration operator*(const double & scal, const CDuration & ddr)
8      { return (ddr * scal); }
9
10      CDuration operator-(const CDuration & ddr , const CDuration & dr)
11      {
12         CDuration dur(ddr);
13         dur.year -= dr.year;  dur.month  -= dr.month ; dur.day    -= dr.day;
14         dur.hour -= dr.hour;  dur.minute -= dr.minute; dur.second -= dr.second; dur.timestep -= dr.timestep;
15         return (dur);
16      }
17
18      CDuration operator+(const CDuration & ddr , const CDuration & dr)
19      {
20         CDuration dur(ddr);
21         dur.year += dr.year;  dur.month  += dr.month ; dur.day    += dr.day;
22         dur.hour += dr.hour;  dur.minute += dr.minute; dur.second += dr.second; dur.timestep += dr.timestep;
23         return (dur);
24      }
25
26      CDuration operator*(const CDuration & ddr , const double & scal)
27      {
28         CDuration dur(ddr);
29         dur.year *= scal;  dur.month  *= scal; dur.day    *= scal;
30         dur.hour *= scal;  dur.minute *= scal; dur.second *= scal; dur.timestep *= scal;
31         return (dur);
32      }
33
34      CDuration operator-(const CDuration & ddr)
35      {
36         CDuration dur(ddr);
37         dur.year   = -dur.year;
38         dur.month  = -dur.month;
39         dur.day    = -dur.day;
40         dur.hour   = -dur.hour;
41         dur.minute = -dur.minute;
42         dur.second = -dur.second;
43         dur.second = -dur.timestep;
44         return (dur);
45      }
46
47      //-----------------------------------------------------------------
48
49      CDate operator+(const CDate & dt, const CDuration & dr)
50      {
51         CDuration drr (dr);
52         int year = 0, month = 0, day = 0, hour = 0, minute = 0, second = 0;
53         const CCalendar & c = dt.getRelCalendar();
54         
55         drr.timestep=0 ;       
56         drr=drr+dr.timestep*dt.getRelCalendar().getTimeStep() ;
57         
58         drr.resolve(dt.getRelCalendar());
59
60         // Ajustement des minutes par rapport aux secondes.
61         second += dt.getSecond() + drr.second;
62         if (second <  0) { minute --; second += c.getMinuteLength(); }
63         if (second >= c.getMinuteLength()) { minute ++; second -= c.getMinuteLength(); }
64
65         // Ajustement des heures en fonction des minutes.
66         minute += dt.getMinute() + drr.minute;
67         if (minute < 0) { hour --; minute += c.getHourLength(); }
68         if (minute >= c.getHourLength()) { hour ++; minute -= c.getHourLength(); }
69
70         // Ajustement des jours en fonction des heures.
71         hour += dt.getHour() + drr.hour;
72         if (hour <  0) { drr.day --; hour += c.getDayLength(); }
73         if (hour >= c.getDayLength()) { drr.day ++; hour -= c.getDayLength(); }
74
75         // Ajustement des mois en fonction des jours.
76         CDate dtt(dt); 
77         drr.day+=dtt.getDay()-1 ;
78         dtt.setDay(1) ;         
79
80         if ( drr.day >= 0 )
81         {
82           for(; c.getMonthLength(dtt) <= drr.day; dtt.addMonth (1))
83           { drr.day -=  c.getMonthLength(dtt); drr.month += 1 ; }
84
85           day = drr.day+1 ;
86         }
87         else 
88         {
89           dtt.addMonth(-1) ;
90           for(; c.getMonthLength(dtt) < -drr.day; dtt.addMonth (-1))
91           { drr.day+=c.getMonthLength(dtt) ; drr.month-=1 ; }
92           day=c.getMonthLength(dtt)+drr.day+1 ;
93         }
94           
95/*
96         if (day <  0) { drr.month --; day += c.getMonthLength(dtt); }
97         if (day > c.getMonthLength(dtt)) { drr.month ++; day -= c.getMonthLength(dtt); } // << ProblÚme ici
98         if (day == 0){ day = c.getMonthLength(dtt); drr.month --; }
99*/         
100         drr.resolve(dt.getRelCalendar());
101
102         // Ajustement des années en fonction des mois.
103         month += dt.getMonth() + drr.month;
104         if (month <  0) { drr.year --; month += c.getYearLength(); }
105         if (month >  c.getYearLength()) { drr.year ++; month -= c.getYearLength(); }
106         if (month == 0){ month = c.getYearLength(); drr.year--; }
107
108         year += dt.getYear() + drr.year;
109
110         return (CDate(dt.getRelCalendar(), year, month, day, hour, minute, second));
111      }
112
113      CDate operator-(const CDate & dt, const CDuration & dr) { return (dt + (-dr)); }
114
115      //-----------------------------------------------------------------
116
117      CDuration operator-(const CDate & dt0, const CDate & dt1)
118      {
119         // TODO :: Vérifier que les deux dates (dt0 et dt1) ont une référence vers le même calendrier.
120         CDuration dur =
121         { dt0.getYear() - dt1.getYear(), dt0.getMonth()  - dt1.getMonth() , dt0.getDay()    - dt1.getDay(),
122           dt0.getHour() - dt1.getHour(), dt0.getMinute() - dt1.getMinute(), dt0.getSecond() - dt1.getSecond() };
123         return (dur.resolve(dt0.getRelCalendar()));
124      }
125
126      //-----------------------------------------------------------------
127
128      /// Les opérateurs de comparaison. (Non testés pour le moment)
129      bool operator==(const CDate& dt0, const CDate& dt1)
130      {
131         // TODO :: Vérifier que les deux dates (dt0 et dt1) ont une référence vers le même calendrier.
132         return ((dt0.getYear() == dt1.getYear()) && (dt0.getMonth()  == dt1.getMonth())  && (dt1.getDay()    == dt0.getDay()) &&
133                 (dt0.getHour() == dt1.getHour()) && (dt0.getMinute() == dt1.getMinute()) && (dt1.getSecond() == dt0.getSecond()));
134      }
135
136      bool operator< (const CDate& dt0, const CDate& dt1)
137      {
138         // TODO :: Vérifier que les deux dates (dt0 et dt1) ont une référence vers le même calendrier.
139         if   (dt0.getYear()  < dt1.getYear())
140         { 
141            return true;
142         }
143         else if (dt0.getYear() == dt1.getYear())
144         { 
145            if   (dt0.getMonth()  < dt1.getMonth())
146            {
147               return true;
148            }
149            else if (dt0.getMonth() == dt1.getMonth())
150            {
151               if   (dt0.getDay()  < dt1.getDay())
152               {
153                   return true;
154               }
155               else if (dt0.getDay() == dt1.getDay())
156               {
157                  if    (dt0.getHour()  < dt1.getHour())
158                  {
159                     return true;
160                  }
161                  else if (dt0.getHour() == dt1.getHour())
162                  { 
163                     if   (dt0.getMinute()  < dt1.getMinute())
164                     {
165                        return true;
166                     }
167                     else if (dt0.getMinute() == dt1.getMinute())
168                     {
169                        if (dt0.getSecond() < dt1.getSecond())
170                           return true;
171                     }
172                  }
173               }
174            }
175         }
176         return false;
177      }
178
179      //-----------------------------------------------------------------
180
181      bool operator!=(const CDate & dt0, const CDate & dt1){ return !(dt1 == dt0); }
182      bool operator> (const CDate & dt0, const CDate & dt1){ return (dt1 < dt0); }
183      bool operator>=(const CDate & dt0, const CDate & dt1){ return ((dt0 > dt1) || (dt1 == dt0)); }
184      bool operator<=(const CDate & dt0, const CDate & dt1){ return ((dt0 < dt1) || (dt1 == dt0)); }
185
186      ///----------------------------------------------------------------
187
188} // namespace xios
189
190
191
192
193
194
Note: See TracBrowser for help on using the repository browser.