RE: [ferret_users] Number of cells

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

Subject: RE: [ferret_users] Number of cells

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

Date: Tue, 19 Mar 2019 21:13:53 +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=ObsCuy9y; 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=ObsCuy9y; 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=ObsCuy9y; 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.228 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=ObsCuy9y; spf=pass (google.com: domain of joshab1995@xxxxxxxxxxx designates 2a01:111:f400:fe4c::808 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 :in-reply-to:references:date:thread-index:thread-topic:subject:cc:to :from:sender:dkim-signature; bh=FTQA/RfOxBoDClbxZkjrPe3xdDrqoQWRkx/Moe4xjuM=; b=dLqijMCWOp2lQqp66UiDiWPrPArR7B9RWTI1+Yf36W12200BfUwMMLwzWQrwDV2nMH Ly90yQlbkZBlfSgAJESFSGL3v6QTOk/WCGw5iXM5449vMvxj0hgbHdRUweAK9WNow/Cm n5DMvTIjVQ91xE+LADCkAQDzHFGBRhfTFcQrPCiVyT0Egwo8h+JBAn2j6SkxHhrc1noG xkd1y3KO1IsxDPvrvHP6O8tIze5dveXwTia+sfmusplzUIak3n+9vJu0cJ1I10V6QQug OIq7oUeG15Emi5l4INWIG9OpWRTYl5I8UUg5O2LC95HSLWvyPrztqPiqvyGXqTmHudQw 6/5Q==

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 :in-reply-to:references:date:thread-index:thread-topic:subject:cc:to :from:dkim-signature; bh=FTQA/RfOxBoDClbxZkjrPe3xdDrqoQWRkx/Moe4xjuM=; b=sMbYwglBDmrkE01Gt4sHr9TBxZmLKA3Q+Qqfaf053/Qf4ycJN1Rwx+15vFzhFyevkX KgOPkmJOK0XMeVTNa5Owp2kMgcSQwKTuKzBqam1sotgKudCE4jmcijUMvkMQ6+c18+e2 1TIhsWiIpQWnqC+rUUZdXecUraNnOc0pjPS1MVys2BP1ln2pA8SXyrskorepwhKpwwwC 3Y8wHg/GT/QwgLMK3b8ZteJKx+GSxpUGt/2Zr6KbxzY0CANmlKSm/VTJmOwqrdMjdWXg FV/SWZXTWEAFc0+qGlSBOQQIkx7xlsBZHW1cZLcmD8HfTxvppHwWsaqGZoFH0d55cXBt P2Qg==

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:cc:to :from:dkim-signature; bh=FTQA/RfOxBoDClbxZkjrPe3xdDrqoQWRkx/Moe4xjuM=; b=K5LLSrCxdmdlfac8byZ19m883JNAF+h+GAJw0Zgw9+tpJSAdS5QM6i9PoHLKIG8Oga IdWVRZ4XiCLxkL1AwUsbPdc4KViMT+sz7+4VrCCfYmQd2xO0ZAqLRte6qTLJ6hK0moQT 7mfLltqXKrJOIzJ4ly6JNmC8vH7Rn25YRvd4nAmxyI+dXmhAJ1WKKQsUePLrtPHAgkIu t6yWZGc8QXQlEwrraFWW/+glRjg7RpwBTG2EVBebdLkjhvA5w9bATJDRy9G0n6uhjPW2 6VUT8gyaLEaHNtGNAxooGlV13zXDdM7jPIgziPH5KagWfQsVDgH+61t8PHVeA4PQQX4g V1HA==

Arc-message-signature: i=1; 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:cc:to :from:dkim-signature; bh=FTQA/RfOxBoDClbxZkjrPe3xdDrqoQWRkx/Moe4xjuM=; b=aRvgNcETZXGydBc63Z2ClIoDXcoCLlCvbDYFR2M3NjQxAUHWztyYrFHvZ41/7GDs4y O6qVDdo32dSRGRjr6iSl9jEZ9Ug9A/giEUOZP4fM7r71n1tXHkhv49N6MiZMkGbwSOls YVSABTIMWgYMCgBC0zT/BSVz0jNmRFWxTO2vEL8YfC2BRiu9YQG8P6GBJoYMfmDOioEv wglM+NZY5tG6el+FaTM9U37707oF61b5zRMhviYz6Ht4zOoV9CPs/hky3BnXmp+76Jvy Ip8/bJn4P4XCtS2eWXfsc7HsBOzkAbTZX0LhFsarTeH/CjUXGZY1pd9XMJiX9TA3YfAw pvPg==

Arc-seal: i=4; a=rsa-sha256; t=1553030047; cv=pass; d=google.com; s=arc-20160816; b=gmGupYktEHBKJedPPhWXj6wizc6fjEnSWERxOLyB8J00I0CPxogx5vedKVKB1XY3Nm WvhbUd7g/qIheWsB68g00BVijdwhshLiLCvKCt2g7B1V5ne/w85UiITyvfpBx8uGzyZc yR497ueWaNYmSv5pf6WzYBC2+fgXcxU36WuerOrEfVRIrocVs3hcEl7d+45UP7LLau3P TVtIbKgrSOBGMe9IuTxMv9LSJZK8+/kRe9IqhSwFyApIUY4qWug1obE7hEX+jBU7J9tu qk9ZD/dyEe4wKaLxkZoLBow8vLmGfyL26+s4pf4WhcrBAcVdGjfPvuRacMrAGarLsaEz l+Iw==

Arc-seal: i=3; a=rsa-sha256; t=1553030046; cv=pass; d=google.com; s=arc-20160816; b=zYonsKZR/0O/bB6i7H5UeR+ELMMBrYzEs+1JwGLFbSOefRouKy+e4iYBzRBSbaYWAI sORgFP4srxAHOAYweQ80CEl4S0A1QmGF243GII34TyX+JqJ2PODuJKttBjFxE13WKWbJ 9xfCrXNnqh4NmAcCixtFcsHDiRUeaoqN+C08iFlDhzk/EToiIr9VXqcL5gFR6ferN0xV Rbwa46Me9D+QBZbkqEobtjKQJYWbjQCU0MmVbWAtws/ubrQFNMrSsorw1VxyKg1yieyG kC4K4n9JGMWE3tPh0VDkgmU3vffj1E3PtNQhl5jAdLuEudJMre1Ih+u0kP5h0TSU/yfY iWjQ==

Arc-seal: i=2; a=rsa-sha256; t=1553030044; cv=pass; d=google.com; s=arc-20160816; b=TcZ1U7CviK0LBR1cXwGiFldF8dBeUpGg2/2MMDazaRmgtJMA0Xv3GbeVy4Gf2Z5bpY 7x55Adfzh+LWCNZ4EAlqsInWpZ55GuHx3Bu65gWRk+WpFHqfZCVeGTy2cNYFw0XE3GVO ybXjbFaVKnnbogmLA2auOSU8DRDIDIdrisYJHwfnhzr5rgll4aWLFwTxQBeSChQ77gA4 l80SmYTM3L4LUAA7mxCrA+d9bVdL/JbISToQME1eqtCKSE/hy1tLxs1iUKqJHbDoZbxG NDAYeya6qbnyQCwrGWIy75NhKXer1fWG6IsA2xkQVNB9qW+Jx0Wb9XkxMPs2F5alsZla uMuA==

Arc-seal: i=1; a=rsa-sha256; t=1553030036; cv=none; d=google.com; s=arc-20160816; b=k44QbOQ73nWLcHmKHjEH3j7EpoukqLLsXEyE901zWI/Vgi9as/GvciNG8faHuHC8pv o7Gda3IRso8/T5ZMNsGIBVWiN4H9Nrdk6GHUGahqZtMfy+eRS+Dtm6tqRrJOsG+uHlJq syotYkEqaTngrp6685A5/OGuKYrlTbl7FWRC9zJG//v2QRlvfZBmNs4F5oJYg2tojbzA DSiYYHJW86IYWwsiMkilrcCYuD7mkLwqqO2o99B3Bg99NDyB+9g3e/Cgo8Ahp/Cw0ukb qPEJWiyEzhZ9DUIetIQn77aL1o+KWifUumc7GtusCRXo7LCmooGYrkCW9bLd1nZLrf4+ uAqw==

Cc: "ferret_users@xxxxxxxx" <ferret_users@xxxxxxxx>

Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa.gov; s=google; h=sender:from:to:cc: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=FTQA/RfOxBoDClbxZkjrPe3xdDrqoQWRkx/Moe4xjuM=; b=OzjYqpEnrdXEg5eLGnsLarMcNynm2Jz8JsR2JdriR4kueLlUbSYeU/YbqiBA2WnR0R ma7cCaJQhNFzSlG3o7gUmSwsRMFQ+RGydFGxMcWy5mHvP6u7k09aBDoA/M/cnmPq34Gc dlPPX7SQzSxLSElGPdGVhUrARTy8L0h8iKyoYRTNUWvZDrm/o5LXX41C+ZbuCWyMuS8E BeD8EZXK3QJN1aI2FlLVa1p+k8A0s37Fe4AEkaAuxqdJFlpywQnIEcWzlt3zTdhBqj8W UCU9jgzpqtr0OXNhcFzJLR7ttCr2wxNiOOZrlRnBjbn/h8F0kXa2+sfKeMfQcYHzSwwe gjGQ==

In-reply-to: <CAELLC7XSDMFMBzhOSSK3RSyDErzD0d61GeN3k8JuTGUh=1Xdhg@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: <DM6PR03MB43484C9338F48F0B9AADCB3DC5400@DM6PR03MB4348.namprd03.prod.outlook.com>,<CAELLC7XSDMFMBzhOSSK3RSyDErzD0d61GeN3k8JuTGUh=1Xdhg@mail.gmail.com>

Sender: owner-ferret_users@xxxxxxxx

Thread-index: AQHU3n9UeHnFIaOQ50OUeaWLkPNsIKYTX5WAgAAUX7c=

Thread-topic: [ferret_users] Number of cells

Hi Karl,

Thanks for this shortcut, hadn’t realised this was the case.

Unfortunately I’m still getting an unusual trend for my data – leaning towards the idea of the issue being in the data.

Cheers,

Josh

Sent from Mail for Windows 10

From: Karl Smith - NOAA Affiliate <karl.smith@xxxxxxxx>
Sent: Tuesday, March 19, 2019 7:59:41 PM
To: JoshAB1995@xxxxxxxxxxx
Cc: ferret_users@xxxxxxxx
Subject: Re: [ferret_users] Number of cells

Leave off the else in the first let. The "implicit else value" is "undefined" (you have not specified what the value should be, so makes sense it should be "undefined").

LET mymask = IF myval LT mymax THEN 1

Karl

On Tue, Mar 19, 2019 at 11:28 AM JoshAB1995@xxxxxxxxxxx <JoshAB1995@xxxxxxxxxxx> wrote:

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,

J

Sent from Mail for Windows 10

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."

Follow-Ups:

Re: [ferret_users] Number of cells
- From: Karl Smith - NOAA Affiliate