[ferret_users] Re: Extract time excursion

To: ferret users <ferret_users@xxxxxxxx>

Subject: [ferret_users] Re: Extract time excursion

From: "Fiedler, Russell (O&A, Hobart)" <Russell.Fiedler@xxxxxxxx>

Date: Tue, 19 Jan 2021 06:44:28 +0000

Accept-language: en-AU, en-US

Arc-authentication-results: i=5; mx.google.com; dkim=pass header.i=@csiro.au header.s=dkim header.b=Z6mnnj76; dkim=pass header.i=@CSIROAU.onmicrosoft.com header.s=selector1-CSIROAU-onmicrosoft-com header.b=bCkkPcRj; arc=pass (i=3 dkim=pass dkdomain=csiro.au dkim=pass dkdomain=csiroau.onmicrosoft.com); spf=fail (google.com: domain of prvs=6469f4a1c=russell.fiedler@xxxxxxxx does not designate 209.85.219.69 as permitted sender) smtp.mailfrom="prvs=6469f4a1c=russell.fiedler@xxxxxxxx"

Arc-authentication-results: i=4; mx.google.com; dkim=pass header.i=@csiro.au header.s=dkim header.b=Z6mnnj76; dkim=pass header.i=@CSIROAU.onmicrosoft.com header.s=selector1-CSIROAU-onmicrosoft-com header.b=bCkkPcRj; arc=pass (i=3 dkim=pass dkdomain=csiro.au dkim=pass dkdomain=csiroau.onmicrosoft.com); spf=fail (google.com: domain of prvs=6469f4a1c=russell.fiedler@xxxxxxxx does not designate 209.85.219.69 as permitted sender) smtp.mailfrom="prvs=6469f4a1c=russell.fiedler@xxxxxxxx"

Arc-authentication-results: i=3; mx.google.com; dkim=pass header.i=@csiro.au header.s=dkim header.b=Z6mnnj76; dkim=pass header.i=@CSIROAU.onmicrosoft.com header.s=selector1-CSIROAU-onmicrosoft-com header.b=bCkkPcRj; arc=pass (i=2 spf=pass spfdomain=csiro.au dkim=pass dkdomain=csiro.au dkim=pass dkdomain=csiroau.onmicrosoft.com dmarc=pass fromdomain=csiro.au); spf=fail (google.com: domain of prvs=6469f4a1c=russell.fiedler@xxxxxxxx does not designate 65.153.203.169 as permitted sender) smtp.mailfrom="prvs=6469f4a1c=Russell.Fiedler@xxxxxxxx"

Arc-authentication-results: i=2; mx.google.com; dkim=pass header.i=@csiro.au header.s=dkim header.b=Z6mnnj76; dkim=pass header.i=@CSIROAU.onmicrosoft.com header.s=selector1-CSIROAU-onmicrosoft-com header.b=bCkkPcRj; arc=pass (i=1 spf=pass spfdomain=csiro.au dkim=pass dkdomain=csiro.au dmarc=pass fromdomain=csiro.au); spf=pass (google.com: domain of prvs=6469f4a1c=russell.fiedler@xxxxxxxx designates 2405:b000:e00:257::7:41 as permitted sender) smtp.mailfrom="prvs=6469f4a1c=Russell.Fiedler@xxxxxxxx"; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=csiro.au

Arc-authentication-results: i=1; mx.microsoft.com 1; spf=pass smtp.mailfrom=csiro.au; dmarc=pass action=none header.from=csiro.au; dkim=pass header.d=csiro.au; arc=none

