1 | mth = ['Jan', 'Feb', 'Mar', 'Apr', 'May', 'Jun', 'Jul', 'Aug', 'Sep', 'Oct', 'Nov', 'Dec', 'DJF', 'MAM', 'JJA', 'SON'] |
---|
2 | CASE STRMID(stddev_mth, 0, 2) OF |
---|
3 | '01': strd = [1] |
---|
4 | '02': strd = [2] |
---|
5 | '03': strd = [3] |
---|
6 | '04': strd = [4] |
---|
7 | '05': strd = [5] |
---|
8 | '06': strd = [6] |
---|
9 | '07': strd = [7] |
---|
10 | '08': strd = [8] |
---|
11 | '09': strd = [9] |
---|
12 | '10': strd = [10] |
---|
13 | '11': strd = [11] |
---|
14 | '12': strd = [12] |
---|
15 | ELSE : CASE STRMID(stddev_mth, 0, 3) OF |
---|
16 | 'DJF': BEGIN & strd = [1, 2, 12] & ntxt = 'DJF' & END |
---|
17 | 'MAM': BEGIN & strd = [3, 4, 5] & ntxt = 'MAM' & END |
---|
18 | 'JJA': BEGIN & strd = [6, 7, 8] & ntxt = 'JJA' & END |
---|
19 | 'SON': BEGIN & strd = [9, 10, 11] & ntxt = 'SON' & END |
---|
20 | ELSE: CASE STRMID(stddev_mth, 0, 6) OF |
---|
21 | 'JASOND' : BEGIN & strd = [7, 8, 9, 10, 11, 12] & ntxt = 'JASOND' & END |
---|
22 | 'JFMAMJ' : BEGIN & strd = [1, 2, 3, 4, 5, 6] & ntxt = 'JFMAMJ' & END |
---|
23 | ELSE : BEGIN & strd = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12] & ntxt = 'Annual' & END |
---|
24 | ENDCASE |
---|
25 | ENDCASE |
---|
26 | ENDCASE |
---|
27 | |
---|
28 | nmth = (size(strd))[1] |
---|
29 | njpt = n_elements((lindgen(jpt))[0:*:12]) |
---|
30 | IF debug_w THEN print, ' nmth,njpt', nmth, njpt |
---|
31 | |
---|
32 | idxm = lonarr(nmth, njpt) |
---|
33 | |
---|
34 | FOR imth = 0, nmth-1 DO BEGIN |
---|
35 | idxm (imth, *)= (lindgen(jpt))[strd(imth)-1:*:12] |
---|
36 | ENDFOR |
---|
37 | |
---|
38 | IF nmth EQ 1 THEN ntxt = mth[strd-1] |
---|
39 | |
---|
40 | print, ' Warning: monthly computation and average for following period: ', ntxt |
---|