【脚本编辑】请大家帮忙指出代码错误

dnh

想求PDO选取的十年高值年与低值年400hPa的风场差异，自己写的gs，由于不是大气专业，代码全是网上查，自己写的，不知道是否正确，希望大佬帮忙看一下，指出错误或者不足，十分感谢。

'reinit'
'sdfopen F:\School\NCEPNCAR\wind\data\uwnd.mon.mean.nc'
'sdfopen F:\School\NCEPNCAR\wind\data\vwnd.mon.mean.nc'
'set grid off'
'set grads off'
'set lon 60 110'
'set lat 20 50'
'set lev 400'
'set mpdset country1'

'set rgb 16 0 120 255'
'set rgb 17 0 160 255'
'set rgb 18 0 200 200'
'set rgb 19 0 210 140'
'set rgb 20 0 220 0'
'set rgb 21 160 230 50'
'set rgb 22 230 220 50'
'set rgb 23 230 175 45'
'set rgb 24 240 130 40'
'set rgb 25 240 80 40'
'set rgb 26 240 20 40'

*季风期PDO低值年风场
'set grid off'
'set grads off'
'define uave1973=ave(uwnd,t=306,t=309)'
'define uave1974=ave(uwnd,t=318,t=321)'
'define uave1975=ave(uwnd,t=330,t=333)'
'define uave1990=ave(uwnd,t=510,t=513)'
'define uave1991=ave(uwnd,t=522,t=525)'
'define uave1999=ave(uwnd,t=618,t=621)'
'define uave2000=ave(uwnd,t=630,t=633)'
'define uave2001=ave(uwnd,t=642,t=645)'
'define uave2008=ave(uwnd,t=726,t=729)'
'define uave2009=ave(uwnd,t=738,t=741)'
'define uaveyl=(uave1973+uave1974+uave1975+uave1990+uave1991+uave1999+uave2000+uave2001+uave2008+uave2009)/10'
'define vave1973=ave(vwnd.2,t=306,t=309)'
'define vave1974=ave(vwnd.2,t=318,t=321)'
'define vave1975=ave(vwnd.2,t=330,t=333)'
'define vave1990=ave(vwnd.2,t=510,t=513)'
'define vave1991=ave(vwnd.2,t=522,t=525)'
'define vave1999=ave(vwnd.2,t=618,t=621)'
'define vave2000=ave(vwnd.2,t=630,t=633)'
'define vave2001=ave(vwnd.2,t=642,t=645)'
'define vave2008=ave(vwnd.2,t=726,t=729)'
'define vave2009=ave(vwnd.2,t=738,t=741)'
'define vaveyl=(vave1973+vave1974+vave1975+vave1990+vave1991+vave1999+vave2000+vave2001+vave2008+vave2009)/10'

*季风期PDO高值年风场
'set grid off'
'set grads off'
'define uave1981=ave(uwnd,t=402,t=405)'
'define uave1983=ave(uwnd,t=426,t=429)'
'define uave1984=ave(uwnd,t=438,t=441)'
'define uave1986=ave(uwnd,t=462,t=465)'
'define uave1987=ave(uwnd,t=474,t=477)'
'define uave1992=ave(uwnd,t=534,t=537)'
'define uave1993=ave(uwnd,t=546,t=549)'
'define uave1995=ave(uwnd,t=570,t=573)'
'define uave1997=ave(uwnd,t=594,t=597)'
'define uave2003=ave(uwnd,t=666,t=669)'
'define uaveyh=(uave1981+uave1983+uave1984+uave1986+uave1987+uave1992+uave1993+uave1995+uave1997+uave2003)/10'
'define vave1981=ave(vwnd.2,t=402,t=405)'
'define vave1983=ave(vwnd.2,t=426,t=429)'
'define vave1984=ave(vwnd.2,t=438,t=441)'
'define vave1986=ave(vwnd.2,t=462,t=465)'
'define vave1987=ave(vwnd.2,t=474,t=477)'
'define vave1992=ave(vwnd.2,t=534,t=537)'
'define vave1993=ave(vwnd.2,t=546,t=549)'
'define vave1995=ave(vwnd.2,t=570,t=573)'
'define vave1997=ave(vwnd.2,t=594,t=597)'
'define vave2003=ave(vwnd.2,t=666,t=669)'
'define vaveyh=(vave1981+vave1983+vave1984+vave1986+vave1987+vave1992+vave1993+vave1995+vave1997+vave2003)/10'

