Re: [ferret_users] vertical weighted averagng

[Abhishek Savita <abhisheknagar2@xxxxxxxxx>]

Hi... all

         I have also faced same problem earlier. I have checked  it analytically as well as with the matlab also. Actully ferret gives the weighted average if we give like var[x=@ave,y=@ave,z=@ave,.....]. But If anybody interested to compute normal average then in place of giving X, Y etc we should give corresponding index number as var[i=@ave,j=@ave,.....]. Where i and j are the index number corresponding Lon and lat.

Please check it I hope it will work.

with regards
Abhishek

On Fri, Mar 24, 2017 at 9:55 AM, Ryo Furue <furue@xxxxxxxxxx> wrote:

Hi Robert and Martin,

yes? list/k=1 temp[x=@ave,y=@ave]

...

      -1.788

but I'd like to have area weighed averaged value. How should I told ferret that x and y is long and lat for averaging computations,

With temp[x=@ave,y=@ave], Ferret does compute the area weighted average if the horizontal axes are longitude and latitude.  If the axes aren't longitude and latitude, Ferret computes x=@ave and y=@ave separately (I think), which may not coincide with the correct area-weighted average.

I don't know if it's possible to compute the area of each gridbox within Ferret. You may have to supply that information from outside. Then,

  let area_temp = temp * boxarea

  let area_average_temp = area_temp[x=@sum,y=@sum] / boxarea[x=@sum,y=@sum]

Ryo

--

Regards
Abhishek Savita

Research Scholar (Earth System Science Technology)

Center For Oceans, Rivers, Atmosphere & Land Science Technology
Indian Institute of Technology, Kharagpur
+91-8609704619