[gradsusr] Converting to regular horizontal grid spacing.
Roger Rodrigues Torres
torres.fisico at gmail.com
Tue Jul 8 07:51:17 EDT 2014
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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140708/c5dab6d8/attachment-0001.html
More information about the gradsusr
mailing list