Arc-message-signature: i=5; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=list-archive:list-help:list-post:list-id:mailing-list:precedence :message-id:mime-version:content-language:accept-language :in-reply-to:references:date:thread-index:thread-topic:subject:to :from:ironport-phdr:ironport-sdr:sender:dkim-signature; bh=UT0oC0WUN1TJyw7Bz2MnAStVs1fdy8x0mytzr9nX6Zc=; b=NW/QIJ85hcQtDuCF0I7UbKb9d6S/bOhsHR/4ZQcTpWhyXiDD9EC07cHdYIXgmaOg3k McL9U5WtZqEgHnwuR39GRJqePu0Hxfom2n7BlW+gOr5tX5cpUH0+uKD+kA1XI7HwwdsH ThTXh5DO0XmKUFrxCmvH/ylf+upizkvieyemW0sdgVIyyYObnRjuSpRVvTujvCCm1K9m D0SUGQt2QwcOR4CJC4hiBx8p9oMo5D/w6LeoUv0rgCX1VGYK3H/UjiuBgRHIhi9Y0ePd RbFRRi/ge3Gz9piwJyVf4+lCruj+1pbtCK+8alWKYSKZhgzVmJaUIZWzh4JoTWfqvRzG LbQg==

Arc-message-signature: i=4; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=message-id:mime-version:content-language:accept-language :in-reply-to:references:date:thread-index:thread-topic:subject:to :from:dkim-signature:ironport-phdr:ironport-sdr:dkim-signature; bh=UT0oC0WUN1TJyw7Bz2MnAStVs1fdy8x0mytzr9nX6Zc=; b=DLHPqxLeeNbddZNYFDjCVP6WQfY0FIu/UM4Xti96TjJWoAA+NCwQ+qgKcaZOr/vE3P xywS9lN5jEu0qP2V6iFHWuIJoQUzQP07co7n/xeo0OFQlTfWd0w+WW9BqUuuIloVu/XR f+RndfJLAR4iAZQL7I1GKCFvS5k4q/ufcW57ghe6cFa6RcmGK6iewpbq3OduBwVtXEpZ GB7XVs5AuTFxHjfi4wepNevUtC92IpUcdO7H+AKRxTwfGL/EtxqcQPFDLC2f3WkpQg5r F/WFjFTcIF+j/IfS/yso3c4X8K0fSit7zK3dnXZ/nF3zn2PVQPNaWKzpJnl66OrgR6oo Ryqw==

Arc-message-signature: i=3; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=mime-version:content-language:accept-language:in-reply-to :references:message-id:date:thread-index:thread-topic:subject:to :from:dkim-signature:ironport-phdr:ironport-sdr:dkim-signature; bh=UT0oC0WUN1TJyw7Bz2MnAStVs1fdy8x0mytzr9nX6Zc=; b=YALqKczkQPKL92Gm/oR7IhP54YAQX5ngBIbHxkqgNG2G4c0Mf5LrG27pCgNu3LwlXE CHk5eePZIQjxbMljvwvjzb8FTCzTwT9U0DaBywe6LS0aVxT2M8rNrEsBW9ah4Ppb/GPw FOcAHXHK01WWEiZK+b4NfcUUIHqGTeh73vHfDRKXE4FzpetzWXFKv3cWN4MqUTLU34O/ Uea3XC0TPGeq9cExPPQbk7kIl/hgFIhVoOytamYr65kui4HMxMtdYyPRj18nXcjvk8qB dlE31SIMOLlhbCXX+JMg4s/tWscMJDfBVShaoOXNaKkz/Jea7eKhy9H1gL1aDvwguqi7 oWKA==

