[gradsusr] Converting to regular horizontal grid spacing.
Jose Augusto Paixão Veiga
veiga.uea at gmail.com
Tue Jul 8 14:50:05 EDT 2014
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
>> <https://twitter.com/intent/follow?screen_name=andrewrevering&source=email&variant=2.0&xd_token=ceae5c40173e5>
>> 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
>>>> <https://twitter.com/intent/follow?screen_name=andrewrevering&source=email&variant=2.0&xd_token=ceae5c40173e5>
>>>> 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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140708/9a573ba4/attachment-0001.html
More information about the gradsusr
mailing list