[ncl-talk] problem while gnerating surface temperature time series

Bill Ladwig ladwig at ucar.edu
Thu Jul 28 12:15:36 MDT 2016


Hi Muhammad,

If I'm understanding what you're trying to accomplish in your script, I
think you want something like this:

f = addfile("./wrfout_d04_2009-07-10_00:00:00.nc","r")

tsk = wrf_user_getvar(f, "TSK", -1) ; Extract all times for TSK

loc = wrf_user_ll_to_ij(f, 36.85693, 2.53115, True)

x = loc(0) - 1
y = loc(1) - 1

time_series = tsk(:, y, x)  # Time series for TSK at location x,y

You don't want to use the 'ndtooned' routine since this routine just
flattens out a multidimensional array in to a single dimensional array.
You want to slice along the time dimension.

Hope this helps,

Bill



On Thu, Jul 28, 2016 at 10:34 AM, Muhammad Omer Mughal <
m.mughal1 at postgrad.curtin.edu.au> wrote:

> Hi
>
> I am trying to generate the time series for surface temperature but I am
> not successful. Till now I have reached upto the attached script. I would
> appreciate your help and support.
>
>
>
> Muhammad Omer Mughal
> MSc BSc Mechanical Engineering
> PhD  Research Scholar
> Remote Sensing and Satellite Research Group
> Department of Imaging and Applied Physics
> Curtin University
>
> Curtin University
> Tel | +61 8 9266 7962
> Fax | +61 8 9266 2377
> Mobile | 0470 237 525
>
> Email | m.mughal1 at postgrad.curtin.edu.au <m.lynch at curtin.edu.au>
> Web | http://curtin.edu.au
>
> Curtin University is a trademark of Curtin University of Technology.
> CRICOS Provider Code 00301J (WA), 02637B (NSW)
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160728/7aa40109/attachment.html 


More information about the ncl-talk mailing list