Re: data set extracted using NCL and read_wrf_nc slightly varies

I am not familiar with read_wrf_nc, so hopefully the wrfhelp folks
can help there.

However, can you explain what you mean about the data being slightly
different? Do you mean the values themselves are different?

Perhaps you are not extracting the same subregion? When you use
this syntax:

  u = a->U(12,0:31,49,49)
  v = a->V((12,0:31,49,49)
  w=a->W(12,0:31,49,49)

these are index values. NCL is 0-based indexing, so "U(12,...)" is
actually the 13th value of the rightmost dimension.

--Mary

--Mary

On Wed, 25 Jun 2008, Basit Khan wrote:

> Hi there,
>
> I am forwarding this question to both ncl_talk and wrfhelp, because i think
> it is related to both forums, i don't know which one is more appropriate
> though.
>
> When I extract data (vertical profile or time series) using NCL, it is
> slightly different form the data extracted using WRF utility 'read_wrf_nc',
> for the same time periods and vertical levels. Could you please advice why
> there is difference between two extracted datasets
>
> For NCL i am using the following to extract data
> ----------------------
> a = addfile("wrfoutfile","r")
> u = a->U(12,0:31,49,49)
> v = a->V((12,0:31,49,49)
> w=a->W(12,0:31,49,49)
>
> print(u)
> print(v)
> print(w)
> -----------
>
>
> Regards.
>
> --
> ---
> Basit A. Khan
> PhD Candidate
> Centre for Atmospheric Research
> Department of Geography
> University of Canterbury Te Whare Wananga o Waitaha Private Bag 4800
> Christchurch 8020, New Zealand
> Tel: (643) 364 2987 Ext. 7912
> Fax: (643) 364 2907
> bak41_at_student.canterbury.ac.nz
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jun 25 2008 - 11:57:47 MDT