[gradsusr] Determine tcorr.
syairah atiqah atiq
syairahatiqah at yahoo.com
Tue Jun 2 02:52:59 EDT 2015
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-DError from TCORR: Error getting gridsOperation Error: Error from tcorr function Error ocurred at column 1DEFINE error: Invalid expression.
Thank in advance
Atiqah
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150602/74b8bb29/attachment.html
More information about the gradsusr
mailing list