Re: [ferret_users] fill the bottom value of z-coordinate data

To: Changshui Xia <xiacs@xxxxxxxxxx>

Subject: Re: [ferret_users] fill the bottom value of z-coordinate data

From: "Jian (Tony) Ma" <tonyj.ma@xxxxxxxxx>

Date: Mon, 16 Apr 2018 12:45:22 +0800

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

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

Arc-authentication-results: i=2; mx.google.com; dkim=pass header.i=@gmail.com header.s=20161025 header.b=YtkvhkDt; 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 tonyj.ma@xxxxxxxxx does not designate 65.127.216.226 as permitted sender) smtp.mailfrom=tonyj.ma@xxxxxxxxx

Arc-authentication-results: i=1; mx.google.com; dkim=pass header.i=@gmail.com header.s=20161025 header.b=YtkvhkDt; spf=pass (google.com: domain of tonyj.ma@xxxxxxxxx designates 209.85.220.41 as permitted sender) smtp.mailfrom=tonyj.ma@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 :message-id:to:references:content-transfer-encoding:cc:date :in-reply-to:subject:mime-version:from:arc-authentication-results :arc-message-signature:arc-authentication-results :arc-message-signature:arc-authentication-results :arc-message-signature:sender:dkim-signature :arc-authentication-results; bh=fDW5MFugPH9zzkPIH6WcNNqNZI+tVUMChJcbVG0W+IU=; b=haTwMh18q/P5945gjC39e9+Ae78EmyiBS1OmSMdJgTCEwOQDPI5tMgmxYTDvwHcnSd 8chIPG+nGlrACcCW+IHjszvyd6PV7RTBtqPNHsV3aajqT1df381rGI+6V+aqXhinTUZ4 5rK1HeaxwhzxQTTIyr3XlrKyrTjAGZyavMjO2wtuMfMEa7n05mc6sEXGNBOygIOep+iW 2BaxErPvHRTzCkbgb6Lw/jU4madGuo7LTCLO7CE6zyoTQsvoVMuVmIpbpSh5khGbaIA8 Vsq8W7ANARiqJcetfM9C9IJIcXmCIzkq6qrtDw6sAPbcQz6GdXfSfGIUTPT/gAXICgjh T7lg==

Arc-message-signature: i=3; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=message-id:to:references:content-transfer-encoding:cc:date :in-reply-to:subject:mime-version:from:dkim-signature :arc-authentication-results:arc-message-signature :arc-authentication-results:arc-message-signature :arc-authentication-results; bh=fDW5MFugPH9zzkPIH6WcNNqNZI+tVUMChJcbVG0W+IU=; b=UZmq5YGePqOxJB6UTsQrDbew6rNvO2dgKMEN7PHfJLQHWyiYELGl3UqD41QKbv7kHp yuVvrRST3Mf6b35w813CnasZ0307KzQCQkcpBJH7b45o+J+CS8HbRJjFL1M67lEkHUj0 sHbssW7HzdJ6YAN9HnV/QuM2JrHrIspfjrBB8aQFI7DTgQDoChaEQpw177Z09+BrR4U2 lCcNcFEmcb97LYM0tog/6T2Is0NkATIJiHUJwHYrrupTk/PJG8ZBtnJUYzh1St5jL6hs folYex7dfRtCq7oXmALYP0E9iLnEPUzhRdLQ0pGhRdAxcU1h/nGr3o7NRkfr1rq3HyjT 6E6A==

