Re: [ferret_users] Re: bin data files in ferret

["ESTANISLAO GAVILAN PASCUAL-AHU" <e.gavilan@xxxxxxxxxx>]

Hi community,

Sorry, I did not reply before. I tried to use the E- aggregation but it did not work because it requires the type of the file and I am not sure how to set that in the script. Currently, I was trying load the file with a loop following this script but I having an error which I am not sure how to solve. This is the script 

define axis/X=1:3060:1: xc
define axis/Y=1:510:1 yc
define grid/X=xc/Y=yc ecco
set region/x=1021:1530
use time_iterr.nc ! this is a file with the different nsteps 
let Lend=`time_mo,return=lsize`
repeat/L=1:`Lend`:1 \
(file/var=ssh/grid=ecco/format=stream/type=r4/swap ETAN.000000`time_mo`.data)

If I do this it says "variable unknown or not in dataset: time_mo". Interestingly, if I remove the "." of data, it works but obviously ferret cannot find the file.

Kind regards,

Estanis

>-----Original Messages-----
> From: "Ansley C. Manke" <ansley.b.manke@xxxxxxxx>
> Sent Time: Wednesday, November 6, 2019
> To: ferret_users@xxxxxxxx
> Cc: 
> Subject: [ferret_users] Re: bin data files in ferret
> 
> Hi,
> 
> I am working only part time so the Ferret List is the best place to 
> write.  I just wrote an answer to the list suggesting swapping X and Y, 
> so it's good to hear you found that solution.  I'm writing back to the 
> List here, because use of aggregations is a good discussion for the users
> 
> It's true that the T aggregation works only with NetCDF files, or 
> OPeNDAP datasets.  I suppose you could try opening the files, defining 
> an ensemble aggregation, which would give you data on a X-Y-E grids, and 
> then reshape that E-aggregation to a grid in X-Y-T with the appropriate 
> T axis.  Or, if you don't mind making copies of the files, when reading 
> the individual .bin files, define each grid with a single-point Time 
> axis, and write to a netCDF file that you could aggregate in time later on.
> 
> 
> On 11/5/2019 12:33 AM, ESTANISLAO GAVILAN PASCUAL-AHU wrote:
> > Hi Ansley,
> >
> > sorry, I am contacting you to your job mail instead of the user ferret 
> > list. I sent a email a few days ago but I did not have any reply. I am 
> > trying to analyse data from MITgcm model with Ferret. I was able to 
> > open the bin files using these commands.
> >
> > define axis/X=1:3060:1: xc
> > define axis/Y=1:510:1 yc
> > define grid/X=xc/Y=yc ecco
> > set region/x=1021:1530
> > file/var=bathy/grid=ecco/format=stream/type=r4/swap 
> > GEBCO_510x6x510_ver06_dig.bin
> > set variable/bad=0 bathy
> > set window/aspect=1
> > SHADE/palette=rainbow  bathy
> >
> > These commands works perfectly fine to open the I/O data from MITgcm. 
> > However, I would like to open a string of T dependent files. I am not 
> > sure if you are familiar with this model but the model creates an 
> > output data in each output frequency. Therefore, at the end of the 
> > simulation we have X amount of files from the same variable 
> > (SALT.00040.data,SALT.00080.data,SALT.00120.data). I tried to use the 
> > aggregate.jnl function, but it seems that function only works with 
> > netcdf files.. Thanks in advance for all your help
> >
> >
> > Kind regards,
> >
> > Estanislao
> >