[gradsusr] Dear All,

Andrew Friedman andfried at gmail.com
Tue May 26 05:32:57 EDT 2015


Hi Atiqah,

I believe you want to use the function tcorr for temporal correlation instead of score:

http://www.iges.org/grads/gadoc/gradfunctcorr.html

Best,
Andrew

On May 26, 2015, at 1:12 AM, syairah atiqah atiq <syairahatiqah at yahoo.com> wrote:

> 
> Hi Andrew, 
> 
> Thank you for the reply. Okay, I attempt to find the relationship for two variables (anomalies sea ice cover and normalized rainfall) at different point of latitude and longitude (antarctic and maritime continent) in GrADs, but the answer was not exactly what I want. I know there was something wrong with my script and I do need help with that.Plus,  I really new to this software. Below was the script I wrote:
> 
> 'sdfopen erai_ci_1979_2014_anom.nc '
> 'sdfopen prate.sfc.mon.mean.nc'
> 'q file 1'
> 'q file 2'
> 'nv1=re2(icec.1,0.25,0.25,ba)'
> 'nv2=re2(prate.2,0.25,0.25,ba)'
> 'set dfile 1'
> 'set lon 0 360'
> 'set lat -90 -60'
> 'set dfile 2'
> 'set lon 90 165'
> 'set lat -12 20'
> 'climprat=ave(nv2*24*60*60*30,t=1,t=432,12)'
> 'anom=nv2*24*60*60*30-climprat'
> 'stDv=sqrt(avg(pow(climprat-nv2,2)),t=1,t=432,12)'
> 'Norm=anom/stDv'
> 'set gxout shaded'
> 'set t 1 12'
>  'set grads off'
>  'set grid off'
>  'd scorr(nv1,anom,lon=0,lon=360,lat=-90,lat=20)'
> 
> However, the answer shows was a value for each month. So in may case, I have 12 values appear in my screen. Hope you can help me.Kindly appreciate it.
> 
> Thank you in advance
> 
> Atiqah
> 
> 
> 
> 
> 
> On Tuesday, May 26, 2015 12:36 AM, Andrew Friedman <andfried at gmail.com> wrote:
> 
> 
> Hi Atiqah,
> 
> Can you please post specifically what you tried and what error you received?
> 
> Andrew
> 
> On May 25, 2015, at 4:44 AM, syairah atiqah atiq <syairahatiqah at yahoo.com> wrote:
> 
> > I am Atiqah. I attempt to find the relationship of two variables by doing the spatial correlation of sea ice cover in Antarctica and rainfall at Maritime Continent. The command I'm using was scorr but the result error. Can someone here help me in writing the script. I really new in this software.
> > 
> > Kindly help me. Thank you in advance.
> > 
> > best regards,
> > Atiqah
> 
> >  
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> 
> 
> 




More information about the gradsusr mailing list