Re: 回复： how to get seasonal time series?

[Katherin Kullgren]

This is what i used to extract a seasonal time series for JJA from a  
monthly time series for 1950 - 2000
This example is precipiation, so you would sum  the values.  If you  
were doing it with Temperature, or somthing else, you would average  
the value.  Simply change sum(***) to ave(***)
Hope this helps.


* This creates the JJA season of Prate from Jun1950/Aug1950 - Jun1999/ 
Aug1999 (50 years)

'set gxout fwrite'
'set fwrite jja_prate.data'
imap = 6
jmap = 8
while (imap <= 594)
'd sum(prate,t='imap',t='jmap')'
imap = imap + 12
jmap = jmap + 12
endwhile

'disable fwrite'










On Nov 30, 2006, at 11:36 PM, dongdong wrote:

> Hi Manik,
> Thank you for your reply!  But, what I need is a JJA series (range  
> from 1958 to 2000), not a climatological variable. I am sorry if i  
> did't make myself clearly.
>
> I found a letter simily to my question, and Charles Seman suggested  
> to write out the results to a GrADS fwrite file. I tried, but the  
> fwrite variable isn't the same as what I defined.
>
> I attached my files. These descriptions are:
> form.vint.WaterFlux.gs    To define waterflux and write out to  
> GRADS fwrite file
> fuq.vint.ctl                      To open the fwrite file in grads
> define.fuq.gif                   horizontal water flux which I  
> drawed after defining it
> fwirte.fuq.gif                    horizontal water flux but form  
> fwite file
> grads_process.txt           commands I issued in grads
>
> Can anyone help me finding the problems? Thanks in advance.
>
> PS. my work here is checking whether my fwrite variables are the  
> same as what i defined. I haven't calculated JJA mean here. I will  
> do that after solving this problem.
>
> Best Regards
> dongdong
>
>
> ----- 原始邮件 ----
> 发件人： Manik Bali <manik at cdac.in>
> 收件人： dongdong <happydong365 at YAHOO.COM.CN>
> 已发送： 2006/12/1(周五), 上午7:14:15
> 主题： Re: how to get seasonal time series?
>
>
> Hi
> You can get the seasonal cycle using.
> Say you have 600 months of data of sst
> 'set t 1 12'
> 'define sstclim = ave(sst, t+0, t=600, 12)'
> 'modify sstclim seasonal'
>
> sstclim then becomes a seasonal variable having 12 time steps  
> ( jan...dev)
> Later on you can take average of JJA as ave(sstclim,t=6,t=8)
> Hope this helps
> Manik
>
>
>
>
> On Thu, Nov 30, 2006, dongdong <happydong365 at YAHOO.COM.CN> said:
>
>> Hello friends,
>>
>> Can you tell me how to get seasonal series from a monthly data  
>> set? I mean, The variables I have are monthly mean, from 1957.9 to  
>> 2002.8.  How can I extract JJA datasets, and get a summer mean  
>> time series from 1958 to 2000?
>>
>> Any suggestion will be appreciated! Thank you !
>>
>>
>>
>> ___________________________________________________________
>> 锟脚伙拷锟斤拷锟斤拷锟斤拷锟 -3.5G锟斤拷锟斤拷 
>> 锟 20M锟斤拷锟斤拷
>> http://cn.mail.yahoo.com/
>>
>
> --
>
>
>       	
> ___________________________________________________________
> 抢注雅虎免费邮箱-3.5G容量，20M附件！
> http://cn.mail.yahoo.com
> <Thanks.rar>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20061201/b3ed930f/attachment.html