[gradsusr] Converting to regular horizontal grid spacing.
Jose Augusto Paixão Veiga
veiga.uea at gmail.com
Wed Jul 9 19:22:02 EDT 2014
Thank you Andrew for your patience to help me with this,
but, I tried again to use the CDO software, as Kishore suggest me, and I
solve my the problem with the re-griding.
I tried to use lterp command to regrid the data. However, after running the
gs program I just received the message *segmentation fault*.
Thank again
José Augusto P. Veiga,
======================================
Coordenador do curso de Meteorologia da UEA
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, Ramal 4333
Cell phone : (92) 8196 7122
Skype: veiga_j.a.p.
-----------------------------------------------------------------------------
CV: http://lattes.cnpq.br/4027612512091565
URL:http://scientificmet.wordpress.com/
-----------------------------------------------------------------------------
On Tue, Jul 8, 2014 at 3:06 PM, Andrew Friedman <andfried at berkeley.edu>
wrote:
> The GrADS function ‘lterp’ is similar to OpenGrADS ‘re’, but doesn't have
> all the same features.
> http://www.iges.org/grads/gadoc/gradfunclterp.html
>
> On Jul 8, 2014, at 11:55 AM, Andrew Revering <andy at f5data.com> wrote:
>
> > You may need the OpenGrADS distribution.
> >
> >
> >
> > Andrew Revering - Follow @andrewrevering on Twitter
> > Convective Development, Inc.
> > PO Box 511, Cambridge, MN 55008
> > http://www.convectivedevelopment.com/
> >
> >
> > On Tue, Jul 8, 2014 at 1:50 PM, Jose Augusto Paixão Veiga <
> veiga.uea at gmail.com> wrote:
> > Hi Andrew and Roger,
> >
> > thank you by relevant hints and patience with my problem.
> >
> > I tried to do you've suggest me. However, after giving the command
> 'define precnew = re(prec,0.2)' in my gs file I get the message bellow:
> >
> >
> > Syntax Error: Invalid Operand
> > 're' not a variable or function name
> > Error ocurred at column 1
> > DEFINE error: Invalid expression.
> >
> > What should be happening?
> >
> >
> >
> >
> >
> > José Augusto P. Veiga,
> >
> > ======================================
> > Coordenador do curso de Meteorologia da UEA
> > 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, Ramal 4333
> > Cell phone : (92) 8196 7122
> > Skype: veiga_j.a.p.
> >
> -----------------------------------------------------------------------------
> > CV: http://lattes.cnpq.br/4027612512091565
> > URL:http://scientificmet.wordpress.com/
> >
> -----------------------------------------------------------------------------
> >
> >
> >
> >
> >
> > On Tue, Jul 8, 2014 at 7:51 AM, Roger Rodrigues Torres <
> torres.fisico at gmail.com> wrote:
> > Dear Jose,
> > as Andrew said, the RE function on Grads will perfectly do what you need.
> > 'define precnew = re(prec,0.5)'
> > In this way, the regrid function will perform a bilinear interpolation
> to the new grid,and it will not create undefined values.
> > After you have defined the new variable in the new grid (in this
> example, precnew), you can save the output in a netcdf file using the
> following comand:
> > 'set sdfwrite output.nc
> > 'sdfwrite precnew'
> >
> > Let me know if this work.
> >
> > Cheers,
> > Roger Torres
> >
> > Dr. Roger Rodrigues Torres
> > Professor Auxiliar I
> > Universidade Federal de Itajubá - Unifei
> > Instituto de Recursos Naturais - IRN
> > Curriculum Lattes: http://lattes.cnpq.br/5641165392004404
> > Webpage: http://rtorres.webnode.com.br
> >
> >
> > 2014-07-08 2:30 GMT-03:00 Andrew Revering <andy at f5data.com>:
> >
> > It looks like your irregular grid is in a variable 'prec'.
> >
> > So if you were to do:
> >
> > 'define precnew = re(prec,0.5)'
> >
> > you would have it interpolated to a 0.5 degree grid. You can change the
> 0.5 to whatever spacing you think is appropriate.
> >
> > I haven't done anything like this before, so I'm not sure if the re()
> function is going to give you an error if it's not a regular grid, but it
> might work. :-)
> >
> > I'm also not sure what will happen when you get to regions where your
> grid has sparse or missing data. That's where you might want to experiment
> with the degree in your new grid spacing. A larger number will probably
> allow it to interpolate farther out, but a very small value may return a
> bunch of missing data points...I'm not sure.
> >
> > Others who know more about this can hopefully chime in.
> >
> >
> >
> > Andrew Revering - Follow @andrewrevering on Twitter
> > Convective Development, Inc.
> > PO Box 511, Cambridge, MN 55008
> > http://www.convectivedevelopment.com/
> >
> >
> > On Mon, Jul 7, 2014 at 11:41 PM, Jose Augusto Paixão Veiga <
> veiga.uea at gmail.com> wrote:
> > Hi Andrew,
> >
> > I presume the command 'define newgrid = re(oldgrid,1.0)' you've suggest
> me can be work in GrAds. If so, old grid would be a function or a specific
> value? It is not clear yet for me.
> >
> > Yes, I am trying to use regular grid date instead not regular.
> >
> > My original ctl file is bellow. But I need have my dx and dy equally
> spaced.
> >
> > I have total precipitation.
> >
> > Thanks for your e-mail.
> >
> > dset ^sum-prec-OCT-APR-KF-d01.bin
> > undef 1.e30
> > title OUTPUT FROM WRF V3.5 MODEL
> > xdef 103 levels
> > 262.8432
> > 263.5718
> > 264.3003
> > 265.0289
> > 265.7574
> > 266.4860
> > 267.2146
> > 267.9431
> > 268.6717
> > 269.4003
> > 270.1288
> > 270.8574
> > 271.5860
> > 272.3145
> > 273.0431
> > 273.7717
> > 274.5002
> > 275.2288
> > 275.9574
> > 276.6859
> > 277.4145
> > 278.1431
> > 278.8716
> > 279.6002
> > 280.3287
> > 281.0573
> > 281.7859
> > 282.5144
> > 283.2430
> > 283.9716
> > 284.7001
> > 285.4287
> > 286.1573
> > 286.8858
> > 287.6144
> > 288.3430
> > 289.0715
> > 289.8001
> > 290.5287
> > 291.2572
> > 291.9858
> > 292.7144
> > 293.4429
> > 294.1715
> > 294.9001
> > 295.6286
> > 296.3572
> > 297.0858
> > 297.8143
> > 298.5429
> > 299.2714
> > 300.0000
> > 300.7286
> > 301.4571
> > 302.1857
> > 302.9142
> > 303.6428
> > 304.3714
> > 305.0999
> > 305.8285
> > 306.5571
> > 307.2856
> > 308.0142
> > 308.7428
> > 309.4713
> > 310.1999
> > 310.9285
> > 311.6570
> > 312.3856
> > 313.1142
> > 313.8427
> > 314.5713
> > 315.2999
> > 316.0284
> > 316.7570
> > 317.4856
> > 318.2141
> > 318.9427
> > 319.6713
> > 320.3998
> > 321.1284
> > 321.8569
> > 322.5855
> > 323.3141
> > 324.0426
> > 324.7712
> > 325.4998
> > 326.2283
> > 326.9569
> > 327.6855
> > 328.4140
> > 329.1426
> > 329.8712
> > 330.5997
> > 331.3283
> > 332.0569
> > 332.7854
> > 333.5140
> > 334.2426
> > 334.9711
> > 335.6997
> > 336.4283
> > 337.1568
> > ydef 89 levels
> > -34.71973
> > -34.11872
> > -33.51342
> > -32.90385
> > -32.29005
> > -31.67207
> > -31.04995
> > -30.42374
> > -29.79348
> > -29.15923
> > -28.52103
> > -27.87895
> > -27.23304
> > -26.58335
> > -25.92998
> > -25.27294
> > -24.61234
> > -23.94821
> > -23.28067
> > -22.60976
> > -21.93555
> > -21.25815
> > -20.57760
> > -19.89401
> > -19.20746
> > -18.51802
> > -17.82580
> > -17.13087
> > -16.43335
> > -15.73330
> > -15.03085
> > -14.32607
> > -13.61906
> > -12.90993
> > -12.19879
> > -11.48574
> > -10.77088
> > -10.05432
> > -9.336166
> > -8.616524
> > -7.895515
> > -7.173248
> > -6.449829
> > -5.725380
> > -5.000008
> > -4.273834
> > -3.546967
> > -2.819534
> > -2.091644
> > -1.363419
> > -0.6349640
> > 9.3582153E-02
> > 0.8221283
> > 1.550529
> > 2.278679
> > 3.006462
> > 3.733765
> > 4.460464
> > 5.186440
> > 5.911591
> > 6.635788
> > 7.358917
> > 8.080879
> > 8.801544
> > 9.520813
> > 10.23856
> > 10.95470
> > 11.66911
> > 12.38168
> > 13.09231
> > 13.80090
> > 14.50734
> > 15.21154
> > 15.91339
> > 16.61279
> > 17.30966
> > 18.00390
> > 18.69541
> > 19.38412
> > 20.06992
> > 20.75273
> > 21.43247
> > 22.10907
> > 22.78243
> > 23.45250
> > 24.11917
> > 24.78239
> > 25.44209
> > 26.09818
> > zdef 1 linear 1 1
> > tdef 213 linear 12Z01OCT2011 1440MN
> > VARS 1
> > prec 1 0 ACUMULATED TOTAL PRECIPITATION (CUMULUS + GRID SCALE) (mm)
> > ENDVARS
> >
> >
> >
> >
> > José Augusto P. Veiga,
> >
> > ======================================
> > Coordenador do curso de Meteorologia da UEA
> > 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, Ramal 4333
> > Cell phone : (92) 8196 7122
> > Skype: veiga_j.a.p.
> >
> -----------------------------------------------------------------------------
> > CV: http://lattes.cnpq.br/4027612512091565
> > URL:http://scientificmet.wordpress.com/
> >
> -----------------------------------------------------------------------------
> >
> >
> >
> >
> >
> > On Tue, Jul 8, 2014 at 12:19 AM, Andrew Revering <andy at f5data.com>
> wrote:
> > This is a real guess since I've never worked with an irregular grid in
> GrADS, but couldn't you just regrid to a certain degree grid and it would
> interpolate to a regular grid?
> >
> > IE
> >
> > 'define newgrid = re(oldgrid,1.0)'
> >
> > would regrid it to a 1.0° grid??
> >
> > While we're on the topic, how would you import/display an irregular grid
> in GrADS? Lets say you have surface temperatures to plot?
> >
> > Andrew Revering - Follow @andrewrevering on Twitter
> > Convective Development, Inc.
> > PO Box 511, Cambridge, MN 55008
> > http://www.convectivedevelopment.com/
> >
> >
> > On Mon, Jul 7, 2014 at 11:13 PM, Jose Augusto Paixão Veiga <
> veiga.uea at gmail.com> wrote:
> > Dear Kishore,
> >
> > I tried this (cdo -f nc import_binary my.ctl output.nc), however I just
> receive the message bellow:
> >
> > cdo import_binary (Abort): I/O error reading record!
> >
> >
> > Do you have any idea what does it mean?
> >
> > Thanks in advance.
> >
> >
> >
> >
> > On Mon, Jul 7, 2014 at 11:35 PM, Kishore Ragi <kishoreragi at gmail.com>
> wrote:
> > cdo -f nc import_binary xxxxxxx.ctl yyyyyyy.nc
> >
> >
> >
> >
> > José Augusto P. Veiga,
> >
> > ======================================
> > Coordenador do curso de Meteorologia da UEA
> > 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, Ramal 4333
> > Cell phone : (92) 8196 7122
> > Skype: veiga_j.a.p.
> >
> -----------------------------------------------------------------------------
> > CV: http://lattes.cnpq.br/4027612512091565
> > URL:http://scientificmet.wordpress.com/
> >
> -----------------------------------------------------------------------------
> >
> >
> >
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> _______________________________________________
> 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/20140709/7a4ea416/attachment-0001.html
More information about the gradsusr
mailing list