[gradsusr] Data request error: Invalid grid coordinates

[Susana Cardoso Pereira]

Dear all,

I've been struggling with this issue for the last days without any glimpse
of solution.

I'm trying to calculate the difference between
two fields from different models with different grids and similar
resolutions (see ctl info and code below).

One is a netcdf file from E-OBS ECAD and the other one a grib file from WRF
(previous a netcdf converted to grib with the UPP tool and then used the
grib2ctl and gribmap to create the ctl and the idx files, respectively).

Each time I tried to calculated the difference between both fields I've got
the following message:

Data Request Error: Invalid grid coordinates
   World coordinates convert to non-integer  grid coordinates


   Variable = apcpsfc.2  Dimension = 3


   Error ocurred at column 1


DEFINE error:  Invalid expression.


I'd tried to work around with regrid but still with the same results.
Any help will be much appreciated.

Tanks in advance,
Susana Cardoso



Code
-------------------------------------------------------------------------------
sdfopen eobs_rr.nc
-----> daily precip files from 1950 to 2012)
set t 1
define rrdec=sum(rr,time=1dec2009,time=31dec2009)
open wrfprs_d01_744.ctl
 -------> month precip for dec 2009
set lat 34 48
set lon -16 -2
set t 1
define d=apcpsfc.2-rrdec

EOBS q ctlinfo
--------------------------------------------------
dset eobs_rr.nc
title


undef -9.99e+33


dtype netcdf


xdef 464 linear -40.375 0.25


ydef 201 linear 25.375 0.25


zdef 1 linear 0 1


tdef 23011 linear 00Z01JAN1950 1440mn


vars 1


rr=>rr  0  t,y,x  rainfall


endvars


WRF description file
------------------------------------------------------------------------------------------------------------------------
dset ^wrfprs_d01.744
index ^wrfprs_d01.744.idx
undef 9.999E+20
title wrfprs_d01.744
*  produced by grib2ctl v0.9.12.6
dtype grib 255
pdef 59 72 lccr 33.253000 -17.031000 1 1 41.500000 41.500000 -9.300000
25000 25000
xdef 59 linear -17.031000 0.271773636616826
ydef 72 linear 33.253000 0.227272727272727
tdef 2 linear 16Z10dec2009 512hr
*  z has 10 levels, for prs
zdef 10 levels
1000 875 750 625 500 375 250 125 30 2
vars 6
ACPCPsfc  0 63,1,0  ** surface Convective precipitation [kg/m^2]
APCPsfc  0 61,1,0  ** surface Total precipitation [kg/m^2]
NCPCPsfc  0 62,1,0  ** surface Large scale precipitation [kg/m^2]
PRATEsfc  0 59,1,0  ** surface Precipitation rate [kg/m^2/s]
REFD2mb  0 211,100,2 ** 2 mb Derived radar reflectivity [dbZ]
RWMRprs 10 170,100,0 ** (profile) Rain water mixing ratio [kg/kg]
ENDVARS
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130408/4f9bb1f3/attachment-0003.html