[ferret_users] selecting data

[Nitin Patil <nitinpatil85@xxxxxxxxx>]

Dear Ferret users,

I am using aod data which is ranged from 0 to >1.5 in the netcdf file (time range is daily for season l=122). I am trying to average the JJAS for specific region but I have to exclude data which is >0.8 in the file.

For that I have tried the following:

-------------------------------------------------------------------------------------

let aod_cne = AOD[d=1,l=1:122]

let data_cne_aod = if aod_cne gt 0.8 then 0 else 1

let xcne = {77,78,78,81,81,82,85,85,88,88,86,86,87,87,86,86,82,82,83,83,84,84,83,83,79,79,78,77,77}

let ycne = {31,31,30,30,29,29,28,28,27,24,24,24,23,20,20,19,19,20,20,21,21,24,24,25,25,27,27,27,27}

let cne = pt_in_poly(data_cne_aod,xcne, ycne)

let cne_data_2000 = if cne eq 1 then cne*AOD[d=1,l=1:122@ave]

------------------------------------------------------------------------------------- 

Whether it is correct or Making some mistake?

--

Regards,

Nitin Patil