How to expand the compression files to standard lon by lat grids

From: Song Feng (songfeng AT XXXXXX)
Date: Tue Jan 28 2003 - 13:41:43 MST

  • Next message: ENOMOTO Takeshi: "writing characters"

    Hello:

    To save space in the netCDF file, it is always desired to eliminate points from data arrays that are invariably missing. I am trying to read those kind of compression data file. The following are the information about the data showed by "ncdump":

    dimensions:
            x = 76576 ;
            tstep = UNLIMITED ; // (607 currently)
    variables:
            float lon(x) ;
                    lon:long_name = "Longitude" ;
                    lon:units = "degrees_east" ;
            float lat(x) ;
                    lat:long_name = "Latitude" ;
                    lat:units = "degrees_north" ;
            double time(tstep) ;
                    time:long_name = "time" ;
                    time:units = "months since 1950-01" ;
                    time:calendar = "gregorian" ;
            int timestep(tstep) ;
                    timestep:long_name = "timesteps" ;
                    timestep:units = "timesteps since 1950-01" ;
                    timestep:tstep_month = 1 ;
            float Qs(tstep, x) ;
                    Qs:long_name = "Average Monthly Surface Runoff Rate" ;
                    Qs:units = "kg/m2/s" ;

    The lat(x) varies from 53 to 25N and lon(x) varies form -125 to -67W and the data resolution is 0.125 longitude by 0.125 latitude. I am try to expand the two dimension variable Qs(tstep, x) to standard three dimensions variable qs(time, longitude, latitude) with the invariable missing grid points (Oceans or lakes area) filled by miss value (such as 999.9) because the expanded 3D data is much easier to handle. I wonder if you will be so kind to offer me some help or some suggestions.

    Thanks a lot.

    Song Fang

    _______________________________________________
    ncl-talk mailing list
    ncl-talk AT ucar.edu
    http://mailman.ucar.edu/mailman/listinfo/ncl-talk



    This archive was generated by hypermail 2b29 : Wed Jan 29 2003 - 09:56:54 MST