[Thread Prev][Thread Next][Index]
Re: [ferret_users] ensemble or compare data from two time series
- To: ferret_users@xxxxxxxx
- Subject: Re: [ferret_users] ensemble or compare data from two time series
- From: Ansley Manke <a.c.manke@xxxxxxxxx>
- Date: Wed, 24 Nov 2021 15:23:11 -0800
- Arc-authentication-results: i=4; mx.google.com; dkim=pass header.i=@gmail.com header.s=20210112 header.b=pE1osAmD; arc=pass (i=2 dkim=pass dkdomain=gmail.com); spf=pass (google.com: domain of a.c.manke@xxxxxxxxx designates 209.85.160.197 as permitted sender) smtp.mailfrom=a.c.manke@xxxxxxxxx
- Arc-authentication-results: i=3; mx.google.com; dkim=pass header.i=@gmail.com header.s=20210112 header.b=pE1osAmD; arc=pass (i=2 dkim=pass dkdomain=gmail.com); spf=pass (google.com: domain of a.c.manke@xxxxxxxxx designates 209.85.160.197 as permitted sender) smtp.mailfrom=a.c.manke@xxxxxxxxx
- Arc-authentication-results: i=2; mx.google.com; dkim=pass header.i=@gmail.com header.s=20210112 header.b=pE1osAmD; arc=pass (i=1 spf=pass spfdomain=gmail.com dkim=pass dkdomain=gmail.com dmarc=pass fromdomain=gmail.com); spf=softfail (google.com: domain of transitioning a.c.manke@xxxxxxxxx does not designate 65.153.203.166 as permitted sender) smtp.mailfrom=a.c.manke@xxxxxxxxx
- Arc-authentication-results: i=1; mx.google.com; dkim=pass header.i=@gmail.com header.s=20210112 header.b=pE1osAmD; spf=pass (google.com: domain of a.c.manke@xxxxxxxxx designates 209.85.220.41 as permitted sender) smtp.mailfrom=a.c.manke@xxxxxxxxx; dmarc=pass (p=NONE sp=QUARANTINE dis=NONE) header.from=gmail.com
- Arc-message-signature: i=4; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=list-archive:list-help:list-post:list-id:mailing-list:precedence :in-reply-to:from:references:to:content-language:subject:user-agent :mime-version:date:message-id:sender:dkim-signature; bh=GBcjKdix33oDljsYRCQYWElCjopnhugn58Ty8S6PKjc=; b=sFxWfgOxOKCj8ONt1ffTJPFHdhshnlL3zLwO9rAR/r9i1j09Nc3NrTmmRgiupPrjN2 ffpdxDQYPPwqdqw03RIvNWe5w2pUb/+RqzNyL8mRw5jU3VJZ+emGn2aeJRj/sXB0KjxC V2bYNOp/YhJBEWQMMUA4oflpkD/qIqGVtjIi/Jc81Qg4mSdyxQkIetNgr7EsZr0wGFAU PqyIiHuEbdkxpL7wl4rMKB130PEciP22v1k+tN+XIGu36C87Jw17CCwVne3EEmrR3aLF jVWQTuwhNfEhM/m0NXC6Z4r7OhJ8v9CSmTCHMvqrOU5T1fO6hIu2cZN5rO+X9mApYsgI g3sw==
- Arc-message-signature: i=3; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=in-reply-to:from:references:to:content-language:subject:user-agent :mime-version:date:message-id:dkim-signature; bh=GBcjKdix33oDljsYRCQYWElCjopnhugn58Ty8S6PKjc=; b=Tm/LsUZHYIebfFfF6Ohbgg8eqvQSkAhGJdBiSkIbWxnZ/PFYEfj0GywmyoICRyI/Nc y9OnDEk+6mzifxkDXLLPJQgqJe02rcmq+vPtE+AJh9C153Z3FGhEwDd43G69WmAhY6Ih /8b3jDoMg6Y5IT1WMmfphZmGrA+d4P0OWBaGuEgEJnYthH1QTAQYg9MKErhTdN8FoeaL fCZ+Km9qRb6b0dBxgU4mOIJBfdLVEr8+5UgCYk10ETtdm7qlXMu7j+N3zUu6yd/4cLt1 Hib+gfmuCqiiPd9WdpjRrGX1xiCH1Pelryn/aiOeubJEqtWkwduibgFfajFidw02EN6Q pOaQ==
- Arc-message-signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=in-reply-to:from:references:to:content-language:subject:user-agent :mime-version:date:message-id:dkim-signature; bh=GBcjKdix33oDljsYRCQYWElCjopnhugn58Ty8S6PKjc=; b=Y75a75i2Gklknn/Uwp8LzVqwXSl938U1isa7QeWB6HndlgZ9HWYzlxx3dQGIJkFW9u kJIJo4/YpW1oAdB20ZaCYtVjTa1Q4BAMPf8CLBzPGKCYY6q9PabH5BnP0uvOwST5WPhp m6bnH6TonoqE/hpD5D7KKNyclkIYO8/jWumW4MSsh7OBOz/OyOrExiF+RApT9GqmyUPF /j6pw5mNlqeEMOVCxXDkaoxl9LCZ88v6dzRnWGK43Wd0dos9uFmjLdDX0tmVHjwQBXmV yS728WVPaFcq699ijwrJWwqarnQ3F3x0fA+VIqubrVp6rfF/ri7aJ78ewfpBGY1vGopz 0uGg==
- Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=in-reply-to:from:references:to:content-language:subject:user-agent :mime-version:date:message-id:dkim-signature; bh=GBcjKdix33oDljsYRCQYWElCjopnhugn58Ty8S6PKjc=; b=0PNrXYZT0xSN4DeY9yooE0SFgxeT8jrBeLcIjct0konJMfQtTCF3yZTWzf+cRmA28t C2z3ax/CeICetxeEy9n2bYbGVpicRyCMPrWCJVyw1CCbPTfTg6HkTp3LOfe6JvfZ+G2e bfWgDAgKJtabwTM1HiUm+PMjKvGFoxVGi/S2hRj/9uT9D4yWtd3pU0x8N6jPM9HNT1fI WqtMP/j/TC40GjD6o5iH0IJfemssGm6SJRxI4/NYKoloPtQf+OcHTJyKN4htxUbOLKGY vr6P+6DPIfBMT6uiHIGv4JZevxZFtj+zzh3u1q2yYY3Ue+LM/R+vXJLg6Avwblk9C8zu 06tw==
- Arc-seal: i=4; a=rsa-sha256; t=1637796196; cv=pass; d=google.com; s=arc-20160816; b=05wHL9KPfC+DDrsn3juWGZVLXtxUu9n7FxAYb/ztlSsw4Ij5bg6RE5eKtIs2ETgtY4 pLvH4cjX5OnftTEifPLo1wa3FjKlt4f6VxMwAMw7XCUgvWvVWfl/uQHgJNzuc4A1Rtmv zgx3N2olcz+D/cYjdO1fATXAk4pu8NiNwE5h/ciWm8pVMzAQ7BzAaJ+S3ck8U312eAJI gAIL86H+argO3CVXeKJF5SZhxqChh87PbwFEh9FekOd2gwbqREQYkrc4VVG1iQ/Cei1+ 1MkH2R5oVND4Q7o/LnbOYZhbvC3wEmxJnARzzq+oUavfMYaXxY8GhEnV3JEgCrfdu45f 8yTA==
- Arc-seal: i=3; a=rsa-sha256; t=1637796195; cv=pass; d=google.com; s=arc-20160816; b=t0V82zahewX3xqg2l+83PbHP/lQxhlFC+QutzLu6rcIwfddWCVBVaPLtXVXjTKClVa 5Mt+XHFRE06IevM3jojC13D4RTQ3TqlmR0Yyol+Qo3rmRxh15wom0JNqJm8l2y7+VxLC gHQCBaN+mH5CBEyLPtxYJgu1ieIBUixXBSG7upFicXfbKxgKX6TKJudsX+8n1bOJhNzE eiSwlc6tH3NzA2S4eysCGEgNOykP/JMr7rzKQ+cOssQrVezfNi6mVWyDOjA7X5LdaGbt +wZGzAavUsQg0dRxkgFDzCqMPRdougQd3FIq5747c0GMPYm3pToNkNMjKGLPZfPY6svC 20iw==
- Arc-seal: i=2; a=rsa-sha256; t=1637796192; cv=pass; d=google.com; s=arc-20160816; b=OwUjQz4iBPuxdWx67Pli8XDSXaW8jetIzFUTJ+Jl7ABuR80Fq6FYSfTZFZ2briBaRK OAKrYqy9wVpZa33z4DYanRK0h0tNyT/7rrsxDBmCkXvD/1tCqVBsSmus0F0utj1tpdD7 KzMtJYU911OGqhBgVAOEOyIOaYWylYBHoHJxU334qXajIlih5HKyJZ9WQYihXEwwmVB6 RNQCaxqu3YsXSxWBsXfAAPYQt+gk103gwvTZdSdCfDcyYHXW2S/b6IQI7y1TglVrp26p 18ddi8PSOboRmggHZeDNTxtHet7zUFC5aP4dem2LiNYUP8Gb1bNr8padSgmgSaHQjBAx SXgQ==
- Arc-seal: i=1; a=rsa-sha256; t=1637796191; cv=none; d=google.com; s=arc-20160816; b=f73l74/gV6nvfyG/65QNIFyxMnzsy2xFzKFlJWuYKN/jw/ku2VSBSzSs7DtP7qVPFV PZ6LAPm10uKlEisjDvwfoGHm2pMHAQmj9JoBds6Vcq6ue0S75bca/WfyLdvWBy3R0zOe fhfYkXc524t2rsWoH959mrX2w0jb9mBy89bRA35BVEn3gVtWQNH/v1wAPWtPMbVGeCeB 9vrZWigfxp11b36zLumvRkqgH2zRQKoHlS6HEb3+fBcNV+Ej8Q7YfCanERx43/UZWj6+ xZC9xtTkcwZt+foxxzFFYyq2Qr/jB+QOPrZ6fBhmnuggUB56ZJbLv/PIhfpSzdE7HZO3 B0Ig==
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa.gov; s=google; h=sender:message-id:date:mime-version:user-agent:subject :content-language:to:references:from:in-reply-to:x-original-sender :x-original-authentication-results:precedence:mailing-list:list-id :list-post:list-help:list-archive; bh=GBcjKdix33oDljsYRCQYWElCjopnhugn58Ty8S6PKjc=; b=H+g1by2Plwdfb2R4eHnacU6UHeMxc3MhbMwJSIAcMWM9QgrtExNxWgpgT+cQY2iyno pdJbeqMml38GB/1q0Uaw0KhCsHJsof83cqKQIJFuD1Q/DgN2NvrEVE8fMGLfCFW8LFlc 7nZeHbtsW/FFKhs1KO+WkPHvLb4+BVkGGXzlGy+SM59P9x57fwdiyfgPJB6BUz1CU5Cf FEQ3La1ecrZasiKmNCpCz+IGBOIWAZltVQbS9uxdcF8rqaRfQOkJ6qGdaP8K8FvbCZNP fpCVFADCfHkcUBkE3Pbrk7MwEVHJQd6hprwxYC3xxqXs8GewTBa34lnno8bf+XFPSy+i yy+A==
- In-reply-to: <CAKE4FWqUQPj5WmRwNZDxHj2X2pMdUDFMo_70JdVepXJMJUVrFw@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>
- Mailing-list: list ferret_users@xxxxxxxx; contact ferret_users+owners@xxxxxxxx
- References: <CAKE4FWqUQPj5WmRwNZDxHj2X2pMdUDFMo_70JdVepXJMJUVrFw@mail.gmail.com>
- Sender: owner-ferret_users@xxxxxxxx
- User-agent: Mozilla/5.0 (Windows NT 10.0; Win64; x64; rv:91.0) Gecko/20100101 Thunderbird/91.3.2
Hi,
Your datasets do not share a time range, so this is not a typical
ensemble.
Usually, you can define an ensemble when the datasets do not have
the same space/time grids, by first defining a regridding
operation so that all datasets have the same space/time grid as
the first dataset in your list. On the page,
https://ferret.pmel.noaa.gov/Ferret/documentation/users-guide/commands-reference/DEFINE#_define_data_agg
scroll down in that page to find "Ensemble Examples". Example 3
shows how to define an ensemble for data on two different grids,
including for this data a modulo regridding operation to define
variables with consistent grids. The variable names must also
match, so that example also shows how to rename the file variable
within the Ferret session, and then define the new regridded
variables using the desired names.
I am not sure what commands you tried using @AVE. If you did a
time average, then the resulting grid no longer depends on time,
so there should not be an error message about time axes. Look at
this example:
! Make a file containing the variables in
the first dataset, averaged in time.
! An ensemble needs one file with the grid that is going to be
used in the
! ensemble; not all of the grids can come from LET definitions.
yes? use dataset_1
yes? let/like=precip p_tave = precip [t=@ave]
yes? save/file=data1.nc p_tave
yes? cancel data 1
! Now open this data and define variables
in the other dataset(s) so that
! they can be included in an ensemble.
yes? use data1.nc
yes? use dataset_2
yes? let/d=2/like=precip p_tave = precip[t=@ave,gxy=p_tave
[d=1]]
yes? ensemble time_ave = 1,2
! The ensemble dataset should contain the variable p_tave with a
grid in X-Y-E
yes? show data time_ave
On 11/18/2021 3:22 PM, Xiaoyu Bai
wrote:
Dear ferreters,
I am trying to ensemble RCPs (2005-2099) and control runs
(1949-2005) together but keep getting T axis differs error.
*** NOTE: Variable "PRECIPITATION" excluded from aggregate:
T axis differs among members
**ERROR: error defining aggregate dataset: No valid datasets
or datasets share no variables that can be aggregated.
It makes sense for pyferret to think the time axis are
different. But comparing the @AVE value from different time
series should also make sense to pyferret. I am wondering how
to solve this problem.
By the way, I already tried to read data in separately and
do @AVE first then compare. It still showed the same
error because the data still come from different time series.
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
Ferret
Privacy Policy | Disclaimer | Accessibility Statement