[ferret_users] RMSE (Spatial... I think?)

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

Subject: [ferret_users] RMSE (Spatial... I think?)

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

Date: Sun, 21 Apr 2019 21:14:50 +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=kMQnEWII; 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=kMQnEWII; 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=kMQnEWII; 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=kMQnEWII; spf=pass (google.com: domain of joshab1995@xxxxxxxxxxx designates 2a01:111:f400:fe45::822 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=ZbYVu9Wv9gqDbcDHOfEOs0GqM7m5NfMPFkCbFySxQ6k=; b=swyJQBGa1GRpH7N7I4jJIwrmBoOMjaMIRe0D79IZutA518u5coGijWSECjPr+QFPCs XE7+EnA4d9e9tmDoICC7ZzC303TEDMWnBJx+VPEAtlKGESu81JGcW09efZrGm8TGacA9 k9I54TNBihF+b2u+TT+E0ccbn6dxBr+FpZ6tsrfJzORj0ShspqlgBosk+Py9QM1GNKG+ Ya+FadcEgBMQdm8z28a6zGPcI+ciA63cbt2T+HASMMYXpJaO2J5tpOrAlCFS0Pq6sxZp P5Q+fT0hvsT9EWzosYH5tgI6jFiZLaFM7KE380jdQdL5HANfRSu/9eP2feN/P8BZuEcd MkjA==

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=ZbYVu9Wv9gqDbcDHOfEOs0GqM7m5NfMPFkCbFySxQ6k=; b=WioxMdsOZ+F9DNXVWWpa15q36fYYnLNyxvmYvSwIPDfdCQy5rPTSozdIfl3k5TKNJK 1ZVKXaYTQE2yqljfhRSd98mUB3+0SlCWOkRJwR/18sh5b0GFOSyqUBb/ZuCv4G0nNsvc +JaMyQAySs09C6+Zkk9m7XucAcrzXw+P+8zHiDf/pmkeVJ2aohgz3HHIVNjd4j6tPMQd aaKk3ziWNjCTw6YZS/jpd/WIBpZfvR6XwMA7C3sI7axORLWg4ar5sGoHu/+vJqznn70w 75WnAAZG7AlcpU8kaagA+kyoRCpfOklBziQkE6eH0T1xLRrGhBY8mtJ6MiftFDNtrKbr Wryw==

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=ZbYVu9Wv9gqDbcDHOfEOs0GqM7m5NfMPFkCbFySxQ6k=; b=nluIcKxNgSQK9pkY7ZnH7YvnmQKCGudqWi8KZt/8m/nSj6dqwmAFCGENvL9idKnUD1 BkBcOghyZf4hrb/1bDhz7PkhXUj2t1AqandZthFwLfvZVX+ZLBs448ByzKsnV2RGpXjF CxuWMFUbFb28Pd3eCZU2cPrCP/kq3qWPFJTE5JkVWtN/hADpcQTYQJAHEK8VM47CNRIF kgFP2o8TSBqxOHTdpG3Ti5/Wt8YF/aELzJxAsFqKCZ62Mrxil/N8DeIOvTw4uFdcpf7N 8wK9xb5oGdndu5WnNkai8cKuOByRV9BkllosC2dcf68Dhu9yRb3ZlxMp5w3dlzm16w+6 orfw==

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=ZbYVu9Wv9gqDbcDHOfEOs0GqM7m5NfMPFkCbFySxQ6k=; b=dGxwN5h6NgcvN/jyVeqMn2l0jG139vn9Epr+MNMFZpuwEwRd9s3pMJaYyj9c8P5N3k H4NivHnNbj1MQwcdmKjT9yJuNy6cvIuxuc4y+IcScjntykgJqpJhI7ExOSNRuQrfp3U4 ZfgftxNV1Q6oFEoqM61ur6ojXJAtWM+iZoks+lClXpTqWO04NAiM4t10490yMx1RW3i5 2OLGzLCvYvYrPm3GQSE+eNuyuPDTm7Le/+ct4adnaMkA4LVUq4zOiifikf/uWqX9y5H2 4oxQysSaC1LYbszby403Ms6Qg0tVQAm4n/Rb0jXzZCWozsDG5SjkKStHApuN7XHGg1mW uc9w==

Arc-seal: i=4; a=rsa-sha256; t=1555881301; cv=pass; d=google.com; s=arc-20160816; b=kyf5yB5OkA26k3RVPEJg6FbXrLZtOfh650A5PfwxxeVnjmDjvBGWi7UyqGYIhnn3i8 NiCeVOArmD/ffljD+k7JswKLtSVn6XrHcQ/tdkkAbUxPUB0YKSg0LAQjhsh75nWvj7dn F7ddvBtK0kfHPUt82K9Ao+a5/F0IQpWebh2lM7uMHcZT7WDqxOJ4LgD6Qz5TtjEnzgP6 g91z9VupN+yBEhLcM0GJr0QdQTmKIj7eY8OYQTJn9o5TofVDHn2my9yarxLUwAQKhcwT iZft6Yhe1Sz8EfpQglqi++9BvfvSOGgXlmQN1CCzGkARNak371XGDcyOAUD/f4Tn14ZG 9WOg==

Arc-seal: i=3; a=rsa-sha256; t=1555881300; cv=pass; d=google.com; s=arc-20160816; b=gGKcBbRL2axq058UeJXLP3vfRfuD4BQUKFa1UEdMcTgJ3IbBJR4rRkBcBltfZfOPfq mZiAVOsDVQ30Ta8+cA2AAmh4Tr/sjY2bIArBYufvjVxmvs/CzrRmJlsANReyvrULMnYo v602e5s5cblNLwxfgHeprC+okUKm4OJWa9GPpJJcKhuVW2Akv55ytu6YvbvJRUSu/4sc Kw0yJCcpqM0oILbEMuiYZ3lVfUBacZGF2+mm8N0PC2OvlIe5G7ykNH/qEqaPRc7SaiBH 2zO4i+BkycVzuMQYQJdemusWrz7F2zJ7VXhr+gcdvFbbAp+TCZfB/C7FIKQAFWBrhfpW fIKQ==

Arc-seal: i=2; a=rsa-sha256; t=1555881298; cv=pass; d=google.com; s=arc-20160816; b=TtawQJ4BvVY3SU6Sa9nVQdG22ECGijQzXGw/evQqFB+WM+mhu1+mOOBVRU/dom5q8A rA190LTwGjyTLCul+N33URFrC77Rr94nIYFmRVGkm7CZgL+yZNVxoKFzpWK8zcuM2zRF CUEU+0Pwf/H6Albmw9gBoASGfAAOvukcopcs7jIDDOTw1dpkcn5Cgt8KJVQ8FI/NGQRV kSrL1z9ML/YIsEVvdt7PRP8o2PUBjvMk4Kq39d53rcA8wVHCdJCEahWHqgdRiC8SIZY6 6lyZstktL+e7GWNOYYuVfcmE9OGXsDKnVPUbCUKvSbLg4g5+1wsFOeGYHptQ0eiIukys yomg==

Arc-seal: i=1; a=rsa-sha256; t=1555881291; cv=none; d=google.com; s=arc-20160816; b=TW79wO+RLfIV5RV+N3R4/i7lPfbn2coMYKv9md+MdW6uykJ/pUFc45sKVx0o9lC88y lx/AEUMS+rCJHgs3vEedhoE+Yl2llvL2xX/dScp0irGOtqmGgtQC7Pfo80V8FpEjNb8g j6n3NhS1Q6yMuUcXxRJIPqsr59LPb4KFW6PMc165kDkClvz+QFQEsoNVr6WOECNjGwZz s8JrZYOazaDCnhboLrjfMAIvEYaAcezR5P6shIs0KPX/3DDuu96c3HLPRWQEUiGzSaL+ C6vNqgU34YdH+fZSCDAO0K+UtwI7YFoR5GYMNHcbYILCGt/5pJuT5SuoWie0lORZl8WM 4DeA==

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=ZbYVu9Wv9gqDbcDHOfEOs0GqM7m5NfMPFkCbFySxQ6k=; b=f2OhDRhrPwgiigYEC7eqmbJ/x9g4XYSZfPh1o/WNxiY8v2F45E8NWk2vWLDBNrPcKV Qiw/qz7TewntTB9OaZi1pnKpZ6CsOHyZjfh0m9oXDkQxBt3MoQ+KYqquI1BPCVFS5wNn CATTwgNySJKOwYHD5ju3C6sPcw35HuScEDaOjUXtYSGANP7DVBrDhpdkvpLNBFCtL1+k rjKsWKOq01QY7N2NBr6RnSvTqN6eljuq4UpPqAAvaa1rSVgValmcQHZeXrzXEacccHaZ x0oq6mac5L+WzIjtLAvGYDm9NaeRx3683Ia80r3jatQA9mofFC1lVUBiKCxf+M2GY2Im aWpA==

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: AQHU+IdBCSNkIruPCUmkM/aww3676w==

Thread-topic: RMSE (Spatial... I think?)

Hi guys,

I’ve been trawling through the mail archives for a solution to this and have been following what others have done but still not sure about my calculation of RMSE.

I want to explore the RMSE at different depths between modelled and observed temperature. Thus far I’ve come up with the following to generate my RMSE values:

yes? use "./obs_temp.nc"

yes? use "./model_temp.nc"

yes? let diff = ocn_temp - regrid[d=1]

yes? let diffvar = diff[x=@var,y=@var]

yes? let diffstd = diffvar^0.5

yes? list diffstd !this shows the values at different depths.

The datasets both have only 1 year and 1 timestep, hence why L has been left out above (not great I know!)..

Is this the correct method for calculating RMSE – values seem a bit funny?

Cheers,

Josh

Sent from Mail for Windows 10