Dear Ferret Users,
I am using the following script to find the days when I have more than 70% data covered in the required area. I am able to do this correctly. I want to write those days data into a saparate file when I have more than 70% grid points available. Can anybody give me a clue to do that.
yes? go mask
set da avhrr_v5data.cdf
l dum if x[g=temp] gt 330.8 and x[g=temp] lt 331.2 and y[g=temp] gt 37.4 and y[g=temp] lt 38.4 and quality_level ge 4 then temp/temp
l np if x[g=temp] gt 330.8 and x[g=temp] lt 331.2 and y[g=temp] gt 37.4 and y[g=temp] lt 38.4 then x[g=temp]*y[g=temp]
l dum1 np/np
l per 100*(dum[x=@sum,y=@sum]/dum1[x=@sum,y=@sum])
l dum2 if per ge 70 and quality_level ge 4 then dum*temp
save/file=kk.nc dum2
LISTing to
file kk.nc
!repeat/l=1:729 (set win 1; shade dum2;set win 2; shade temp; pause)
yes? sh da
currently SET data sets:
1> ./avhrr_v5data.cdf (default)
name title I J K L
TEMP SEA_SURFACE_TEMPERATURE[D=1]-27 1:49 1:49 ... 1:729
QUALITY_LEVEL
SST measurement quality 1:49
1:49 ... 1:729
Thanks and Regards,
Sriram