[gradsusr] Converting to regular horizontal grid spacing.
Andrew Revering
andy at f5data.com
Tue Jul 8 14:55:37 EDT 2014
You may need the OpenGrADS distribution.
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 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
>>> <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
>>
>>
>
> _______________________________________________
> 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/e3fcbf4c/attachment-0001.html
More information about the gradsusr
mailing list