transect

From: Mateus Teixeira <mtex2k3_at_nyahnyahspammersnyahnyah>
Date: Wed, 21 Jun 2006 10:32:40 -0300

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
Received on Wed Jun 21 2006 - 07:32:40 MDT

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