[gradsusr] Plot wind using regression

[Andrew Friedman]

Hi,

What issue are you having? What specific line is causing the problem?

Andrew

> On Jan 13, 2017, at 9:12 AM, syairah atiqah atiq <syairahatiqah at yahoo.com> wrote:
> 
> Hello everyone.
> 
> I do need help to plot wing using regression. Below are my script:
> 
> * For HIP
> * *******
>     'open E:\IMPORTANT\Desktop\Latest\hip\sie.-2.ctl'
>     'set t 1 6'
>     'a=sie'
>     'close 1'
> *
>     'open E:\IMPORTANT\Desktop\Latest\wind\hp.mjjw.-2.ctl'
> *
>    'set grid off'
>    'set grads off'
> *
>    'set lon 20 115'
>    'set lat -75 0'
> *
> * For zonal wind
> * ==============
>    'define coeff1 = tregr(a,u,t=1,t=6)'
>    'define aavg = ave(a,t=1,t=6)'
>    'define mavg1 = ave(u,t=1,t=6)'
> *
>    'define reg1 = coeff1*(a-aavg)+mavg1'
> *
> * For meridional wind
> * ===================
>    'define coeff2 = tregr(a,v,t=1,t=6)'
>    'define aavg = ave(a,t=1,t=6)'
>    'define mavg2 = ave(v,t=1,t=6)'
> *
>    'define reg2 = coeff1*(a-aavg)+mavg2'
> *
> *  Anomalies plot
> * ===============
>    'set arrscl 0.6 6'
>    'set arrlab on'
>    'd skip(reg1,2,2);reg2'
> *
>    'set map 1 1 12'
>    'draw map'
> *
>    'printim C:\Users\atiqah\Desktop\reg.hp.uv.png white'  
> *
> 
> Thank you in advance
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr