RE: [ncl-talk] transect

From: Debasish Pai Mazumder <debasish_at_nyahnyahspammersnyahnyah>
Date: Wed, 21 Jun 2006 09:20:13 -0800

Dear Mateus da Silva Teixeira,

Please give this try,

trans =

according to example3 in

Good luck


-----Original Message-----
From: [] On Behalf
Of Mateus Teixeira
Sent: Wednesday, June 21, 2006 5:33 AM
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 =


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]
            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!

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"

  ; Abre arquivo de dados
  arquivo = addfile("./","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 =

Novidade no Yahoo! Mail: receba alertas de novas mensagens no seu celular.
Registre seu aparelho agora!

ncl-talk mailing list

ncl-talk mailing list
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:25 MDT