Hi,
I don't know what you are intending with the combination of the units of a @DIN definite integral and the @RSUM running sum, but there is a big difference in how those transformations are computed. (See the Information about Transformations: http://ferret.pmel.noaa.gov/
Ferret/documentation/users- ). When the grid is in known units on the surface of the earth, the results of the integral are returned in data * square meters. The running sum is returned as an unweighted sum of the data in each grid cell, so those are just in the units of the heat content variable. The @DIN, in units*square meter will overwhelm the calculation.guide/variables-xpressions/ XPRESSIONS#_GeneralInfo -Ansley
On 12/1/2016 1:15 PM, binny ra wrote:
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!my script is looking like thisDear Ferret UsersI am quite new to ferret and I am performing normalisation of ocean heat content after calculating the running sum and zonal integration.
!!!!!!!! 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-no v-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/calenda r=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-no v-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/calenda r=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