[gradsusr] help with regrid
Nick Heath
nkh09 at fsu.edu
Tue Nov 29 09:22:25 EST 2011
Hey,
I think I have found the problem in your script. Make these changes:
'reinit'
* Data to be regrided into the defalt data
'open uvtwg_ncep.ctl'
* Default data
'open uvtwg_ecmwf.ctl'
'set dfile 2'
'set x 1 240'
'set y 1 121'
'variavel=lterp(u.1,u.2)'
'set gxout fwrite'
'set fwrite uvtwg_new.bin'
'd variavel'
'disable fwrite'
I believe that will work.
-Nick
----- Original Message -----
From: Jose Augusto Paixão Veiga <veiga.uea at gmail.com>
Date: Tuesday, November 29, 2011 2:14 pm
Subject: Re: [gradsusr] help with regrid
To: GrADS Users Forum <gradsusr at gradsusr.org>
> 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
>
>
Nick Heath
Graduate Student
Florida State University
Department of Earth, Ocean, and Atmospheric Science
Love Bldg. Rm 311
More information about the gradsusr
mailing list