Re: [ferret_users] Spatial average from multiple files;

[Karnan C <karnanc6@xxxxxxxxx>]

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.ncdownloaded 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 theimage, 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), ledthe empty regions, though other files having values.
>>
>> I should take theaverage 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