<div dir="ltr">Dear all,<div><br></div><div>I've been struggling with this issue for the last days without any glimpse of solution.</div><div><br></div><div>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).</div>
<div><br></div><div>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).</div>
<div><br></div><div>Each time I tried to calculated the difference between both fields I've got the following message:</div><div><br></div><div><div>Data Request Error: Invalid grid coordinates</div>
<div> World coordinates convert to non-integer grid coordinates </div><div> Variable = apcpsfc.2 Dimension = 3 </div>
<div> Error ocurred at column 1 </div><div>DEFINE error: Invalid expression.</div></div><div>
<br></div><div><br></div><div>I'd tried to work around with regrid but still with the same results.</div><div>Any help will be much appreciated.</div><div><br></div><div>Tanks in advance,</div>
<div>Susana Cardoso</div><div><br></div><div><br></div><div><br></div><div>Code </div><div>-------------------------------------------------------------------------------</div><div>
sdfopen <a href="http://eobs_rr.nc" target="_blank">eobs_rr.nc</a> -----> daily precip files from 1950 to 2012)</div><div>set t 1</div><div>define rrdec=sum(rr,time=1dec2009,time=31dec2009)</div>
<div>open wrfprs_d01_744.ctl -------> month precip for dec 2009</div><div>set lat 34 48</div><div>set lon -16 -2</div><div>set t 1</div><div>define d=apcpsfc.2-rrdec</div>
<div><br></div><div>EOBS q ctlinfo</div><div>--------------------------------------------------</div><div><div>dset <a href="http://eobs_rr.nc" target="_blank">eobs_rr.nc</a></div><div>title </div>
<div>undef -9.99e+33 </div><div>dtype netcdf </div>
<div>xdef 464 linear -40.375 0.25 </div><div>ydef 201 linear 25.375 0.25 </div>
<div>zdef 1 linear 0 1 </div><div>tdef 23011 linear 00Z01JAN1950 1440mn </div>
<div>vars 1 </div><div>rr=>rr 0 t,y,x rainfall </div>
<div>endvars</div><div><br></div><div><br></div></div><div>WRF description file</div><div>------------------------------------------------------------------------------------------------------------------------</div>
<div><div>dset ^wrfprs_d01.744</div><div>index ^wrfprs_d01.744.idx</div><div>undef 9.999E+20</div><div>title wrfprs_d01.744</div><div>* produced by grib2ctl v0.9.12.6</div><div>dtype grib 255</div><div>pdef 59 72 lccr 33.253000 -17.031000 1 1 41.500000 41.500000 -9.300000 25000 25000</div>
<div>xdef 59 linear -17.031000 0.271773636616826</div><div>ydef 72 linear 33.253000 0.227272727272727</div><div>tdef 2 linear 16Z10dec2009 512hr</div><div>* z has 10 levels, for prs</div><div>zdef 10 levels</div><div>1000 875 750 625 500 375 250 125 30 2</div>
<div>vars 6</div><div>ACPCPsfc 0 63,1,0 ** surface Convective precipitation [kg/m^2]</div><div>APCPsfc 0 61,1,0 ** surface Total precipitation [kg/m^2]</div><div>NCPCPsfc 0 62,1,0 ** surface Large scale precipitation [kg/m^2]</div>
<div>PRATEsfc 0 59,1,0 ** surface Precipitation rate [kg/m^2/s]</div><div>REFD2mb 0 211,100,2 ** 2 mb Derived radar reflectivity [dbZ]</div><div>RWMRprs 10 170,100,0 ** (profile) Rain water mixing ratio [kg/kg]</div><div>
ENDVARS</div><div><br></div><div><br></div></div><div><br></div></div>