[gradsusr] help with regrid

Jose Augusto Paixão Veiga veiga.uea at gmail.com
Tue Nov 29 09:12:36 EST 2011


Dear Nick I changed things in the script and in CTL as you suggest me, but the problem persist.

My Script:

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


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

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

'd variavel'

'disable fwrite'


The CTL modified

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 0 99  zonal wind            (U)  [m/sec]   
endvars

Any other hint?


Thanks for help.


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/
-----------------------------------------------------------------------------

On Nov 29, 2011, at 9:48 AM, Nick Heath wrote:

> José Augusto P. Veiga,
> 
> You can remove the 'set gxout shaded' when you are writing to a file.  In the control file for the new data you created trying making this change:
> 
> 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 0 99  zonal wind            (U)  [m/sec]  
> endvars
> 
> In the past, I have found that in the variable description I needed 0 99 instead of what was originally in there.  I'm not sure if that will fix it, but it's worth a shot.
> Good luck.
> 
> -Nick
> 
> 
> ----- Original Message -----
> From: Jose Augusto Paixão Veiga <veiga.uea at gmail.com>
> Date: Tuesday, November 29, 2011 1:36 pm
> Subject: [gradsusr] help with regrid
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> 
>> 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/
>> --------------------------------------------------------------------
>> ---------
>> 
>> 
>> 
>> 
> 
> 
> 
> 
> 
> 
> 
> Nick Heath
> Graduate Student
> Florida State University
> Department of Earth, Ocean, and Atmospheric Science
> Love Bldg. Rm 311
> 
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20111129/3e8328a1/attachment-0003.html 


More information about the gradsusr mailing list