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

[dongdong]

Thank you Katherin, what's more, The fwrite variables do not consistent with what I defined. Perhaps I should change the title.

I wrote the defined water flux out and reopened it in grads. Surprisingly, the figures are not identical.

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 any one help me finding out the problem?
Thank you very much! Any suggestion will be greatly appreciated!
Best Regrads

dongdong

----- 原始邮件 ---- 
发件人： Katherin Kullgren <kullgren at CRCES.ORG> 
收件人： GRADSUSR at LIST.CINECA.IT 
已发送： 2006/12/1(周五), 下午10:00:10 
主题： Re: 回复： how to get seasonal time series? 



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>


	

	
		
___________________________________________________________ 
Mp3疯狂搜-新歌热歌高速下
http://music.yahoo.com.cn/?source=mail_mailbox_footer
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Thanks.rar
Type: application/octet-stream
Size: 99663 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20061202/b469ba99/attachment.obj