[ferret_users] Re: normalization query

On Fri, Dec 2, 2016 at 8:15 AM, binny ra <binny2020srb@xxxxxxxxx> wrote:

Dear Ferret Users

I am quite new to ferret and I am performing normalisation of ocean heat content after calculating the running sum and zonal integration.
my script is looking like this

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! anomaly files !!!!!!!!!!!!!!!!!!!!!!!!!!!!
use all_anm.nc      !global (x=0E:0W,y=62.5S:62.5N)
use indi_anm.nc    !indian ocean (20E:120E, 62.5S:30N)

sh da
     currently SET data sets:
    1> ./all_anm.nc
name     title                               I         J         K         L         M         N
OHCG5_ANM_ALL
          OHCG5-OHCG5_CLI[GT=OHCG5@ASN]    1:360     1:126     ...       1:142     ...       ...

    2> ./indi_anm.nc
name     title                            I         J         K         L         M         N
OHCG5_ANM_IN
          OHCG5-OHCG5_CLI[GT=OHCG5@ASN]    1:100     1:93      ...       1:142     ...       ...

sh gr ohcg5_anm_all[d=1]
    GRID GAN1
name       axis                       # pts   start                  end                 subset
LON       LONGITUDE          360mr   0.5E                 0.5W                full
LAT28_153 LATITUDE           126 r   62.5S                62.5N               full
normal    Z
TIME      TIME               142 i   15-JAN-2005 00:00    15-OCT-2016 00:00   full
normal    E
normal    F

GRID GFA1
name       axis              # pts   start                 end                 subset
LON21_120 LONGITUDE          100mr   20.5E                119.5E              full
LAT28_1531_93 LATITUDE        93 r    62.5S                29.5N               full
normal    Z
TIME      TIME               142 i   15-JAN-2005 00:00    15-OCT-2016 00:00   full
normal    E
normal    F
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! zonal integration and running sum !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

let ohcg_g5s=ohcg5_anm_all[d=1,x=@din,y=@rsum]     ((!! value will come in 10^15 ))
let ohcg_in5s=ohcg5_anm_in[d=2,x=@din,y=@rsum]

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! yearly average !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
define axis/edges/t/units=months/t0=15-dec-2004/t=15-jan-2005:15-nov-2016:1/calendar=360_days tmonthly
let ohcc_gd=ohcg_g5s[gt=tmonthly@asn]

define axis/edges/t/units=years/t=15-jan-2005:15-oct-2016:1/calendar=360_days t_yearly
let ohcc_g5=ohcc_gd[gt=t_yearly@ave]

define axis/edges/t/units=months/t0=15-dec-2004/t=15-jan-2005:15-nov-2016:1/calendar=360_days tmonthly
let ohcc_ind=ohcg_in5s[gt=tmonthly@asn]

define axis/edges/t/units=years/t=15-jan-2005:15-oct-2016:1/calendar=360_days t_yearly
let ohcc_in5=ohcc_ind[gt=t_yearly@ave]

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! difference of two years !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

let ohcc_gg5s=ohcc_g5[l=12]-ohcc_g5[l=1]
let ohcc_ii5s=ohcc_in5[l=12]-ohcc_in5[l=1]

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! normalisation by global value !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

let ohcc_noi=ohcc_ii5s/ohcc_gg5s

list ohcc_noi is showing 1 at every y value

if i try like this way such as
let ohcc_noi=ohcc_ii5s[y=62.5S:30N,x=20E:120E]/ohcc_gg5s[y=62.5S:30N,x=20E:120E]

then also the same thing happens and i got a straight line plot.

this is the list output
list ohcc_noi
29.5N / 93:   ....
28.5N / 92:   ....
27.5N / 91:   ....
26.5N / 90:   ....
25.5N / 89:   ....
24.5N / 88: 1.000
23.5N / 87: 1.000
22.5N / 86: 1.000
21.5N / 85: 1.000
20.5N / 84: 1.000
19.5N / 83: 1.000
18.5N / 82: 1.000
17.5N / 81: 1.000
16.5N / 80: 1.000
15.5N / 79: 1.000
14.5N / 78: 1.000
13.5N / 77: 1.000
12.5N / 76: 1.000
11.5N / 75: 1.000
10.5N / 74: 1.000
9.5N / 73: 1.000
8.5N / 72: 1.000
7.5N / 71: 1.000
6.5N / 70: 1.000
5.5N / 69: 1.000
4.5N / 68: 1.000
3.5N / 67: 1.000
2.5N / 66: 1.000
1.5N / 65: 1.000
0.5N / 64: 1.000
0.5S / 63: 1.000
1.5S / 62: 1.000
2.5S / 61: 1.000
3.5S / 60: 1.000
4.5S / 59: 1.000
5.5S / 58: 1.000
6.5S / 57: 1.000
7.5S / 56: 1.000
8.5S / 55: 1.000
9.5S / 54: 1.000
10.5S / 53: 1.000
11.5S / 52: 1.000
12.5S / 51: 1.000
13.5S / 50: 1.000
14.5S / 49: 1.000
15.5S / 48: 1.000
16.5S / 47: 1.000
17.5S / 46: 1.000
18.5S / 45: 1.000
19.5S / 44: 1.000
20.5S / 43: 1.000
21.5S / 42: 1.000
22.5S / 41: 1.000
23.5S / 40: 1.000
24.5S / 39: 1.000
25.5S / 38: 1.000
26.5S / 37: 1.000
27.5S / 36: 1.000
28.5S / 35: 1.000
29.5S / 34: 1.000
30.5S / 33: 1.000
31.5S / 32: 1.000
32.5S / 31: 1.000
33.5S / 30: 1.000
34.5S / 29: 1.000
35.5S / 28: 1.000
36.5S / 27: 1.000
37.5S / 26: 1.000
38.5S / 25: 1.000
39.5S / 24: 1.000
40.5S / 23: 1.000
41.5S / 22: 1.000
42.5S / 21: 1.000
43.5S / 20: 1.000
44.5S / 19: 1.000
45.5S / 18: 1.000
46.5S / 17: 1.000
47.5S / 16: 1.000
48.5S / 15: 1.000
49.5S / 14: 1.000
50.5S / 13: 1.000
51.5S / 12: 1.000
52.5S / 11: 1.000
53.5S / 10: 1.000
54.5S / 9: 1.000
55.5S / 8: 1.000
56.5S / 7: 1.000
57.5S / 6: 1.000
58.5S / 5: 1.000
59.5S / 4: 1.000
60.5S / 3: 1.000
61.5S / 2: 1.000
62.5S / 1: 1.000

kindly help me out i tried in every way but not getting any success. i hope ferreters are having any clue to this problem

regards, Binny