|
Hi, You'll need to subscribe to the Ferret Users List before you can
send questions to the list. I have gone ahead and added your email
to the list, but you will need to write a new message to
ferret_users@xxxxxxxx . You should have already got one or two
email messages saying you have subscribed. I have not studied your question in detail but I see one or two things for you to try. If you want your xax and yax to represent longitude and latitude, use the /units qualifier on the DEFINE AXIS command, /units=degrees_east and /units=degrees_north. Other than that I don't see a problem with what you have done, but instead of RESHAPE, you might instead do a regridding in just the X and Y directions, using the @ASN regridding transformation. Ansley -------- Forwarded Message --------
From ferret_users-owner@xxxxxxxxxxxxx Mon Jan 15 23:46:00 2018 Received: from mail-pl0-f69.google.com (mail-pl0-f69.google.com [209.85.160.69]) by eclipse.pmel.noaa.gov (8.14.4/8.14.4) with ESMTP id w0G7k0MA034504 for <ferret_users@xxxxxxxxxxxxxxxxxxxxx>; Mon, 15 Jan 2018 23:46:00 -0800 Received: by mail-pl0-f69.google.com with SMTP id m1sf5474400pls.20 for <ferret_users@xxxxxxxxxxxxxxxxxxxxx>; Mon, 15 Jan 2018 23:46:00 -0800 (PST) ARC-Seal: i=4; a=rsa-sha256; t=1516088760; cv=pass; d=google.com; s=arc-20160816; b=ID1iq6BgpmBO2r+Tq7wUYfs1VNEsCpgvHhJtB7OS+vAIITWyBa0bz30xkQChNrPfBy QlGdUM4HtuRYIAt8OKVRokn0kKZ1jEyIoZHYAaWVPhaTQgX3WXC6fmdgeGg2AdcLPeKf BWn1/mp8nhQMCaPOSEV7PSTLGXcokbDx9gYLgA882BSPkd0FFdB+Cn65cqPdtdjZm7iH zvCCkxmphEEVMOrGFqMd2E5MU/SBuugCrUDTEGv8v6q5Pa9NfihvRMX1seHf6XjTHftY GA5yYWm4jktBNDCGSlE3A0YBdjjLYUxVONa3xWhBm3nE28SJfUwZg5P9CJ4pUQ0IN8+b XMyw== 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 :content-language:content-transfer-encoding:mime-version:user-agent :date:subject:from:to:arc-authentication-results :arc-message-signature:arc-authentication-results :arc-message-signature:message-id:arc-authentication-results :arc-message-signature:dkim-signature:arc-authentication-results; bh=ttU0QeIUuT+NdOXxg4fuHzX2ve63Ph0zFEJ3cwSc/fI=; b=1J/pa6oZiWR09nernvXICsfJR5Bo7JhbJUFFZZ21brpkQ9SJzMnchAQe4vrmsOrXTr LztgQ+Aq/e2ikdE7TE3aJDzcOtTSl9e5jyA+xeFX76NYSAjNNdy/joz7nPV4eARieZpp Pfxf3qLCXgCJCX8ZY/V6JLMsBUtQjphNIs2mM2HMvoII27J3+/bxDQqeSKQ0SRPqYMHm xbfHLX4TfN38EwMrqG0iB7qQ8vhfiwq5jGXD2Pc2l0dcwgkQqWEAm6213rv6VpKVEuje uwMr1oACKhAUlpJx4y3I96En7sQkspJ0FrVlaFKW+z3683hfOlqiM/ALkRV89IsQrypu nxXQ== ARC-Authentication-Results: i=4; mx.google.com; dkim=pass header.i=@gmail.com header.s=20161025 header.b=r9SE+rFn; arc=pass (i=2 dkim=pass dkdomain=gmail.com); spf=softfail (google.com: domain of transitioning mgocean@xxxxxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=mgocean@xxxxxxxxx; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=gmail.com DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa-gov.20150623.gappssmtp.com; s=20150623; h=message-id:to:from:subject:date:user-agent:mime-version :content-transfer-encoding:content-language:x-original-sender :x-original-authentication-results:precedence:mailing-list:list-id :list-post:list-help:list-archive; bh=ttU0QeIUuT+NdOXxg4fuHzX2ve63Ph0zFEJ3cwSc/fI=; b=Yr+pKUSaeQGK5LAOONsro4fU8Bsuee35Yuw2F+4FaWXadKdg3a69+TEzpuD6gnMlO6 qhBa8kIYVv5ezrlCUm/D831RKIjopcNQweEzB6QTWtiNWZnG9MVTHqT+bks4zoYhZoLN rNN5q7BXaT7wQlQEhPCYifmHWXQBzjOxrcDhYyQM5XKNf2n7DPcAfl+j/xw9jid8mebY d8OWknl/CCP9WBdrueuTmptkMlWGeD76ICbYj0e2b0Y6dI3WmXQ3yTtZy6/bKi1r3nCx MR7XInqVXWR0DA0B9ONWYywWHd94JAnhC7yrIKhnVUX84NeRfzCHcMN4HbufQ1MXCb+1 7lWw== X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20161025; h=x-gm-message-state:message-id:to:from:subject:date:user-agent :mime-version:content-transfer-encoding:content-language :x-original-sender:x-original-authentication-results:precedence :mailing-list:list-id:list-post:list-help:list-archive; bh=ttU0QeIUuT+NdOXxg4fuHzX2ve63Ph0zFEJ3cwSc/fI=; b=gG73K4zhb0Gey7x9SPIwd/nohjhHVDZ8k71DdsJb6XBN6I3lRqPGqrB23YR0H58E12 2n9gBSefDLIA2Fec7s9aJFu+5f4PP1WZ004/d56nWtefzLdW54JS67UNMRCvRZ28rIv1 x9LWIHQEozMttX6PDuvHYOm0dhzBB5j5ZWjU83bsq5W2RvrJroOks+n3R6QgMzzQkNVE HCqWeu3MMll3uM1U7KNVLTlJBCaUzz2DdwWJYM0T6hXKUZEpmVbyZmRxUKCpcAd7ArTU MGRpmkKwXnbdGTuqzikMrPWMsMxpksswWe/b880S2KdC9s1dcRnvmyVuZG8z6CvMzRqS llMA== X-Gm-Message-State: AKwxyteP48QWk/0odsd4Hf80c8r8vZzCHeJUhddQvXd5IxPy4eKSAQM2 nBhXIY4Qd6tq9u/AhsCPsePXwWps X-Google-Smtp-Source: ACJfBosfBttMht17WWhG/qd3aJ6v3yvLxLXyUTwgyBe6hF4hPtVgAuHvKCEoZUaxrlOW3LUxN2kfJg== X-Received: by 10.99.49.80 with SMTP id x77mr1605685pgx.52.1516088760622; Mon, 15 Jan 2018 23:46:00 -0800 (PST) X-BeenThere: ferret_users@xxxxxxxx Received: by 10.159.204.146 with SMTP id t18ls299949plo.10.gmail; Mon, 15 Jan 2018 23:46:00 -0800 (PST) X-Received: by 10.99.177.76 with SMTP id g12mr7148322pgp.269.1516088760278; Mon, 15 Jan 2018 23:46:00 -0800 (PST) ARC-Seal: i=3; a=rsa-sha256; t=1516088760; cv=pass; d=google.com; s=arc-20160816; b=dDHKY/dTVTvxQqYWix9X9EcBUUSbwrfeZsaPP4dtKwafO3ZL3o7Z2T+jthHPmVnr70 JhsPvaxAKDn0q8FmOUl9LgN4+i3ocbhuwJDDzDtH1T+KhfUjVZ6ThmAd9dIIOtxFqO10 VRNniSzPFZJ34Bzia8Ef1ieEBrPX6OY1/Bj42BDLO4oJebSFwWVLCtjN6z8oygJGRLSC knFmGaOa3GPvgadvja46qeLGTqD44cmXjl6YsHloFxVttgfdvc1gCPKc0Gq0+4j8vgCt UE2JmKIg61ck9csjo+6PnW9a5oLBFxAx+CFk+okGFpeGaiiSi6BFxOL0FigKSd1doSCj X2wg== ARC-Message-Signature: i=3; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=content-language:content-transfer-encoding:mime-version:user-agent :date:subject:from:to:dkim-signature:arc-authentication-results :arc-message-signature:arc-authentication-results :arc-message-signature:message-id:arc-authentication-results; bh=ttU0QeIUuT+NdOXxg4fuHzX2ve63Ph0zFEJ3cwSc/fI=; b=NV8KbAR1VXoC/EBlYq8H+u/761fZUBxqvpX2EgMStKo5aBoqjfCIVQafpBs02tHZ+A /0J5MmjLNs9nT91PqzPrDhQH4sTqVZrsDP6h2QY305YWCfwwJgox5gRKGozE4ChiCShd f/KZUXWfyispIv7G4Padlc7frlzPUOgpY9mG5EVDcpZd9V0Io76Lo8cFy0npX4utSGit M35SIUS0WD9WWG45AEUwJo9y+lHvrPmoMbRrz4vTUbxKvCQxbVZQZHM3Jt/FWRxkhYBu rNh+GTcLhfcsSC6W8ZCbuAGANtC7jIvwvvOydxblJaHePNUaaWtwpNaPzo7UmCBZTelx R3VQ== ARC-Authentication-Results: i=3; mx.google.com; dkim=pass header.i=@gmail.com header.s=20161025 header.b=r9SE+rFn; arc=pass (i=2 dkim=pass dkdomain=gmail.com); spf=softfail (google.com: domain of transitioning mgocean@xxxxxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=mgocean@xxxxxxxxx; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=gmail.com Received: from mx.us.etp.fireeyegov.com (mx.usgov.email.fireeyecloud.com. [52.222.5.176]) by mx.google.com with ESMTPS id m61si1421277plb.141.2018.01.15.23.46.00 for <ferret_users@xxxxxxxx> (version=TLS1_2 cipher=ECDHE-RSA-AES128-GCM-SHA256 bits=128/128); Mon, 15 Jan 2018 23:46:00 -0800 (PST) Received-SPF: softfail (google.com: domain of transitioning mgocean@xxxxxxxxx does not designate 52.222.5.176 as permitted sender) client-ip=52.222.5.176; Received: from smf-etp-govcloud-emps-prod-06.fireeyecloud.com ([10.87.36.250]) by USGOVA04 with SMTP id bLwVeunE48g65bLwWetfzg; Tue, 16 Jan 2018 07:46:00 +0000 Message-ID: <579a15a2-5a5a-5d2f-04e9-e81e1902019d@xxxxxxxxx> Received: from localhost.localdomain (localhost [127.0.0.1]) by smf-etp-govcloud-emps-prod-06.fireeyecloud.com (Postfix) with SMTP id 3zLMj35Vvfz2VKQp for <ferret_users@xxxxxxxx>; Tue, 16 Jan 2018 07:45:59 +0000 (UTC) Received: from mx.us.etp.fireeyegov.com (unknown [10.87.38.253]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (No client certificate requested) by smf-etp-govcloud-emps-prod-06.fireeyecloud.com (Postfix) with ESMTPS id 3zLMj33m7Tz2VKQk for <ferret_users@xxxxxxxx>; Tue, 16 Jan 2018 07:45:59 +0000 (UTC) Received: from mail-qk0-f198.google.com ([209.85.220.198]) by USGOVA03 with SMTP id bLwUede9vOMWabLwVehOvu; Tue, 16 Jan 2018 07:45:59 +0000 X-FE-ETP-SENDER-IP: 209.85.220.198 X-FE-ETP-CONNECTING-IP: 209.85.220.198 Received: by mail-qk0-f198.google.com with SMTP id v195so7441912qka.10 for <ferret_users@xxxxxxxx>; Mon, 15 Jan 2018 23:45:59 -0800 (PST) X-Received: by 10.237.49.37 with SMTP id 34mr50471629qtg.36.1516088758241; Mon, 15 Jan 2018 23:45:58 -0800 (PST) X-Received: by 10.237.49.37 with SMTP id 34mr50471620qtg.36.1516088758099; Mon, 15 Jan 2018 23:45:58 -0800 (PST) ARC-Seal: i=2; a=rsa-sha256; t=1516088758; cv=pass; d=google.com; s=arc-20160816; b=BvaLH67AhH8tKtilMiBOtW1SkTMmcqulh+0Cavxj9ffaRGodOxhvcOPRTukh/FjR5N 7NLZJ/t1jNQCl9eMAteMW7HGi2BvKUa+fBL/3e7jBz668+EJyX9KgB6InHeylzJ6efo0 NpvPGTbQowYty5n0zZdTgnoyPz5ghz20/mL7Q+NrPgYupj4mYJ5rPkV7yzhjzXWOtw/G NRCdCq1hihgbtGAzj8A1FoxCZUzUXwbfSm5vsTs31dngpX/Lt9BMkSmlv9ppY/K5LUts Q+KM0ws1NWu7xcxLxZmsQUvPYLoekR7YyQtSRdZY+MlmBtg8S2911qJ3xBkZbggTmrDj jFng== ARC-Message-Signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=content-language:content-transfer-encoding:mime-version:user-agent :date:message-id:subject:from:to:dkim-signature :arc-authentication-results:arc-message-signature :arc-authentication-results; bh=ttU0QeIUuT+NdOXxg4fuHzX2ve63Ph0zFEJ3cwSc/fI=; b=K2sGKcYL6hItfPqvHSaToymx/nQbjs6GuKPM78zSZjphSVjv2zd7Gi0KHmzUsB+ge8 AF7sUyidTWe+ljPJb07Ofn+uqysxfEQfvqeqXP7mGVkwuTNwHnKcfFIb/WvNbK0MGQ9l JCk7rXp69AqHniPA/B8VD6mZTUqG/VB158H8Z0l+a1DA6wXDFH6SmBrM45K/vktYFDcg BYAUxXs/oqIu3hlm3s6Zfn3SrFtT3BGTHN5+qxH+CT6y3v8Q12lv60YZfC82pFzN+PgT 4vIKn5BEdmZYV4GpgUj8Gblni8txmULAEbCRtG/tNxenry5QMvWFnjsstiup3BE8VIZS lgQA== ARC-Authentication-Results: i=2; mx.google.com; dkim=pass header.i=@gmail.com header.s=20161025 header.b=r9SE+rFn; 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 mgocean@xxxxxxxxx does not designate 65.127.216.229 as permitted sender) smtp.mailfrom=mgocean@xxxxxxxxx Received: from ipss1-host81.is.centurylink.net (ipss1-host81.is.centurylink.net. [65.127.216.81]) by mx.google.com with ESMTPS id 3si1769435qkb.289.2018.01.15.23.45.57 for <ferret_users@xxxxxxxx> (version=TLS1_2 cipher=ECDHE-RSA-AES256-GCM-SHA384 bits=256/256); Mon, 15 Jan 2018 23:45:58 -0800 (PST) Received-SPF: softfail (google.com: domain of transitioning mgocean@xxxxxxxxx does not designate 65.127.216.229 as permitted sender) client-ip=65.127.216.229; X-Qwest-Status: hEUBLG9ErOsrnKeZ Received: from ipss1-host229.is.centurylink.net (ipss1-tls04-vpn7 [65.127.216.229]) by ipss1-host81.is.centurylink.net (Postfix) with ESMTP id B55E44BCDE05 for <ferret_users@xxxxxxxx>; Tue, 16 Jan 2018 07:45:54 +0000 (UTC) X-Qwest-Status: hEUBLG9ErOsrnKeZ Received: from mail-wm0-f70.google.com (unknown [74.125.82.70]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client CN "smtp.gmail.com", Issuer "Google Internet Authority G2" (verified OK)) by ipss1-host229.is.centurylink.net (Postfix) with ESMTPS id 0B89A1061132 for <ferret_users@xxxxxxxx>; Tue, 16 Jan 2018 07:45:54 +0000 (UTC) Received: by mail-wm0-f70.google.com with SMTP id g16so1671284wmg.6 for <ferret_users@xxxxxxxx>; Mon, 15 Jan 2018 23:45:54 -0800 (PST) X-Received: by 10.80.196.3 with SMTP id v3mr12467013edf.232.1516088753166; Mon, 15 Jan 2018 23:45:53 -0800 (PST) X-Received: by 10.80.196.3 with SMTP id v3mr12466997edf.232.1516088752984; Mon, 15 Jan 2018 23:45:52 -0800 (PST) ARC-Seal: i=1; a=rsa-sha256; t=1516088752; cv=none; d=google.com; s=arc-20160816; b=MOWQk2GuITztNVOfh67NppoQ4kCFBM01MUgmvX2z206zyXLW30AZ3A3GNlPTxoM/A/ xohZ2H5trR2BJxK2oMXqlt6ca1uk788VPdScU5ik7RWUx34qu8z/Q83yT5pfbRvYVnQK +gABYx6cq5YyryDwkAZl0YlnMl30ctG94PJ/EsBbDBIMrPauflOUEbknKKgv+/7A12Nn eJ+p0Q1m8Dnqq+aT7KI7Vupg0g6DIxh9REKAQcr/ajjfBerIwPY6DUriXEG63yTNJqzK pmUpa13Y8jjDpIYKvlCsZFmbpS9WuX/gp+go/NX8L/QhZ16yi1UuBWKLsJNDnCT7lJNZ SLcg== ARC-Message-Signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=content-language:content-transfer-encoding:mime-version:user-agent :date:message-id:subject:from:to:dkim-signature :arc-authentication-results; bh=ttU0QeIUuT+NdOXxg4fuHzX2ve63Ph0zFEJ3cwSc/fI=; b=vSIap0Een0Zmc50vanBVGwswQF2UY4PDJAYqKZVBXQpsgQXKBD1WhGzissdqwvaUcn iP4Kq3qoAIJILnZD6nEDOynyIOdd8ilG3ALrmw8b2b9ZDWSV72FHq2Fp0JXp6ZoGS63c zk+A/Qcnkw0FOYnbZQYv1Dt4ST+aKujVOWVUBpQIlwH7AN9DR+X54ny0ckvy2M3nHkui UxC2R9d1DWtHJkb+srhQabyAiB5c1CXEQ+V9Q/QgEjwyY/zubo1LVOysOtzRGtkEh7Cc HPbrwhWI04ZQw6bP3aiZHT+8nFsRVb52VD9JI68z3I0Ly3eRC8q3gy/e9w7jtwqUnXu5 JIdQ== ARC-Authentication-Results: i=1; mx.google.com; dkim=pass header.i=@gmail.com header.s=20161025 header.b=r9SE+rFn; spf=pass (google.com: domain of mgocean@xxxxxxxxx designates 209.85.220.41 as permitted sender) smtp.mailfrom=mgocean@xxxxxxxxx; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=gmail.com Received: from mail-sor-f41.google.com (mail-sor-f41.google.com. [209.85.220.41]) by mx.google.com with SMTPS id t23sor92946edb.49.2018.01.15.23.45.52 for <ferret_users@xxxxxxxx> (Google Transport Security); Mon, 15 Jan 2018 23:45:52 -0800 (PST) Received-SPF: pass (google.com: domain of mgocean@xxxxxxxxx designates 209.85.220.41 as permitted sender) client-ip=209.85.220.41; X-Received: by 10.80.186.114 with SMTP id 47mr28890706eds.252.1516088752351; Mon, 15 Jan 2018 23:45:52 -0800 (PST) Received: from [95.183.158.170] ([95.183.158.170]) by smtp.gmail.com with ESMTPSA id y5sm939801edj.5.2018.01.15.23.45.50 for <ferret_users@xxxxxxxx> (version=TLS1_2 cipher=ECDHE-RSA-AES128-GCM-SHA256 bits=128/128); Mon, 15 Jan 2018 23:45:51 -0800 (PST) To: ferret_users@xxxxxxxx From: Murat Gunduz <mgocean@xxxxxxxxx> Subject: RESHAPE function - a bug? Date: Tue, 16 Jan 2018 10:45:49 +0300 User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:52.0) Gecko/20100101 Thunderbird/52.4.0 MIME-Version: 1.0 Content-Type: text/plain; charset=utf-8; format=flowed Content-Transfer-Encoding: 8bit Content-Language: en-US X-Gm-Spam: 0 X-Gm-Phishy: 0 X-ums-g: true X-Original-Sender: mgocean@xxxxxxxxx X-Original-Authentication-Results: mx.google.com; dkim=pass header.i=@gmail.com header.s=20161025 header.b=r9SE+rFn; arc=pass (i=2 dkim=pass dkdomain=gmail.com); spf=softfail (google.com: domain of transitioning mgocean@xxxxxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=mgocean@xxxxxxxxx; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=gmail.com Precedence: list Mailing-list: list ferret_users@xxxxxxxx; contact ferret_users+owners@xxxxxxxx List-ID: <ferret_users.noaa.gov> X-Google-Group-Id: 776276095559 List-Post: <https://groups.google.com/a/noaa.gov/group/ferret_users/post>, <mailto:ferret_users@xxxxxxxx> List-Help: <https://support.google.com/a/noaa.gov/bin/topic.py?topic=25838>, <mailto:ferret_users+help@xxxxxxxx> List-Archive: <https://groups.google.com/a/noaa.gov/group/ferret_users/> Dear Ferret users, I am using, FERRET v7.2 (optimized) version. I have a netcdf file in which a variable named "VOMECRTY" is exist. The x and y dimension of the variable are in "i,j" coordinate (i.e. 1,2,3 etc), i want to reshape the variable to "lon,lat" coordinate. I simply define new lon,lat dimension from available variables in the netcdf file and generate an output grid (i.e. out). let lon = nav_lon[j=1,d=1] let lat = nav_lat[i=1,d=1] define axis/from_data/name=xax/x lon define axis/from_data/name=yax/y lat let out = x[gx=xax]*0+y[gy=yax]*0+z[gz=VOMECRTY[d=1]]*0+t[gt=VOMECRTY[d=1]]*0 However, when I reshape the variable to new grid; let VOMECRTY_new = reshape(VOMECRTY[d=1],out) it gives wrong reshaped variable. It look like output variable (VOMECRTY_new) is on new lon,lat, time and depth grid, yes? show grid VOMECRTY_new    GRID (G005)  name      axis             # pts  start end                subset  XAX      X                  11 r  28.935 29.248             (axis from grid-changing func)  YAX      Y                  13 r  40.991 41.266             (axis from grid-changing func)  DEPTHV   DEPTH (m)          60 i- 1.002 2300.5             (axis from grid-changing func)  TIME_COUNTER TIME          4416 r  01-JUL-2010 00:30 31-DEC-2010 23:30  (axis from grid-changing func) but when i try to shade new variable, it gives below error; shade VOMECRTY_new[k=1,l=1]  **ERROR: illegal limits: Z axis limits don't match variable def'n: VOMECRTY_NEW          Are you giving explicit limits when SET REGION would suffice ? shade VOMECRTY_new[k=1,l=1] Command file, command group, or REPEAT execution aborted But, when I try following command without [k=1,l=1], it perfectly plot what i want; shade VOMECRTY_new It seems that ferret do not reshape the time and depth axis? I would be very happy, if you could help on this. Thank you very much, Murat |