[ferret_users] Re: How to save NetCDF data on curvilinear grid?

To: ferret_users@xxxxxxxx

Subject: [ferret_users] Re: How to save NetCDF data on curvilinear grid?

From: Szymon Roziewski <szymon.roziewski@xxxxxxxxx>

Date: Thu, 03 Mar 2011 10:58:32 +0100

Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:in-reply-to:references:date :message-id:subject:from:to:content-type; bh=c4plzUMNQRaP46Tk60vcghLV75lxnB+iBTje575n+to=; b=mv94TWbv4m3KhQaQyWDZgfK3N6rHdVq3GHS65sPR/Wq73U8UlMOZyUjMg/lb4CCuYM zUAlEWGWE8xoK+SLMFtoddNm3uNn6B8e8CeYxSCyJKRYm0E8Fi+ssgcPl+wbm6bVDp1h 88rLOL6150OcEfCfoRD0+f4cKwFQY6uQHdls4=

Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:date:message-id:subject:from:to :content-type; b=DCifAZN24xiVxY+V09bw/DxvI5itlUEyjkcM+PL5rMPUURpGm4HsVvLmHffrZKXGHP lnyt+jDFyK8ghm/eweFviTGZvmiYRZ0fkmADdykbqduW9HAmYmkEL2zGB2qWXHQOvQ7H S7YcVOSL3n90J/aPS/CL1273OiqgWraN9BrHo=

In-reply-to: <AANLkTi=T1A4669RfNkbb8mNndmy=stcyMetm2VYsc2EV@xxxxxxxxxxxxxx>

References: <AANLkTi=T1A4669RfNkbb8mNndmy=stcyMetm2VYsc2EV@xxxxxxxxxxxxxx>

Sender: owner-ferret_users@xxxxxxxx

2011/3/2 Szymon Roziewski <szymon.roziewski@xxxxxxxxx>

Hi Ferreters!
I'm trying to save some data in 660 colums and 750 rows matrix which is defined on curvilinear grid but I don't know how to do.

This is my ferret script:

DEFINE AXIS/x=1:660:1 xaxis
DEFINE AXIS/y=1:750:1 yaxis
DEFINE GRID/x=xaxis/y=yaxis gridlonlat
FILE/VARIABLES=rose/COLUMNS=660/GRID=gridlonlat warnTopography660x750_1space.dat
FILE/VARIABLES=clon_in/COLUMNS=660/GRID=gridlonlat LonIcmWarn660x750.dat
FILE/VARIABLES=clat_in/COLUMNS=660/GRID=gridlonlat LatIcmWarn660x750.dat
set var/title="LONGITUDE" clon_in[d=2]
set var/title="LATITUDE" clat_in[d=3]

save/clobber/file=Warn_in_ICM_Topography.nc rose[d=1],clon_in[d=2],clat_in[d=3]

First of all I read the data on some abstract grid (660x750) i.e. rose (main data), clon_in (matrix of curvilinear longitudes), clat_in (matrix of curvilinear latitudes).
Afterwards, I would like to save the rose data on curvilinear grid e.g. rose(clon_in, clat_in) as NetCDF file.
How could I manage to do that?

Kind regards,
Szymon Roziewski