[ferret_users] Number of cells

To: "ferret_users@xxxxxxxx" <ferret_users@xxxxxxxx>

Subject: [ferret_users] Number of cells

From: "JoshAB1995@xxxxxxxxxxx" <JoshAB1995@xxxxxxxxxxx>

Date: Tue, 19 Mar 2019 18:26:39 +0000

Accept-language: en-GB, en-US

Arc-authentication-results: i=4; mx.google.com; dkim=pass header.i=@outlook.com header.s=selector1 header.b=kDE4hc2C; arc=pass (i=2 dkim=pass dkdomain=outlook.com); spf=softfail (google.com: domain of transitioning joshab1995@xxxxxxxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=joshab1995@xxxxxxxxxxx; dmarc=pass (p=NONE sp=QUARANTINE dis=NONE) header.from=outlook.com

Arc-authentication-results: i=3; mx.google.com; dkim=pass header.i=@outlook.com header.s=selector1 header.b=kDE4hc2C; arc=pass (i=2 dkim=pass dkdomain=outlook.com); spf=softfail (google.com: domain of transitioning joshab1995@xxxxxxxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=joshab1995@xxxxxxxxxxx; dmarc=pass (p=NONE sp=QUARANTINE dis=NONE) header.from=outlook.com

Arc-authentication-results: i=2; mx.google.com; dkim=pass header.i=@outlook.com header.s=selector1 header.b=kDE4hc2C; arc=pass (i=1 spf=pass spfdomain=outlook.com dkim=pass dkdomain=outlook.com dmarc=pass fromdomain=outlook.com); spf=softfail (google.com: domain of transitioning joshab1995@xxxxxxxxxxx does not designate 65.127.216.226 as permitted sender) smtp.mailfrom=JoshAB1995@xxxxxxxxxxx

Arc-authentication-results: i=1; mx.google.com; dkim=pass header.i=@outlook.com header.s=selector1 header.b=kDE4hc2C; spf=pass (google.com: domain of joshab1995@xxxxxxxxxxx designates 40.92.14.96 as permitted sender) smtp.mailfrom=JoshAB1995@xxxxxxxxxxx; dmarc=pass (p=NONE sp=QUARANTINE dis=NONE) header.from=outlook.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 :message-id:mime-version:content-language:accept-language:date :thread-index:thread-topic:subject:to:from:sender:dkim-signature; bh=dZy227b35i3ctNzvr2zbL6vUKSlPjDs8xvnp4UAdiR4=; b=RQeJTDCsBuLD0TB+jYjAd0SPkEq3vIRhgQZ1i00KUL3q/ugu4SUX8Lp0DTvyYGLHEA roKk+3WxQWADvKwp7PClkR4PMNgxNo4LGQinA7pjDKfh5krQXfu952VfyDh4EuSyPnHF cpAwQ++RfyIm1cz2fes+WCEBEXwsPEr59J7iHdX4x43raVgyY/Kaj4r42JtO/Rk6K0R7 3GdONK9UyGyVxxCzB2EN+n13yr6MgaXfP6F3SI5hTdUNwG+XmMpfGG84UP9C2WaEPIFs XfVLoxGeuhbx2q+zg+dJK5Rgf4Spfn0TkJWFC+aWGfdQxqsDoD+CyWg7xXZ9F/jToHHK f7tQ==

Arc-message-signature: i=3; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=message-id:mime-version:content-language:accept-language:date :thread-index:thread-topic:subject:to:from:dkim-signature; bh=dZy227b35i3ctNzvr2zbL6vUKSlPjDs8xvnp4UAdiR4=; b=aj92F92oGdLzB00pelD79U0YvkvwGFlfR+RrcgZVbDhpCS1t1S1I6BoNM1jjq4GPfE lRs2GOFPUq/cKIXdsUekz7jb4pOdj3QlqOwj8n4NZbIg42zMNsY3TvenB9kshJeJQGRg BieMh3gBvFG+1aLlcWUAnorOScF17xJgtVYldn7UOxxLEY1H/oKeqUb5cnNhhJWHHxn+ Co8hfl/WMLRLiWyI6RW6XhJ7e1clYBl86gFUKuBJ+tBzk+TdZTcBMZZox/xswC/9Ut6q Qr7hoO9VZJ1bqedlDI59QINqhqbGGItBCx2zneCmYUkPg1gqp7uMG2SEjlU/ox5ORSB2 4sEg==

Arc-message-signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=mime-version:content-language:accept-language:message-id:date :thread-index:thread-topic:subject:to:from:dkim-signature; bh=dZy227b35i3ctNzvr2zbL6vUKSlPjDs8xvnp4UAdiR4=; b=MgiIStbqa2G9qz43i6yUbnEn0DOxdHTMdMhXYTZoPSHh/ccufatnktYlju7kS3T6Yi nXICRe+DBA3QCvkWvvbuEeQcnV80GCSQthb4bG03O2VcyrEHiMeuCi5nic2C4BTajc+r Zjierli6C43+cAL9zCXCa+hQKETETaW1V4YQoRfV3qHMB26uTdlhHdaS/Wr6ia1nxoQc 1drID0FFkRK7ArJYfIv+yfbw+adZ9DGWPuuU2RECLfl63N86YDpZtyyTyW3THAdVCV34 tTdmoiSkQUNaSxrOW5Fkzw+l3am7QC95jTWD+n9rsk1OcBNvJEKkdtpupfzhzMqB3vgI 8f7Q==

Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=mime-version:content-language:accept-language:message-id:date :thread-index:thread-topic:subject:to:from:dkim-signature; bh=dZy227b35i3ctNzvr2zbL6vUKSlPjDs8xvnp4UAdiR4=; b=fHDvcoJAvWOKGEsv/k7jxX35NYkW1qlzsbBieYUWggpxkVv0IvXXVXVIwtELF6bJeQ EvAIxsYULrZWNI74azbEl0Rd7PlWTsGG7QjjJOeILoG7I+njXpBkj0aXOTDnJUPD1Hvr LFAYCdLw+qtTI6Nb84MqNdroaz1OKxsharkxu1SCwcpBG9mGJ+hYMOab7zsEPc31c27u rbAZBn/2L1PRFKwIGKxDfhzF5NTJEQBVEij+y9Y9Zi8A9g8iHn1noYQYTTfqg4kUaA0I wUMYEaR5K9Ggoy9xkMSf09NVbDxlWwEN1e7zaaoX/uBKHibrJIAzDyHKhoDcNvrpQKhH V+WQ==

Arc-seal: i=4; a=rsa-sha256; t=1553020011; cv=pass; d=google.com; s=arc-20160816; b=PrA/5aUwrHEFe0quQrkf0Mu6vnICXEiQJoaZKkLwQd9QM5yZGWOqqZ6ajHi2w/LJ8N v/m1RSZPrD/tSktKvabVXWWIDQguLJqfrm0abRyPN1bSplSJaC5qiE7SNPcezp2Nawip yQxBQT7LdJltHrLQu1F6ZgjXRuvBGVhaGJZuDxRv1Yzc8teof6U2jLstDT0ShtjTSbWZ Qrm16LsO6y5utb6zno0zlyZPAq82gPuJifB6cBwnxk/CWO2n8lhxI/FqtIWxb+cslvku nHk/JnXg8pZIEgttPqR9yF4pPKN8/9OcePAPvL4Fr8kEK+2yQjQLRNqEBlvza7hL8w0g 6yzA==

Arc-seal: i=3; a=rsa-sha256; t=1553020010; cv=pass; d=google.com; s=arc-20160816; b=SzeNhgdvSFf6vHd0zL4Moro4KWCVqwxSmSRBr9ZCkTMBFXsZj1/Rz8qbNjSx2Fr3Kf APADRY8hiKeLkHwfXWUuy1yal3VHAcVQoVI0cVcCjeFQQywl7KZeYN+FbsCZyDRvuLD8 LWKoMdNPCWH8DBQGdKOAW4xHYGAGAH2I05LE5JH0lg+CZVXgrg8+XQJKlsNZW62PT6YT nskYsMieDsyhOCuP9Z2oRg0yBoQkyPGAXuVzN9oMqFvIRNFuto+qz5Z7tKTDM5+w/oE0 viNSJjwUb1qukpfGUyjiKYwQEiQs8q+8dKm1R12IGR4su4ADlkF1OZErUpDtRmP3VnsP 3aiw==

Arc-seal: i=2; a=rsa-sha256; t=1553020007; cv=pass; d=google.com; s=arc-20160816; b=Vf9497phli94A4jDisJ0bs3SYybhxv5VPPXBgmIAjvCfXRD5VROJjx3ZbIumJqIL+d tdRk3aTlaFRojvYT9AVTxqjobNd8TpuVqx4ehYrutu3nQkd+erpWTDQGt1QwWCvcNGeO 6yYRmMgbK3GG3REXI2pg3Dbba217Yd3lG5lAviu/1sv6kMkG6Y+Xu3e5QlcY2+9b8LMq Q560gIDNLbfiGbH8vtWZ+hU2NRxwGA8Rj7TP0bwHfeHvftouDHMOaokjPVXRU8K45t/J 8qLm0ovth/St2tSaNCgCpX/DDOtoUmJbrxHK+oK+osAWcC6xpwUF5HDddttonnduGpWD nauA==

Arc-seal: i=1; a=rsa-sha256; t=1553020000; cv=none; d=google.com; s=arc-20160816; b=x0nuYaza5+jrfL8GcvGEvUgxQCX/rl2+XEaJ0mInN0zraNlGRwOaq27fUkOFjxl6TA Xk55XTkljI1RV+5T2KpikerV4G8LUvz6HScBNRtQUTOa5FPD/dhw1k15g7KH6GSHV0sH PcHpe7aUJ0R1PrKYUHCKJF3xQNlCKBeuEPMMhL6AXQhUrlmevmrL0MxfeAiYsUKB4O+T yeYxo3GMJKbX94fiXhkm9/hgCIY1N9qb/H6eLuNq/mAOSaLxEYbDjQtBE6lW8XTaMVUq YHPjg2xbEu0uR9TclSwtQBqtDi9agCjCZOnBCNLXKj118ShQI5Xnqd9NGgjrHGvj1UDX C6aw==

Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa.gov; s=google; h=sender:from:to:subject:thread-topic:thread-index:date :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=dZy227b35i3ctNzvr2zbL6vUKSlPjDs8xvnp4UAdiR4=; b=FED47pHXLprniMnobk23hY8wMjM16dHeg0RzIM/Ksu9AsgOkWiw7v0mDPuDq3P2R5p mn7k0+odOcNOEKfJVpIWHefR4hlfqJG7n/4G+aJ+EpLlKaXT0txvu8pboW28OGoBBsDr Df93yjYpPcJQjQ5bN6vbIqcwE0yHC3EthXU5vPgTdMFIRNNn2k3BQgkxTZxKgWrAayfY cis+o5phKDUjS1zQUOzJPogO48pTawEXSZWoUvVtZ2+jPmFY3wguLYyz/jLnnGO7ABgk LdxPbxbOjirOPipU2ZOdftvns/pGy5Nk9uWh/2vv1IFAS2TU1ZYp316tFmKEw8gNGKTD mceQ==

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

Sender: owner-ferret_users@xxxxxxxx

Thread-index: AQHU3n9UeHnFIaOQ50OUeaWLkPNsIA==

Thread-topic: Number of cells

Dear Ferreters,

I’ve been stuck for a considerable length of time trying to calculate the number of grid cells that contain a desired value, which I could then multiply by some value to gain a volume.

For the values I’m after I’ve created the following mask:

Let X = if x LT x then 1 else 0 !This sets all of the values that are within the range I’m after as ‘1’

Set var/bad = 0 X !To remove the values from my dataset that fall outside of the desired range

Plot X[x=@ngd,y=@ngd,k=@ngd] !Hoping that this plots the total number of grid cells that are not ‘missing’ against time

From this plot I get an unexpected result (a decrease when it should be increasing!!!)... Have been beginning to think that there is an issue with the data.

So I was wondering if there was anything obvious that anyone can see I’m doing incorrectly in the masking? I’m not sure if my ‘plot’ is adding the values correctly – the best I could think to do was for in each axis (lat,lon,depth) to select the ‘good values’ or to ‘sum’ my values...

Any help is greatly appreciated,

Sent from Mail for Windows 10