[gradsusr] problem in calculation of precipitation bias

[pyarimohan maharana]

Dear all,

i am new to this forum.  I am trying to calculate the precipitation bias
(regcm output -CRU ) and when i am running the script, it is ending with
segmentation fault.....
 this is the script


'reinit'

'set display color white'

'c'

'enable print test2.meta'

'open /export/data/RegCM_Pyari/pyari/RegCM/NNRP2run/1988/SRF.1988060100.ctl'

'open /export/data/RegCM_Pyari/pyari/RegCM/NNRP2run/1988/SRF.1988070100.ctl'

'open /export/data/RegCM_Pyari/pyari/RegCM/NNRP2run/1988/SRF.1988080100.ctl'

'open /export/data/RegCM_Pyari/pyari/RegCM/NNRP2run/1988/SRF.1988090100.ctl'

'sdfopen /export/data/GPCP.PRECIP.NC'

'define t1=ave(tpr.1,t=1,t=240)'

'define t2=ave(tpr.2,t=241,t=488)'

'define t3=ave(tpr.3,t=489,t=736)'

'define t4=ave(tpr.4,t=737,t=976)'

'define T=(t1+t2+t3+t4)/4.0'

'define P=ave(precip.5,t=90,t=93)'

'set grid off'

'set grads off'

'set mpdset hires'

'set lat 5 40'

'set lon 65 98'

'set gxout shaded'

'set clevs -5 -4 -3 -2 -1 0 1 2 3 5 8 11 14 17 20'

'set ccols 3 14 4 11 5 13 0 10 7 12 8 2 6 15'

'd T-P'

'set gxout contour'

'set cint 1'

'd T-P'

'run cbar.gs'

'close 5'

'close 4'

'close 3'

'close 2'

'close 1'

 'print'

'disable print'







can anyone help?

thanks.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101124/df3278f2/attachment-0003.html