Re: Missing in daily observation.

The issue is not clear to me.

IDEALLY: TIME(NTIM), X(NTIM) NTIM=1440
                            TIME(TIME), X(TIME) ; netCDF
                                                            
                            TIME(0) = 0
                            TIME(1) = 1
                            TIME(59) = 59
                            TIME(60) = 100
                            TIME(61) = 101
                            TIME(1439) = 2359

REALITY: time(ntim), x(ntim) (say) ntim=500

                            time(0) = 8
                            time(1) = 9
                            time(2) =23

   TIME = new (1440, "integer", "No_FillValue")
   TIME!0 = "TIME"
   X = new (1440, "float" , 1e20)
   X!0 = "TIME"
where
    TIME(0) = 0
    TIME(1) = 1
    TIME(59) = 59
    TIME(60) = 100
    TIME(61) = 101
     TIME(1439) = 2359

     TIME&TIME = TIME ; create coordinate array

Maybe
     X(time) = x

Wang, Yi wrote:
> Hi Folks:
>
> I am not sure if there is a good way to deal with this situation.
>
> I have lots of daily observation at 1-minute interval. However, among them,
> there are lots of missed observation time slices.
>
> The time of observation is set up as this:
>
> 0:59, 100:159, 200:259, .......... 2300:2359 (total of 1440 for one day)
>
> How do I locate (detect) the missing observation period and set them to the
> missing value (Fill_Value)?
>
> I can do a loop over each daily file. However, is there any better approach
> than the loop?
>
> Thanks,
>
> Yi Wang, PhD
> Pacific Northwest National Lab
>
>

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Feb 15 2008 - 13:14:49 MST