预报时效计算

柿子柿子柿子

被封在家，无聊摆弄GrADS，出图的时候想不起来怎么查询当前时间并将其设置为标题，于是搞出了这个计算预报时间的东西。

1. * 是否闰年。是的话返回1，不是的话返回0；然后在计算2月的日数的时候直接加就好
   
2. function isLeapYear(y)
   
3.     if ( math_mod(y, 4) != 0)
   
4.         return 0
   
5.     else
   
6.         if (math_mod(y, 100) = 0)
   
7.             return 0
   
8.         else
   
9.             if (math_mod(y, 400) = 0)
   
10.                 return 1
    
11.             endif
    
12.         endif
    
13.         return 1
    
14.     endif
    
15. 
    
16. * 计算本月多少天。y为年份，m为月份
    
17. function getMonthDays(y, m)
    
18.     i = 1
    
19.     while (i<=12)
    
20.         mon.i = 31
    
21.         i = i + 1
    
22.     endwhile
    
23.     mon.2 = 28 + isLeapYear(y)
    
24.     mon.4 = 30
    
25.     mon.6 = 30
    
26.     mon.9 = 30
    
27.     mon.11 = 30
    
28.     m = m + 0
    
29.     return mon.m
    
30. 
    
31. * 计算时间。sdate为起报时间YYYYMMDDHH；period为预报时效，单位为小时
    
32. function calcDate(sdate, period)
    
33.     y = substr(sdate, 1, 4)
    
34.     m = substr(sdate, 5, 2)
    
35.     d = substr(sdate, 7, 2)
    
36.     h = substr(sdate, 9, 2)
    
37.     y = y + 0
    
38.     m = m + 0
    
39.     d = d + 0
    
40.     h = h + period
    
41.     while (h>=24)
    
42.         h = h - 24
    
43.         d = d + 1
    
44.         if (d > getMonthDays(y,m))
    
45.             d = 1
    
46.             m = m + 1
    
47.             if (m>12)
    
48.                 m = 1
    
49.                 y = y + 1
    
50.             endif
    
51.         endif
    
52.     endwhile
    
53.    
    
54.     ret = y
    
55.     if (m<10)
    
56.         ret = ret % '0' % m
    
57.     else
    
58.         ret = ret %  m
    
59.     endif
    
60.     if (d<10)
    
61.         ret = ret % '0' % d
    
62.     else
    
63.         ret = ret %  d
    
64.     endif
    
65.     if (h<10)
    
66.         ret = ret % '0' % h
    
67.     else
    
68.         ret = ret %  h
    
69.     endif
    
70.     return ret

复制代码

测试内容如下

1. say calcDate('2020072412', 2)
   
2. say calcDate('2020073112', 24)
   
3. say calcDate('2020022812', 24)
   
4. say calcDate('2020022812', 48)
   
5. say calcDate('2020022812', 72)
   
6. say calcDate('2021022812', 24)
   
7. say calcDate('2021022812', 48)
   
8. say calcDate('2021022812', 72)

复制代码

测试结果为

1. 2020072414
   
2. 2020080112
   
3. 2020022912
   
4. 2020030112
   
5. 2020030212
   
6. 2021030112
   
7. 2021030212
   
8. 2021030312

复制代码

（目前看着暂时没发现问题……如果用着的时候发现问题再说）