| 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: BEGIN & strd = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12] & ntxt = 'Annual' & END |
|---|
| 21 | ENDCASE |
|---|
| 22 | ENDCASE |
|---|
| 23 | |
|---|
| 24 | nmth = (size(strd))[1] |
|---|
| 25 | njpt = n_elements((lindgen(jpt))[0:*:12]) |
|---|
| 26 | IF debug_w THEN print, ' nmth,njpt', nmth, njpt |
|---|
| 27 | |
|---|
| 28 | idxm = lonarr(nmth, njpt) |
|---|
| 29 | |
|---|
| 30 | FOR imth = 0, nmth-1 DO BEGIN |
|---|
| 31 | idxm (imth, *)= (lindgen(jpt))[strd(imth)-1:*:12] |
|---|
| 32 | ENDFOR |
|---|
| 33 | |
|---|
| 34 | IF nmth EQ 1 THEN ntxt = mth[strd-1] |
|---|
| 35 | |
|---|
| 36 | print, ' Warning: correlation computed and average for following period: ', ntxt |
|---|
