incompatible grids

Ranyére Nóbrega ranyere.nobrega at SIPAM.GOV.BR
Wed Mar 12 19:02:42 EDT 2008


Use regrid function.
define example = re(var,1,ba)
for example, were 1 is grid resolution.


Ted Pollock escreveu:
> Hi all,
>
> I have two global grids of different resolutions and I'd like to
> compare them but I'm having problems since the two sets of grid points
> do not line up. (I get an incompatible grid error if I try to take the
> difference between the two.) The higher resolution grid is netcdf
> format and the other is not.
>
> I thought the easiest way to do this would be to use the pdef entry
> and its interpolations to map the non-netcdf grid onto the higher
> resolution netcdf grid. Unfortunately I'm not sure how to get pdef to
> work with the mercator projection.
>
> I'm also open to mapping the netcdf grid onto the lower resolution
> grid using a ctl file if thats possible but I think that might also
> require a pdef entry. (I'm still using gradsv1.9b4.)
>
> Any help on this would be much appreciated. (Below I've included parts
> of the ctl file and ncdump from the two datasets.)
>
> Thanks in advance
> -Ted
>
> CTL FILE:
>
> XDEF 96 LINEAR  0.0   3.75
> *
> YDEF 80 GAUSR30 1
> *
>
> NCDUMP:
> dimensions:
>          lon = 192 ;
>          lat = 94 ;
>          time = UNLIMITED ; // (722 currently)
>
> variables:
>          float lat(lat) ;
>                  lat:units = "degrees_north" ;
>                  lat:actual_range = 88.542f, -88.542f ;
>                  lat:long_name = "Latitude" ;
>          float lon(lon) ;
>                  lon:units = "degrees_east" ;
>                  lon:long_name = "Longitude" ;
>                  lon:actual_range = 0.f, 358.125f ;
>          double time(time) ;
>


--
__________________________________________________________
Ranyére Silva Nóbrega
Analista Intelectual
Divisão de Meteorologia e Climatologia - SIPAM
Centro Técnico e Operacional de Porto Velho



More information about the gradsusr mailing list