*季风期PDO高低值年风场差异
'set grid off'
'set grads off'
'define uaveyhl=uaveyh-uaveyl'
'define vaveyhl=vaveyh-vaveyl'
'set gxout shaded'
'set clevs 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2'
'set ccols 16 17 18 19 20 21 22 23 24 25 26'
'set csmooth on'
'd mag(uaveyhl,vaveyhl)'
'set gxout vector'
*'d uaveyhl;vaveyhl'
'drawvc uaveyhl vaveyhl 10.0 0.95 0.6 2'
'cbarn 1 0 5.6 0.2'
*'draw title 10years of highest-lowest PDO values'
'printim F:\School\NCEPNCAR\PDOwind\季风期PDO高低值年风场差异77.tif white'
'c'

*非季风期PDO低值年风场
'define nsuwnd1973=(ave(uwnd,t=301,t=305)+ave(uwnd,t=310,t=312))/2'
'define nsuwnd1974=(ave(uwnd,t=313,t=317)+ave(uwnd,t=322,t=324))/2'
'define nsuwnd1975=(ave(uwnd,t=325,t=329)+ave(uwnd,t=334,t=336))/2'
'define nsuwnd1990=(ave(uwnd,t=505,t=509)+ave(uwnd,t=514,t=516))/2'
'define nsuwnd1991=(ave(uwnd,t=517,t=521)+ave(uwnd,t=526,t=528))/2'
'define nsuwnd1999=(ave(uwnd,t=613,t=617)+ave(uwnd,t=622,t=624))/2'
'define nsuwnd2000=(ave(uwnd,t=625,t=629)+ave(uwnd,t=634,t=636))/2'
'define nsuwnd2001=(ave(uwnd,t=637,t=641)+ave(uwnd,t=646,t=648))/2'
'define nsuwnd2008=(ave(uwnd,t=721,t=725)+ave(uwnd,t=730,t=732))/2'
'define nsuwnd2009=(ave(uwnd,t=733,t=737)+ave(uwnd,t=742,t=744))/2'
'define uaveylns0=(nsuwnd1973+nsuwnd1974+nsuwnd1975+nsuwnd1990+nsuwnd1991+nsuwnd1999+nsuwnd2000+nsuwnd2001+nsuwnd2008+nsuwnd2009)/10'
'define nsvave1973=(ave(vwnd.2,t=301,t=305)+ave(vwnd.2,t=310,t=312))/2'
'define nsvave1974=(ave(vwnd.2,t=313,t=317)+ave(vwnd.2,t=322,t=324))/2'
'define nsvave1975=(ave(vwnd.2,t=325,t=329)+ave(vwnd.2,t=334,t=336))/2'
'define nsvave1990=(ave(vwnd.2,t=505,t=509)+ave(vwnd.2,t=514,t=516))/2'
'define nsvave1991=(ave(vwnd.2,t=517,t=521)+ave(vwnd.2,t=526,t=528))/2'
'define nsvave1999=(ave(vwnd.2,t=613,t=617)+ave(vwnd.2,t=622,t=624))/2'
'define nsvave2000=(ave(vwnd.2,t=625,t=629)+ave(vwnd.2,t=634,t=636))/2'
'define nsvave2001=(ave(vwnd.2,t=637,t=641)+ave(vwnd.2,t=646,t=648))/2'
'define nsvave2008=(ave(vwnd.2,t=721,t=725)+ave(vwnd.2,t=730,t=732))/2'
'define nsvave2009=(ave(vwnd.2,t=733,t=737)+ave(vwnd.2,t=742,t=744))/2'
'define vaveylns0=(nsvave1973+nsvave1974+nsvave1975+nsvave1990+nsvave1991+nsvave1999+nsvave2000+nsvave2001+nsvave2008+nsvave2009)/10'

