Re: wind speed

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Wed Sep 15 2010 - 12:26:41 MDT

Use addfiles to read across all the files
http://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml

    u = fu[:]->uwnd
    v = fv[:]->vwnd

if type short
http://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml

    u = short2flt(fu[:]->uwnd)
    v = short2flt(fv[:]->vwnd)

    w = sqrt(u^2 + v^2)
    copy_VarMeta(u, w)
    w@long_name = "wind speed"

On 9/15/10 1:58 AM, Vishali P. wrote:
> Dear NCL users
>
> I am using 4x Daily NCEP/DOE AMIP-II Reanalysis wind data. The data
> files are separated by years, and, 2 files per year (one for component u
> and another for v).
>
> I would like to calculate wind speed, for each time step, from 1998 to
> 2004 with the formula
>
> Windi = sqrt(Ui^2+Vi^2)
>
> Then, create an output ( .nc ) of wind speed with all 4
> dimensions(time,lev,lat,lon).
>
>
>
> Thanks in advance
>
> Vp
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 15 11:26:46 2010

This archive was generated by hypermail 2.1.8 : Thu Sep 16 2010 - 11:05:24 MDT