Re: [ferret_users] open bin mitgcm output

To: ferret_users@xxxxxxxx
Subject: Re: [ferret_users] open bin mitgcm output
From: "Ansley C. Manke" <ansley.b.manke@xxxxxxxx>
Date: Tue, 5 Nov 2019 11:10:57 -0800
In-reply-to: <13604938.1035e.16e24cd297d.Coremail.e.gavilan@hhu.edu.cn>
List-archive: <https://groups.google.com/a/noaa.gov/group/ferret_users/>
List-help: <https://support.google.com/a/noaa.gov/bin/topic.py?topic=25838>, <mailto:ferret_users+help@noaa.gov>
List-id: <ferret_users.noaa.gov>
List-post: <https://groups.google.com/a/noaa.gov/group/ferret_users/post>, <mailto:ferret_users@noaa.gov>
Mailing-list: list ferret_users@xxxxxxxx; contact ferret_users+owners@xxxxxxxx
References: <13604938.1035e.16e24cd297d.Coremail.e.gavilan@hhu.edu.cn>
Sender: owner-ferret_users@xxxxxxxx
User-agent: Mozilla/5.0 (Windows NT 10.0; WOW64; rv:68.0) Gecko/20100101 Thunderbird/68.2.0

Hi Estanislao,

I'm not familiar with this file type, but here are some ideas for you in trying to read it with Ferret.

* check that the data type of r4 translates to values that make sense, by just reading a few values

DEFINE AXIS/X=1:10:1 xaxis

DEFINE GRID/X=xaxis xgrid

FILE/VAR=depth/GRID=xgrid/FORMAT=UNF/type=r4 Depth.data

LIST depth

Are the values correct? If not try FORMAT=stream or other data types

* Ferret grids are all in the order XYZTEF, or some subset of those, so that X varies fastest. If the file has Y varying fastest, then you'd need to define the grid with

DEFINE AXIS/X=1:3060:1 xaxis

DEFINE AXIS/Y=1:510:1 yaxis

DEFINE GRID/X=xaxis/Y=yaxis xyecco

...

and once you've read the data define a new variable using a TRANSPOSE function to reorder the data to a correct XY grid.

On 10/31/2019 7:30 PM, ESTANISLAO GAVILAN PASCUAL-AHU wrote:

Hi Community,

I would like to open a data output file from ECCO2 in MITgcm. I have tried the following commands but the conversion is wrong

DEFINE AXIS/X=1:510:1 x510

DEFINE AXIS/Y=1:3060:1 y3060

DEFINE GRID/X=x510/Y=y3060 xyecco

FILE/VAR=depth/GRID=xyecco/FORMAT=UNF/type=r4 Depth.data

I have also tried

FILE/VAR=depth/GRID=xyecco/FORMAT=stream/type=r4 Depth.data

I know I can use the matlab tool called readbin but I would like to learn to do it with FERRET.

Kind regards,

Estanislao

References:
- [ferret_users] open bin mitgcm output
  - From: ESTANISLAO GAVILAN PASCUAL-AHU

Previous by thread: [ferret_users] open bin mitgcm output
Next by thread: [ferret_users] installing pyferret on public Linux machine

[Thread Prev][Thread Next][Index]