tcorr

Hyacinth Cyprain Nnamchi hc.nnamchi at LASG.IAP.AC.CN
Thu Oct 29 00:06:57 EDT 2009


 Dear Abbas,
Before you display tcorr, you have to reset x, y, and t. I have make some changes to your code. I hope it helps or you can look at again. Also see the documentation page: http://www.iges.org/grads/gadoc/gradfunctcorr.html

Good luck,
hyacinth

*****************Correlation***************
   'reinit'
'sdfopen D:/Monsoon/DATA/Vwind1000-100mb-1Apr-31Oct2000-4timeDaily.nc'
'sdfopen D:/Monsoon/DATA/Hgt850mb-1Apr-30Sep2000-4timeDaily.nc'

    'set x 1'
   'set y  1'
    'set z 2'   *(please you have to confirm that your lev 850 is z=2!)
    'set t 245 612'
    'define p = aave(vwnd.1,lon=59,lon=63,lat=30,lat=36)'    *indicated .1 for the first var


    'set gxout contour'
    'set csmooth on linear'
    'set grid off'
    'set grads off'
    'set font 1'
    'set cint 0.1'
    'set clskip 1'
    'set xlevs 41 43 45 47 49 51 53 55 57 59 61 63 65 67 69 71 73 75'
    'set ylint 2'
    'set ylevs 23 25 27 29 31 33 35 37 39 41 43 45 47'
    'set cthick 12'
    'set xlopts 1 6 0.15'
    'set ylopts 1 6 0.15'
    'set clopts 1 6 0.13'
    'set mpdset hires'
    'set map 1 1 15'
    'set ccolor 1 1 6'

    'set lat 22 45'
    'set lon 40 76'
    'set z 2'
  'set t 1'
    'd tcorr(p,hgt.2,t=245,t=612)'
    'draw title Correlation'
    'printim D:/Monsoon/Correlation850hPa.gif x1200 y1100 white'
say






--



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20091029/d4c1b3c5/attachment.html 


More information about the gradsusr mailing list