Undefined variable

Hi users,

My script below say:

fatal:Variable (nome) is undefined
fatal:Execute: Error occurred at or near line 42 in file ex1.ncl

I can't understand why it's occur because before I was using the same
script to create figures.

Can anyone help me?

Thanks,

Guilherme.

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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/ut_string.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/time_axis_labels.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"

begin

lin = 105120
col = 14

  ctrl =
asciiread("single_point-output-ctrl.txt",(/lin,col/),"float")
  var1 =
asciiread("single_point-output-vmax60.txt",(/lin,col/),"float")
  var2 =
asciiread("single_point-output-vmax120.txt",(/lin,col/),"float")
  var3 =
asciiread("single_point-output-vmax240.txt",(/lin,col/),"float")

  do coluna = 1,3

  if ( coluna .eq. 1) ncoluna = 10 end if
  if ( coluna .eq. 2) ncoluna = 11 end if
  if ( coluna .eq. 3) ncoluna = 13 end if

  if (ncoluna .eq. 4) then
     titulo = "Net Primary Productivity (NPP)"
     eixoy = "NPP (tonC ha~S~-1~N~)"
     nome = "NPP"
  end if

  if (ncoluna .eq. 7) then
     titulo = "Net Ecosystem Exchange (NEE)"
     eixoy = "NEE (tonC ha~S~-1~N~)"
     nome = "NEE"
  end if

  if (ncoluna .eq. 5) then
     titulo = "Heterotrofic Respiration (RH)"
     eixoy = "RH (tonC ha~S~-1~N~ )"
     nome = "RH"
  end if

  print("Variavel --> " + nome)

nt1=0
nt2=743

  x = ctrl(nt1:nt2,0)
  yctrl = ctrl(nt1:nt2,ncoluna)
  yvar1 = var1(nt1:nt2,ncoluna)
  yvar2 = var2(nt1:nt2,ncoluna)
  yvar3 = var3(nt1:nt2,ncoluna)

  data1 = new((/4,dimsizes(x)/),float )

  data1(0,:) = yctrl
  data1(1,:) = yvar1
  data1(2,:) = yvar2
  data1(3,:) = yvar3

  wks_type = "ps" ; "pdf" "eps" "epsi" "ps" "x11"
  wks_type@wkOrientation = "portrait"
; wks_type@wkPaperSize = "A4"
  wks_type@wkPaperWidthF = 10.5 ; in inches
  wks_type@wkPaperHeightF = 14.0 ; in inches
  wks = gsn_open_wks(wks_type,"FIG.")

  res = True
; res@gsnMaximize = True
  res@tiMainString = titulo
  res@gsnLeftString = ""
  res@gsnCenterString = ""
  res@gsnRightString = ""
  res@tiXAxisString = "Month"
  res@tiYAxisString = eixoy
  res@tmYLMinorOn = True ; desabilita (False) o
minortick do eixo y
  res@gsnYRefLine = 0.0 ; linha de referência no
valor zero
; res@trYMinF = -0.040 ; valor mínimo do eixo y

; res@trYMaxF = 0.060 ; valor máximo do eixo y

; res@tmYLMode = "Explicit"
; res@tmYLLabels =
(/"-0.030","-0.020","-0.010","0.0","0.010","0.020","0.030","0.040","0.050","0.060","0.070"/)
; res@tmYLValues =
(/-0.030,-0.020,-0.010,0.0,0.010,0.020,0.030,0.040,0.050,0.060,0.070/)

  res@pmLegendDisplayMode = "Always" ; liga a legenda
  res@xyExplicitLegendLabels =
(/"CTRL","VMAX60","VMAX120","VMAX240"/)
  res@xyLineThicknesses = (/4.0,2.0,2.0,2.0/) ; espessura das
linhas
  res@xyLineColors = (/"black","red","green","purple"/)

  res@xyDashPattern = 0 ; força todas as
linhas para estilo sólida
  res@xyMarkLineMode = "MarkLines"
; res@xyMarkers = (/12,11,16,3/)
  res@xyMarkers = (/1,1,1,1/)
  res@xyMarkerColors = (/"black","red","green","purple"/) ;
3 different colors
  res@gsnFrame = False ;
avança (True) ou nao (False) as figuras
  res@pmLegendParallelPosF = 0.87 ;
desloca legenda em x
  res@pmLegendOrthogonalPosF = -1.17
 ; desloca legenda em y
  res@pmLegendWidthF = 0.3 ;
largura da linha da legenda
  res@pmLegendHeightF = 0.15 ;
altura entre as palavras da legenda
  res@lgBoxMinorExtentF = 0.15 ;
Shorten the legend lines
  res@lgLabelFontHeightF = 0.015 ;
altera o tamanho da fonte da legenda
  res@lgPerimOn = False ;
True = liga caixa da legenda
  res@lgItemOrder = (/3,2,1,0/) ;
Reorder the legends
  res@tmXBLabelAngleF = 0.0 ;
rotaciona os valores do eixo x
  res@tmXBLabelConstantSpacingF = 0.90 ;
espaçamento horizontal das palavras do eixo x
  res@tmXBFormat = "f" ;
remove os pontos decimais do eixo x
  res@tmXBLabelStride = 1
  res@txFontHeightF = 0.01 ;
modifica texto esquerdo, central e direito
  res@tmYROn = True ;
retira o tickmark da direita da figura
  res@tmXTOn = True ;
retira o tickmark da esquerda da figura
  res@tmXBLabelFontHeightF = 0.02 ;
tamanho da fonte do eixo y
; res@tmXBMode = "Explicit"

; res@tmXBValues = (/1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20, \
;
21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40, \
;
41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60/)

; res@tmXBLabels = (/"Jan~C~2001"," "," "," "," "," "," "," "," "," ", "
", " " , \
; "Jan~C~2002"," "," "," "," "," "," "," "," "," ", "
", " " , \
; "Jan~C~2003"," "," "," "," "," "," "," "," "," ", "
", " " , \
; "Jan~C~2004"," "," "," "," "," "," "," "," "," ", "
", " " , \
; "Jan~C~2005"," "," "," "," "," "," "," "," "," ", "
", " "/)

  plot = gsn_csm_xy(wks,x,data1,res)

  system("convert -trim FIG."+nome+".ps FIG."+nome+".png ; rm -f
FIG."+nome+".ps *~ ")

; delete(wks) ; apaga o pdf

  end do

end

-- 
*Guilherme Martins*
http://guilhermemartins.hd1.com.br/

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat May 5 07:59:29 2012