[gradsusr] Converting to regular horizontal grid spacing.
Andrew Friedman
andfried at berkeley.edu
Tue Jul 8 15:06:01 EDT 2014
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
More information about the gradsusr
mailing list