Arc-message-signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=to:references:message-id:content-transfer-encoding:cc:date :in-reply-to:subject:mime-version:from:dkim-signature :arc-authentication-results:arc-message-signature :arc-authentication-results; bh=fDW5MFugPH9zzkPIH6WcNNqNZI+tVUMChJcbVG0W+IU=; b=gF0o1j7bMPuHnX8nCCCo88UKVcl4lOEHDUDDFhpRd/1VU8c3AnE+muhZNDR7tP3lk1 QcAWxHHfzq0YHZB8N19WnaV93aP1gPB/7QNxwOUcuYia9BIQ21SzvT/heC/bCBlK4eg/ r4biUiFGTP78tEPUFdJo5IERsg/HpzkO1RhZFgg83E85u6U+q3f4NXexZU3bsjqc5PsY +NVntoFIVuzzyP0WY/qoWG4x336gSZVL4aY6UDXvzWHpuWyo8/sH0Cnaqyt55aAzIIfD LK6Zr/WJjHrX8pKTnXqRFVI+VbkpFcitp5qFYTbAJOiT5a2ipagUN5JnoMUQ3zJqlSWL IIcg==

Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=to:references:message-id:content-transfer-encoding:cc:date :in-reply-to:subject:mime-version:from:dkim-signature :arc-authentication-results; bh=fDW5MFugPH9zzkPIH6WcNNqNZI+tVUMChJcbVG0W+IU=; b=J0+pj0EmPuCMan03Ds/58Upkkxxd4riM6HBncejqbeIw5t2LJX1sSYXbflUWtjVH19 1YSTa/9NDqRtDXJwU5VsQpwIEAmUZNp/XQ1NDYZNxlzawh6vlD213Ciad/ZW8liLYNL7 ir9pC+3/vRTg1hEfDnzP9vbeJ7HeIb7az3Y3RBhEwAGb1NYpXgh51rajKfnQsHpypQjR udZiYhThzcMyx5MA2s6LsWWMe+uMa8y6PEd6nET2xBtziLGPovDl5yaMV8CAClyLe7Bo YDk8FQtHmKqarI2sRatF0yP53tDHuCAqk1ecJgT97bN/p98XcTrhBLDOgf8FpIFE08Mw DzOQ==

Arc-seal: i=4; a=rsa-sha256; t=1523853958; cv=pass; d=google.com; s=arc-20160816; b=Iu+5ZQqnc9HfC5n6Suywdm9U2A68qXxf7JjiE+rlU7Ap2Rqo4/7YJwKUOoE2OJnDgU aEci7Jvl0byktYsuuFW2/upjoDJ21gT6GnMc6S18RtyAscbjik8IUbahSO0RfBYzutVd kdvffrkZCcVzvguGlaZGCYa6DLGrmeok3VntZj5PWMcSZeBsrxye0Zvex/DO2++bKCAu V4NNxhUL3b+QV29gK4G9Ox1v2Mi7RxGZ8VEJ1ZDTp6B3LrG1K9PaKWXi9t6AQAjSUHYu n5wilkTZqg1yTOtfOdPV3AAqbSwAzYqSmFv2gNLcPVkrkd9uXV8jpk3jjaf/0DkxYiSA Mdhw==

Arc-seal: i=3; a=rsa-sha256; t=1523853956; cv=pass; d=google.com; s=arc-20160816; b=UokEbA77VQVwkQ/i9rKIVX6/Qka02GGOSgCDa3mfRgHlbBt2O1uaYGHHLcUCBX963F BEJ6oaddZBfjvKGTL10VEaw33cC9I+uBN9hs4otLW5QcCzUlEocfRGE6Xr85X+ysZvvR qX2O4GzhOycjlmIshvLngHL+1t9Lt6kEnty2aIXVupB9TRkGMXdMGb6UGQgjuniTTRIV P6uqxXvFkTywDLikFS5RBmX56IETxdYZGEg9BHoLUScpUWUnL1Rk27k2xbAGlBWFJTsQ U8NnbeSpDLpSmpMwH4TlpsbGz/rMCScOBRizJwVPQZiMQWYB9YwpxfrwoNPyisd6sDDh O2GQ==