Arc-message-signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=mime-version:content-language:accept-language:in-reply-to :references:message-id:date:thread-index:thread-topic:subject:to :from:dkim-signature:ironport-phdr:ironport-sdr:dkim-signature; bh=UT0oC0WUN1TJyw7Bz2MnAStVs1fdy8x0mytzr9nX6Zc=; b=NOBfR8nwSEd+EeMvtOD4LltOhnpQsQRj6QpnFEBboaA7h7eNYhtZ6jnc9f3WZ0A/GB QHJO3baoo0JUspuzFUfLDTRJIlYphwf0dRWMY/bNN3zBWgq6+EOIizsStDsyb0mcdkSt Q9ei+jm8o1ShUf7D6PjDaslgOSbMeQLDmaFUoN7XmiNlQbkOLgUiPjvMlO1/7BDHiiIW wM4/kOY+L/cnKjaMIinrh5675IDMTne2x2XMRAzDiE3oixYen0B4WdxHCqBZT1Lllryo tDrxD9drea31sPY/bkWYqHCibL4+SnOI/fzTbeVpgFeTSdbDVUND8xtM7webi1H6TL1i LH8g==

Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=microsoft.com; s=arcselector9901; h=From:Date:Subject:Message-ID:Content-Type:MIME-Version:X-MS-Exchange-SenderADCheck; bh=UT0oC0WUN1TJyw7Bz2MnAStVs1fdy8x0mytzr9nX6Zc=; b=GgRSlGsCWPLmZPuppe2MeHi+7QJwhWWhMA3dNMcIozVEN3ACJLbeBrTsDiL9pl00QlTcNA+FyRUY2J/1GkolD3xbjHlGKVgMQAORIfQGwcHOnajjeUThs3d5/f3sTLFQVzp5EYh3zAJZynB0FuZuZ881duCocKySFDp6pj/aYDnPTLYIhF8WIgscBejJSHc+fxwqbm4HhkyjOl0Ie6ZxNoGPrEyhUywRB3AK61VRuDwxScPxMC8FCQuQ586xctilwRKQ3y/pB2AHkfMCGJ4UAywvOIios6/IAfTVYEaooK/7tgVtXytu/tNYmN0pSfhW+/GYwtHchFoD1dePRuSQ3Q==

Arc-seal: i=5; a=rsa-sha256; t=1611038680; cv=pass; d=google.com; s=arc-20160816; b=ELbVRLtymeD2SkBWCSfPfM+OIuVMLtH1y597MrcJ0YqKVpy33Pa8kVX4PFe/KeKELs E2y9BOixpXIFIdf2LsWFnhzZ9kw6aT5/jOyg+sQQmfNxQRG+tGMjVx2A2fH2oKsNnqWE EohUsQVZfjpKZhNdXDhEbCR8Nk99qzL1kfudus8wbXXxcvV3ENdbRucuSpNMfpp7tnWP MOnbVnH+EI9+M4AVa64e5+UwO3lHKFwYIhBFkcX9MoJ37aR9v0T4n4WBNgZOojlFxwoA CDDqEAmpMAbOKW+qqcfCwucw4/2iB6wpnJKNI6d0i6QcHOb7swhaptq6PBzwny8Uudb3 Ak2A==

Arc-seal: i=4; a=rsa-sha256; t=1611038679; cv=pass; d=google.com; s=arc-20160816; b=coh8VJlD3qt4qo8cI6jLhzuoHrwcTiS4r53QV++D6ABg7WBadHwae73+7wHGdgKNqj jUEpppBY72uKsoSzEjLvc9FyyIiCNgxK5iACs3jVKrvBj7gIluuhayxhlvPrUHQv6POV VVDpHy4y/Hp7w37NAQMLgkIXZ6tkob9K+qpmAExjNTOVcCIvI46W5xzx6sb1dlCoSysE XpiFEO2iV2VLoD9TxjMULYxRtaZqIHrirzoi5F5hXZGDYcGBXq+hLYv2eyKELoEpEK8G ndei4tj7WuM8AAxRdJ5FeRA58JkzmbdYU8Q1FTaJu/Ur2l1p2ORdumo0O0limZtzQdd7 MUmQ==

Arc-seal: i=3; a=rsa-sha256; t=1611038677; cv=pass; d=google.com; s=arc-20160816; b=BWCbswpbM+ps6T33y2Iqap9KmdpG6MJyFri8hynCforQkMa+tLpZj1WKuh5TDRcjXa LlF5Cpxbu84+pYR5voDLZcwxVkqkaQKnP4yXl2gTC8HZYuhCWrrg8yNOW96n9RiyxfpL jxZHg49bYXTHBs3TgqqZZxihydlI+PuNBAzzXi16FLWMETNu9VF2U4ljuOJdNRXtZT2C xhw1+qzESyL8MzAzGLRrSDb41UJSYWK3ojFSKlFBK6YwEE97FA6KwJvZHCYI8977UZ1d JlLKNjl9HU6WzQ7cSGHj4duPpKoGm8y6NGq2BNMX2heOYpv0zBPDxpRamz9NZfyAJt26 cXBA==

