compprocchains: Compare EUMETSAT OSISAF products from two different...

In steingod/R-osisaf: Useful functions for handling EUMETSAT Ocean and Sea Ice SAF data

Description

Compare OSISAF products from two different processing chains. This function takes two directory names as input, it collects filenames from these directories, checks that the file is HDF and that it is contained within both directories, before reading it. Data read are compared by plotting a time series of the mean differences at the active graphics device.

Usage

compprocchains(chain1,chain2,from=ISOdate(1996,1,1),to=Sys.time(),area=NULL) 

Arguments

chain1	Full path to first directory to read
chain2	Full path to second directory to read
from	ISOdate specification of oldest file to process
to	ISOdate specification of newest file to process
area	Perl regexp to handle directories with several products (e.g. sea ice)

Details

Two directories are given as input to the function. Within each directory all products fulfilling the time and area specification is read. For each product a mean value is calculated and the difference in the mean values for the same product within each directory is generated.

Value

Returns a data.frame containing a time identification, the difference between the two products, and the mean value for each of the two products.

time	Time specification as generated by ISOdatetime
data	Difference between the two processing chains for each time specification

Note

This function is probably of most interest to the development team when performing software updates.

This function has sometimes proven buggy, requiring a restart of R. The most plausible reason is a memory bug when reading HDF5 files, but this is yet undetected. Results shown are correct.

Author(s)

Øystein Godøy ([e-mail])

See Also

readosisaf

Examples

    # Compare operational and pre-operational processing
    compprocchains("/data/osisaf/output/ice/multi/conc/",
    "/data/osisaf_it/output/ice/multi/conc/",
    from=ISOdate(2006,4,1),to=ISOdate(2006,6,30),area="sh")

steingod/R-osisaf documentation built on May 30, 2019, 2:32 p.m.