Arc-seal: i=2; a=rsa-sha256; t=1523853941; cv=pass; d=google.com; s=arc-20160816; b=VM1p07ZLVWVEfSIyskMpRTOmdbJM0K/pvz+8MZYDZNpQbrRmCAzWzgnSIyHClTlP9D RGXB7xETF4U+O/ZVB/o6bZldk5mhKFA2maYiE38rDuvRovKZ2n9IsLUIBfzrKSF3uF3O wYMOzeOqF3siHuM1Bqyayo43tEp6MDCW2gkR3Pnirqb1Ea43fLe0q5vChLec0fw2O8Zt 8UAv0+rTBqaIO+nbwzvt6UOosNBnOoGSTM6qhohakFFOHCn+I9GxN6nh7f9E8F6Dq1iZ 0LwtjZIfUKO6Tt7jK+HxqFFZg0JC414MGDBNu6tAPTuh9rzuMY6wIFQaIpD5DoiIj0pW 3iqw==

Arc-seal: i=1; a=rsa-sha256; t=1523853932; cv=none; d=google.com; s=arc-20160816; b=B/G/vfV/P73+qyu+JYHf4VqKQ9XbVZKHCp5BFRHrkFPTDQUlmkSzAGU6XpFg9EUVFG NpXijqiNBGdeVyay19qpJQk+XJ3fMqRqkflRcEzLspd3NRwyekjSkXd68eFMFIHgmsOM QpazepJ29hmeo68Hbkymq9Np3XEJeZF65qxcsEx3x5lnQEf5Mns6rl+9UXRGkwLqsqzy 4IY4CkEFaAUiW6PyqxNIq3N9zGeQ0UdAv/K24wGEEO+y7kIxsh/m77pNnuwoT/HZb7ZN wrL8c6capMWpWyiZA3kaW3L2v0l/UmiSGnvUeb5eUwxKRqoPXifrTFgbcWt8SUhnD0Gm xXSw==

Cc: ferret_users@xxxxxxxx, 舒启 <shuqi@xxxxxxxxxx>

Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa.gov; s=google; h=sender:from:mime-version:subject:in-reply-to:date:cc :content-transfer-encoding:references:to:message-id :x-original-sender:x-original-authentication-results:precedence :mailing-list:list-id:list-post:list-help:list-archive; bh=fDW5MFugPH9zzkPIH6WcNNqNZI+tVUMChJcbVG0W+IU=; b=JsKgcytkoTN3yvn4ZKtyvO6zh3PWrlg8EKAclU62MicpWUNWdTFSrhvTP8E1gv+Dz6 t295uW3+5aHYfMclqwvpxEJ/Oh+LxFTuufo1abEeT7EdV9HozZi5eLyL3pvDDc2+eyG/ wljnfKoJfsLa+6pW3BSwquIJk70XOkRLHnw8j12DlmF34qvANDZ+IlqUK60k9NEfUBDR JG5CxURRPzqhbDBc20rQwxQEcCNKOK9A+FF0vz1hxvVGYTpCF6E0QEy2OlIW4YXzIeWo VYP9+ZkKqCXag+IuB+bjtUizqrA1YAffk1yziHZlUHqE7NYdy5gKvhgPWF9PLKawIGnk vGZQ==

In-reply-to: <8FC0047FE39D4223B3599833882AD69A@adminPC>

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: <8FC0047FE39D4223B3599833882AD69A@adminPC>

Sender: owner-ferret_users@xxxxxxxx

Dear Changshui,

Pls see the attached method and sample.

Attachment: Xia.zip
Description: Zip archive

Basically, you just need to find out the bottom level of your data, set data at the other levels to be missing, and then do vertical average. You can get your bottom level data.

Cheers,

Tony

Sent from my iPhone

On Apr 16, 2018, at 9:54 AM, Changshui Xia <xiacs@xxxxxxxxxx> wrote:

Hello ferreters,

I want to draw the fill/contour figure of botom level temperature from a z-coordinate netcdf file data (MOM5 output ) . That means the level below this bottom level is sea bed, value is missing_value in netcdf file.

Does anyone have any idea how to achieve this?

Thank you in advance.

Changshui