[gradsusr] Hello to all users of Grads

Mon Apr 1 21:52:09 EDT 2013

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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130402/f6567dc7/attachment-0003.html 