Arpege grid

Hi,
Anybody willing to share an example script to open Arpege grid data (
http://www.ncl.ucar.edu/Document/Graphics/contour_grids.shtml#Arpege) and
regrid into regular lat/lon. A sample nc header is given below.
Thanks
siji

netcdf zg_6h_19991201_19991231_3D {
dimensions:
        nb_points = 21692 ;
        gaussian_grid = 1 ;
        lat_circle = 120 ;
        time = UNLIMITED ; // (124 currently)
        plev = 16 ;
variables:
        int plev(plev) ;
                plev:standard_name = "Air pressure" ;
                plev:units = "hPa" ;
        double lon(gaussian_grid, nb_points) ;
                lon:standard_name = "longitude" ;
                lon:units = "degrees_east" ;
        double lat(gaussian_grid, nb_points) ;
                lat:standard_name = "latitude" ;
                lat:units = "degrees_north" ;
        int mask(gaussian_grid, nb_points) ;
                mask:standard_name = "Land/Sea mask" ;
                mask:units = "1 for land, 0 for ocean" ;
        double orog(gaussian_grid, nb_points) ;
                orog:standard_name = "surface_altitude" ;
                orog:units = "m" ;
        int nlon(lat_circle) ;
                nlon:standard_name = "number of longitudes" ;
        int time(time) ;
                time:axis = "t" ;
                time:calendar = "gregorian" ;
                time:standard_name = "Time axis" ;
                time:time_origin = "1999-11-01 12:00:00" ;
                time:units = "hours since 1999-11-01 12:00:00" ;
        int date(time) ;
                date:standard_name = "current date (YYYYMMDD)" ;
        int datesec(time) ;
                datesec:standard_name = "current seconds of current date" ;
        float zg(time, plev, gaussian_grid, nb_points) ;
                zg:standard_name = "Geopotential height" ;
                zg:units = "m" ;
                zg:vairmer_name = "PLGEO000" ;

// global attributes:
                :title = "SAFO output" ;
                :history = "Creation date : 2006-08-12 20:54:00" ;
                :subtitle = "Cerfacs Netcdf conversion" ;
                :Conventions = "CF-1.0" ;
                :institution = "Meteo-France" ;
                :source = "Arpege Climat" ;
                :contact = "Global Change Team (Cerfacs)" ;
                :reference = "Deque, M., C. Dreveton, A. Braun and D.
Cariolle, 1994: the climate version of Arpege/IFS : a contribution to the
french community climate modelling. Climate Dyn., 10, 249-266." ;
                :troncature = 106 ;
                :grid_type = 2 ;
                :pole_longitude = 0.209f ;
                :pole_sine_latitude = 0.643f ;
                :stretch_factor = 3.f ;
                :regular_longitude_nb = 240 ;
                :model_levels = 31 ;
}

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Aug 23 2006 - 09:09:48 MDT