Re: [ferret_users] define axis / grid from variables

["Paulo B. Oliveira" <pboliveira@xxxxxxx>]

Hi Marina,

Try the attached script and see if you get the same output.

You may wish to look at the curv_to_rect function
( http://ferret.pmel.noaa.gov/Ferret/documentation/users-guide/appendix-a-external-functions/CURV-TO-RECT?searchterm=curv )
if are going to combine your data with data from other sources.

Saludos,

Paulo


On Mon, 2015-04-20 at 08:54 +0000, Marina Chifflet wrote:

Hi Folks

 

I apologize for the easy question, but I can’t find the right syntax and I’m becoming crazy….

 

I have files with model outputs on a i/j/k grid (not geo-ref)

I=1:236 / j=1:150 / k=1:72

(example of the diatoms file “1997_4.3v4_DIAT_march.nc” on ftp://ftp.azti.es/resto/mchifflet/)

On other hand, I have the latitudes and longitudes in a separated file, with the same i/j grid: ll5km.nc (attached file)

LAT and LON are defined on the 2D grid i,j

 

How can I create axis/grid from the variables in ll5km.nc ?

To be able to save the output variables of the model on the right grid.

 

I tried a lot of syntaxes but can’t find the right one.

 

Because each of the variables LAT and LON are defined on a 2D i,j grid, I cannot use :

define axis/x/from_data/name=longulf lon[d=1]

 

Do you know the syntax to do this?

Thanks a lot.

Regards,

______________________________________

Dr Marina Chifflet

AZTI - Tecnalia / Marine research Division
Herrera Kaia Portualde z/g
20110 Pasaia (Gipuzkoa) - Spain
Tel: (+34) 667 174 495
e-mail: mchifflet@xxxxxxx
http://www.researchgate.net/profile/Marina_Chifflet
http://www.azti.es

 

Attachment: diat_map.png
Description: PNG image

 ! NOAA/PMEL TMAP
 ! FERRET v6.93 (PyFerret 1.1.0)
 ! Linux 2.6.35.14-106.fc14.i686 - 12/05/14
 ! 20-Apr-15 11:15     

use ll5km.nc
use 1997_4.3v4_DIAT_march.nc

shade/k=1/l=1 diat,lon[d=1],lat[d=1]
go land_detail

frame/file=diat_map.png

quit