grads中计算水平螺旋度和垂直螺旋度的gs程序

《grads中计算水平螺旋度和垂直螺旋度的gs程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、GrADS中计算水平螺旋度和垂直螺旋度的gs程序水平螺旋度的gs文件:'reinit''clear'prefix.1='u';prefix.2='v'pi=3.1415926'sdfopenf:hyxdatacepuwnd.2003.nc''sdfopenf:hyxdatacepvwnd.2003.nc'levstr='850700600500400300''setlon80110''setlat3050''settime06z20jul''definecr=uwnd'*****************

2、*********************************************************************************i=1while(i<=9)j=1while(j<=13)uave.i.j=0.0vave.i.j=0.0k=1while(k<=5)'setlev'subwrd(levstr,k)'undefinefld''definefld=uwnd.1''qdefvalfld'32+j''48+inwval=subwrd(result,3)uave.i.j=uave.i.j

3、+nwval'undefinefld''definefld=vwnd.2''qdefvalfld'32+j''48+inwval=subwrd(result,3)vave.i.j=vave.i.j+nwval*say'nwval='nwval'vave.'i'.'j'='vave.i.jk=k+1endwhileuave.i.j=uave.i.j/5vave.i.j=vave.i.j/5if(vave.i.j!=0)qq.i.j=math_atan(math_abs(vave.i.j/uave.i.j))elseif(ua

4、ve.i.j>0)qq.i.j=0endifif(uave.i.j<0)qq.i.j=piendifif(uave.i.j=0)say'chaoshasarrised!'endifendifif(uave.i.j>=0&vave.i.j>=0);qq.i.j=qq.i.j-40.0;endif;if(uave.i.j<0&vave.i.j>=0);qq.i.j=180-40.0-qq.i.j;endif;if(uave.i.j<=0&vave.i.j<=0);qq.i.j=180-40.0+qq.i.j;endif;if(

5、uave.i.j>0&vave.i.j<=0);qq.i.j=270-40.0+qq.i.j;endif;*if(uave.i.j>=0&vave.i.j>=0);qq.i.j=qq.i.j-40.0*3.1415926/180.0;endif;*if(uave.i.j<0&vave.i.j>=0);qq.i.j=(180-40.0)*3.1415926/180.0-qq.i.j;endif;*if(uave.i.j<=0&vave.i.j<=0);qq.i.j=(180-40.0)*3.1415926/180.0+qq.

6、i.j;endif;*if(uave.i.j>0&vave.i.j<=0);qq.i.j=(270-40.0)*3.1415926/180.0+qq.i.j;endif;uave.i.j=math_sqrt(uave.i.j*uave.i.j+vave.i.j*vave.i.j)*math_cos(qq.i.j*3.1415926/180.0)*0.75vave.i.j=math_sqrt(uave.i.j*uave.i.j+vave.i.j*vave.i.j)*math_sin(qq.i.j*3.1415926/180.

7、0)*0.75*uave.i.j=math_sqrt(uave.i.j*uave.i.j+vave.i.j*vave.i.j)*math_cos(qq.i.j)*0.75*vave.i.j=math_sqrt(uave.i.j*uave.i.j+vave.i.j*vave.i.j)*math_sin(qq.i.j)*0.75j=j+1endwhilei=i+1endwhile***************************************************************************

8、********************************************************************************************************************************************************