GrADS中计算水平螺旋度和垂直螺旋度的gs程序.doc

水平螺旋度的gs文件:'reinit''clear'='u';='v'pi='sdfopenf:\hyx\data\ncep\''sdfopenf:\hyx\data\ncep\'levstr='850700600500400300''setlon80110''setlat3050''settime06z20jul''definecr=uwnd'**************************************************************************************************i=1while(i<=9)j=1while(j<=13)===1while(k<=5)'setlev'subwrd(levstr,k)'undefinefld''definefld=''qdefvalfld'32+j''48+inwval=subwrd(result,3)=+nwval'undefinefld''definefld=''qdefvalfld'32+j''48+inwval=subwrd(result,3)=+nwval*say'nwval='nwval'vave.'i'.'j'='=k+=.j=(!=0)=math_atan(math_abs(.j))elseif(>0)=0endifif(<0)=piendifif(=0)say'chaoshasarrised!'endifendifif(>=0&>=0);=-;endif;if(<0&>=0);=180--;endif;if(<=0&<=0);=180-+;endif;if(>0&<=0);=270-+;endif;*if(>=0&>=0);=-*;endif;*if(<0&>=0);=(180-)*-;endif;*if(<=0&<=0);=(180-)*+;endif;*if(>0&<=0);=(270-)*+;endif;=math_sqrt(