*非季风期PDO高值年风场
'set grid off'
'set grads off'
'define nsuwnd1981=(ave(uwnd,t=397,t=401)+ave(uwnd,t=406,t=408))/2'
'define nsuwnd1983=(ave(uwnd,t=421,t=425)+ave(uwnd,t=430,t=432))/2'
'define nsuwnd1984=(ave(uwnd,t=433,t=437)+ave(uwnd,t=442,t=444))/2'
'define nsuwnd1986=(ave(uwnd,t=457,t=461)+ave(uwnd,t=466,t=468))/2'
'define nsuwnd1987=(ave(uwnd,t=469,t=473)+ave(uwnd,t=478,t=480))/2'
'define nsuwnd1992=(ave(uwnd,t=529,t=533)+ave(uwnd,t=538,t=540))/2'
'define nsuwnd1993=(ave(uwnd,t=541,t=545)+ave(uwnd,t=550,t=552))/2'
'define nsuwnd1995=(ave(uwnd,t=565,t=569)+ave(uwnd,t=574,t=576))/2'
'define nsuwnd1997=(ave(uwnd,t=589,t=593)+ave(uwnd,t=598,t=600))/2'
'define nsuwnd2003=(ave(uwnd,t=661,t=665)+ave(uwnd,t=670,t=672))/2'
'define uaveyhns0=(nsuwnd1981+nsuwnd1983+nsuwnd1984+nsuwnd1986+nsuwnd1987+nsuwnd1992+nsuwnd1993+nsuwnd1995+nsuwnd2003+nsuwnd1997)/10'
'define nsvave1981=(ave(vwnd.2,t=397,t=401)+ave(vwnd.2,t=406,t=408))/2'
'define nsvave1983=(ave(vwnd.2,t=421,t=425)+ave(vwnd.2,t=430,t=432))/2'
'define nsvave1984=(ave(vwnd.2,t=433,t=437)+ave(vwnd.2,t=442,t=444))/2'
'define nsvave1986=(ave(vwnd.2,t=457,t=461)+ave(vwnd.2,t=466,t=468))/2'
'define nsvave1987=(ave(vwnd.2,t=469,t=473)+ave(vwnd.2,t=478,t=480))/2'
'define nsvave1992=(ave(vwnd.2,t=529,t=533)+ave(vwnd.2,t=538,t=540))/2'
'define nsvave1993=(ave(vwnd.2,t=541,t=545)+ave(vwnd.2,t=550,t=552))/2'
'define nsvave1995=(ave(vwnd.2,t=565,t=569)+ave(vwnd.2,t=574,t=576))/2'
'define nsvave1997=(ave(vwnd.2,t=589,t=593)+ave(vwnd.2,t=598,t=600))/2'
'define nsvave2003=(ave(vwnd.2,t=661,t=665)+ave(vwnd.2,t=670,t=672))/2'
'define vaveyhns0=(nsvave1981+nsvave1983+nsvave1984+nsvave1986+nsvave1987+nsvave1992+nsvave1993+nsvave1995+nsvave1997+nsvave2003)/10'

*非季风期PDO高低值年风场差异
'set grid off'
'set grads off'
'define uaveyhlns=uaveyhns0-uaveylns0'
'define vaveyhlns=vaveyhns0-vaveylns0'
'set gxout shaded'
'set clevs 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2'
'set ccols 16 17 18 19 20 21 22 23 24 25 26'
'set csmooth on'
'd mag(uaveyhlns,vaveyhlns)'
'set gxout vector'
*'d uaveyhlns;vaveyhlns'
'drawvc uaveyhlns vaveyhlns 10 0.95 0.6 2'
'cbarn 1 0 5.6 0.2'
*'draw title 10years of highest-lowest PDO values'
'printim F:\School\NCEPNCAR\PDOwind\非季风期PDO高低值年风场差异87.tif white'
'c'

dnh

版主，我不是故意发这个主题的，选择脚本编辑时一直让我选择主题分类，发送不出去，选择这个主题就能发送出去。谁能帮我改下主题。。