Re: [ferret_users] Regridding from a regular grid to natural coordinates defined by along/cross isotherms, isobars, bathymetry, streamlines etc.

["'Drijfhout, S.S. (Sybren)' via _OAR PMEL Ferret Users" <ferret_users@xxxxxxxx>]

Thank you Ryo,

 

I already figured a script that is doing something like you suggest, but it is quite slow and not elegant, including concatenating a few hundred single curves to one file and having to repeat this more than 40 times with an insecure remote connection. Maybe you’re right and should I combine it with a fortran script.

 

Sybren

 

Van: Ryo Furue <furue@xxxxxxxxxx>
Datum: dinsdag, 4 mei 2021 om 18:19
Aan: Drijfhout, S.S. (Sybren) <s.s.drijfhout@xxxxx>
CC: ferret_users@xxxxxxxx <ferret_users@xxxxxxxx>
Onderwerp: Re: [ferret_users] Regridding from a regular grid to natural coordinates defined by along/cross isotherms, isobars, bathymetry, streamlines etc.

Hi Sybren,

 

In theory I have to map xv(i,j) to xvnew(I,jstar)  with jstar=j-j1

Ferret wants me to define a new y-axis but does not allow me to define ystar by a let command, and using j and ja gives errors because of recursive problems.  

Has anybody experience with regridding to something like streamline-coordinates and/or does anybody know how to perform such a regridding operation in ferret?

 

I do that using Fortran. I didn't even consider using Ferret. . . . But I guess you can at least define xstar(i) and ystar(i) for a single curve.  Then you can get the values along the curve using the SAMPLEXY(var, xstar, ystar) function.  Save the result and move on to the next curve. . . . Finally, you combine these "layers" in the ystar direction.

 

In my Fortran code, I just shift the array in the j direction by j1.

 

Regards,

Ryo

 

 

 

Sybren