PRO wind2xycoordinates,wdir,x,y,xwind,ywind

;wdir; degree

; transform xy coordinates to wind angle to  coordiante system

wdirtemp=270.0-wdir
;print,wdirtemp
rad=wdirtemp*3.1415926/180.0


; find a sign mistake in the coordinate transfoamtion in the code of
; /davis/s2/wang/d2wwg/code/fortran/footprintpaper/code/integral-flux.pro
; maybe not, becsaue the different definition of dx
; 10/14/2003, maybe right now.

xwind=x*cos(rad)+y*sin(rad)
ywind=-x*sin(rad)+y*cos(rad)

return
end