[gradsusr] Hello to all users of Grads

Wed Apr 3 06:20:11 EDT 2013

I got here a good two Scripts, one for writing data on the Grads command
window, and the second save as an image file. This first script is very
good and I use it, the second script has a problem and  i wanted to know if it
possible to set up a loop to automatically make the result a separate line
When I want to save the result in a file pictures?.
I'm pretty new in Grads, and the question came to for me to learn.

avi.

2013/4/2 Muhammad Yunus Ahmad Mazuki <ukm.yunus at gmail.com>

> Well, I have been wondering, whether you want to save in text form or
> picture form? If you use 'draw string' it will save as pictures, and just
> overlapping on each other unless you set different coordinate for the
> string drawing.
>
>
> On Mon, Apr 1, 2013 at 6:25 PM, avi kojokro <kojokroavi at gmail.com> wrote:
>
>> The problem was in two lines,in this command
>> while (time <= end_time)
>> and this command
>> time=start_time
>>
>> Now there is another problem, the result comes out from grads is signed
>> in one line, Is there a way to loop the result that it will record each
>> year separately?.I add the image file
>>
>>
>> 2013/4/1 Rupak Rajbhandari <rupak.rajbhandari at gmail.com>
>>
>>> It says error occured at line 22.
>>> your code should be
>>> 'set t 'sart_time
>>>
>>> -rupak
>>>
>>>
>>>
>>> On Mon, Apr 1, 2013 at 3:05 PM, avi kojokro <kojokroavi at gmail.com>wrote:
>>>
>>>> To Thomas Robinson Thanks for the example, I tried the script in this
>>>> way as you showing me :
>>>>
>>>>
>>>>
>>>> 'reinit'
>>>> 'set display color white'
>>>> f = 'http://goldsmr2.sci.gsfc.nasa.gov:80/dods/MATMNXSLV'
>>>> 'sdfopen 'f
>>>>
>>>>
>>>> start_time=3
>>>> end_time=12*34
>>>>
>>>> 'set parea 0.7 10.5 0.6 8'
>>>> 'set lat 32.5'
>>>> 'set lon 35'
>>>>
>>>> set t 'start_time'
>>>> while (time <=end_time
>>>> 'set t  'time
>>>>
>>>>
>>>> 'd t850-273.16'
>>>> cc=subwrd(result,4)
>>>> 'q time'
>>>> ctime = subwrd(result,3)
>>>> chour = substr(ctime,9,15)
>>>> 'draw string 0.5 8 'chour'='cc
>>>>
>>>>
>>>> time=time+12
>>>> endwhile
>>>>
>>>>
>>>> 'printim avrg.png x1000 y700'
>>>> 'end'
>>>>
>>>>
>>>> But I get an error , And I attach a screenshot of the error , What did
>>>> I do wrong?.
>>>>
>>>>
>>>> 2013/4/1 Thomas Robinson <ter at hawaii.edu>
>>>>
>>>>> Aloha Avi,
>>>>>
>>>>> It looks like your code for each time step is the same.  You can loop
>>>>> using while/endwhile.  Try this and change the variable end_time to
>>>>> whatever you want it to be:
>>>>>
>>>>>  'reinit'
>>>>> 'set display color white'
>>>>> f = 'http://goldsmr2.sci.gsfc.nasa.gov:80/dods/MATMNXSLV<http://goldsmr2.sci.gsfc.nasa.gov/dods/MATMNXSLV>
>>>>> '
>>>>> 'sdfopen 'f
>>>>>
>>>>> ** Set the start and end times (month number)
>>>>> start_time=3
>>>>> end_time=12*34
>>>>>
>>>>>
>>>>> 'set parea 0.7 10.5 0.6 8'
>>>>> 'set lat 32.5'
>>>>> 'set lon 35'
>>>>>
>>>>> ** Set the time to the start time
>>>>> time=start_time
>>>>> ** Loop through the times
>>>>> while (time <=end_time
>>>>> **Set time
>>>>> 'set t ' time
>>>>>
>>>>> 'd t850-273.16'
>>>>> cc=subwrd(result,4)
>>>>> 'q time'
>>>>> ctime = subwrd(result,3)
>>>>> chour = substr(ctime,9,15)
>>>>> 'draw string 0.5 8 'chour'='cc
>>>>>
>>>>> ** Set the time to the next year
>>>>> time=time+12
>>>>>
>>>>> endwhile
>>>>>
>>>>>
>>>>>
>>>>> On Sun, Mar 31, 2013 at 8:26 PM, avi kojokro <kojokroavi at gmail.com>wrote:
>>>>>
>>>>>>
>>>>>> Thank you.
>>>>>> I did not phrase the question correctly, I think that the example I
>>>>>> gave is not good. I will try to explain myself again.
>>>>>> First I must say that I am working with windows not Linux.
>>>>>> I want to get data from a server that already gives me a monthly
>>>>>> averages, here's a better example of what i means:
>>>>>>
>>>>>>
>>>>>> 'reinit'
>>>>>> 'set display color white'
>>>>>> f = 'http://goldsmr2.sci.gsfc.nasa.gov:80/dods/MATMNXSLV'
>>>>>> 'sdfopen 'f
>>>>>>
>>>>>>
>>>>>>
>>>>>> 'set parea 0.7 10.5 0.6 8'
>>>>>> 'set lat 32.5'
>>>>>> 'set lon 35'
>>>>>>
>>>>>> 'set t 3 '
>>>>>> 'd t850-273.16'
>>>>>> cc=subwrd(result,4)
>>>>>> 'q time'
>>>>>> ctime = subwrd(result,3)
>>>>>> chour = substr(ctime,9,15)
>>>>>> 'draw string 0.5 8 'chour'='cc
>>>>>>
>>>>>>
>>>>>> 'set t 15 '
>>>>>> 'd t850-273.16'
>>>>>> cc=subwrd(result,4)
>>>>>> 'q time'
>>>>>> ctime = subwrd(result,3)
>>>>>> chour = substr(ctime,9,15)
>>>>>> 'draw string 0.5 7.7 'chour'='cc
>>>>>>
>>>>>>
>>>>>> 'set t 27 '
>>>>>> 'd t850-273.16'
>>>>>> cc=subwrd(result,4)
>>>>>> 'q time'
>>>>>> ctime = subwrd(result,3)
>>>>>> chour = substr(ctime,9,15)
>>>>>> 'draw string 0.5 7.7 'chour'='cc
>>>>>>
>>>>>> How do I create a loop to do it automatically, the server provides
>>>>>> data from 1979 to 2013, and I want to get only the data of March.
>>>>>>
>>>>>> I hope I have explained myself better.
>>>>>> Regards
>>>>>> avi.
>>>>>>
>>>>>>
>>>>>>
>>>>>> 2013/4/1 Muhammad Yunus Ahmad Mazuki <ukm.yunus at gmail.com>
>>>>>>
>>>>>>> The average function itself has what you needed. ave(var1,t=1,t=360,
>>>>>>> *12*), this means it will average values from t=1,t=13,t=25 ...
>>>>>>> until t=360.Example of point value extraction.
>>>>>>>
>>>>>>> 'sdfopen foo.nc'
>>>>>>> lonvalue=90
>>>>>>> latvalue=0
>>>>>>> levvalue=850
>>>>>>> tvalue=1
>>>>>>> 'set lon 'lonvalue
>>>>>>> 'set lat 'latvalue
>>>>>>> 'set t 'tvalue
>>>>>>> 'set lev 'levvalue
>>>>>>> 'd variable1'
>>>>>>> catch=sublin(result,4)
>>>>>>> 'set gxout fwrite'
>>>>>>> 'set fwrite point.dat'
>>>>>>> 'd 'catch
>>>>>>> 'disable fwrite'
>>>>>>>
>>>>>>> This is a basic one, modify it as your need requires. I'm not sure
>>>>>>> about excel, but this will create a binary file, which can be opened by
>>>>>>> GrADS (using a controller file) and Scilab (which I use for further
>>>>>>> analysis).
>>>>>>>
>>>>>>> Yunus.
>>>>>>>
>>>>>>>
>>>>>>> On Sun, Mar 31, 2013 at 2:25 PM, avi kojokro <kojokroavi at gmail.com>wrote:
>>>>>>>
>>>>>>>> I'm trying to do a script of averages like this example :
>>>>>>>>
>>>>>>>>
>>>>>>>> 'set t 1 '
>>>>>>>> 'define varav=ave(t850-273.16,t=1,t=1,1)'
>>>>>>>> 'd varav'
>>>>>>>> cc=subwrd(result,4)
>>>>>>>> 'draw string 0.5 8 jan='cc
>>>>>>>>
>>>>>>>> I want my script calculates every 12 months, how do I create a loop
>>>>>>>> script that calculates every 12 month average of the same month.
>>>>>>>> I will give an example of what I mean: I want to begin to calculate average
>>>>>>>> of January 1979 and then later to calculate directly the January
>>>>>>>> 1980 and so on.
>>>>>>>> Is it possible that the result will write directly to the Excel
>>>>>>>> file, then how to write the command to it.
>>>>>>>>
>>>>>>>> Please note that I first started Grads.
>>>>>>>>
>>>>>>>> Thank you so much.
>>>>>>>>
>>>>>>>> avi
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> _______________________________________________
>>>>>>>> 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
>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>> --
>>>>> Tom Robinson
>>>>> President - Graduate Student Organization
>>>>> Graduate Student - Department of Meteorology
>>>>> 732-718-2323
>>>>>
>>>>> _______________________________________________
>>>>> 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/20130403/686c7100/attachment-0003.html 