Re: Arpege grid to regular lat/lon?

Hello Aideen

I assume it looks like the grid description below.

There is no NCL built-in function for this type of regrid
operation on the ARPEGE grid.

The grid has 'nodes' (nb_points) so I don't think the SCRIP
regrid software will work.

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 ;
}

Aideen Foley wrote:
> Hi everyone,
> I am trying to regrid an Arpege grid to regular lat/lon.
> I searched the NCL site and found that someone else had asked the same
> thing a couple of years ago, and someone else had replied saying they
> would be putting a script on the NCL site to do just this.
> I cannot find this script! I have found arpege_1.ncl, the script for
> plotting Arpege data, but this is not really what I need... Does anyone
> have an example of a script for regridding Arpege data, like an
> rcm2rgrid for Arpege?
> Many thanks,
> Aideen Foley
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Dec 16 2008 - 16:06:08 MST