[Thread Prev][Thread Next][Index]
Re: [ferret_users] Spatial average from multiple files;
- To: Russ Fiedler <russell.fiedler@xxxxxxxx>
- Subject: Re: [ferret_users] Spatial average from multiple files;
- From: Karnan C <karnanc6@xxxxxxxxx>
- Date: Tue, 20 Jun 2017 10:50:53 +0530
- Arc-authentication-results: i=2; mx.google.com; dkim=pass header.i=@gmail.com header.b=G98mc8y0; spf=pass (google.com: domain of karnanc6@xxxxxxxxx designates 2a00:1450:400c:c09::22e as permitted sender) smtp.mailfrom=karnanc6@xxxxxxxxx; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=gmail.com
- Arc-authentication-results: i=1; mx.google.com; dkim=pass header.i=@gmail.com header.b=G98mc8y0; spf=pass (google.com: domain of karnanc6@xxxxxxxxx designates 2a00:1450:400c:c09::22e as permitted sender) smtp.mailfrom=karnanc6@xxxxxxxxx; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=gmail.com
- Arc-message-signature: i=2; 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-transfer-encoding:cc:to:subject:message-id:date:from :references:in-reply-to:mime-version:arc-authentication-results :arc-message-signature:dkim-signature:arc-authentication-results; bh=8Gy7Bq/ULW+xNtkavngbmaXCfyuLNGYrVFm/OqeQNAc=; b=v98iia02xu1YHt99+AV0MKp3O+Mt6Z1kIKBAJxhqsjASfltUuscWg3a8z+yGnzUjDu vxSbJmRkFIKnvJ4Sn645E0xG59Uq5HVqM0lUZe6ha/1ZDpzEEy1weZ4r+RgdCSNSNIUJ +w/dRI5J2gzWIySmFE5F4HJoqtCK5jCfRRaVzNCZycIG5uolkPI15TH0kYA+8KWZaea4 YqRNW3xo8mpfEsF1n1BxnEGqnFOGqeEiH/MU97z0uU9axDItktrJ2O78PejJf4mQKs0k xWB8beHJMutKBmYSEBW/Xq+mwZJ1zJ7kVWCA6+d01fHRW0TA66dyWEGdb4bOG1HnUUBK e6LA==
- Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=content-transfer-encoding:cc:to:subject:message-id:date:from :references:in-reply-to:mime-version:dkim-signature :arc-authentication-results; bh=8Gy7Bq/ULW+xNtkavngbmaXCfyuLNGYrVFm/OqeQNAc=; b=XMdOCtF6OYiDweemPTIS3FdeY0j7ldYZHl457pJB5GhnqiWq2lYicCmUGQHRd5ESts X+aOxdcX0FY1JO/a7rYJRZPEwbpYD+bdjDmg2kzw+bd3XArLpCY+SCh+k0Ca9Zl9lLBm ApxjVbXuN7nV2ex25gANQVgNVUgzoHj5WjTXuFWFFg6fK8Xw/lPwTulDD0DgN3VZL/Uk PqXf/dcnJqu4JEuuVgKCp/jk1xT5h2lKqUPZAgJuVNyMlcblpvYIKTekNjKdKc+cnAos 7IKdFLknUQ0NBFLCV3HhTUwVdKboicegVhxCvwTWxnmJyQn2lYlDt1p9qpq9dn44MQGE 2Qmg==
- Arc-seal: i=2; a=rsa-sha256; t=1497936055; cv=pass; d=google.com; s=arc-20160816; b=a48j/DPTok2wFDMBwHa9IZlwTGDyFXK0R35vdJ2yKLfnOAiA9a99GS0vdF3PCMilKD L1KOgn/w+i1WyZeIT5aR3IqB9YkO+f210d0arXMHSFJSExhY81z0Ro9wo41sfmL4nU2q eKbcFGfSFE9X3+40OfsoGxYLU/6PFmSFqNc/C9mldofK78nRHIka6Kj2bopZymwn15b1 oS2je+7rmj9VVSfAb1hjmWBK6h1ns/1FDy0fj5z81AVy6S1EBg3lSWzYKXhe/1cioqBq o+Z9TL2MKD4jl+zEnaRY8OMD+ERQCVxdCCnZt2V5O1TrIDwHnaU6WOTtx+w1ilwYzTo0 GEAQ==
- Arc-seal: i=1; a=rsa-sha256; t=1497936054; cv=none; d=google.com; s=arc-20160816; b=ILrIauly1ywMUSkM/K3YTJm+43lG8kQtzfm30JwycQnljRu9c33oa52AILllLH+Hun 319cT46s0aqWioNe9UhEjXWtKGVFxxFeMQg0wsw8CmHDnhx9Qbvm34H6scxvv66jS+Qt sOlGL9vb4HzJweGmyNwfOU0qTQe3wU/kFee3az3Jeid7AdhCow+BTIX9D9/xe6fdNytW rc2gYcX1RJ2LjI3ESZRcmjZNYk8Ppnh+LcT6Kf6N1RrZ+bW8k6M+FgkQmS5SIsmBtLYF 9KK0UI8K0lR9MNI9itwm7F7ZZTst0zp1sciepz3W01XSX/vu/1sX92Wx/bBToUtBDAhZ ec+g==
- Cc: ferret_users@xxxxxxxx
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa-gov.20150623.gappssmtp.com; s=20150623; h=mime-version:in-reply-to:references:from:date:message-id:subject:to :cc:content-transfer-encoding:x-original-sender :x-original-authentication-results:precedence:mailing-list:list-id :list-post:list-help:list-archive; bh=8Gy7Bq/ULW+xNtkavngbmaXCfyuLNGYrVFm/OqeQNAc=; b=Dvwq1hdTm+7/B8EjCepRd63tBncj9dzXXsByqvKY0nG92Ou/beCLFdc5GvWLQoprOn UPPPwKpUSsHZFE5pHToOK5v9TVpgO7udrbXiujYt+0LvkPIPCR/PdkmxfxqnO1DLhY4v rLoNCioASw8QcHFYxcfNiQ9w8dnLv3bt7xtERUXdxoAhdL4IJrcYP2no8zrOOyaNwynd fWPxwz+Ubaq6hBt3bErZLW8PcIrcECC/pXqZqWWpNCk5IfwrOB6a20aLnk00Antljk0F CgIAdvurHNPdON9JLlR9OK82IQNTxrCN0FHG2k8yZhM1EYy9dBnVWwNFFUjEiU1skZRG DZoA==
- In-reply-to: <aac88d67-fbbf-4326-519c-3960edbbeb89@csiro.au>
- 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: <CAFYXeb_NKcc8SyPKgLrAK7jfGxmNHJCjsG5YDTynRQsRwJ39Og@mail.gmail.com> <aac88d67-fbbf-4326-519c-3960edbbeb89@csiro.au>
- Sender: owner-ferret_users@xxxxxxxx
Thank you Dr Russ,
The suggested ensemble (DEFINE/DATA) method is working well.
yes? use 1.nc , 2.nc , 3.nc , 4.nc
yes? DEFINE DATA/AGGR/E/TITLE="kd" KD490_aggr= 1, 2, 3, 4
yes? shade MODISA_L3M_KD_2014_KD_490[d=KD490_aggr,m=@ave]
Here the fill value (-32767) also took role in deriving average,
though my actual data range is 0.01 to 5.
Then I used the following scripts to change the fill value for each file
yes? set var/bad=0 MODISA_L3m_KD_2014_Kd_490
Then I followed above method (ensemble).
Now I got the spatial average.
Thank you very much.
2017-06-16 11:43 GMT+05:30 Russ Fiedler <russell.fiedler@xxxxxxxx>:
>
>
> Hi,
>
> Try making an ensemble. This will make virtual dataset using the E/M dimension
>
> ENSEMBLE/title="KD490 Ensemble" kd490_ens=1,2,3,4
> show data kd490_ens
>
> shade MODISA_L3M_KD_2014_KD_490[d=kd490_ens,m=@ave]
>
>
> Look up DEFINE/DATA in the documentation for more info.
>
> Russ
>
>
> On 16/06/17 15:24, Karnan C wrote:
>>
>> Hi,
>>
>> I tried to derive spatial average of KD490 from four different files (1.nc, 2.nc, 3.nc and 4.nc downloaded from giovanni database)
>>
>> \\I used the following stpes
>>
>> yes? use 1.nc , 2.nc , 3.nc , 4.nc
>> yes? sh da
>> currently SET data sets:
>> 1> ./1.nc
>> name title I J K L M N
>> MODISA_L3M_KD_2014_KD_490
>> Diffuse attenuation coefficient 1:156 1:96 ... ... ... ...
>>
>> 2> ./2.nc
>> name title I J K L M N
>> MODISA_L3M_KD_2014_KD_490
>> Diffuse attenuation coefficient 1:156 1:96 ... ... ... ...
>>
>> 3> ./3.nc
>> name title I J K L M N
>> MODISA_L3M_KD_2014_KD_490
>> Diffuse attenuation coefficient 1:156 1:96 ... ... ... ...
>>
>> 4> ./4.nc (default)
>> name title I J K L M N
>> MODISA_L3M_KD_2014_KD_490
>> Diffuse attenuation coefficient 1:156 1:96 ... ... ... ...
>>
>>
>> yes? fill MODISA_L3M_KD_2014_KD_490[d=1]
>> yes? frame/file=Figure1.gif //Figure1 in the attachement
>> yes? fill MODISA_L3M_KD_2014_KD_490[d=2]
>> yes? frame/file=Figure2.gif //Figure2 in the attachement
>> yes? fill MODISA_L3M_KD_2014_KD_490[d=3]
>> yes? frame/file=Figure3.gif //Figure3 in the attachement
>> yes? fill MODISA_L3M_KD_2014_KD_490[d=4]
>> yes? frame/file=Figure4.gif //Figure4 in the attachement
>>
>> yes? fill (MODISA_L3M_KD_2014_KD_490[d=1]+MODISA_L3M_KD_2014_KD_490[d=2]+MODISA_L3M_KD_2014_KD_490[d=3]+MODISA_L3M_KD_2014_KD_490[d=4])/4
>>
>> yes? frame/file=Figure5.gif //Figure5 in the attachement
>>
>>
>> //In the image, most of the regions were empty while taking the average (the method read from ferret user group).
>> //It was because of the missing data in any of the file (eg. 4.nc), led the empty regions, though other files having values.
>>
>> I should take the average from these files.
>> How can I solve this?
>> Is it possible to replace the 'no data (....)' as 'zero (0)'?
>> or Is it possible to feed K and L axis with some values like L value = 1,2,3 and 4 for each .nc files, respectively? May it be useful to take average like
>> "yes? fill MODISA_L3M_KD_2014_KD_490 [l=1:4@AVE]".
>>
>> The sample data files attached here.
>>
>> Thanks.
>> --
>> C. Karnan, CSIR - SRF
>> National Institute of Oceanography (CSIR)
>> India.
>
>
--
C. Karnan, CSIR - SRF
Ph. D Research Scholar
National Institute of Oceanography (CSIR)
Regional Centre, Dr Salim Ali Road
Kochi - 682 018.
India.
Mobile No. +91 9633032345; +91 9941811608
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
Ferret
Privacy Policy | Disclaimer | Accessibility Statement