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
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jun 21 2006 - 11:20:13 MDT
This archive was generated by hypermail 2.2.0 : Wed Jun 21 2006 - 17:13:41 MDT