[gradsusr] Determine tcorr.

[Andrew Friedman]

Hi Atiqah,

I think one problem is that zICnew varies in space and time, but the first
argument of tcorr can only vary in time.

Andrew

On Tue, Jun 2, 2015 at 8:52 AM, syairah atiqah atiq <syairahatiqah at yahoo.com
> wrote:

> Can anyone help me with the below scripts. I get an error when I run it.
>
> 'reinit'
> *
> * sea ice cover;
> * -------------
> *
>   'sdfopen D:\atiqah\GrADs\Real\Data\montly\SeaIceCover_19792014.nc'
> *
> *
>   'set dfile 1'
>   'set lon 0 360'
>   'set lat -90 -55'
> *
>   'icMJJAS=ave(ci*100,t=5,t=417,12)'
> *
>   'clim5=ave(ci*100,t=5,t=413,12)'
>   'clim6=ave(ci*100,t=6,t=414,12)'
>   'clim7=ave(ci*100,t=7,t=415,12)'
>   'clim8=ave(ci*100,t=8,t=416,12)'
>   'clim9=ave(ci*100,t=9,t=417,12)'
> *
>   'climavg=clim5+clim6+clim7+clim8+clim9/5'
> *
>   'anomIC=icMJJAS-climavg'
> *
>
> 'avg2=(ave(pow(clim5-ci,2),t=5,t=413,12)+ave(pow(clim6-ci,2),t=6,t=414,12)+ave(pow(clim7-ci,2),t=7,t=415,12)+ave(pow(clim8-ci,2),t=8,t=416,12)+ave(pow(clim9-ci,2),t=9,t=417,12))/5'
> *
>   'stDv=sqrt(avg2)'
> *
>   'zIC=anomIC/stDv'
> *
> * MSLP MH;
> * ========
> *
>   'sdfopen D:\atiqah\GrADs\Real\Data\montly\MSLP_ST_TP_1979_2014.nc'
> *
>   'set dfile 2'
>   'set lon 40E 90E'
>   'set lat 25S 35S'
> *
>   'mslMJJAS=ave(msl,t=5,t=417,12)'
> *
>   'clim5=ave(msl,t=5,t=413,12)'
>   'clim6=ave(msl,t=6,t=414,12)'
>   'clim7=ave(msl,t=7,t=415,12)'
>   'clim8=ave(msl,t=8,t=416,12)'
>   'clim9=ave(msl,t=9,t=417,12)'
> *
>   'climavg=clim5+clim6+clim7+clim8+clim9/5'
> *
>   'anomMH=mslMJJAS-climavg'
> *
> *
>
> 'avg2=(ave(pow(clim5-msl,2),t=5,t=413,12)+ave(pow(clim6-msl,2),t=6,t=414,12)+ave(pow(clim7-msl,2),t=7,t=415,12)+ave(pow(clim8-msl,2),t=8,t=416,12)+ave(pow(clim9-msl,2),t=9,t=417,12))/5'
> *
>   'stDv=sqrt(avg2)'
> *
> *
>   'zMH=anomMH/stDv'
> *
> * set same grid for both file;
> * ==========================
> *
>   'zICnew=re2(zIC,0.25,0.25,ba)'
>   'zMHnew=re2(zMH,0.25,0.25,ba)'
> *
> * Correlation;
> * ===========
> *
>  'set t 1 12'
>  'set lon 80 200'
>  'set lat -90 28'
>  'define corr=tcorr(zICnew,zMHnew,t=5,t=417)'
>  'd corr'
>
> Error from TCORR:  1st arg must be 0-D
> Error from TCORR:  Error getting grids
> Operation Error:  Error from tcorr function
>   Error ocurred at column 1
> DEFINE error:  Invalid expression.
>
> Thank in advance
>
> Atiqah
>
>
> _______________________________________________
> 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/20150603/cf3cf726/attachment-0001.html