Re: [ferret_users] normalization query

["Ansley C. Manke" <ansley.b.manke@xxxxxxxx>]

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-guide/variables-xpressions/XPRESSIONS#_GeneralInfo). 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.

-Ansley

On 12/1/2016 1:15 PM, binny ra 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