Re: [ferret_users] post processing (interpolating) MOM4 (CM2.1) output with CURV_TO

Re: [ferret_users] post processing (interpolating) MOM4 (CM2.1) output with CURV_TO_RECT

To: Malte Stuecker <stuecker@xxxxxxxxxxxxxxxx>

Subject: Re: [ferret_users] post processing (interpolating) MOM4 (CM2.1) output with CURV_TO_RECT

From: Karl Smith - NOAA Affiliate <karl.smith@xxxxxxxx>

Date: Wed, 2 Nov 2016 05:56:04 -0700

Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa-gov.20150623.gappssmtp.com; s=20150623; h=mime-version:in-reply-to:references:from:date:message-id:subject:to :cc:x-original-sender:x-original-authentication-results:precedence :mailing-list:list-id:list-post:list-help:list-archive :list-unsubscribe; bh=C+KFDuFOVsQxUeD1ypcxJKKeIBz7kdyVPkKJ/Kq5KEk=; b=ZjtmNfTZ0uW8c+VHHHgrq1iO/8bCUSzoMt71Qo/K7iinG+ZiSIoWempgpU6Xu0sYlh BAQcN53Kb9kWG8EUPCiL61hkiVCGnp6IQwTaD7xixMBM4Sp43mOINd0apjgQq18k1ISQ H9uZeuywpKH5c1Rgn65sNlxST6mjIOHiEjN0JaJr5cr0pEAX50RY+8LTlNTs9z878j4C IkGbIsoXtXrXxOcK2wcC8jDhUF/9G7eNG4voq7bQ71V2VBJfC8ZT3UfC4lbuaF+VWByW SDzTyHAlNUWm+CvSMOAwiUlqTAUUCMX2wTXIgQFYFRag/pM3XXyHe8Ogf+/WpsHRevHy ChXQ==

In-reply-to: <CAOi8gWhfDmfcgUHFarqKhBkOFOhFk7CYio=SCTTXziiK2cyMRg@mail.gmail.com>

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>

List-unsubscribe: <mailto:googlegroups-manage+809428289204+unsubscribe@googlegroups.com>, <https://groups.google.com/a/noaa.gov/group/ferret_users/subscribe>

Mailing-list: list ferret_users@xxxxxxxx; contact ferret_users+owners@xxxxxxxx

References: <CAOi8gWhfDmfcgUHFarqKhBkOFOhFk7CYio=SCTTXziiK2cyMRg@mail.gmail.com>

Sender: owner-ferret_users@xxxxxxxx

Hi Malte,

Is this on a Mac or Linux system? I assume you are using Ferret and not PyFerret.

Regards,

Karl

On Tue, Nov 1, 2016 at 5:39 PM, Malte Stuecker <stuecker@xxxxxxxxxxxxxxxx> wrote:

Dear mailing list,

I am currently postprocessing MOM4 (CM2.1) output following this example:

http://ferret.pmel.noaa.gov/Ferret/documentation/users-guide/appendix-a-external-functions/CURV-TO-RECT

Having the following variable on the curvilinear grid:

name title I J K L
CELL_AREA
cell area 1:360 1:200 ... ...
COSROT COSROT,SINROT points east 1:360 1:200 ... ...
GEOLAT latitude 1:360 1:200 ... ...
GEOLON longitude 1:360 1:200 ... ...
SINROT -SINROT,COSROT points north 1:360 1:200 ... ...
CN ice concentration 1:360 1:200 1:5 1:900

I followed the instructions in the provided example as:

yes? let lonin = geolon[d=1]
yes? let latin = geolat[d=1]
yes? define axis/x=0:359:1/modulo/units=degrees xax
yes? define axis/y=50:90:1/units=degrees yax
yes? let lonlatout = y[gy=yax] + x[gx=xax]
yes? let map = curv_to_rect_map(lonin,latin,lonlatout,10)
yes? save/clobber/file=curv_map_NH_1deg.nc map

yes? cancel var/all
use curv_map_NH_1deg.nc

2> ./curv_map_NH_1deg.nc (default)
name title I J K L
MAP CURV_TO_RECT_MAP(LONIN,LATIN,LO 1:360 1:41 1:4 1:3

However the generated map matrix has zero values everywhere. I assume this is not correct? I get the same result when using a smaller radius (e.g. 3 degrees) in the interpolation function:

yes? let map = curv_to_rect_map(lonin,latin,lonlatout,3)

I am wondering what I am doing wrong here.

Thank you very much in advance!

Karl M. Smith, Ph.D.
JISAO Univ. Wash. and PMEL NOAA
"The contents of this message are mine personally and do
not necessarily reflect any position of the Government
or the National Oceanic and Atmospheric Administration."