[gradsusr] help with regrid

[Jose Augusto Paixão Veiga]

Dear all,

I am trying to regrid NCEP (2.5 x 2.5, one time, one level) data into a ECMWF (1.5x1.5, one time, one level) data, with lterp command, issuing to compare both. To do this job, I used the script bellow:

'reinit'
* Data to be regrided into the defalt data
'open uvtwg_ncep.ctl'

* Default data
'open uvtwg_ecmwf.ctl'

'variavel=lterp(u.1,u.2)'
'd variavel'

'set gxout fwrite'
'set fwrite uvtwg_new.bin'

'set gxout shaded'
'd variavel'

'disable fwrite' 



After running the above script I construct the ctl file to receive the uvtwg_new.bin. And so far so good. However, when I opened GrAds and displayed the variable recorded it saw as it was shifted in relation to the original data.

I have attached both figures and included bellow the CTL files that I am using as a hope to some help.


Original CTL (NCEP)

dset   uvtwg_ncep.bin
undef  9.999E+20 
title  NCEP 2.5x2.5 graus 
xdef 144  linear     0.0 2.5
ydef  73  linear   -90.0 2.5
zdef   1  linear     1   1 
tdef   1  linear  JAN1979 1mo 
vars  5
u 1  33,100  zonal wind            (U)  [m/sec]   
v 1  34,100  meridional wind       (V)  [m/sec]   
t 1  11,100  temperature           (T)      [C]   
w 1  39,100  vertical velocity     (w)  [Pa/s]    
g 1  40,100  geopotencial          (w)  [m²/s²]

Defalt CTL (ECMWF)

dset   uvtwg_ecmwf.bin
undef  9.999E+20 
title  ERA INTERIN 1.5x1.5 grau 
xdef 240 linear     0.0 1.5
ydef 121 linear   -90.0 1.5
zdef   1 linear     1   1 
tdef   1 linear  JAN1979 1mo 
vars  5
u 1  33,100  zonal wind            (U)  [m/sec]   
v 1  34,100  meridional wind       (V)  [m/sec]   
t 1  11,100  temperature           (T)      [C]   
w 1  39,100  vertical velocity     (w)  [Pa/s]    
g 1  40,100  geopotencial          (w)  [m²/s²]   
endvars

CTL for the regrided data (uvtwg_new.bin)

dset   uvtwg_new.bin
undef  9.999E+20 
title  NCEP 1.5x1.5 graus 
xdef 240  linear     0.0 1.5
ydef 121  linear   -90.0 1.5
zdef   1  linear     1   1 
tdef   1  linear  JAN1979 1mo 
vars  1
u 1  33,100  zonal wind            (U)  [m/sec]   
endvars


Thanks in advance,


José Augusto P. Veiga, 

======================================
Universidade do Estado do Amazonas 
Departamento de Meteorologia
Escola Superior de Tecnologia (EST)
-----------------------------------------------------------------------------
Av. Darcy Vargas, 1200, Manaus-AM Brasil
Work phone:    (92) 3878 4333
Skype: veiga_j.a.p.
-----------------------------------------------------------------------------
CV: http://lattes.cnpq.br/4027612512091565
URL:http://scientificmet.wordpress.com/
http://nmet.wordpress.com/
-----------------------------------------------------------------------------



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20111129/21cb8216/attachment-0003.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: original.gif
Type: image/gif
Size: 43692 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20111129/21cb8216/attachment-0006.gif 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: regrided.gif
Type: image/gif
Size: 43851 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20111129/21cb8216/attachment-0007.gif