Arc-seal: i=2; a=rsa-sha256; t=1611038674; cv=pass; d=google.com; s=arc-20160816; b=rUzl3axJJQOIEikHxysLk4QWZZhO2n78uV9NP69b/8GJpV5r3k/T0RhZCDeHhOi3lw yuOvlt/LGrs43nWXOMfxp5P58E8BVvqq8v7fM9gfc/ufKEVe4z5maE0OfzUSYZOJVcNs XHhKhCzs/Yqx9PSSVP6GXhK9jIcfNZlSMQ8FOW+sjD3xUCWh80cdOKtIG3cZ0Glvoh/8 TPxk/h9KbgY2Go7zXoN7f5Bn7zQ4BngRYNVewcpcqb83PbrM9l8u3Qbf4ysUuG7SFqkf E297AFUr6uPiVQbT6zmqxJvRFOHF6dibF9+JGxp81yYy8ru6JNw4WQYvZkKIFS8wJBI1 yvmg==

Arc-seal: i=1; a=rsa-sha256; s=arcselector9901; d=microsoft.com; cv=none; b=lwyhHR2DEB4iq2qcoT+HH9TcvwWimbCeMGK6NIRyHJheBxlXHl0gbMFRTui0X8WcaXDo1aTjiU376i/ZfwqTVKncg2MA5EcFk9VWksYPc2x4lblErC3I1DDDkzM+0lYGG9/WCz7gg63T8Oi9ZlXFIQJxdP5JiOON1SGuSf//6RuIxdrBiFsacQPqPRKg7ymyi7jPOKZuIAfttq3C+XrCz3rIPsPMv/CTD/jiS3yQLiWCte+TH+6YZ9MvChB0JX9jnlZbmk3Tr5DDbqwldRyCNfh7mgch+R0p8i7Jwtqd3DWW8AOIYMckdtmIqdKVlXlLhqOM37GN1NY6JtQPF5+4bA==

Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa.gov; s=google; h=sender:ironport-sdr:ironport-phdr:from:to:subject:thread-topic :thread-index:date:references:in-reply-to:accept-language :content-language:mime-version:message-id:x-original-sender :x-original-authentication-results:precedence:mailing-list:list-id :list-post:list-help:list-archive; bh=UT0oC0WUN1TJyw7Bz2MnAStVs1fdy8x0mytzr9nX6Zc=; b=K9xkCqPtKJ4/sCXoiKQKoz1vQaUfaIYagOaAwumfbC9MJmdBXygi0moU1D5p3OX+PH U78zcevZ0wpcxGSD1EOC4FFyuVIP9+OHCSX5dmUzq+aDIyO8bm5OjwRUw/sRjiNL3zJp MiN6sZZnbGxuXYYuNiWiuS6bDQ5girfj4nDcFNCY2y9aJeLHWm7X62lwr4DUOFwbNYRG x+dTfKZ1BvNlS+ZR3hcrmA4WsUVgRjOITpJpduEfJ/GUaWxMriSD7Js0TYCMuyLdv2cx Jn9xxuYYcnCT+TblK93VsNXFfvvOLozcUNyJHRFUjgipff7LQjR97PkWm4Y11Lgh8tiM 723Q==

In-reply-to: <933804681.672195.1611003827873.JavaMail.zimbra@lsce.ipsl.fr>

