[ncl-talk] Need Help

[Jack Ogaja]

Assuming that ‘everything’ else in your code is correct, the snippet(of your code) below simply explains your concern;

****
  do it = 0,ntimes-1	        ;Loop for the time: it= starting time
   time = it
    	…
          ; Calculate Wind Speed from Vectors
        	  spd = (u_plane*u_plane + v_plane*v_plane)^(0.5)
         	 windspd0(it)=spd(x,y) ;<=== instantaneous[compare to “dir0” below]

          ; Wind direction at 0 Meters
       		r2d = 45.0/atan(1.0)            ; conversion factor (radians to degrees)
       		dir = atan2(u_plane,v_plane) * r2d + 180
       		dir0 = dir(x,y)        ;<=== is always overwritten[compare to “windspd0(it)” above]
 ****
i.e. your directions will always be the last direction calculated in the loop.

Jack.

> On 20 Oct 2015, at 08:56, isakhar sakhar isakhar <isakhar.inside13 at gmail.com> wrote:
> 
> Dear NCL Team
> 
> I use the script in attachment to find Wind Speed and Wind Direction in a location. I can get the wind speed but I get odd result for wind direction because all same. The output also in the attachment.
> 
> Regards,
> Soares
> <output.txt><wrf_script.ncl>_______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk