Re: transect

From: Mateus Teixeira <mtex2k3_at_nyahnyahspammersnyahnyah>
Date: Wed, 21 Jun 2006 17:04:56 -0300

Thanks Debasish, I did not pay attention to this!

Now, I'm having another problem! The NCEP reanalysis gives their
variables, in case of temperature, as

air ( time, level, lat, lon)

where lon = (/90., 87.5, 85., 82.5, ... , -82.5, -85., -87.5, -90/). So,
using this as it comes, the function linint2_point gives the following
warning messages:

warning:linint2_points: xi and yi must be monotonically increasing
warning:linint2_points: xi and yi must be monotonically increasing
warning:linint2_points: xi and yi must be monotonically increasing
warning:linint2_points: xi and yi must be monotonically increasing
warning:linint2_points: xi and yi must be monotonically increasing
warning:linint2_points: xi and yi must be monotonically increasing
warning:linint2_points: xi and yi must be monotonically increasing

Is there a way to invert this dimension in my array directly or should I
perform a DO loop to this task?

Thank you!

Mateus da Silva Teixeira

Debasish Pai Mazumder a écrit :
> Dear Mateus da Silva Teixeira,
>
> Please give this try,
>
> trans =
> linint2_points(t_ar&lon_T,t_ar&lat_T,t_ar,True,dist_at_gclon,dist_at_gclat,2)
>
> according to example3 in
> http://www.ncl.ucar.edu/Document/Functions/Built-in/linint2_points.shtml
>
> Good luck
>
> Debasish
>
>
> -----Original Message-----
> From: ncl-talk-bounces_at_ucar.edu [mailto:ncl-talk-bounces_at_ucar.edu] On Behalf
> Of Mateus Teixeira
> Sent: Wednesday, June 21, 2006 5:33 AM
> To: ncl-talk_at_ucar.edu
> Subject: transect
>
> Dear NCL users,
>
> I'm new in NCL language, so I'm having some problems to do things that
> seems to be easy!
> I'm trying to create a transect, following the example script that exist
> in NCL website, but some is going wrong with linint2_points(), that
> gives the following error message:
>
> fatal:linint2_points: The rightmost dimensions of fi must be nyi x nxi,
> where nyi and nxi are the lengths of yi and xi respectively
> fatal:Execute: Error occurred at or near line 30 in file transec.ncl
>
> This happens with
>
> trans =
> linint2_points(t_ar&lat_T,t_ar&lon_T,t_ar,True,dist_at_gclon,dist_at_gclat,2)
>
> where
>
> Variable: t_ar
> Type: float
> Total Size: 714816 bytes
> 178704 values
> Number of Dimensions: 3
> Dimensions and sizes: [z_T | 17] x [lat_T | 73] x [lon_T | 144]
> Coordinates:
> lat_T: [90..-90]
> lon_T: [ 0..357.5]
>
> I do not understand this message. I think that is all ok, but I get this
> error message!
> Please, see the entire script below!
>
> Thanks
> Mateus da Silva Teixeira
>
>
> Since that I'm learning NCL, I doing the transect by steps, so the
> script below is only a part of the Transect script present in website.
> The script :
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>
> begin
> ; Abre arquivo de dados
> arquivo = addfile("./T_Marco2006.nc","r")
> t = arquivo->air(0,:,:,:)
>
> t_ar = t*t_at_scale_factor+t_at_add_offset
> t_ar!0 = "z_T"
> t_ar!1 = "lat_T"
> t_ar!2 = "lon_T"
> t_ar&lat_T = t&lat
> t_ar&lon_T = t&lon
>
> lat1 = -35.
> lon1 = -60.
> lat2 = -25.
> lon2 = -48.
> npts = 100 ; numero de pontos no transect
>
> dist = gc_latlon(lat1,lon1, lat2,lon2, npts,2)
> pts = ispan(0,npts-1,1)*1.0
>
> printVarSummary(t_ar) ; informacoes de t_ar
>
> ; interpola os dados ao "grande circulo"
> trans =
> linint2_points(t_ar&lat_T,t_ar&lon_T,t_ar,True,dist_at_gclon,dist_at_gclat,2)
> end
>
>
>
> _______________________________________________________
> Novidade no Yahoo! Mail: receba alertas de novas mensagens no seu celular.
> Registre seu aparelho agora!
> http://br.mobile.yahoo.com/mailalertas/
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
>

                
_______________________________________________________
Abra sua conta no Yahoo! Mail: 1GB de espaço, alertas de e-mail no celular e anti-spam realmente eficaz.
http://mail.yahoo.com.br/
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jun 21 2006 - 14:04:56 MDT

This archive was generated by hypermail 2.2.0 : Wed Jun 21 2006 - 17:13:25 MDT