Ironport-phdr: 9a23:hYho9h8NRXquov9uRHGN82YQeigqvan1NQcJ650hzqhDabmn44+7ZRWN4/JkhRnGWpuIo/5Hiu+DtafmVCRA5Juaq3kNfdRKUANNksQZmQEsQavnQU32JfLndWo2ScJFUlI2+XayOA5aH9itL1HXq2e5uDgVHBi3PAFpJ+PzT4jVicn/1+2795DJJQtSgz/oaqhzNhLwpBiCuw==

Ironport-sdr: DG5KYAbruimdGRQX067+7Qv+DrC8NsSOC/qVETUjc0SN5sHOla/k7BtaMvhmh/X0TPLsLkJwM7 ED3VtN+LA7dg==

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: <933804681.672195.1611003827873.JavaMail.zimbra@lsce.ipsl.fr>

Sender: owner-ferret_users@xxxxxxxx

Thread-index: KeH7LplLXIKbA0fYGyTGM8PvVHIWnunlqg/S

Thread-topic: Extract time excursion

Hi Patrick,

The problem here is the conflict between a regular year axis and the irregular gregorian calendar which only repeats every 400 years.

yes? define axis/t=1998:2009:1/units=years/t0=01-JAN-0000 tyear

yes? sho axis tyear

name axis # pts start end

TYEAR TIME 12 r 31-DEC-1997 12:21 31-DEC-2008 04:22

T0 = 01-JAN-0000

Axis span (to cell edges) = 12

! Show 400 year cycle

yes? define axis/t=1600:2000:1/units=years/t0=01-JAN-0000 tyear

Replacing definition of axis TYEAR

yes? sho axis tyear

name axis # pts start end

TYEAR TIME 401 r 01-JAN-1600 00:00 01-JAN-2000 00:00

T0 = 01-JAN-0000

Axis span (to cell edges) = 401

In your first example there is no calendar implied. In the second a standard/gregorian calendar is assumed. If you had the /CAL=NOLEAP qualifier in your axis definition

all would be good.

yes? say `a,return=calendar`

!-> MESSAGE/CONTINUE GREGORIAN

yes? define axis/t=1998:2009:1/units=years/t0=01-JAN-0000/cal=noleap tyear

Replacing definition of axis TYEAR

yes? sho axis tyear

name axis # pts start end

TYEAR TIME 12 r 01-JAN-1998 00:00 01-JAN-2009 00:00

T0 = 01-JAN-0000

CALENDAR = NOLEAP

Axis span (to cell edges) = 12

Ah, I just saw that Ryo replied with something similar.

Cheers,

Russ

From: owner-ferret_users@xxxxxxxx <owner-ferret_users@xxxxxxxx> on behalf of Patrick Brockmann <patrick.brockmann@xxxxxxxxxxxx>
Sent: Tuesday, 19 January 2021 8:03 AM
To: ferret users <ferret_users@xxxxxxxx>
Subject: [ferret_users] Extract time excursion

Hi ferreters,

I would like to extract at ISO 8601 format

start date and end date of a netcdf file.

Here are simplified lines that illustrate a problem of precision.

yes? define axis/t=1998:2009:1/units=years tyear

yes? let a=t[gt=tyear]
yes? say `a, return=tstart`
!-> MESSAGE/CONTINUE 1998
1998

---> ok

yes? define axis/t=1998:2009:1/units=years/t0=01-JAN-0000 tyear
yes? let a=t[gt=tyear]
yes? say `a, return=tstart`
!-> MESSAGE/CONTINUE 1997-12-31 12:21:36
1997-12-31 12:21:36

---> why this 1997-12-31 ???

Attached a netcdf file where the problem exists.

yes? use toto.nc

yes? let a=t[gt=var]

yes? say `a, return=tstart`
!-> MESSAGE/CONTINUE 31-DEC-1949 21:00:00
31-DEC-1949 21:00:00

PS: tested with last 7.6 release

Regards

Patrick

--
Data Analytics and Visualization Engineer / Project Manager
LSCE/IPSL, CEA-CNRS-UVSQ laboratory
LSCE - Climate and Environment Sciences Laboratory
IPSL - Institut Pierre Simon Laplace
--