Re: [ferret_users] Regridding the point locations data into spatial manner

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

Hi,
Here I think a good way to treat the data in three datasets as one variable will be to call a function which concatenates them.

use AD01_temp.nc

use AD02_temp.nc

use AD03_temp.nc

def axis/x=69:69:1/units=degrees xlong

def axis/y=12:18:3/units=degrees ylat

! This puts the data onto an abstract y axis (with y=1,2,3,...)
let temp_3y = YCAT( AD03_temp[d=3], AD03_temp[d=2], AD03_temp[d=1])

! define a variable on your grid and use RESHAPE to put your data onto that grid.

let var_on_grid = x[gx=xlong] + y[gy=ylat]
let temp_grd_xy = RESHAPE(temp_3y, var_on_grid)

show grid temp_grd_xy

The way I wrote that, this will put AD03_temp[d=3] at y=12, AD03_temp[d=2] at y=15, and so on.  If this is not the correct order, change the order of concatenation.  Now, if you use the variable temp_grd_xy, Ferret will read the data from each file as needed.

Ansley

On 2/14/2016 10:00 PM, Nanda Kishore Reddy wrote:

Hi,

         I have a 3 locations of data in which it contains different latitudes (i.e. uniformly spaced) but longitude, depth and time axis are fixed for all the locations. How can I regrid these point locations data into spatial manner.

can data/all; can var/all; set me/size=300

use AD01_temp.nc

use AD02_temp.nc

use AD03_temp.nc

def axis/x=69:69:1/units=degrees xlong

def axis/y=12:18:3/units=degrees ylat

def grid/x=xlong/y=ylat grd_xy

let sp_grd0 = AD03_temp[d=3,g=grd_xy@asn]

let sp_grd1= AD02_temp[d=2,g=grd_xy@asn]

let sp_grd2 = AD01_temp[d=1,g=grd_xy@asn]

Thanks & Regards

Nanda Kishore reddy. B