[MITgcm-support] ECCO regrid interpolation
ESTANISLAO GAVILAN PASCUAL-AHU
e.gavilan at hhu.edu.cn
Fri Jan 17 00:51:17 EST 2020
Hi Ryan,
thanks so much for you reply. I will take a look. I did not know about this tool.
Kind regards,
Estanislao
> -----Original Messages-----
> From: mitgcm-support-request at mitgcm.org
> Sent Time: Thursday, January 16, 2020
> To: mitgcm-support at mitgcm.org
> Cc:
> Subject: MITgcm-support Digest, Vol 199, Issue 13
>
> Send MITgcm-support mailing list submissions to
> mitgcm-support at mitgcm.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> or, via email, send a message with subject or body 'help' to
> mitgcm-support-request at mitgcm.org
>
> You can reach the person managing the list at
> mitgcm-support-owner at mitgcm.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of MITgcm-support digest..."
>
>
> Today's Topics:
>
> 1. Re: ECCO regrid interpolation (ESTANISLAO GAVILAN PASCUAL-AHU)
> 2. Re: ECCO regrid interpolation (Martin Losch)
> 3. Re: ECCO regrid interpolation (Ryan Abernathey)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Wed, 15 Jan 2020 15:28:28 +0800 (GMT+08:00)
> From: "ESTANISLAO GAVILAN PASCUAL-AHU" <e.gavilan at hhu.edu.cn>
> To: mitgcm-support at mitgcm.org
> Subject: Re: [MITgcm-support] ECCO regrid interpolation
> Message-ID: <5c0d2ce2.11c6.16fa81a856a.Coremail.e.gavilan at hhu.edu.cn>
> Content-Type: text/plain; charset=us-ascii
>
>
> Hi Martin,
>
> That what I thought at the beginning. However, the grid since of ECCO2 changes, so I was reluctant to follow that approach.
>
> Kind regards,
>
> Estanislao
>
> > -----Original Messages-----
> > From: mitgcm-support-request at mitgcm.org
> > Sent Time: Wednesday, January 15, 2020
> > To: mitgcm-support at mitgcm.org
> > Cc:
> > Subject: MITgcm-support Digest, Vol 199, Issue 12
> >
> > Send MITgcm-support mailing list submissions to
> > mitgcm-support at mitgcm.org
> >
> > To subscribe or unsubscribe via the World Wide Web, visit
> > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> > or, via email, send a message with subject or body 'help' to
> > mitgcm-support-request at mitgcm.org
> >
> > You can reach the person managing the list at
> > mitgcm-support-owner at mitgcm.org
> >
> > When replying, please edit your Subject line so it is more specific
> > than "Re: Contents of MITgcm-support digest..."
> >
> >
> > Today's Topics:
> >
> > 1. ECCO regrid interpolation (ESTANISLAO GAVILAN PASCUAL-AHU)
> > 2. Re: ECCO regrid interpolation (Martin Losch)
> >
> >
> > ----------------------------------------------------------------------
> >
> > Message: 1
> > Date: Tue, 14 Jan 2020 09:53:11 +0800 (GMT+08:00)
> > From: "ESTANISLAO GAVILAN PASCUAL-AHU" <e.gavilan at hhu.edu.cn>
> > To: mitgcm-support at mitgcm.org
> > Subject: [MITgcm-support] ECCO regrid interpolation
> > Message-ID: <73997e45.db5.16fa1c13616.Coremail.e.gavilan at hhu.edu.cn>
> > Content-Type: text/plain; charset="gbk"
> >
> > Dear community,
> >
> >
> > I was trying to regridded the sea surface elevation data to the u and v points in order to compute the normal transport across each cell. At the beginning I thought I could use griddata from MATLAB or python. However, since my grid is curvilinear I think I cannot do such thing. Does MITgcm have any kind of tool for regridding this type of data? Thanks in advance for your help
> >
> >
> > Kind regards,
> >
> >
> > Estanislao
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200114/3f0a8587/attachment-0001.html>
> >
> > ------------------------------
> >
> > Message: 2
> > Date: Tue, 14 Jan 2020 11:23:12 +0100
> > From: Martin Losch <Martin.Losch at awi.de>
> > To: MITgcm Support <mitgcm-support at mitgcm.org>
> > Subject: Re: [MITgcm-support] ECCO regrid interpolation
> > Message-ID: <BB97573C-E6FD-4A1A-8BEA-4890AAE003F0 at awi.de>
> > Content-Type: text/plain; charset="utf-8"
> >
> > Hi Estanislao,
> >
> > wouldn?t the natural C-grid way of computing EtaN at U points be
> > EtaNatU(i,j) = 0.5*(EtaN(i-1,j)+EtaN(i,j))
> > and similarly at V points
> > EtaNatV(i,j) = 0.5*(EtaN(i,j-1)+EtaN(i,j))
> > ?
> >
> > Martin
> >
> > > On 14. Jan 2020, at 02:53, ESTANISLAO GAVILAN PASCUAL-AHU <e.gavilan at hhu.edu.cn> wrote:
> > >
> > > Dear community,
> > >
> > > I was trying to regridded the sea surface elevation data to the u and v points in order to compute the normal transport across each cell. At the beginning I thought I could use griddata from MATLAB or python. However, since my grid is curvilinear I think I cannot do such thing. Does MITgcm have any kind of tool for regridding this type of data? Thanks in advance for your help
> > >
> > > Kind regards,
> > >
> > > Estanislao
> > > _______________________________________________
> > > MITgcm-support mailing list
> > > MITgcm-support at mitgcm.org
> > > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> >
> >
> >
> > ------------------------------
> >
> > Subject: Digest Footer
> >
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> >
> >
> > ------------------------------
> >
> > End of MITgcm-support Digest, Vol 199, Issue 12
> > ***********************************************
>
>
>
> ------------------------------
>
> Message: 2
> Date: Wed, 15 Jan 2020 09:40:37 +0100
> From: Martin Losch <Martin.Losch at awi.de>
> To: MITgcm Support <mitgcm-support at mitgcm.org>
> Subject: Re: [MITgcm-support] ECCO regrid interpolation
> Message-ID: <4BC69496-D432-48C7-B213-DA06820849D5 at awi.de>
> Content-Type: text/plain; charset="us-ascii"
>
> Hi again,
>
> yes, ECCO2 uses a curvilinear grid, but the cell-to-cell variations are very small.
>
> In general, the internal MITgcm averaging and differencing works in the same way (i.e. the gradient of eta at a U-point is actually something like (etaN(i,j)-etaN(i-1,j))*recip_dxC), so I think you are safe doing it like that.
>
> Martin
>
>
> > On 15. Jan 2020, at 08:28, ESTANISLAO GAVILAN PASCUAL-AHU <e.gavilan at hhu.edu.cn> wrote:
> >
> >
> > Hi Martin,
> >
> > That what I thought at the beginning. However, the grid since of ECCO2 changes, so I was reluctant to follow that approach.
> >
> > Kind regards,
> >
> > Estanislao
> >
> >>
> >> Dear community,
> >>
> >>
> >> I was trying to regridded the sea surface elevation data to the u and v points in order to compute the normal transport across each cell. At the beginning I thought I could use griddata from MATLAB or python. However, since my grid is curvilinear I think I cannot do such thing. Does MITgcm have any kind of tool for regridding this type of data? Thanks in advance for your help
> >>
> >>
> >> Kind regards,
> >>
> >>
> >> Estanislao
> >> -------------- next part --------------
> >> An HTML attachment was scrubbed...
> >> URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200114/3f0a8587/attachment-0001.html>
> >>
> >> ------------------------------
> >>
> >> Message: 2
> >> Date: Tue, 14 Jan 2020 11:23:12 +0100
> >> From: Martin Losch <Martin.Losch at awi.de>
> >> To: MITgcm Support <mitgcm-support at mitgcm.org>
> >> Subject: Re: [MITgcm-support] ECCO regrid interpolation
> >> Message-ID: <BB97573C-E6FD-4A1A-8BEA-4890AAE003F0 at awi.de>
> >> Content-Type: text/plain; charset="utf-8"
> >>
> >> Hi Estanislao,
> >>
> >> wouldn?t the natural C-grid way of computing EtaN at U points be
> >> EtaNatU(i,j) = 0.5*(EtaN(i-1,j)+EtaN(i,j))
> >> and similarly at V points
> >> EtaNatV(i,j) = 0.5*(EtaN(i,j-1)+EtaN(i,j))
> >> ?
> >>
> >> Martin
> >>
> >>> On 14. Jan 2020, at 02:53, ESTANISLAO GAVILAN PASCUAL-AHU <e.gavilan at hhu.edu.cn> wrote:
> >>>
> >>> Dear community,
> >>>
> >>> I was trying to regridded the sea surface elevation data to the u and v points in order to compute the normal transport across each cell. At the beginning I thought I could use griddata from MATLAB or python. However, since my grid is curvilinear I think I cannot do such thing. Does MITgcm have any kind of tool for regridding this type of data? Thanks in advance for your help
> >>>
> >>> Kind regards,
> >>>
> >>> Estanislao
> >>> _______________________________________________
> >>> MITgcm-support mailing list
> >>> MITgcm-support at mitgcm.org
> >>> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> >>
> >>
> >>
> >> ------------------------------
> >>
> >> Subject: Digest Footer
> >>
> >> _______________________________________________
> >> MITgcm-support mailing list
> >> MITgcm-support at mitgcm.org
> >> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> >>
> >>
> >> ------------------------------
> >>
> >> End of MITgcm-support Digest, Vol 199, Issue 12
> >> ***********************************************
> >
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
>
>
> ------------------------------
>
> Message: 3
> Date: Wed, 15 Jan 2020 10:16:14 -0500
> From: Ryan Abernathey <ryan.abernathey at gmail.com>
> To: mitgcm-support at mitgcm.org
> Subject: Re: [MITgcm-support] ECCO regrid interpolation
> Message-ID:
> <CAJAZx5AgcPfdaZkPOu4A29vT00fmL=+NQTPN-5C49Wwg46=ZJQ at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> This sort of grid interpolation / difference application is precisely what
> our software package xgcm is for:
> https://xgcm.readthedocs.io/en/latest/
>
> Here is an example of using it with ECCOv4:
> https://xgcm.readthedocs.io/en/latest/example_eccov4.html
>
> Perhaps this will be helpful to you.
>
> On Wed, Jan 15, 2020 at 7:18 AM Martin Losch <Martin.Losch at awi.de> wrote:
>
> > Hi again,
> >
> > yes, ECCO2 uses a curvilinear grid, but the cell-to-cell variations are
> > very small.
> >
> > In general, the internal MITgcm averaging and differencing works in the
> > same way (i.e. the gradient of eta at a U-point is actually something like
> > (etaN(i,j)-etaN(i-1,j))*recip_dxC), so I think you are safe doing it like
> > that.
> >
> > Martin
> >
> >
> > > On 15. Jan 2020, at 08:28, ESTANISLAO GAVILAN PASCUAL-AHU <
> > e.gavilan at hhu.edu.cn> wrote:
> > >
> > >
> > > Hi Martin,
> > >
> > > That what I thought at the beginning. However, the grid since of ECCO2
> > changes, so I was reluctant to follow that approach.
> > >
> > > Kind regards,
> > >
> > > Estanislao
> > >
> > >>
> > >> Dear community,
> > >>
> > >>
> > >> I was trying to regridded the sea surface elevation data to the u and v
> > points in order to compute the normal transport across each cell. At the
> > beginning I thought I could use griddata from MATLAB or python. However,
> > since my grid is curvilinear I think I cannot do such thing. Does MITgcm
> > have any kind of tool for regridding this type of data? Thanks in advance
> > for your help
> > >>
> > >>
> > >> Kind regards,
> > >>
> > >>
> > >> Estanislao
> > >> -------------- next part --------------
> > >> An HTML attachment was scrubbed...
> > >> URL: <
> > http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200114/3f0a8587/attachment-0001.html
> > >
> > >>
> > >> ------------------------------
> > >>
> > >> Message: 2
> > >> Date: Tue, 14 Jan 2020 11:23:12 +0100
> > >> From: Martin Losch <Martin.Losch at awi.de>
> > >> To: MITgcm Support <mitgcm-support at mitgcm.org>
> > >> Subject: Re: [MITgcm-support] ECCO regrid interpolation
> > >> Message-ID: <BB97573C-E6FD-4A1A-8BEA-4890AAE003F0 at awi.de>
> > >> Content-Type: text/plain; charset="utf-8"
> > >>
> > >> Hi Estanislao,
> > >>
> > >> wouldn?t the natural C-grid way of computing EtaN at U points be
> > >> EtaNatU(i,j) = 0.5*(EtaN(i-1,j)+EtaN(i,j))
> > >> and similarly at V points
> > >> EtaNatV(i,j) = 0.5*(EtaN(i,j-1)+EtaN(i,j))
> > >> ?
> > >>
> > >> Martin
> > >>
> > >>> On 14. Jan 2020, at 02:53, ESTANISLAO GAVILAN PASCUAL-AHU <
> > e.gavilan at hhu.edu.cn> wrote:
> > >>>
> > >>> Dear community,
> > >>>
> > >>> I was trying to regridded the sea surface elevation data to the u and
> > v points in order to compute the normal transport across each cell. At the
> > beginning I thought I could use griddata from MATLAB or python. However,
> > since my grid is curvilinear I think I cannot do such thing. Does MITgcm
> > have any kind of tool for regridding this type of data? Thanks in advance
> > for your help
> > >>>
> > >>> Kind regards,
> > >>>
> > >>> Estanislao
> > >>> _______________________________________________
> > >>> MITgcm-support mailing list
> > >>> MITgcm-support at mitgcm.org
> > >>> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> > >>
> > >>
> > >>
> > >> ------------------------------
> > >>
> > >> Subject: Digest Footer
> > >>
> > >> _______________________________________________
> > >> MITgcm-support mailing list
> > >> MITgcm-support at mitgcm.org
> > >> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> > >>
> > >>
> > >> ------------------------------
> > >>
> > >> End of MITgcm-support Digest, Vol 199, Issue 12
> > >> ***********************************************
> > >
> > > _______________________________________________
> > > MITgcm-support mailing list
> > > MITgcm-support at mitgcm.org
> > > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> >
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> >
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200115/13b71ca8/attachment-0001.html>
>
> ------------------------------
>
> Subject: Digest Footer
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
>
> ------------------------------
>
> End of MITgcm-support Digest, Vol 199, Issue 13
> ***********************************************
More information about the MITgcm-support
mailing list