修改脚本

好好学习1

'reinit'
'c'
'sdfopen E:\ncep\uwnd.sig995.mon.mean.nc'
'sdfopen E:\ncep\vwnd.sig995.mon.mean.nc'
'define slp=hcurl(uwnd.1,vwnd.2)'
'set lon 90 190'
'set lat 0 50'
'set t 1 700'
'define p11=ave(slp,t=151,t=154)'
'define p12=ave(slp,t=163,t=166)'
'define p13=ave(slp,t=175,t=178)'
'define p14=ave(slp,t=187,t=190)'
'define p15=ave(slp,t=199,t=202)'
'define p16=ave(slp,t=211,t=214)'
'define p17=ave(slp,t=223,t=226)'
'define p18=ave(slp,t=235,t=238)'
'define p19=ave(slp,t=247,t=250)'
'define p110=ave(slp,t=259,t=262)'
'define p111=ave(slp,t=271,t=274)'
'define p112=ave(slp,t=283,t=286)'
'define p113=ave(slp,t=295,t=298)'
'define p114=ave(slp,t=307,t=310)'
'define p115=ave(slp,t=319,t=322)'
'define p1=(p11+p12+p13+p14+p15+p16+p17+p18+p19+p110+p111+p112+p113+p114+p115)/15'
'define p21=ave(slp,t=331,t=334)'
'define p22=ave(slp,t=343,t=346)'
'define p23=ave(slp,t=355,t=358)'
'define p24=ave(slp,t=367,t=370)'
'define p25=ave(slp,t=379,t=382)'
'define p26=ave(slp,t=391,t=394)'
'define p27=ave(slp,t=403,t=406)'
'define p28=ave(slp,t=415,t=418)'
'define p29=ave(slp,t=427,t=430)'
'define p210=ave(slp,t=439,t=442)'
'define p211=ave(slp,t=451,t=454)'
'define p212=ave(slp,t=463,t=466)'
'define p2=(p21+p22+p23+p24+p25+p26+p27+p28+p29+p210+p211+p212)/12'
'define dif1=p2-p1'
'define p31=ave(slp,t=475,t=478)'
'define p32=ave(slp,t=487,t=490)'
'define p33=ave(slp,t=499,t=502)'
'define p34=ave(slp,t=511,t=514)'
'define p35=ave(slp,t=523,t=526)'
'define p36=ave(slp,t=535,t=538)'
'define p37=ave(slp,t=547,t=550)'
'define p38=ave(slp,t=559,t=562)'
'define p39=ave(slp,t=571,t=574)'
'define p3=(p31+p32+p33+p34+p35+p36+p37+p38+p39)/9'
'define dif2=p3-p2'
'define p41=ave(slp,t=583,t=586)'
'define p42=ave(slp,t=595,t=598)'
'define p43=ave(slp,t=607,t=610)'
'define p44=ave(slp,t=619,t=622)'
'define p45=ave(slp,t=631,t=634)'
'define p46=ave(slp,t=643,t=646)'
'define p47=ave(slp,t=655,t=658)'
'define p48=ave(slp,t=667,t=670)'
'define p49=ave(slp,t=679,t=682)'
'define p410=ave(slp,t=691,t=693)'
'define p4=(p41+p42+p43+p44+p45+p46+p47+p48+p49+p410)/10'
'define dif3=p4-p3'
'set t 1'

's1=( (p11-p1)*(p11-p1)+(p12-p1)*(p12-p1)+(p13-p1)*(p13-p1)+(p14-p1)*(p14-p1)+(p15-p1)*(p15-p1)+(p16-p1)*(p16-p1)+(p17-p1)*(p17-p1)+(p18-p1)*(p18-p1)+(p19-p1)*(p19-p1)+(p110-p1)*(p110-p1)+(p111-p1)*(p111-p1)+(p112-p1)*(p112-p1)+(p113-p1)*(p113-p1)+(p114-p1)*(p114-p1)+(p115-p1)*(p115-p1) )/15.*14.'
's2=((p21-p2)*(p21-p2)+(p22-p2)*(p22-p2)+(p23-p2)*(p23-p2)+(p24-p2)*(p24-p2)+(p25-p2)*(p25-p2)+(p26-p2)*(p26-p2)+(p27-p2)*(p27-p2)+(p28-p2)*(p28-p2)+(p29-p2)*(p29-p2)+(p210-p2)*(p210-p2)+(p211-p2)*(p211-p2)+(p212-p2)*(p212-p2))/12.*11.'
's3=((p31-p3)*(p31-p3)+(p32-p3)*(p32-p3)+(p33-p3)*(p33-p3)+(p34-p3)*(p34-p3)+(p35-p3)*(p35-p3)+(p36-p3)*(p36-p3)+(p37-p3)*(p37-p3)+(p38-p3)*(p38-p3)+(p39-p3)*(p39-p3))/9.*8.'
's4=((p41-p4)*(p41-p4)+(p42-p4)*(p42-p4)+(p43-p4)*(p43-p4)+(p44-p4)*(p44-p4)+(p45-p4)*(p45-p4)+(p46-p4)*(p46-p4)+(p47-p4)*(p47-p4)+(p48-p4)*(p48-p4)+(p49-p4)*(p49-p4)+(p410-p4)*(p410-p4))/10.*9.'

't12=25.*(p1-p2)/(s1+s2)/(sqrt(1./15.+1./12.) )'
't23=19.*(p2-p3)/(s2+s3)/(sqrt(1./12. +1./9.) )'
't34=17.*(p3-p4)/(s3+s4)/(sqrt(1./9. +1./10.) )'
'ta=1.71'
'tb=1.73'
'tc=1.74'

'set vpage 1 7 5.0 7.3'
'set grads off'
'set grid off'
'set gxout shaded'
'set black -1.71 1.71'
'set rbcols 8 0 14'
'set clevs -0.06 0.06'
'd t12'
'set rbcols auto'
'set gxout contour'
'set cthick 5'
'd dif1'


'draw title (a) Vor-DIF1'
'set vpage 1 7 2.6 4.9'
'set grads off'
'set grid off'
'set gxout shaded'
'set rbcols 8 0 14'
'set clevs -0.06 0.06'
'd t23'
'set rbcols auto'
'set gxout contour'
'set cthick 5'
'd dif2'
'draw title (b) Vor-DIF2'
'set vpage 1 7 0.2 2.5'
'set grads off'
'set grid off'
'set gxout shaded'
'set rbcols 8 0 14'
'set clevs -0.06 0.06'
'd t34'
'set rbcols auto'
'set gxout contour'
'set cthick 5'
'd dif3'
'draw title (c) Vor-DIF3'
'set vpage off'

'enable print E:\ncep\Vor-dif.gmf'
'print'
'disable print'
红色区域可能存在问题，请求大神修改

付亚男

你报错的结果是什么？把所有无关紧要的修饰全部去掉，看看能不能显示出结果。

1099221723

'define slp=hcurl(uwnd.1,vwnd.2)'
直接对三维变量进行hcurl,貌似是不可以的
应该是先循环计算每一个时刻的vor，并存储到dat里，再写ctl读取