[Thread Prev][Thread Next][Index]

Re: [ferret_users] Spatial average from multiple files;




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.



[Thread Prev][Thread Next][Index]