13楼的是风向角度数是从正Y轴开始顺时针算,按照平常的习惯和三角函数里度数从正X轴开始逆时针算,改进为:
if(u(i,j,t,yr)==0.0.and.v(i,j,t,yr)>0.0) dir(i,j,t,yr)=180.0
if(u(i,j,t,yr)==0.0.and.v(i,j,t,yr)<0.0) dir(i,j,t,yr)=0.0
if(u(i,j,t,yr)>0.0.and.v(i,j,t,yr)==0.0) dir(i,j,t,yr)=270.0
if(u(i,j,t,yr)<0.0.and.v(i,j,t,yr)==0.0) dir(i,j,t,yr)=90.0
if(u(i,j,t,yr)==0.0.and.v(i,j,t,yr)==0.0) dir(i,j,t,yr)=999.0
if(u(i,j,t,yr)>0.0.and.v(i,j,t,yr)>0.0) dir(i,j,t,yr)=atan(v(i,j,t,yr)/u(i,j,t,yr))*180/pi
if(u(i,j,t,yr)<0.0.and.v(i,j,t,yr)>0.0) dir(i,j,t,yr)=180.0+atan(v(i,j,t,yr)/u(i,j,t,yr))*180/pi
if(u(i,j,t,yr)<0.0.and.v(i,j,t,yr)<0.0) dir(i,j,t,yr)=180.0+atan(v(i,j,t,yr)/u(i,j,t,yr))*180/pi
if(u(i,j,t,yr)>0.0.and.v(i,j,t,yr)<0.0) dir(i,j,t,yr)=360.0+atan(v(i,j,t,yr)/u(i,j,t,yr))*180/pi |